Research on Green Innovation Efficiency of Listed Chinese Energy Companies Based on Triple Bottom Line

,is paper, based on the panel data, measures the comprehensive efficiency of green innovation, the economic output efficiency of green innovation, the social output efficiency of green innovation, and the environmental output efficiency of green innovation of listed energy companies in China by means of the window DEA model from the perspective of the triple bottom line. Moreover, the impacts of environmental factors and random factors on the green innovation efficiency are analyzed. It is found that firstly, green innovation efficiency of listed Chinese energy companies shows a V-shaped change trend without excluding environmental factors and random factors. Secondly, environmental factors and random factors are the primary factors to promote the green innovation efficiency of energy listed companies in China. Finally, after removing the interference of environmental factors and random factors, comprehensive efficiency of green innovation of the listed energy companies in China is relatively high and the distribution is relatively concentrated, but there is a downward trend during the research period.


Introduction
Energy enterprises play a very important role in China's economic development and assume the historical responsibility of realizing the beneficial development strategy of China as well. Relevant statistical results show that China's energy production and consumption enjoy sustaining growth. In 2019, the total energy production in China is 3.97 billion tons of standard coal, which increases by 5.3%, compared to that in 2018. Meanwhile, total energy consumption in 2019 reaches 4.86 billion tons of standard coal, which increases by 4.7% compared to that in 2018 [1], as is shown in Figure 1.
However, with the increasing of the energy production and consumption in China, some problems related to ecological and environmental governance appear. erefore, 2018 China Energy Development Report points out that China should continuously improve energy efficiency and develop towards high quality after basically meeting the demand of scale [2]. To retain the sustainable and high quality development of the energy industry, it is necessary to provide the enterprise innovation ability as soon as possible, improving the efficiency of green innovation of enterprises. In particular, it is beneficial to transform the energy companies from traditional extensive form to technology-oriented innovative type, allocate the resources, and promote the development of the green economy. Innovative development and green development are the key to transformation of energy enterprises. e green innovation of energy enterprises can break the constraints of resources and environment. erefore, it is necessary to objectively evaluate the innovation efficiency of green technology in Chinese energy enterprises under the condition that the R&D investment of energy enterprises continues to increase and innovation resources are seriously insufficient.
In accordance with the EU's Horizon 2020 Framework Plan, technological innovation should be expressed as social hope, safe, and sustainable [3]. erefore, in the process of sustainable development of enterprises, environmental responsibility and social responsibility are becoming as important as corporate profits. ese three factors constitute the triple bottom line (TLB) sustainable assessment framework and can evaluate the sustainable development of enterprises from a multidimensional perspective [4,5]. From a management perspective, efficiency refers to the ratio between the various inputs and outputs of an organization at a given time. Green innovation efficiency of energy enterprises from the perspective of triple bottom line also belongs to the category of input and output. at is to say, under the condition of limited investment in green innovation, energy enterprises should fulfill their economic responsibility, social responsibility, and environmental responsibility under the triple bottom line principle.
As the main body of sustainable innovation activities, it is necessary to discuss the efficiency of green innovation from the perspective of energy enterprises and from the perspective of the triple bottom line, but there is a lack of relevant research. At present, although some studies focus on the energy industry or the innovation efficiency of the energy industry, they do not focus the sustainable assessment framework and do not distinguish the internal characteristics of the energy industry, such as the green innovation efficiency of energy enterprises in different directions must be significantly different. In addition, it is difficult for existing research methods to effectively solve the problem of efficiency estimation. e window DEA method can distinguish the expected output from the nonexpected output, which is in line with the connotation of green innovation efficiency of energy enterprises.
Based on the above analysis, this paper conducts the research from the following aspects: at first, the efficiency differences of green innovation subjects can be analyzed in a more detailed way by taking listed energy companies as the research object, and the number of samples can be expanded to analyze the microforms of green innovation of energy in a broader scope. Second, from the perspective of the TLB theory, the window DEA model is constructed to comprehensively measure and analyze the economic, social, and environmental efficiency of China's listed energy companies in the process of green innovation, which is helpful to guide the sustainable development of energy companies. Finally, the three-stage DDF-DEA model adopted in this paper can effectively eliminate the impact of environmental factors and random errors on green innovation efficiency of listed energy companies so as to more truly reflect the actual level of green innovation efficiency of listed energy companies. e study is arranged as follows: Section 2 discusses the literature on TBL and the evaluation of green innovation efficiency and puts forward the innovation points. Section 3 illustrates the window DEA three-stage model and the source of the empirical data set in details. Section 4 discusses the main findings of empirical analysis. e concluding section draws a conclusion and puts forward the direction of future study.

Literature Review
In recent years, the development of the energy industry has aroused people's concern about the ecological environment and pollution. At the same time, more and more countries and regions are paying attention to green development and sustainable innovation. Relevant studies mainly focus on the following aspects.

Concept of Green Innovation.
At present, studies on the concept and connotation of green innovation are still in the exploration stage. e concept of green innovation originated from "sustainable development" and was first used in 1980 by the International Union for Conservation of Nature and Natural Resources in its World Conservation Strategy Report. Subsequently, "sustainable development" was defined as the deeds meeting current needs without compromising the ability of future generations to meet their own   2 Complexity needs [6]. On this basis, the literatures related to sustainable innovation [7][8][9][10], ecological innovation [11][12][13][14][15], green innovation [16][17][18][19][20][21], and environmental innovation [22][23][24][25] are also gradually increasing. Since 2005, the international use of the concept of green innovation and related studies have increased rapidly. Driessen and Hillebrand argue that green innovation should not aim at reducing the environmental burden but should aim at producing significant environmental benefits [26]. Chen et al. believe that color innovation is equal to hardware or software innovation related to green products or green processes, including those related to energy conservation, pollution prevention, waste recycling, green product design, and enterprise environmental management technology [17]. According to Dai, the definition and connotation of "green innovation" in the academic community are unified, and "green" tends to be abused because of its over-conceptualization and vague connotation [27]. Liu sorted out the contents and basic conclusions of green innovation research and analyzed the rationality of the conclusions from the perspective of classification definition [28].
ere are also some studies that comprehensively evaluate the efficiency of green innovation through spatial econometrics [30] and entropy evaluation [42]. In terms of influencing factors of green innovation efficiency, it mainly focuses on macroenvironmental factors and microfactors. e macroenvironment mainly includes the institutional environment [29,43], natural resources [44,45], related industries [46], and related international trade relations [47][48][49]. At the microlevel, internal factors related to enterprises mainly include enterprise costs [50,51], enterprise goals, and social responsibilities [52,53]. It can be seen that empirical research on the efficiency of green innovation is varied according to different research problems.

Triple Bottom Line eory.
In 1997, Elkington first proposed the concept of the triple bottom line. He believed that in the field of responsibility, CSR can be divided into economic responsibility, environmental responsibility, and social responsibility [5,54]. e triple bottom line provides a framework for measuring business performance and has achieved economic, social, and environmental success [55]. At present, the triple bottom line has become an important part of enterprise and management decision and has been paid more and more attention in practice. e authors of [56][57][58] studied the importance of different enterprises or governments on profit, environmental, or social responsibility from the perspective of practical framework of sustainability. It shows that the triple bottom line is a practical research method to measure the enterprise's sustainable performance [59][60][61].
ere are also some relevant studies that combine the triple bottom line theory with accounting theory to study the company's equity [62], annual report disclosure [63,64], input-output [65], and other factors.
In short, the triple bottom line always takes economic, social, and environmental values as the main line, which provides a practical framework for the study of sustainable development.

Research Design
is paper mainly studies the energy green innovation efficiency of listed companies in China. In view of data availability, this paper selected the listed companies in Shanghai Stock Exchange and Shenzhen Stock Exchange as the research object.
A total of qualified 49 listed energy companies are screened, including Shenzhen Energy, Tunghsu Azure, Victory Stock, Nanjing Public, and Guangdong Electric Power Development. e research time range is from 2016 to 2018, in order to ensure the timeliness of the research conclusions.

Model Construction
3.1.1. Traditional DEA Model. As a common method of efficiency measurement, DEA is widely used, but its shortcomings are also criticized. e shortcomings of the traditional DEA methods are reflected in the following aspects: first, the model assumes that decision units are homogeneous. However, the environment of each decisionmaking unit is different, so the decision-making unit cannot be completely homogeneous, and its efficiency is more or less affected by environmental factors and random factors. It is difficult for traditional DEA to eliminate the influence of the above factors, and the accuracy of efficiency value is negatively reflected. Second, the DEA model analyzes the efficiency value by constructing the frontier surface, which is different in different periods, resulting in the non-interperiod comparability of DEA efficiency [66].
In order to solve the above problems, scholars at home and abroad have revised the traditional DEA model appropriately. e three-stage DEA model proposed by Fried et al. well eliminates the influence of environmental factors and random factors on decision-making units, thus ensuring the homogeneity of decision-making units [67]. In order to make interpersonal comparison of DEA efficiency, Charnes et al. proposed a window DEA model, which can be used to analyze the evolution law of the efficiency value of the decision-making unit over time [68].
In this paper, panel data are used to study the green innovation efficiency of listed energy companies in China from the perspective of TBL. In order to eliminate the impact of environmental factors and random factors on efficiency measurement and analyze the dynamic evolution characteristics and influencing factors of green innovation efficiency, a three-stage window DEA model is constructed.

ree-Stage Window DEA Model.
e three-stage DEA model proposed by Fried et al. (2002) is mainly aimed at cross-sectional data. If panel data are used, the three-stage DEA model must be modified. In this context, the threestage window DEA model constructed in this paper is proposed. Its core idea is as follows: (1) e efficiency is evaluated by means of window DEA. (2) e relaxation variables are decomposed by the panel SFA model. (3) e panel management inefficiencies separation formula is used to adjust the decision unit. (4) e panel Tobit model is used to analyze the influencing factors of DEA efficiency of homogeneous decision-making unit. e details are listed as follows.
(1) e First Stage. In the first stage, the DEA-BBC model is adopted to measure the efficiency. Since the efficiency values of different fronts cannot be compared, this paper proposes to adopt window DEA to measure the dynamic efficiency of the decision-making unit. e window DEA model compares each decision-making unit with the frontier formed by fixed-width windows, so as to achieve the comparability of the efficiency of all decision-making units in time series.
(2) e Second Stage. Relaxation variables can reflect the initial low efficiency. In the second stage, this paper mainly focuses on relaxation variables and considers that relaxation variables are composed of environmental factors, management inefficiency, and statistical noise. By means of SFA regression, the relaxation variables are decomposed into the above three effects. Panel data model can better control individual heterogeneity, increase freedom, provide more information, have better statistical characteristics, and more accurate estimation results. erefore, in the second stage, panel SFA regression is adopted in this paper: S nit � f Z nit ;β nit + ] nit + μ nit , i � 1, 2, . . . , I, n � 1, 2, . . . , N. (1) where S nit is the relaxation value of the n th input in i th decision-making unit at the t phase. Z nit is the environment variable, and β nit is the coefficient of the environment variable. v nit + μ nit is the mixed error term, v nit represents random interference, and μ nit represents management inefficiency. It is generally assumed that the random disturbance term v∼N (0, σ 2 v ) obeys the standard normal distribution, and the management inefficiency term μ∼N + (0, σ 2 μ ) obeys the normal distribution truncated at zero.
According to the relationship between management inefficiency and time, the panel SFA can be divided into time-invariant model and time-varying recession model. e time-invariant model assumes that management inefficiency does not change with time and is expressed as μ it � μ i . e corresponding formula for separating management inefficiency has changed as well: where (3) e time-varying recession model assumes that the management inefficiency changes with time. Battese and Coelli (1992) [69] where T is the last period of decision-making units and η is recession coefficient), and the corresponding management inefficiency separation formula is as follows: where e time-varying recession equation is used to verify the estimated η. If η is not significant and is not 0, then the timeinvariant model is used [70].
(3) e ird Stage. e adjusted input-output variables were analyzed by the DEA model again. Similarly, window DEA was used in this paper to obtain the efficiency value of the decision-making unit after eliminating environmental factors and random factors. At this time, the efficiency value conforms to the assumption of DEA homogeneity and is relatively objective and accurate.

Further Analysis.
In order to make the constructed three-stage window DEA model in this paper easy to understand, the model is further analyzed in combination with Figure 2. As is shown in Figure 2, in the first stage, the window DEA model was used to calculate the green innovation efficiency of listed energy companies in China. e obtained DEA efficiency value did not conform to the assumption of homogeneity of decision units, so the efficiency value was not accurate. According to the research of Cummins et al. [66], there are three factors that affect the initial efficiency value, namely, environmental factors, management factors, and random factors. Due to the existence of environmental factors and random factors, decision units are not homogeneous, so it is necessary to eliminate environmental factors and random factors. In the second stage, with the help of SFA regression estimation parameters, environmental factors and random factors are 4 Complexity eliminated to obtain the adjusted input variables (or output variables). After the homogenization process in the second stage, the green innovation efficiency of listed Chinese energy companies in the third stage is relatively objective, and in theory, the efficiency value at this time is only affected by the internal management factors of the decision unit.

Index of Innovation Input.
As for innovation input, most scholars select indicators from the perspective of "human resources and property." For example, in their study on the efficiency of green innovation in China's heavily polluted industries, Fang et al. [71] focused on the selection of indicators from the two aspects of manpower and financial input, which are measured by the number of R&D personnel and R&D expenditure, respectively. Referring to the previous research literature, this paper selects the innovation input index of the three aspects of manpower input, financial input, and material input. For the manpower input, the index selected is the number of R&D personnel; for the financial input, the index selected is R&D expenditure; and for the material input, the index selected is fixed capital investment, and the fixed asset balance at the end of the year is taken as its proxy variable.

Index of Innovation Output.
Under the background of TBL, the green innovation output of listed energy companies in China includes economic efficiency output, social efficiency output, and environmental efficiency output. e economic benefit output of listed energy companies is measured by annual operating revenue, the social benefit output of listed energy companies is measured by patent quality and annual public welfare donation, and the environmental benefit output is measured by whether there is a positive environmental governance program.

Environmental Indicators of Innovation.
e selection of innovative environmental indicators has an impact on the efficiency of decision-making unit (DMU), but there are factors that DMU cannot control as is shown below: (1) Level of economic development: its impact on green innovation efficiency of listed Chinese energy companies is mainly reflected in two aspects. First, in regions with a high level of economic development, local governments pay more attention to scientific and technological innovation and will try their best to provide high-tech energy enterprises with a good research and policy environment, promote the transformation of scientific research results, maintain a high level of patent quality, and improve the social benefits of green innovation of listed energy companies. Second, in regions with a high level of economic development, the concept of green environment is more popular, which will improve the threshold of environmental access for enterprises, actively promote listed companies to improve the environmental governance program, and improve the environmental benefits of green innovation for listed energy companies. Per capita GDP better reflects the level of economic development, and this paper measures the level of regional economic development by per capita GDP. (2) Openness: its impact on green innovation efficiency of listed energy companies in China is mainly reflected in two aspects. First, in recent years, the Chinese government has vigorously advocated Complexity 5 scientific research to go abroad and carried out cooperative research with internationally renowned research institutions and enterprises, so as to enhance the scientific and technological innovation capacity of Chinese enterprises. Regions with a high degree of openness to the outside world are conducive to local enterprises to actively carry out foreign scientific research cooperation, promote the transformation of scientific research achievements of listed companies, and improve the social benefits of green innovation of listed energy companies. Second, in regions with high openness, the market economy is relatively active and the demand for energy products is relatively high, which helps us to improve the economic benefits of green innovation of listed energy companies. According to the conventional practice of referring to the literature at home and abroad, the level of opening up is measured by the proportion of total import and export volume to GDP, that is, the degree of dependence on foreign trade. (3) Government support: its impact on green innovation efficiency is mainly reflected in two aspects.
On the one hand, the inherent externality of green innovation of listed energy companies may bring high R&D costs and risks to enterprise innovation, and appropriate government intervention and financial subsidies can make up for such defects and improve the social benefits of green innovation of listed energy companies. On the other hand, excessive government intervention may cause "crowding-out effect," which is not conducive to technological innovation of enterprises and reduces social benefits of green innovation of listed energy companies. Government subsidies can directly measure the government's support for listed energy companies. In order to eliminate the influence of scale, the ratio of government subsidies to total assets is used as a proxy indicator of government support.
e above definitions and data sources are summarized in Table 1.

Isotropy Test of Input-Output Indicators.
In the window DEA model, the output and input indicators should have the same change trend; that is, there is no input index increase in the output index decreases, or there is no input index decrease with the input index increase. e Pearson correlation coefficient test method was used to investigate the isotropy of input-output indicators. e test results are shown in Table 2. As shown in the table, the correlation coefficients of the three input variables and four output variables selected in this paper are all positive, and all basically pass the significance test. It can be seen that the innovative input-output index of listed energy companies in China can be empirically studied with the three-stage window DEA model.

Window DEA Efficiency Measure in the First Stage.
DEA-solver Pro5.0 software was used to calculate the green innovation efficiency (including comprehensive efficiency of green innovation, economic output efficiency of green innovation, social output efficiency of green innovation, and environmental output efficiency of green innovation) of 49 listed Chinese energy companies from 2016 to 2018 based on input orientation, assuming variable scale reward of DUM and setting window width d � 3. Table 3       Complexity and c statistics for model 3 is 0.999, which means c statistics of three models are still in the range 0∼1, and all pass the 1% level of significance test. It shows that the energy green innovation efficiency in the first stage is influenced by environmental factors and random factors mutual influence. Finally, η statistics for model 1 is 0, η statistics for model 2 is 0.005, and η statistics for model 3 is 0.009, which means η statistics pass the significance test. is shows that the management inefficiency term of the labor input equation is constant, so the constant coefficient model should be adopted. Management inefficiency of model 2 and model 3 has typical time-varying characteristics, and the timevarying recession model should be adopted. Overall, the panel SFA model set up in stage 2 is reasonable.

e Impact of Environmental Factors on the Green Innovation Efficiency
(1) Level of Economic Development. e regression coefficient of economic development level of the slack variable of manpower input is negative, and it passes the test of 1% significance level. e regression coefficients of financial input and material input relaxation variables are both positive and pass the test of significance level of 1% and 5%, respectively. e regression results show that the economic development level has different effects on different green innovation input variables. First of all, the impact on human resources investment, the higher the level of regional economic development, the more to absorb outstanding scientific research and innovation talent, so the innovation ability of scientific research team is stronger, making full use of their talents and reducing the redundant investment of scientific research personnel. Second, the higher the level of regional economic development, the more the attention will be paid to technological innovation and the concept of green development, and a large amount of R&D funds will be invested in green innovation research projects. However, research projects are difficult to reflect the results in a short term, and there is a certain lag period, resulting in redundant financial investment. Finally, the impact on material input, economic development level of the region, and the investment scale is relatively large and easy to cause investment redundancy.
(2) Openness. e regression coefficients of external openness to the relaxation variables of manpower input and financial input was negative, and both passed the test of 1% significance level. e regression coefficient of the material input relaxation variable did not pass the test of significance level. It can be seen that the degree of openness to the outside world has a positive impact on green innovation efficiency of listed energy companies in China mainly by influencing human and financial input. e main reason is that the higher the degree of openness to the outside world is, the more the exchanges and cooperation between the listed energy companies and foreign scientific research institutions will be, which can ensure the frontier, practicability, and value of scientific research activities, reduce the redundancy of human resources input and scientific research project establishment, and improve the efficiency of scientific and technological innovation.

(3) Government Supports.
e regression coefficients of government support for the relaxation variables of human input, financial input, and material input are all positive and pass the test of significance levels of 1%, 5%, and 5%, respectively. is shows that in the green innovation link of the listed energy companies in China, government intervention is more obvious than the scope of "appropriate intervention," which leads to redundant investment of researchers, redundant investment of research project funds, and redundant investment of research equipment, which has a negative impact on the green innovation of China's listed energy companies.
Due to the interference of environmental factors and random factors, 49 listed energy companies did not conform to the assumption of homogeneity, resulting in the great defects of the green innovation efficiency estimated in the first stage. In order to eliminate environmental factors and random factors, this paper uses the above management inefficiency separation formula to decompose environmental factors and random factors and adopts the following formula to put the decision unit in the same external environment, so as to achieve the assumption of homogeneity: where X A ni is the input after the adjustment, X ni is the input before the adjustment, [max(f(Z i ; β  Table 5 shows the values of green innovation economic output efficiency, social output efficiency, environmental output efficiency, and comprehensive efficiency of listed energy companies in China. As shown in the table, green innovations of listed energy companies based on the TBL are different and vary greatly.

Analysis on Economic Output Efficiency of Green
Innovation. Figure 3 shows Certainly, this is related to the country's strong promotion of green development. Some interest groups related to listed energy companies have some resistance to this, which will directly affect the economic output of listed energy companies. By region, in 2016, the efficiency value of green innovation economic output of listed energy companies in the eastern region was 0.996, the efficiency value of green  innovation economic output in the central region was 0.932, and the efficiency value of green innovation economic output in the western region was 0.934. In 2017, the output efficiency value of green innovation economy in the eastern region was 0.973, the output efficiency value of green innovation economy in the central region was 0.940, and the output efficiency value of green innovation economy in the western region was 0.923. In 2018, the output efficiency value of green innovation economy in the eastern region was 0.915, the output efficiency value of green innovation economy in the central region was 0.898, and the output efficiency value of green innovation economy in the western region was 0.744, all of which showed a relatively obvious trend of efficiency decline, and the output efficiency of green innovation economy in the western region was the most obvious decline. From the final value of output efficiency of green innovation economy, it shows that "eastern > central > western," which is consistent with the traditional economic pattern. is benefits from the geographical environment and economic development level of eastern regions and central regions. e listed energy companies started their green innovation early and are characterized with high technology and rich innovation talented pool, which is beneficial to energy organization innovation activities of listed companies, and can effectively promote the efficiency of economic output. In addition, governments in the eastern and central regions invest more on R&D equipment and R&D funds, which will also improve the efficiency of enterprises' innovative economic output.

Analysis on Social Output Efficiency of Green
Innovation. Figure 4 shows the social output efficiency of green innovation of China's listed energy companies from 2016 to 2018. At the national level, the social output efficiency value of green innovation in 2016 was 0.778. In 2017, the output efficiency of green innovation society was 0.780, basically unchanged compared with the previous year. In 2018, the output efficiency value of the green innovation economy was 0.685, showing a significant downward trend. e main reason for the decline of green innovation output efficiency in 2018 is that although the number of patents of most listed energy companies has been increasing year by year, the quality has declined, leading to the fact that the output decreases when the input increases year by year. Relevant data show that the licensing rate of invention patents in China has dropped from 56.8 percent in 2015 to 40 percent in 2019, which also reflects the low pass rate and low quality of patents in China.
By region, in 2016, the social output efficiency values of green innovation of listed energy companies in the eastern region, the central region, and the western region were 0.739, 0.772, and 0.940, respectively, which shows the opposite of the traditional economic pattern "eastern > central > western." In 2017, the social output efficiency values of green innovation in the eastern region, the central region, and western region were 0.753, 0.734, and 0.937, respectively. In 2018, the social output efficiency values of green innovation in the eastern region, central region, and western region were 0.685, 0.675, and 0.697, respectively.
Overall, the green innovation social output efficiency in the western region is higher than that in the central region, which is in turn higher than that in the eastern region, while the green innovation output efficiency of the distribution is on the contrary. It can be seen that the listed energy companies in eastern areas are better at creating economic benefits, and the companies in the central areas and western areas are better at creating social benefits. is is because the western region has relatively little input, but its social output is relatively high. From 2016 to 2018, the green innovation social output efficiency of Datong Gas, Hongcheng Water Industry, Tianfu Energy, and other companies is 1, achieving the effectiveness requirement of DEA. ese enterprises basically reach a stable state in terms of social output. In addition, although Yinxing Energy, Chongqing Water, and Chongqing Gas are all reaching the effectiveness requirement of DEA, the social output efficiency of green innovation in most years is 1. erefore, listed energy companies in western China have strong ability in creating social benefits. For the central and eastern regions, firstly, there are more listed energy companies. Individual companies such as Shenzhen Energy, Yuedian Power A, Wanneng Power, Zhangze Power, China Nuclear Power, Datang Power Generation and other companies have large investment, resulting in social output e annual average efficiency is below 0.5. For example, Datang Power is only approximately 0.14, which has great impact on the comprehensive social output in the central eastern areas. In addition, influenced by government innovation policy, it appeared a downward trend in 2018.

Analysis on Environmental Output Efficiency of Green
Innovation. Figure 5  Aiming at these problems, the Ministry of Ecological Environment issued "Ecological Environmental Monitoring Quality Supervision and Inspection Action Plan (2018-2020)" in August 2018, which requires that by 2020, the enterprise needs to constantly improve the quality of ecological environment monitoring data security responsibility system, crack down on inappropriate behavior intervention ecological environment monitoring, effectively curb the ecological environment monitoring institutions and units data falsification problems, and ensure the monitoring data of true, accurate and objective. By region, the output efficiency value of green innovation environment of listed energy companies in the eastern region, the central region and the western region were 0.729, 0.726 and 0.842 in 2016, 2017, and 2018 respectively. e following trend is observed: "western > eastern > central". In 2017, the environmental output efficiency values of green innovation in the eastern region, the central region, and the western region were 0.702, 0.695, and 0.831 respectively. In 2018, the environmental output efficiency values of green innovation in the eastern region, the central region, and the western region were 0.635, 0.672, and 0.664 respectively, all of which showed an obvious trend of efficiency decline, and the most obvious change was in the western region. On the whole, environmental output efficiency of green innovation in the central region is relatively stable, and the environmental output efficiency of green innovation in the eastern and western regions fluctuates significantly. is is mainly due to the relatively small investment in the western region as a whole, which in turn directly affects the efficiency of its environmental output value. e eastern region, with the guide of national policies, constantly introduced effective environmental policies, so that the environmental management plans of the listed companies have positive effects as well, but its overall higher investment makes its final environment output lower than that in the western region. Moreover, there are more listed energy companies in the eastern region, and those companies with lower efficiency values directly affect the overall level of efficiency. erefore, based on the different regional characteristics of the eastern region, central region, and western region, the government should formulate corresponding environmental governance policies for listed energy companies.

Analysis on Comprehensive Efficiency of Green
Innovation. Figure 6  is is mainly affected by the global energy environment as a whole in 2018. At the same time, the macropolicy adjustment in China and increasing economic downward pressure makes the green innovation output efficiency of the listed energy company decline in 2018, which eventually led to a downward trend of the comprehensive efficiency value. However, it should be noticed that the overall comprehensive efficiency value is still high, which, from a certain extent, shows the scale effect of listed energy companies in the field of green innovation development gradually formed, and it is improving as well.
By region, in 2016, the comprehensive efficiency values of green innovation of listed energy companies in the eastern region, the central region, and the western region were 0.978, 0.996, and 1, which shows the "western > central > eastern" pattern. In 2017, the comprehensive efficiency values of green innovation of listed energy companies in the eastern region, the central region, and the western region were 0.993, 0.987, and 0.998. In 2018, the comprehensive efficiency values of green innovation of listed energy companies in the eastern region, the central region, and the western region were 0.943, 0.918, and 0.778. is shows that there is a distinct regional characteristic for the green innovation energy efficiency of the listed energy companies in China. e comprehensive green innovation efficiency value of the western region from 2016 to 2017 is significantly higher than those of central and eastern regions, which shows the effectiveness of pairing assistance policies in the early phase of "the 13 th Five-Year Plan" and the motivation of the environment with low competition pressure. Considering the   actual situation, the number of listed energy companies in the eastern and central regions is much higher than that in the western regions, which also has a great impact on the average value of comprehensive efficiency. However, in 2017-2018, the comprehensive efficiency values of Chongqing Water and Chongqing Gas were only 0.288 and 0.281, resulting in a large fluctuation of the comprehensive efficiency value in the western region. Among the three regions, only the western region had a comprehensive efficiency value of 1 in 2016, which reached the effectiveness requirement of DEA. Although some enterprises in other regions achieved their optimal efficiency, their overall levels were slightly inferior to the western region.
To reflect the space-time evolution path of energy green innovation comprehensive efficiency of listed companies in China, this paper constructs a GE matrix with the horizontal axis representing green innovation comprehensive efficiency changes and the vertical axis representing the initial value of green innovation comprehensive efficiency. e GE matrix divides the green innovation comprehensive efficiency of listed energy companies in China into nine types, as showed in Figure 7. e listed companies are mainly distributed in the four areas of the GE matrix. Among them, listed energy companies such as "Shanxi Gas," "Tunghsu Azure," "Jiangsu Guoxin," "Fuchun Environmental Protection," "Bohai Stock," "Tianhao Environment," "Zhangze Power," "Huajing Power Holdings," and "Huayin Power" are distributed in the range of "high initial efficiency value and decreasing efficiency." "Zhongshan Utilities," "Victory Stock," "Wanneng Power," "Yinxing Power," "Nanjing Utilities," "Solar Energy," "Shouchuang Stock," "Xintai Energy," "Datong Gas," "Mindong Power," "Huaneng International," "Shaoneng Stock," "Xingrong Environment," "Jiantou Energy," and "Shanghai Power" are distributed in the range of "high initial efficiency values and unchanged efficiency," among which the efficiency of "Yinxing Power" decreased at the early stage and then increased. Listed energy companies such as "Guozhong Water Supplies," "Disen Stock," "Datang Power," "Guangzhou Development," and "Capital Environment Protection" are distributed in the range of "high initial efficiency value and increasing efficiency." Listed energy companies such as "State Power," "Hongcheng Water Supplies," "Funeng Stock," "Chuangchun Gas," "Wuhan Stock," "Tianfu Energy," and "Shenneng Stock" are distributed in the range of "low initial efficiency value and increasing efficiency." Five other ranges of the GE matrix are dotted with publicly traded energy companies. "Jiangnan Water Supplies" and "China Nuclear Power" are located in the range of "medium initial efficiency value and decreasing efficiency." "Hubei Energy" is located in the range of "medium initial efficiency value and unchanged efficiency." ree energy companies, namely, "Guangdong Electric Power Development," "Lianmei Stock," and "Jinhong Stock," are located in "medium initial efficiency and decreasing efficiency." ree other energy companies "Chuantou Energy," "Ancai Hi-tech," and "Chongqing Gas" are located in the range of "Low initial value and decreasing efficiency," among which "Chuantou Energy" increases in the early stage and then decreases. In the range of "low initial efficiency and unchanged efficiency," there are three companies, namely,"Yangtze Power," "Shenzhen Gas," and "Chongqing Gas." It can be seen that there are different development environments for listed energy companies in China, leading to different changing characteristics of the comprehensive efficiency of green innovation. erefore, national and local governments should adopt different supporting policies for different enterprises to make full use of government resources.

Conclusions
Based on green innovation input-output data of 49 listed energy companies in China from 2016 to 2018, this paper studied the green innovation efficiency of listed energy companies in China and its influencing factors by constructing a three-stage window DEA model. e main conclusions are as follows.
Firstly, without eliminating environmental factors and random factors, green innovation economic output efficiency, green innovation social output efficiency, green innovation environmental output efficiency, and green innovation comprehensive efficiency of listed energy companies in China show a V-shaped trend.
Secondly, environmental factors and random factors have influenced the green innovation efficiency of listed Chinese energy companies. To be specific, the level of regional economic development reduces the redundancy of manpower input in green innovation of listed energy companies but increases the redundancy of financial input and material input in green innovation of listed energy companies. Openness is helpful to reduce the redundancy of manpower and financial resources in green innovation of listed energy companies, but it has little impact on the redundancy of material resources in green innovation of listed energy companies. e government's support will increase the redundant human, financial, and material input of listed energy companies in green innovation.
Finally, after removing the interference of environmental factors and random factors, green innovation economic output efficiency of the listed energy companies in China presents as "eastern > central > western" pattern which is consistent with the traditional economic pattern. Social output efficiency of green innovation in the western region is higher than that in the central region, which is in turn higher than that in the eastern region. Environmental output efficiency of green innovation in the central region is relatively stable, while that in the eastern and western regions fluctuates obviously. At the national level, the listed energy companies in China have high comprehensive efficiency of green innovation, but there is a downward trend during the research period. According to the distribution of GE matrix, 49 listed energy companies are mainly distributed in four areas, and just a few listed energy companies are scattered in the remaining five areas of GE matrix.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare no conflicts of interest. e funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.