CN117973707B - Carbon emission analysis method and system based on full-element productivity improvement - Google Patents

Carbon emission analysis method and system based on full-element productivity improvement Download PDF

Info

Publication number
CN117973707B
CN117973707B CN202410393923.1A CN202410393923A CN117973707B CN 117973707 B CN117973707 B CN 117973707B CN 202410393923 A CN202410393923 A CN 202410393923A CN 117973707 B CN117973707 B CN 117973707B
Authority
CN
China
Prior art keywords
department
matrix
production
factor
carbon emission
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202410393923.1A
Other languages
Chinese (zh)
Other versions
CN117973707A (en
Inventor
冯翠洋
赵佳龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology Beijing CUMTB
Original Assignee
China University of Mining and Technology Beijing CUMTB
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology Beijing CUMTB filed Critical China University of Mining and Technology Beijing CUMTB
Priority to CN202410393923.1A priority Critical patent/CN117973707B/en
Publication of CN117973707A publication Critical patent/CN117973707A/en
Application granted granted Critical
Publication of CN117973707B publication Critical patent/CN117973707B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The application relates to the technical field of carbon emission, and provides a carbon emission analysis method and system based on full-element productivity improvement. In the analysis method, based on a predetermined classification standard, the acquired multi-region input-output table is subjected to department matching with a carbon emission list, and an HA-IO matrix, a factor decomposition matrix, a multi-MAR weight of each department yield and a multi-MAR weight of each department production factor are determined so as to calculate the elasticity of different department yields to different department technologies; further, determining a corresponding department carbon emission variation when the full factor productivity varies; and determining the carbon emission intensity variation of each department when the productivity of the whole element is changed according to the carbon emission variation of each department so as to further determine the carbon emission key department. By the method, the influence conduction mechanism among all departments is cleared, the impact condition of all departments under the condition of technical change is determined, and key departments or industries for reducing carbon emission are identified.

Description

Carbon emission analysis method and system based on full-element productivity improvement
Technical Field
The application relates to the technical field of carbon emission, in particular to a carbon emission analysis method and system based on full-element productivity improvement.
Background
Carbon peaking and carbon neutralization are important strategic goals for economic and social development, and industry upgrades and technological advances are also important directions of development for modernization, but at present, the conduction mechanisms of specific impact of technological advances on carbon emissions remain unclear.
However, technical knowledge is a public product, because of the technology overflow effect and the network externality, the effect produced by the technological progress of one sector can produce a wide range of nationwide effects on carbon emissions through the value chain, and as the industry development system is gradually perfected, the relationship among the economic sectors is tighter, further amplifying the effect of the technological progress on the total carbon emissions. At present, the influence of technical progress on carbon emission is mainly limited to the influence of technical progress in a single field on the field, and an evaluation system for multi-area, multi-department and all-industry linkage is lacked.
Thus, there is a need to provide a solution to the above-mentioned deficiencies of the prior art.
Disclosure of Invention
The present application is directed to a carbon emission analysis method and system based on full factor productivity enhancement to solve or mitigate the above-mentioned problems of the prior art.
In order to achieve the above object, the present application provides the following technical solutions:
the application provides a carbon emission analysis method based on full-element productivity improvement, which comprises the following steps: step S101, based on a predetermined classification standard, performing department matching on the acquired multi-region input-output table and the carbon emission list;
Step S102, determining an HA-IO matrix, a factor decomposition matrix, multi-MAR weights of the output of each department and multi-MAR weights of the production factors of each department based on the input-output table and the carbon emission list after department matching;
Step S103, calculating the elasticity of the different department yields to the different department technologies according to the multi-MAR weights of the production factors of each department, the HA-IO matrix and the factor decomposition matrix;
Step S104, calculating corresponding department carbon emission variation when the production rate of the whole element is changed according to the elasticity of the production rate of different departments to the technology of different departments;
Step S105, determining the carbon emission intensity variation of each department when the productivity of the whole element is changed according to the carbon emission variation of each department;
And step S106, determining a carbon emission key department according to the carbon emission intensity variation of each department when the full-factor productivity is changed.
Preferably, in step S102, the formula is as follows:
construction of HA-IO matrix ; Wherein/>Respectively represent HA-IO matrix/>Different element matrix,/>Representing a zero matrix;
According to the formula:
determining element matrix Element/>
Wherein,Representation department/>Consumption department/>The total value of the product accounts for department/>The proportion of total value consumed for all departments' products; /(I)Representation department/>Is a price of the product; /(I)Representation department/>To department/>The consumption of the product;
According to the formula:
determining element matrix Element/>
Wherein,Representation department/>Consumption department/>The total value expenditure of the product accounts for department/>The proportion of total revenue; /(I)Representation department/>Accept from department/>The amount of product produced; /(I)Representation department/>Is a price of the product; /(I)Representation department/>Is a yield of (2); representation department/> Price of the product.
According to the formula:
determining element matrix Element/>
Wherein,Representation department/>In the production factor/>The expense on the way accounts for department/>Share of total revenue; /(I)Representing the production factor/>Is a price of (2); /(I)Representation department/>Provided production factor/>Is the number of (3); /(I)Respectively represent departments/>Price of product, department/>Yield of the product;
,/> Representing the number of departments,/> Are all positive integers.
Preferably, in step S102, the formula is as follows:
Constructing a factorial decomposition matrix
Wherein,Representing the production factor/>Revenue distribution to departments/>Share of/>Representation assigned to department/>Production factor/>Is a total value of (2); /(I)Representing the production factor/>Is a total value of (2);
,/> Representing the number of departments,/> Are all positive integers.
Preferably, in step S102, the formula is as follows:
Determining departments Multi-MAT weight of yield of (2)/>
Wherein,Representation department/>The total value of the product accounts for the share of the domestic total value GDP; /(I)Representation department/>Is a total output of (2);,/> Representing the number of departments,/> Are all positive integers;
According to the formula:
Determining production factors Multi-modular weights/>
Wherein,Representing the production factor/>The total value of the product accounts for the total value/>Is a fraction of (2); /(I)Representing the production factor/>Is a total value of (2); /(I),/>Representing the number of departments,/>Are all positive integers.
Preferably, in step S103, the formula is as follows:
Decomposition department A polymma weight of the yield of (a);
Wherein, Representing department/>, as producerThe produced product flows to department/>, which is consumerIs the sum of all paths of (a); /(I)For HA-IO matrix/>Middle/>Line/>Elements of a column; /(I)Representing HA-IO matrix/>Is presented as the freon-tigv inverse matrix/>/>Line/>Elements of a column; /(I)Representing the number of departments;
According to the formula:
Decomposing production factors Is a multi-magnitude weight of (2);
Wherein, Representing the production factor/>To department/>Is the sum of all paths of (a); /(I)Representation matrix/>/>Row of linesElements of a column; /(I)Rich-tigff inverse matrix/>, representing an HA-IO matrix/>Line/>Elements of a column;
,/> Representing the number of departments,/> Are all positive integers.
Preferably, in step S103, the formula is as follows:
Determining input-output covariance;
Wherein, Representation/>And/>Covariance between; /(I)Respectively represent HA-IO matrix/>Is presented as the freon-tigv inverse matrix/>(1 /)Column sum/>A column; /(I)Is an element matrix/>Middle/>Line/>Elements of a column; /(I)Representing HA-IO matrix/>Is presented as the freon-tigv inverse matrix/>Middle/>Line/>Column element,/>Representing HA-IO matrix/>Is of the inverse matrix of Rion-tigvMiddle/>Line/>Elements of a column;
According to the formula:
Determining the elasticity of the polymma weights of different department yields to different producer technologies;
Wherein, Representation department/>Multi-MAT weight of yield of (2)/>To department/>Production technique/>Elasticity of (a); /(I)Representation department/>The alternative elasticity of the input of (a); /(I)Representing HA-IO matrix/>Is presented as the freon-tigv inverse matrix/>(1 /)Column/>And/>Column/>Is a covariance of (2); /(I)Representing HA-IO matrix/>Is presented as the freon-tigv inverse matrix/>(1 /)Column/>And/>Column/>Is a covariance of (2); /(I)Representing the production factor/>Multi-modular weights/>To department/>Production technique/>Elasticity of (a); /(I)Representing the production factor/>Revenue distribution to departments/>Is a fraction of (2); /(I)Representing department/>, as producerThe produced product flows to department/>, which is consumerIs the sum of all paths of (a);
According to the formula:
Determining the elasticity of the polymma weights of the production factors to different producer technologies;
Wherein, Representing the production factor/>Multi-modular weights/>To department/>Production technique/>Elasticity of (a); /(I)Representation department/>A polymma weight of the yield of (a); /(I)Representation represents department/>The alternative elasticity of the input of (a); representation of HA-IO matrix/> Is presented as the freon-tigv inverse matrix/>(1 /)Column/>And/>Column ofIs a covariance of (2); /(I)Representing HA-IO matrix/>Is presented as the freon-tigv inverse matrix/>(1 /)Column/>And/>Column/>Is a covariance of (2); /(I)Representing the production factor/>Multi-modular weights/>To department/>Production technique/>Elasticity of (a); /(I)Representing the production factor/>Revenue distribution to departments/>Is a fraction of (2);
According to the formula:
Determining the elasticity of the price of the production factors to the technologies of different departments;
Wherein, Representing the production factor/>Price/>To department/>Production technique/>Elasticity of/>Representation department/>A polymma weight of the yield of (a); /(I)Representing the production factor/>Multi-modular weights/>To department/>Production technique/>Elasticity of (a);
According to the formula:
determining the elasticity of the yield of different producers to different department technologies;
Wherein, Representation department/>Yield/>To department/>Production technique/>Elasticity of (a); /(I)Representation department/>Multi-MAT weight of yield of (2)/>To department/>Production technique/>Elasticity of (a); /(I)Representing HA-IO matrix/>Is presented as the freon-tigv inverse matrix/>Middle/>Line/>Column element,/>Representing HA-IO matrix/>Is presented as the freon-tigv inverse matrix/>Middle/>Line/>Elements of a column; representing the production factor/> Price versus department/>Production technique/>Elasticity of (a); /(I)Representation department/>A polymma weight of the yield of (a);
,/> Representing the number of departments/> Are all positive integers.
Preferably, in step S104, the formula is as follows:
Determining departments Carbon emission intensity/>
Wherein,For department/>Carbon emission of/(Representation department/>Is a total output of (2); /(I),/>Indicating the number of departments to be reached,Are all positive integers.
Preferably, in step S105, the formula is as follows:
Determining departments Carbon emissions vs department/>Production technique/>Elasticity of (a);
Wherein, Is department/>Carbon emissions of (2); /(I)Representation department/>Yield/>To department/>Production technique/>Elasticity of (a); /(I)Representation department/>Carbon emission intensity of (2);
,/> Representing the number of departments,/> Are all positive integers;
According to the formula:
Determining departments Due to department/>Carbon emission variation/>, caused by variation of full-factor productivity
Wherein,For/>Diagonal matrix of/>Representation department/>Carbon emission of/(For/>Diagonalization of/>Representation department/>A change in the full factor productivity of (2); /(I)Representation department/>Is a production technology of (2).
Preferably, in step S106, the formula is as follows:
Determining due departments Variations in the full element productivity of (a) lead to department/>Sales amount matrix of (2)/>
Wherein,For/>Is a diagonal matrix of (a); /(I)Representation department/>Is a total output of (a); /(I)Representation department/>Is a production technology of (2); ; representation department/> A change in the full factor productivity of (2);
According to the formula:
Determining due departments Variations in the full element productivity of (a) lead to department/>Variation of carbon emission intensity/>
Wherein,Respectively represent departments/>Carbon emissions and department /)Is a total output of (a);
representation department/> Due to department/>Carbon emission variation caused by variation of the full-factor productivity;
,/> Representing the number of departments,/> Are all positive integers.
The embodiment of the application also provides a carbon emission analysis system based on full-element productivity improvement, which comprises the following components: the matching unit is configured to match the acquired multi-region input-output table with the carbon emission list in a department mode based on a predetermined classification standard;
The parameter determining unit is configured to determine an HA-IO matrix, a factor decomposition matrix, a multi-MAR weight of the output of each department and a multi-MAR weight of the production factors of each department based on the input-output table and the carbon emission list after the department matching;
The technology elastic unit is configured to calculate the elasticity of the different department yields to the different department technologies according to the multi-MAR weight of the production factors of each department, the HA-IO matrix and the factor decomposition matrix;
A variation determining unit configured to calculate a corresponding department carbon emission variation when the full-factor productivity varies, based on the elasticity of the different department technologies for the different department yields;
an emission intensity variation unit configured to determine carbon emission intensity variation amounts of respective departments when the full-factor productivity varies, based on the carbon emission variation amounts of the respective departments;
And a key section determining unit configured to determine a carbon emission key section based on the carbon emission intensity variation amounts of the respective sections when the full-factor productivity is changed.
The beneficial effects are that:
according to the carbon emission analysis method based on the full-factor productivity improvement, firstly, based on a predetermined classification standard, the acquired multi-region input-output table and the carbon emission list are subjected to department matching, and based on the input-output table and the carbon emission list after department matching, an HA-IO matrix, a factor decomposition matrix, multi-MAR weights of the output of each department and multi-MAR weights of the production factors of each department are determined; then, calculating the elasticity of the different department yields to the different department technologies according to the multi-MAR weights of the production factors of each department, the HA-IO matrix and the factor decomposition matrix; then, according to the elasticity of the output of different departments to the technology of different departments, calculating the corresponding carbon emission change amount of the departments when the productivity of the whole element changes, and determining the carbon emission intensity change amount of each department when the productivity of the whole element changes according to the carbon emission change amount of each department; finally, the carbon emission key departments are determined according to the variation of the carbon emission intensity of each department when the total biomass productivity is changed.
Therefore, the change of the production technology is described in a full-factor production mode, the limitation of the technology progress in the prior art in a specific industry is effectively avoided, and the technology progress of each department is described by a universal index; meanwhile, through quantitative analysis of the relation between the total element productivity change and the total carbon emission, the impact on the carbon emission of the industrial chain based on the technical change is transmitted among departments, the influence conduction mechanism among the departments is cleared, the impact condition of the departments under the technical change condition is determined, and then the key departments or industries for reducing the carbon emission are identified.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this specification, illustrate embodiments of the application and together with the description serve to explain the application. Wherein:
FIG. 1 is a flow diagram of a carbon emission analysis method based on full factor production enhancement provided in accordance with some embodiments of the present application;
fig. 2 is a schematic diagram of a carbon emission analysis system based on full factor production enhancement according to some embodiments of the present application.
Detailed Description
The application will be described in detail below with reference to the drawings in connection with embodiments. The examples are provided by way of explanation of the application and not limitation of the application. Indeed, it will be apparent to those skilled in the art that modifications and variations can be made in the present application without departing from the scope or spirit of the application. For example, features illustrated or described as part of one embodiment can be used on another embodiment to yield still a further embodiment. All other embodiments, which are derived by a person skilled in the art based on the embodiments of the present application, shall fall within the scope of protection of the embodiments of the present application.
The total element production rate (Total Factor Productivity, abbreviated as TFP) means the total amount of all elements input required for unit production when all production elements are used. In other words, an index for measuring production efficiency, which considers the number and quality of various elements in the production process, includes: comprehensive performance of labor, capital, science and technology, management, and the like. The higher the TFP, the higher the production efficiency, the higher the utilization efficiency of various elements in the production process, and more resources can be used to produce more output. TFP is therefore one of the important indicators of overall economic growth and development.
Aiming at the limitations among technical progress, carbon emission spatial distribution and different economic departments, the application quantitatively describes the technical progress by improving the productivity of all elements, and provides a carbon emission analysis method based on the productivity improvement of all elements.
First, it should be noted that,,/>Representing the number of departments; Are all positive integers. Department/> Department/>Respectively represent the/>Individual department, th/>Individual sector, production factor/>Represents the/>The meaning of the description of the production factors of the individual departments in the present application is not described in detail here.
As shown in fig. 1, the carbon emission analysis method based on the full factor productivity improvement includes:
And step S101, based on a predetermined classification standard, performing department matching on the acquired multi-region input-output table and the carbon emission list.
In the application, the departments of the multi-region input-output table are matched with the data of the carbon emission list by adopting the standard national economic industry classification standard, so that the departments of each region input-output table can be matched with the carbon emission data. Specifically, the departments in the multi-region input-output table are not necessarily consistent with the departments in the carbon emission list, and at this time, the departments in the input-output table and the departments in the carbon emission list need to be disassembled or combined, so that the departments in the input-output table can be consistent with the departments in the carbon emission list. For example, the departments in the input-output table are agriculture, forestry, animal husbandry and fishery (representing one department), but in the carbon emission list, 4 departments of agriculture, forestry, animal husbandry and fishery are required to be combined into agriculture, forestry, animal husbandry and fishery (1 department) in the carbon emission list, or the agriculture, animal husbandry and fishery (1 department) in the input-output table is split into 4 departments of agriculture, forestry, animal husbandry and fishery, so that the departments in the input-output table can be determined to be consistent with the departments in the carbon emission list through department matching.
Step S102, determining an HA-IO matrix, a factor decomposition matrix, multi-MAR weights of the output of each department and multi-MAR weights of the production factors of each department based on the input-output table and the carbon emission list after department matching.
Specifically, the formula is as follows:
construction of HA-IO matrix ; Wherein/>Respectively represent HA-IO matrix/>Different element matrix,/>Representing a zero matrix.
Further, according to the formula:
determining element matrix Element/>; Wherein/>Representation department/>Consumption department/>The total value of the product accounts for department/>The proportion of total value consumed for all departments' products; /(I)Representation department/>Is a price of the product; /(I)Representation department/>Opposite departmentsIs a product consumption amount.
According to the formula:
determining element matrix Element/>; Wherein/>Representation department/>Consumption department/>The total value expenditure of the product accounts for department/>The proportion of total revenue; /(I)Representation department/>Accept from department/>The amount of product produced; /(I)Representation department/>The value of the total output of the input-output table is the total output vector of the input-output table/>An element; /(I)Representation department/>Is a price of the product; /(I)Representation department/>Price of the product; /(I)Representation department/>Is a product of the above process.
According to the formula:
determining element matrix Element/>; Wherein/>Representation department/>In the production factor/>The expense on the way accounts for department/>Share of total revenue; /(I)Representing the production factor/>Is a price of (2); /(I)Representation department/>Provided production factor/>Is the number of (3); Respectively represent departments/> Total output of department/>Price of product, department/>Yield of the product.
In the application, the following formula is adopted:
Constructing a factorial decomposition matrix ; Wherein/>Representing the production factor/>Revenue distribution to departments/>Is used in the proportion of (a),Representation assigned to department/>Production factor/>Is a total value of (2); /(I)Representing the production factor/>Is a total value of (2).
In the application, the following formula is adopted:
Determining departments Multi-MAT weight of yield of (2)/>; Wherein/>Representation department/>The total value of the product accounts for the share of the domestic total value GDP; /(I)Representation department/>Is a total yield of (2).
According to the formula:
Determining production factors Multi-modular weights/>; Wherein/>Representing the production factor/>The total value of the product accounts for the total value/>Is a fraction of (2); /(I)Representing the production factor/>Is a total value of (2).
Step S103, calculating the elasticity of the different department yields to the different department technologies according to the multi-MAR weights of the production factors of the departments, the HA-IO matrix and the factor decomposition matrix.
According to the method, firstly, the elasticity of the multi-MAR weight of the production factors of each department, the multi-MAR weight of the production factors of each department and the price of the production factors of each department to the technologies of different departments is calculated according to the multi-MAR weight of the production factors of each department, the HA-IO matrix and the factor decomposition matrix; and then determining the elasticity of the output of different departments to the technologies of different departments according to the calculated multi-MAR weight of the output of each department, the multi-MAR weight of the production factors of each department and the elasticity of the price of the production factors to the technologies of different departments.
Specifically, first, according to the formula:
Opposite departments Multi-MAT weight of yield of (2)/>/>, Of production factorsMulti-modular weights/>And decomposing. Wherein,Representing department/>, as producerThe produced product flows to department/>, which is consumerA sum of all paths (flow paths); /(I)For HA-IO matrix/>Middle/>Line/>Elements of a column; /(I)Representing HA-IO matrix/>Is presented as the freon-tigv inverse matrix/>/>Line/>Elements of a column; /(I)Representing the production factor/>To department/>Is the sum of all paths of (a); /(I)Representation matrix/>/>Line/>Elements of a column; /(I)Rich-tigff inverse matrix/>, representing an HA-IO matrix/>Line/>Column elements.
And according to the formula:
determining input-output covariance; wherein, Representation/>And/>Is used to determine the covariance of (1),Respectively represent HA-IO matrix/>Is presented as the freon-tigv inverse matrix/>(1 /)Column sum/>A column; /(I)Is an element matrixMiddle/>Line/>Elements of a column; /(I)Representing HA-IO matrix/>Is presented as the freon-tigv inverse matrix/>Middle/>Line/>Column element,/>Representing HA-IO matrix/>Is presented as the freon-tigv inverse matrix/>Middle/>Line/>Column elements.
Then, according to the formula:
determining the elasticity of the polymma weights of different department yields to different producer technologies; wherein, Representation department/>Multi-MAT weight of yield of (2)/>To department/>Production technique/>Elasticity of (a); /(I)Representation department/>The alternative elasticity of the input of (a); representing HA-IO matrix/> Is presented as the freon-tigv inverse matrix/>(1 /)Column/>And/>Column/>Is a covariance of (2); /(I)Representing HA-IO matrix/>Is presented as the freon-tigv inverse matrix/>(1 /)Column/>And/>Column ofIs a covariance of (2); /(I)Representing the production factor/>Multi-modular weights/>To department/>Production technique/>Elasticity of (a); /(I)Representing the production factor/>Revenue distribution to departments/>Is a fraction of (2); /(I)Representing department/>, as producerThe produced product flows to department/>, which is consumerIs a sum of all paths of the (c).
According to the formula:
Determining the elasticity of the polymma weights of the production factors to different producer technologies; wherein, Representing the production factor/>Multi-modular weights/>To department/>Production technique/>Elasticity of (a); /(I)Representation department/>Is a multi-magnitude weight of (2); /(I)Representation represents department/>The alternative elasticity of the input of (a); /(I)Representation of HA-IO matrix/>Is of the inverse matrix of Rion-tigv(1 /)Column/>And/>Column/>Is a covariance of (2); /(I)Representing HA-IO matrix/>Is presented as the freon-tigv inverse matrix/>(1 /)Column/>And/>Column/>Is a covariance of (2); /(I)Representing the production factor/>Multi-modular weights/>To department/>Production technique/>Elasticity of (a); /(I)Representing the production factor/>Revenue distribution to departments/>Is a fraction of (a).
According to the formula:
Determining the elasticity of the price of the production factors to the technologies of different departments; wherein, Representing the production factor/>Price/>To department/>Production technique/>Elasticity of/>Representation department/>A polymma weight of the yield of (a); /(I)Representing the production factor/>Multi-modular weights/>To department/>Technique/>Is a spring of the same type.
Finally, according to the formula:
determining the elasticity of the yield of different producers to different department technologies; wherein, Representation department/>Yield/>To department/>Production technique/>Elasticity of (a); /(I)Representation department/>Multi-modular weights/>To department/>Production technique/>Elasticity of (a); /(I)Representing HA-IO matrix/>Is presented as the freon-tigv inverse matrix/>Middle/>Line/>Column element,/>Representing HA-IO matrix/>Is presented as the freon-tigv inverse matrix/>Middle/>Line/>Column element;/>Representing the production factor/>Price versus department/>Is produced by the production technology of (2)Elasticity of (a); /(I)Representation department/>Is a multi-magnitude weight of (a).
Step S104, calculating and determining corresponding department carbon emission variation when the productivity of the whole element is changed according to the elasticity of the output of different departments to the technologies of different departments.
In the application, firstly, the elasticity of the production of different departments to the technology of different departments is determined, and then the elasticity of the carbon emission of the department to the production technology of other departments is determined by the carbon emission intensity of the departments; finally, the department carbon emission variation of other departments due to the variation of the department full element productivity is determined by the flexibility of the department carbon emission to the production technology of other departments.
Specifically, the formula is as follows:
Determining departments Carbon emission intensity/>; Wherein/>For department/>Carbon emission of/(Representation department/>Is a total yield of (2).
According to the formula:
Determining departments Carbon emissions vs department/>Production technique/>Elasticity of (a);
Wherein, therein Is department/>Carbon emissions of (2); /(I)Representation department/>Yield/>To department/>Is produced by the production technology of (2)Elasticity of (a); /(I)Representation department/>Carbon emission intensity of (c).
Further, according to the formula:
Determining departments Due to department/>Carbon emission variation/>, caused by variation of full-factor productivity; Wherein/>Is thatDiagonal matrix of/>Representation department/>Carbon emission of/(For/>Diagonalization of/>Representation departmentA change in the full factor productivity of (2); /(I)Representation department/>Is a production technology of (2).
Step S105, determining the carbon emission intensity variation of each department when the productivity of the whole element is changed according to the carbon emission variation of each department.
Specifically, first, according to the formula:
Determining due departments Variations in the full element productivity of (a) lead to department/>Sales amount matrix of (2)/>; Wherein,For/>Is a diagonal matrix of (a); /(I)Representation department/>Is a total output of (a); /(I)Representation department/>Is a production technology of (2); /(I); Representation department/>Variation of the full factor production rate of (a).
Then, according to the formula:
/>
Determining due departments Variations in the full element productivity of (a) lead to department/>Variation of carbon emission intensity/>; Wherein/>Respectively represent departments/>Carbon emissions and department /)Is a total output of (a); /(I)Representation department/>Due to department/>Carbon emission variation caused by variation of the full-factor productivity.
And step S106, determining a carbon emission key department according to the carbon emission intensity variation of each department when the full-factor productivity is changed.
And finally, sorting the carbon emission intensity variation of all departments, and selecting the departments corresponding to the carbon emission intensity variation larger than the set emission threshold from the sorting result as carbon emission key departments.
Therefore, the change of the production technology is described in a full-factor production mode, the limitation of the technology progress in the prior art in a specific industry is effectively avoided, and the technology progress of each department is described by a universal index; meanwhile, through quantitative analysis of the relation between the total element productivity change and the total carbon emission, the impact on the carbon emission of the industrial chain based on the technical change is transmitted among departments, the influence conduction mechanism among the departments is cleared, the impact condition of the departments under the technical change condition is determined, and then the key departments or industries for reducing the carbon emission are identified.
As shown in fig. 2, an embodiment of the present application further provides a carbon emission analysis system based on full factor productivity improvement, including:
The matching unit is configured to match the acquired multi-region input-output table with the carbon emission list in a department mode based on a predetermined classification standard;
The parameter determining unit is configured to determine an HA-IO matrix, a factor decomposition matrix, a multi-MAR weight of the output of each department and a multi-MAR weight of the production factors of each department based on the input-output table and the carbon emission list after the department matching;
The technology elastic unit is configured to calculate the elasticity of the different department yields to the different department technologies according to the multi-MAR weight of the production factors of each department, the HA-IO matrix and the factor decomposition matrix;
A variation determining unit configured to calculate a corresponding department carbon emission variation when the full-factor productivity varies, based on the elasticity of the different department technologies for the different department yields;
an emission intensity variation unit configured to determine carbon emission intensity variation amounts of respective departments when the full-factor productivity varies, based on the carbon emission variation amounts of the respective departments;
and a key section determining unit configured to determine a carbon emission key section based on the carbon emission intensity variation of the full-factor productivity variation individual section.
The carbon emission analysis system based on the full-element productivity improvement provided by the embodiment of the application can realize the steps and the processes of the carbon emission analysis method based on the full-element productivity improvement described in any embodiment, and achieve the same technical effects, and are not described in detail herein.
The above description is only of the preferred embodiments of the present application and is not intended to limit the present application, but various modifications and variations can be made to the present application by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the protection scope of the present application.

Claims (2)

1. A carbon emission analysis method based on full factor productivity enhancement, comprising:
step S101, based on a predetermined classification standard, performing department matching on the acquired multi-region input-output table and the carbon emission list;
Step S102, determining an HA-IO matrix, a factor decomposition matrix, multi-MAR weights of the output of each department and multi-MAR weights of the production factors of each department based on the input-output table and the carbon emission list after department matching; wherein, according to the formula:
Constructing an HA-IO matrix omega; wherein, omega CN、ΩNN、ΩNF respectively represents different element matrixes in the HA-IO matrix omega, and O represents a zero matrix;
According to the formula:
Determining an element omega ci in an element matrix omega CN; wherein omega ci represents the proportion of the total value of the products of the department c consuming the department i to the total value of the products of the department c consuming all departments; p i represents the price of the product of department i; c ci represents the consumption of the product of department i by department c;
According to the formula:
Determining an element omega ji in an element matrix omega NN; wherein ω ji represents the proportion of the total value expenditure of the products of department j consuming department i to the total income of department j; x ji represents the number of products that department j accepts from department i produces; p j represents the price of the product of department j; y j represents the yield of department j; p i represents the price of department i product;
According to the formula:
Determining an element omega jf in an element matrix omega NF; wherein ω jf represents the fraction of the expenditure of department j on production factor f to the total revenue of department j; w f represents the price of the production factor f; l jf represents the number of production factors f provided by division j; p j、yj represents the price of the department j product, and the yield of the department j product, respectively;
i. c, j, f E N, N represents the number of departments, i, c, j, f, N are positive integers;
According to the formula:
Constructing a factor decomposition matrix phi cf; where phi cf represents the share of revenue of the production factor f allocated to the division c and Tv cf represents the total value of the production factor f allocated to the division c; tv f represents the total value of the production factor f;
According to the formula:
determining a polymma weight lambda j of the yield of the department j; wherein lambda j represents the share of the total value of the department j products to the domestic total value GDP; s j represents the total yield of department j;
According to the formula:
Determining a multi-MAR weight lambada f of a production factor f; wherein Λ f represents the share of the total value of the production factor f to the domestic total value GDP; tv f represents the total value of the production factor f;
step S103, calculating the elasticity of the different department yields to the different department technologies according to the multi-MAR weights of the production factors of each department, the HA-IO matrix and the factor decomposition matrix; wherein, according to the formula:
Decomposing the polymma weight of the yield of the department j; in the method, in the process of the invention, Representing the sum of all paths that the product produced by section i, which is the producer, flows to section c, which is the consumer; omega cj is the element of row c and column j in the HA-IO matrix omega; phi ji represents the elements of the j-th row and i-th column in the freon-tigv inverse matrix phi of the HA-IO matrix omega; n represents the number of departments;
According to the formula:
decomposing the polymma weight of the production factor f; in the method, in the process of the invention, Representing the sum of all paths from the production factor f to the department c; omega cj represents an element of row c and column j in matrix omega CN; phi jf represents the element of the j-th row and f-th column in the freon-tigff inverse matrix phi of the HA-IO matrix;
According to the formula:
Determining input-output covariance; where Cov Ω(j) represents the covariance between ψ (k) and ψ (l); the k-th column and the l-th column of the Rion-tigv inverse matrix ψ of the HA-IO matrix Ω are respectively represented by ψ (k)、Ψ(l); omega ji is the element of the j-th row i column in the element matrix omega NN; phi ik represents the elements of the ith row and k column in the freon-tigv inverse matrix psi of the HA-IO matrix omega, phi il represents the elements of the ith row and l column in the freon-tigv inverse matrix psi of the HA-IO matrix omega;
According to the formula:
determining the elasticity of the polymma weights of different department yields to different producer technologies; in the method, in the process of the invention, The elasticity of the multimma weight λ i representing the yield of division i to the production technique a k of division k; θ j represents the alternative elasticity of the input of department j; Covariance of the kth column ψ (k) and the ith column ψ (i) of the freon-tigv inverse matrix ψ35 and representing the HA-IO matrix Ω; /(I) Covariance of g-th column ψ (g) and i-th column ψ (i) of the freon-tigv inverse matrix ψ35 representing HA-IO matrix Ω; /(I)The elasticity of the multimma weight g representing the production factor g to the production technique a k of division k; phi cg represents the share of revenue from the production factor g allocated to department c; /(I)Representing the sum of all paths that the product produced by section i, which is the producer, flows to section c, which is the consumer;
According to the formula:
determining the elasticity of the polymma weights of the production factors to different producer technologies; in the method, in the process of the invention, The elasticity of the multi-MAT weight f representing the production factor f to the production technique A k of department k; lambda j represents the polymma weight of the yield of department j; θ j represents the alternative elasticity of the input of department j; /(I)Covariance of the kth column ψ (k) and the f column ψ (f) of the freon tigv inverse matrix ψ35 and representing the HA-IO matrix Ω; /(I)Covariance of g-th column ψ (g) and f-th column ψ (f) of the freon-tigv inverse matrix ψ35 representing HA-IO matrix Ω; /(I)The elasticity of the multimma weight g representing the production factor g to the production technique a k of division k; phi cg represents the share of revenue from the production factor g allocated to department c;
According to the formula:
determining the elasticity of the price of the production factors to the technologies of different departments; in the method, in the process of the invention, The elasticity of the price w f representing the production factor f to the production technology a k of division k, lambda k representing the polymma weight of the yield of division k; /(I)The elasticity of the multi-MAT weight f representing the production factor f to the production technique A k of department k;
According to the formula:
determining the elasticity of the yield of different producers to different department technologies; in the method, in the process of the invention, Representing the elasticity of the production y i of sector i to the production technique a k of sector k; /(I)The elasticity of the multimma weight λ i representing the yield of division i to the production technique a k of division k; phi ik represents the elements of the ith row and k column in the freon-tigv inverse matrix psi of the HA-IO matrix omega, phi ig represents the elements of the ith row and g column in the freon-tigv inverse matrix psi of the HA-IO matrix omega; /(I)The flexibility of the production technique A k to the department k, which represents the price of the production factor g; lambda k represents the polymma weight of the yield of division k;
k. l and g are N, and k, l and g are positive integers;
Step S104, calculating corresponding department carbon emission variation when the production rate of the whole element is changed according to the elasticity of the production rate of different departments to the technology of different departments; wherein, according to the formula:
Determining the carbon emission intensity e k of the department k; wherein E k is the carbon emission amount of the department k, and S k represents the total yield of the department k;
Step S105, determining the carbon emission intensity variation of each department when the productivity of the whole element is changed according to the carbon emission variation of each department; wherein, according to the formula:
Determining the elasticity of the carbon emissions of sector i to sector k production technique a k; wherein E i is the carbon emissions of sector i; Representing the elasticity of the production y i of sector i to the production technique a k of sector k; e k denotes the carbon emission intensity of division k;
According to the formula:
Determining a carbon emission variation H i,k of the department i due to the full factor productivity variation of the department k; in the method, in the process of the invention, For the diagonal matrix of E i, E i represents the carbon emissions of division i,/>For diagonalization of Δa k/Ak, Δa k/Ak represents the change in full element productivity for division k; a k represents a production technology of department k;
Step S106, determining carbon emission key departments according to the carbon emission intensity variation of each department when the full-factor productivity is changed; wherein, according to the formula:
Determining a sales change matrix T i,k for division i due to a change in the full-element productivity of division k; in the method, in the process of the invention, A diagonal matrix of S; s i represents the total yield of department i; a k represents a production technology of department k; Δa k/Ak represents the variation in the full element productivity of division k;
According to the formula:
Determining a change Δd i,k in the carbon emission intensity of division i due to a change in the full factor productivity of division k; wherein E i、yi represents the carbon emission amount of the department i and the total output amount of the department i, respectively;
H i,k represents the carbon emission variation of division i due to the full factor productivity variation of division k.
2. A carbon emission analysis system based on full factor production enhancement, comprising
The matching unit is configured to match the acquired multi-region input-output table with the carbon emission list in a department mode based on a predetermined classification standard;
The parameter determining unit is configured to determine an HA-IO matrix, a factor decomposition matrix, a multi-MAR weight of the output of each department and a multi-MAR weight of the production factors of each department based on the input-output table and the carbon emission list after the department matching; wherein, according to the formula:
Constructing an HA-IO matrix omega; wherein, omega CN、ΩNN、ΩNF respectively represents different element matrixes in the HA-IO matrix omega, and O represents a zero matrix;
According to the formula:
Determining an element omega ci in an element matrix omega CN; wherein omega ci represents the proportion of the total value of the products of the department c consuming the department i to the total value of the products of the department c consuming all departments; p i represents the price of the product of department i; c ci represents the consumption of the product of department i by department c;
According to the formula:
Determining an element omega ji in an element matrix omega NN; wherein ω ji represents the proportion of the total value expenditure of the products of department j consuming department i to the total income of department j; x ji represents the number of products that department j accepts from department i produces; p j represents the price of the product of department j; y j represents the yield of department j; p i represents the price of department i product;
According to the formula:
Determining an element omega jf in an element matrix omega NF; wherein ω jf represents the fraction of the expenditure of department j on production factor f to the total revenue of department j; w f represents the price of the production factor f; l jf represents the number of production factors f provided by division j; p j、yj represents the price of the department j product, and the yield of the department j product, respectively;
i. c, j, f E N, N represents the number of departments, i, c, j, f, N are positive integers;
According to the formula:
Constructing a factor decomposition matrix phi cf; where phi cf represents the share of revenue of the production factor f allocated to the division c and Tv cf represents the total value of the production factor f allocated to the division c; tv f represents the total value of the production factor f;
According to the formula:
determining a polymma weight lambda j of the yield of the department j; wherein lambda j represents the share of the total value of the department j products to the domestic total value GDP; s j represents the total yield of department j;
According to the formula:
Determining a multi-MAR weight lambada f of a production factor f; wherein Λ f represents the share of the total value of the production factor f to the domestic total value GDP; tv f represents the total value of the production factor f;
The technology elastic unit is configured to calculate the elasticity of the different department yields to the different department technologies according to the multi-MAR weight of the production factors of each department, the HA-IO matrix and the factor decomposition matrix; wherein, according to the formula:
Decomposing the polymma weight of the yield of the department j; in the method, in the process of the invention, Representing the sum of all paths that the product produced by section i, which is the producer, flows to section c, which is the consumer; omega cj is the element of row c and column j in the HA-IO matrix omega; phi ji represents the elements of the j-th row and i-th column in the freon-tigv inverse matrix phi of the HA-IO matrix omega; n represents the number of departments;
According to the formula:
decomposing the polymma weight of the production factor f; in the method, in the process of the invention, Representing the sum of all paths from the production factor f to the department c; omega cj represents an element of row c and column j in matrix omega CN; phi jf represents the element of the j-th row and f-th column in the freon-tigff inverse matrix phi of the HA-IO matrix;
According to the formula:
Determining input-output covariance; where Cov Ω(j) represents the covariance between ψ (k) and ψ (l); the k-th column and the l-th column of the Rion-tigv inverse matrix ψ of the HA-IO matrix Ω are respectively represented by ψ (k)、Ψ(l); omega ji is the element of the j-th row i column in the element matrix omega NN; phi ik represents the elements of the ith row and k column in the freon-tigv inverse matrix psi of the HA-IO matrix omega, phi il represents the elements of the ith row and l column in the freon-tigv inverse matrix psi of the HA-IO matrix omega;
According to the formula:
determining the elasticity of the polymma weights of different department yields to different producer technologies; in the method, in the process of the invention, The elasticity of the multimma weight λ i representing the yield of division i to the production technique a k of division k; θ j represents the alternative elasticity of the input of department j; Covariance of the kth column ψ (k) and the ith column ψ (i) of the freon-tigv inverse matrix ψ35 and representing the HA-IO matrix Ω; /(I) Covariance of g-th column ψ (g) and i-th column ψ (i) of the freon-tigv inverse matrix ψ35 representing HA-IO matrix Ω; /(I)The elasticity of the multimma weight g representing the production factor g to the production technique a k of division k; phi cg represents the share of revenue from the production factor g allocated to department c; /(I)Representing the sum of all paths that the product produced by section i, which is the producer, flows to section c, which is the consumer;
According to the formula:
determining the elasticity of the polymma weights of the production factors to different producer technologies; in the method, in the process of the invention, The elasticity of the multi-MAT weight f representing the production factor f to the production technique A k of department k; lambda j represents the polymma weight of the yield of department j; θ j represents the alternative elasticity of the input of department j; /(I)Covariance of the kth column ψ (k) and the f column ψ (f) of the freon tigv inverse matrix ψ35 and representing the HA-IO matrix Ω; /(I)Covariance of g-th column ψ (g) and f-th column ψ (f) of the freon-tigv inverse matrix ψ35 representing HA-IO matrix Ω; /(I)The elasticity of the multimma weight g representing the production factor g to the production technique a k of division k; phi cg represents the share of revenue from the production factor g allocated to department c;
According to the formula:
determining the elasticity of the price of the production factors to the technologies of different departments; in the method, in the process of the invention, The elasticity of the price w f representing the production factor f to the production technology a k of division k, lambda k representing the polymma weight of the yield of division k; /(I)The elasticity of the multi-MAT weight f representing the production factor f to the production technique A k of department k;
According to the formula:
determining the elasticity of the yield of different producers to different department technologies; in the method, in the process of the invention, Representing the elasticity of the production y i of sector i to the production technique a k of sector k; /(I)The elasticity of the multimma weight λ i representing the yield of division i to the production technique a k of division k; phi ik represents the elements of the ith row and k column in the freon-tigv inverse matrix psi of the HA-IO matrix omega, phi ig represents the elements of the ith row and g column in the freon-tigv inverse matrix psi of the HA-IO matrix omega; /(I)The flexibility of the production technique A k to the department k, which represents the price of the production factor g; lambda k represents the polymma weight of the yield of division k;
k. l and g are N, and k, l and g are positive integers;
A variation determining unit configured to calculate a corresponding department carbon emission variation when the full-factor productivity varies, based on the elasticity of the different department technologies for the different department yields; wherein, according to the formula:
Determining the carbon emission intensity e k of the department k; wherein E k is the carbon emission amount of the department k, and S k represents the total yield of the department k;
an emission intensity variation unit configured to determine carbon emission intensity variation amounts of respective departments when the full-factor productivity varies, based on the carbon emission variation amounts of the respective departments; wherein, according to the formula:
Determining the elasticity of the carbon emissions of sector i to sector k production technique a k; wherein E i is the carbon emissions of sector i; Representing the elasticity of the production y i of sector i to the production technique a k of sector k; e k denotes the carbon emission intensity of division k;
According to the formula:
Determining a carbon emission variation H i,k of the department i due to the full factor productivity variation of the department k; in the method, in the process of the invention, For the diagonal matrix of E i, E i represents the carbon emissions of division i,/>For diagonalization of Δa k/Ak, Δa k/Ak represents the change in full element productivity for division k; a k represents a production technology of department k;
A key section determining unit configured to determine a carbon emission key section based on the carbon emission intensity variation of each section when the full-factor productivity varies; wherein, according to the formula:
Determining a sales change matrix T i,k for division i due to a change in the full-element productivity of division k; in the method, in the process of the invention, A diagonal matrix of S; s i represents the total yield of department i; a k represents a production technology of department k; Δa k/Ak represents the variation in the full element productivity of division k;
According to the formula:
Determining a change Δd i,k in the carbon emission intensity of division i due to a change in the full factor productivity of division k; wherein E i、yi represents the carbon emission amount of the department i and the total output amount of the department i, respectively;
H i,k represents the carbon emission variation of division i due to the full factor productivity variation of division k.
CN202410393923.1A 2024-04-02 2024-04-02 Carbon emission analysis method and system based on full-element productivity improvement Active CN117973707B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202410393923.1A CN117973707B (en) 2024-04-02 2024-04-02 Carbon emission analysis method and system based on full-element productivity improvement

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202410393923.1A CN117973707B (en) 2024-04-02 2024-04-02 Carbon emission analysis method and system based on full-element productivity improvement

Publications (2)

Publication Number Publication Date
CN117973707A CN117973707A (en) 2024-05-03
CN117973707B true CN117973707B (en) 2024-06-04

Family

ID=90856672

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202410393923.1A Active CN117973707B (en) 2024-04-02 2024-04-02 Carbon emission analysis method and system based on full-element productivity improvement

Country Status (1)

Country Link
CN (1) CN117973707B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20140101469A (en) * 2013-02-08 2014-08-20 주식회사 휴디콤 The system and method for the carbon emission aim management and carbon emission reduction monitoring
CN109615254A (en) * 2018-12-17 2019-04-12 清华-伯克利深圳学院筹备办公室 Measuring method, device, management equipment and the storage medium of carbon emission performance
CN115358561A (en) * 2022-08-16 2022-11-18 北京师范大学 Multi-region input-output evaluation method based on multi-scenario analysis cause analysis
CN115511147A (en) * 2022-07-05 2022-12-23 国网上海市电力公司 Regional carbon emission prediction method
CN116485068A (en) * 2023-03-17 2023-07-25 北京师范大学 Input-output analysis method based on multivariate statistics

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20140101469A (en) * 2013-02-08 2014-08-20 주식회사 휴디콤 The system and method for the carbon emission aim management and carbon emission reduction monitoring
CN109615254A (en) * 2018-12-17 2019-04-12 清华-伯克利深圳学院筹备办公室 Measuring method, device, management equipment and the storage medium of carbon emission performance
CN115511147A (en) * 2022-07-05 2022-12-23 国网上海市电力公司 Regional carbon emission prediction method
CN115358561A (en) * 2022-08-16 2022-11-18 北京师范大学 Multi-region input-output evaluation method based on multi-scenario analysis cause analysis
CN116485068A (en) * 2023-03-17 2023-07-25 北京师范大学 Input-output analysis method based on multivariate statistics

Also Published As

Publication number Publication date
CN117973707A (en) 2024-05-03

Similar Documents

Publication Publication Date Title
Fei et al. How land transfer affects agricultural land use efficiency: Evidence from China’s agricultural sector
Amaza et al. Identification of factors that influence technical efficiency in rice-based production systems in Nigeria
CN106845767A (en) A kind of industry development in science and technology power quantitative estimation method and assessment system
He et al. Multi-dimensional diagnosis model for the sustainable development of regions facing water scarcity problem: A case study for Guangdong, China
Kancs et al. Factor content of bilateral trade: The role of firm heterogeneity and transaction costs
Iliyasu et al. Productivity growth, technical change and efficiency change of the Malaysian cage fish farming: an application of Malmquist Productivity Index approach
CN116882562A (en) Agricultural carbon emission influence factor analysis method and system
Bai et al. Agricultural trade impacts global phosphorus use and partial productivity
CN117973707B (en) Carbon emission analysis method and system based on full-element productivity improvement
Wang et al. The impact of environmental regulation on water resources utilization efficiency
Xin et al. Small farmer’s planting confidence and willingness to pay for leguminous green fertilizer: environmental attributes perspective
Umanath et al. Data envelopment analysis to estimate technical and scale efficiency of farms in Periyar-Vaigai irrigation system of Tamil Nadu, India
Li et al. [Retracted] Analysis Method of Agricultural Total Factor Productivity Based on Stochastic Block Model (SBM) and Machine Learning
CN111768282B (en) Data analysis method, device, equipment and storage medium
Zhao et al. The effect of urbanization on agricultural eco-efficiency and mediation analysis
CN111160635A (en) Regional logistics demand influence factor prediction method based on Relieff algorithm
Deng et al. Identification of key sectors of water resource utilization in China from the perspective of water footprint
Qian et al. RETRACTED ARTICLE: An empirical study on the impact of human capital on cost efficiency: taking the main corn production areas in China as an example
Li et al. Is Conservation Agriculture Adoption a Blessing or a Curse for Famers’ Agricultural Income?
CN112488438A (en) Industrial enterprise comprehensive benefit evaluation method and system combining credit information
Wei et al. Adoption mode of agricultural machinery and food productivity: evidence from China
Wan et al. An estimation framework for economic cost of land use based on artificial neural networks and principal component analysis with R
Dholakia et al. Issues in strategy for export promotion: An inter-industry analysis
CN116757074A (en) Intelligent optimization method for grazing strategies for different soil chemical properties
Qi et al. Trends and driving forces of agricultural carbon emissions: A case study of Anhui, China

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant