CN114881805B - Forestry carbon sink planning method based on carbon sink density - Google Patents
Forestry carbon sink planning method based on carbon sink density Download PDFInfo
- Publication number
- CN114881805B CN114881805B CN202210419419.5A CN202210419419A CN114881805B CN 114881805 B CN114881805 B CN 114881805B CN 202210419419 A CN202210419419 A CN 202210419419A CN 114881805 B CN114881805 B CN 114881805B
- Authority
- CN
- China
- Prior art keywords
- tree
- planting
- carbon sink
- carbon
- density
- 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
Links
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 142
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 141
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000012216 screening Methods 0.000 claims abstract description 9
- 238000005259 measurement Methods 0.000 claims abstract description 6
- 230000002159 abnormal effect Effects 0.000 claims abstract description 4
- 241000894007 species Species 0.000 claims description 78
- 239000002028 Biomass Substances 0.000 claims description 27
- 238000004364 calculation method Methods 0.000 claims description 26
- 230000009919 sequestration Effects 0.000 claims description 13
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 10
- 210000000481 breast Anatomy 0.000 claims description 9
- 241000196324 Embryophyta Species 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 7
- 235000007164 Oryza sativa Nutrition 0.000 claims description 6
- 235000009566 rice Nutrition 0.000 claims description 6
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 5
- 239000001569 carbon dioxide Substances 0.000 claims description 5
- 230000009467 reduction Effects 0.000 claims description 5
- 238000001556 precipitation Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 2
- 238000010276 construction Methods 0.000 claims description 2
- 238000005070 sampling Methods 0.000 claims description 2
- 241000209094 Oryza Species 0.000 claims 3
- 239000000126 substance Substances 0.000 claims 1
- 238000011160 research Methods 0.000 abstract description 3
- 238000007405 data analysis Methods 0.000 abstract description 2
- 238000011835 investigation Methods 0.000 abstract description 2
- 238000011158 quantitative evaluation Methods 0.000 abstract description 2
- 238000007781 pre-processing Methods 0.000 abstract 1
- 241000018650 Pinus massoniana Species 0.000 description 13
- 235000011609 Pinus massoniana Nutrition 0.000 description 13
- 241000219492 Quercus Species 0.000 description 13
- 235000016976 Quercus macrolepis Nutrition 0.000 description 12
- 240000007594 Oryza sativa Species 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 238000007726 management method Methods 0.000 description 3
- 244000050510 Cunninghamia lanceolata Species 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004134 energy conservation Methods 0.000 description 2
- 239000002023 wood Substances 0.000 description 2
- 235000005638 Austrian pine Nutrition 0.000 description 1
- 235000018185 Betula X alpestris Nutrition 0.000 description 1
- 235000018212 Betula X uliginosa Nutrition 0.000 description 1
- 241000218645 Cedrus Species 0.000 description 1
- 244000043261 Hevea brasiliensis Species 0.000 description 1
- 241000218657 Picea Species 0.000 description 1
- 235000008565 Pinus banksiana Nutrition 0.000 description 1
- 244000019397 Pinus jeffreyi Species 0.000 description 1
- 235000013264 Pinus jeffreyi Nutrition 0.000 description 1
- 235000008578 Pinus strobus Nutrition 0.000 description 1
- 235000008585 Pinus thunbergii Nutrition 0.000 description 1
- 235000014030 Podocarpus spicatus Nutrition 0.000 description 1
- 241001149649 Taxus wallichiana var. chinensis Species 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 235000017985 rocky mountain lodgepole pine Nutrition 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/80—Management or planning
- Y02P90/84—Greenhouse gas [GHG] management systems
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Quality & Reliability (AREA)
- Animal Husbandry (AREA)
- Operations Research (AREA)
- Educational Administration (AREA)
- Game Theory and Decision Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Development Economics (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to a carbon sink density-based forestry carbon sink planning method, and belongs to the technical field of forestry carbon sink planning. Firstly, collecting basic information in a planting area, and finishing screening of alternative tree species sets and field measurement of sample trees by combining the collected information; secondly, preprocessing the measured data by adopting an algorithm, and rejecting abnormal data; then, constructing a carbon sink density definition, and calculating the carbon sink density; and finally, performing forestry carbon sink planning based on the carbon sink density to ensure the maximum carbon sink amount of a target area within the operational life. The invention comprises the following steps: based on the national 'double carbon' strategy, the carbon sink density concept definition is constructed to realize the quantitative evaluation of the carbon sink capacity of different tree species within the target operational life; the forestry planning and the forestry carbon sink amount are innovatively combined, and the blank in the prior art is filled; based on theoretical research and actual investigation and by means of actually measured data analysis, accuracy and reliability of forestry carbon sink planning are ensured to the greatest extent.
Description
Technical Field
The invention relates to a carbon sink density-based forestry carbon sink planning method, and belongs to the technical field of forestry carbon sink planning.
Background
Carbon sequestration refers to the process of removing carbon dioxide from the air by artificial means or by natural means. The carbon sink is mainly divided into five types, namely forest carbon sink, grassland carbon sink, farmland carbon sink, ocean carbon sink and artificial carbon sink.
By combining the national 'double-carbon' strategy, the forest carbon sink has a larger development space, and in order to realize the goals of 'carbon peak reaching' and 'carbon neutralization', the tree planting strength is increased while the energy conservation and emission reduction are promoted and the unit GDP energy consumption is reduced, the forest coverage rate and the total forest carbon sink amount are improved through tree planting, and the development of the green economy of China is promoted together with the energy conservation and emission reduction.
At present, the main type of patents related to carbon sink is to calculate the carbon sink amount of forest regions by adopting a technical mode of biomass conversion factor calculation or leaf area calculation, for example, a patent represented by application publication No. CN201510178889.7, an applicant calculates the annual net carbon amount of unit leaf area by collecting and measuring related data, then substitutes the management and cultivation coefficient, the total leaf area amount and the average received light intensity value of unit leaf area into a formula to calculate the carbon sink value of an individual tree, and constructs a carbon sink calculation model of the individual tree in a city by combining crown appearance characteristics, light energy utilization rate and management and cultivation mode on the basis of a gas exchange method to measure the carbon sink capacity of the green space in the city.
The main type of the patents in the aspect of forestry planning is to pursue the maximization of forestry planting economy from the perspective of economy, and simultaneously, to overcome the defects of the existing planting mode, a planting planning or planting method is provided. For example, in a patent represented by application publication No. CN201910936592.0, a patent applicant performs steps such as seed selection, planting and mixed planting on rubber forests, firstly performs regional management on the rubber forests in a grouping mode, and then performs mixed plant and crop planting among rubber trees, so that the land utilization rate is improved, and the economic benefit of forest planting in the same area is improved.
In general, the existing main patents related to carbon sink and forestry planning only start from the perspective of single carbon sink calculation or forestry planning, and do not combine the two to provide a planting plan or method for realizing the forestry carbon sink optimization. Aiming at the defects of the two patents, the concept of carbon sink density is creatively constructed, and the problems of tree species selection, planning and planting and the like are systematically solved under the condition that the maximum carbon collection amount in the target business year is ensured.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: the invention provides a carbon sink density-based forestry carbon sink planning method according to the defects and shortcomings of the prior art and method, and solves the problem that the carbon sink amount and the forestry planning are separated from each other in the current forestry related planning, so that a forestry planning method with optimal carbon sink is lacked.
The technical scheme of the invention is as follows: a forestry carbon sink planning method based on carbon sink density comprises the following steps:
step1: a basic information acquisition link: collecting regional information, climate information, geographic information and internet tree information;
step2: tree species screening and sample tree measurement links: collecting field tree species, screening alternative tree species sets, measuring alternative tree species set sample trees;
step3: and (3) a data processing link: performing data exception processing on the sample tree measurement data acquired at Step2, eliminating abnormal data, and constructing a data set;
step4: and (3) a carbon sink density calculation link: calculating biomass, carbon sink amount, tree planting area and carbon sink density based on the data set constructed by Step 3;
step5: a planting planning link: planning the planting tree species and the planting period thereof under the maximum carbon sink amount target according to the Step4 calculation result and the operation period of the target area, and correspondingly calculating the planting intervals and the planting number of different tree species.
Specifically, the specific implementation steps of Step1 are as follows:
step1.1: collecting regional information: collecting the planting area S of a tree planting area, unit: square meter; width, unit: rice; length, unit: rice; operational age t plan The unit: year;
step1.2: climate information acquisition: obtain planting regional climate information from the internet platform, include: annual average temperature, unit: c, centigrade degree; annual average precipitation, unit: millimeter; annual average sunshine time, unit: h; a climate type;
step1.3: geographic information acquisition: using a positioning device, collecting geographical information of a planting area, the collected information of the planting area comprising: longitude, latitude, altitude, units of altitude: rice;
step1.4: collecting information of the Internet tree species; according to climate and geographic information of planting areas collected in Step1.2 and Step1.3, the tree species suitable for planting under the conditions of the same latitude, the same altitude and the same climate type are collected based on the Internet or an expert judgment method, and the tree species is set as a tree species set K1.
Specifically, the specific implementation steps of Step2 are as follows:
step2.1: collecting field tree species: collecting tree species information near a target planting area in the field by adopting an automatic or manual mode, and setting the tree species information as a tree species set K2;
step2.2: screening an alternative tree species set: combining the tree species sets K1 and K2 collected in Step1.4 and Step2.1, screening out alternative tree species sets K3 and G j Represents different tree species;
K3＝K1∩K2＝{G 1 ,G 2 ,...,G j } (1)
step2.3: measuring a sample tree: sampling the same number m of samples of different tree species and trees corresponding to the tree species with different ages on the spot near a planting area according to the alternative tree species set K3, and measuring detailed information of the samples, wherein the detailed information specifically comprises the following steps: tree species G j Age Ag, height H, diameter D, and crown width W.
Specifically, the specific implementation steps of Step3 are as follows:
step3.1: data clustering: processing and analyzing the tree information measured in field in Step2.3 according to tree species by adopting a clustering analysis method and using a K-means or DBSCAN algorithm, removing abnormal data, and repeatedly measuring and processing the data until the number of the trees of the sample with each age corresponding to each required tree species reaches m; step3.2: and (3) data set construction: constructing the data subjected to noise reduction processing into a data set V1;
v1= { tree number i, tree species Gj, tree age Ag i Height of tree H i Diameter at breast height D i Crown width W i }
Specifically, the specific implementation steps of Step4 are as follows:
step4.1: biomass calculation: calculating biomass M of tree i by using tree height and breast diameter of V1 data set in Step3.2 i ；
M i ＝a×(D i p H i q ) (2)
Wherein a, p and q are regression constants, and can be inquired according to ' the forestry industry standard of the people's republic of China-standing tree biomass model and carbon measurement parameter ' issued by the national forestry bureau according to tree species and breast-height, and the biomass unit is kilogram;
step4.2: and (3) calculating the carbon sink amount: calculating the carbon sink C of the tree i by using the biomass variable quantity and the carbon-containing coefficient i ；
ΔM i ＝M Ag -M Ag ′ (3)
plant＝Ag-Ag′ (4)
Wherein, Δ M i Representing the amount of change in biomass over a period of time, M, of the tree i Ag Represents the current biomass of the tree, M Ag ' represents biomass when the tree is transplanted, and plant represents the planting age of the tree; ag represents the current age of the tree; ag' represents the tree age of the tree in transplanting and planting, and the unit of carbon sink is kilogram;
C i ＝ΔM i ×Tc i (5)
wherein Tc i The carbon-containing coefficient of the tree i can be inquired according to the stand wood biomass model and the carbon metering parameter which are the forestry industry standards of the people's republic of China and are released by the national forestry bureau of China;
step4.3: calculating the planting area of the tree i: calculating the tree crown width of the V1 data set in Step3.2, regarding the tree crown width as a circle, wherein the tree crown width is the diameter of the circle, making a circumscribed square of the circle, the diameter of the circle is equal to the side length of the square, and the area of the square is the tree planting area of the tree in the corresponding planting year;
S i ＝W i 2 (6)
wherein S is i The planting area of the tree corresponding to the planting age limit is expressed in the unit of squareRice, W i Is the crown of tree i;
step4.4: carbon sink density ρ calculation: the carbon sink density is defined as the amount of carbon dioxide absorbed by trees planted in a unit land area in a unit time and is used for measuring the condition that the trees in the area absorb the carbon dioxide within a specific time and range;
wherein the unit of the carbon sink density is kg/(m) 2 ·y)；
Specifically, as shown in formula (7), the carbon sequestration density is determined by dividing the carbon sequestration amount by the planting area by the planting age, and the carbon sequestration density is determined by the unit time and the unit area;
step4.5: according to the data set V1, calculating the carbon sink density of all sample trees of different tree species under different planting years by using the breast height, the tree height and the crown width combined formula (7), and then respectively calculating the corresponding average carbon sink density;
wherein ρ b Representing the carbon sink density of the sample tree;
based on the existing calculation data, a new data set is constructed:
step4.6: and (3) calculating the total carbon sink density f (rho) under the target operational age by combining the data set V2 in Step4.5, wherein f (rho) represents the total value of the carbon sink density of the trees planted within the target operational age, and solving under the condition of the maximum total carbon sink density so that:
wherein, t plan The time limit of the operation is shown,
denotes the G th j Planting age limit of seed and tree species>
Denotes the G th j Average carbon sink density corresponding to planting age of the seed tree.
Calculating planning planting years corresponding to different tree species according to a formula (9), and obtaining a planting planning analysis set V3 of the forestry within a target operation year:
specifically, as shown in formula (9), a function of the average carbon sink density of the trees about the planned planting age of the trees is constructed, the total carbon sink density of the tree species is calculated by multiplying the average carbon sink density of the trees by the planting age of the trees, and different tree species are summed up respectively, so that the total carbon sink density of the forest land is maximized within the operating age, and the maximum total carbon sink amount is further ensured.
Specifically, the specific implementation steps of Step5 are as follows:
step5.1: and (3) planting interval calculation: correspondingly carrying out planting planning according to the analysis set V3 of Step4.6 based on G j The planting age limit of the tree species
The average crown of all samples is calculated to obtain G j Corresponding planting interval of a tree species>
Wherein, W b Representing the crown width of the sample tree;
step5.2: planting planning: according to the length and width of the target planting area collected in Step1.1And planting intervals in Step5.1
Calculate G j Corresponding number of plantable trees in a tree species>
Wherein length is the length of the planting area, width is the width of the planting area,
is G j <xnotran> , [ </xnotran>]Representing rounding;
step5.3: tree planting data set V4 constructed based on Step5.1 and 5.2
The tree detailed planting plan is as follows:
the invention has the following effective effects:
firstly, based on the industry of 'double carbon', the carbon sink density concept is constructed, the carbon sink capacity of different tree species within a target year is quantitatively evaluated, and the tree planting is reasonably planned according to the carbon sink capacity; secondly, the forestry planning and the forestry carbon sink amount are innovatively combined, and the blank in the prior art is filled; thirdly, the accuracy of the carbon sink amount calculation result is ensured to the maximum extent by means of the analysis of the measured data and various data.
Drawings
FIG. 1 is a general flow chart of the present invention;
FIG. 2 is a flow chart of carbon sink density calculation;
FIG. 3 is a schematic diagram of the carbon sink density of different tree species;
FIG. 4 is a Quercus tree algorithm cluster map in an embodiment;
FIG. 5 is a cluster diagram of the fir algorithm in the example;
fig. 6 is a masson pine algorithm cluster map in an embodiment.
Detailed Description
The invention is further described with reference to the following figures and specific examples.
Example 1: as shown in fig. 1-6, a forestry carbon sink planning method based on carbon sink density includes the following steps:
step1.1: a certain area with the area size of 100000 square meters (about 150 mu) needs to be planted with a carbon sink forest, the length of the area is 500 meters, the width of the area is 200 meters, and the operating life of the forest land is 10 years.
Step1.2: the local environment has the annual average temperature of 16 ℃, the annual average precipitation of 1535.6 mm and the sunshine duration of 1586.6 hours, and the region belongs to subtropical monsoon climate.
Step1.3: the latitude and longitude ranges of the planting areas are collected through the positioning equipment, and are 111-20 '5' to 111-54 '39' east longitude, 23-58 '33' to 24-14 '25' north latitude, and the average altitude is 200 meters.
Step1.4: the climate of the area is similar from 23 degrees north latitude to 25 degrees north latitude, and the suitable tree species for planting are collected according to the climate, geography and other information: k1= { cedar, oak, pinus massoniana, birch, taxus chinensis, black pine }.
Step2.1: collecting tree species information near a planting area, wherein common tree species in a range of 10km near the planting area specifically comprise K2= { oak, fir, masson pine and spruce }.
Step2.2: and selecting an alternative tree species set according to the tree species obtained in Step1.3 and Step2.1 through Internet and local actual research, wherein K3= { oak, fir and masson pine }.
Step2.3: selecting the tree age of the transplanted tree to be 5 years according to the operating age limit of the forest land, so that 10 sample trees are extracted from the fir, oak and Chinese red pine with the tree age of 5-15 years in the surrounding area respectively, and the data information of the tree species, the tree age, the breast height, the crown width and the tree height of the tree is measured and recorded.
Step3.1: adopting a DBSCAN algorithm to perform noise reduction treatment on the sample tree data collected in Step2.3 in a tree species division manner, wherein the DBSCAN parameter setting comprises the following steps: radius is 1.5 and the minimum neighborhood point number is 3. As shown in particular in fig. 4-6. And if the noise points exist in the investigated data, cleaning the data, removing the noise points, and performing supplementary measurement until the number m of the trees with different ages of various tree species is 10.
Step3.2: statistical data of tree heights, crown widths and breast diameters of oak, fir and masson pine in different ages are obtained, and the statistical data are respectively shown in tables 2, 3 and 4.
Step4.1: according to the 'PRC forestry industry standard-standing tree biomass model and carbon metering parameters', it is respectively known that the regression constants a, p and q need to adopt different values for trees with the breast diameter more than or equal to 5 cm or less than 5 cm, and the regression constants a, p and q and carbon-containing coefficients in the biomass formulas of oak, fir and masson pine are inquired as follows:
table 1: biomass conversion factor constant and carbon content coefficient corresponding to oak, fir and Chinese red pine
Biomass calculation: the biomass of three trees was calculated by substituting the values of the corresponding regression constants in table 1, and the calculation results are shown in tables 2, 3 and 4:
(1) oak tree:
(2) and (3) fir:
(3) masson pine:
step4.2: carbon sink amount calculation formula: calculating the carbon sink amount (unit: kilogram) according to the biomass;
(1) oak tree: c i ＝0.4802×ΔM i
(2) And (3) fir: c i ＝0.499×ΔM i
(3) Masson pine: c i ＝0.5252×ΔM i
Calculating the planting area of the Step4.3 trees: the tree planting area is calculated by multiplying the crown width by the crown width, and the calculation results are respectively shown in tables 2, 3 and 4:
S i ＝W i 2
table 2: calculation data table of sample trees of oak trees with different ages
Table 3: calculation data table of sample trees of China fir with different ages
Table 4: calculation data table of sample trees of pinus massoniana in different ages
Step4.4: carbon sink density calculation:
(1) oak tree:
(2) and (3) fir:
(3) masson pine:
step4.5: according to the carbon sink densities of different tree ages of different tree species of the sample tree, solving the average carbon sink densities of the sample tree at different tree ages as follows:
table 5: average carbon sink density table (unit: kg/(m) for three species samples with different ages and planting years 2 ·y))
Step4.6: the data of table 5 were substituted into the following calculation:
Step5.1: the planting interval (unit: meter) at the time of planting was determined based on the average crown width at the time of planting 10 years of the fir tree of 15 years of age in table 3.
Step5.2: according to the planting range of 500 × 200, the number of trees which can be planted is calculated to be about:
| tree species | Year limit of planting | Planting space | Number of plants planted |
| Chinese fir wood | For 10 years | 4.2 m | 5593 a plant |
Compared with the prior art, the invention has the following advantages: firstly, based on the national 'double carbon' strategy, the quantitative evaluation of the carbon sink capacity of different tree species within the target operational life is realized by constructing the concept definition of carbon sink density; secondly, the forestry planning and the forestry carbon sink amount are innovatively combined, and the blank in the prior art is filled; thirdly, based on theoretical research and actual investigation and based on actually measured data analysis, accuracy and reliability of forestry carbon sink planning are ensured to the greatest extent; fourthly, the method provided by the text is simple, practical and efficient, and has high application and popularization values.
The above examples are merely illustrative of embodiments of the present invention, which are described in more detail and detail, and should not be construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the inventive concept, which falls within the scope of the present invention.
Claims (9)
1. A forestry carbon sink planning method based on carbon sink density is characterized in that: the method comprises the following steps:
step1: a basic information acquisition link: collecting regional information, climate information, geographic information and Internet tree information;
step2: the tree species screening and sample tree measuring link: collecting field tree species, screening alternative tree species sets, measuring alternative tree species set sample trees;
step3: and a data processing link: performing data exception processing on the sample tree measurement data acquired at Step2, eliminating exception data, and constructing a data set V1;
step4: a carbon sink density calculation link: calculating biomass, carbon sink amount, tree planting area and carbon sink density based on the data set constructed by Step 3;
step5: a planting planning link: planning the planting tree species and the planting age thereof under the condition of the maximum total carbon sink density according to the Step4 calculation result and the operation age of the target area, and correspondingly calculating the planting intervals and the planting quantity of different tree species;
the specific steps of Step4 are as follows:
step4.1: biomass calculation: calculating biomass M of tree i by using the height and the breast diameter of the V1 data set in Step3 i ；
M i ＝a×(D i p H i q ) (2)
Wherein a, p and q are regression constants, and biomass unit is kilogram;
step4.2: and (3) calculating the carbon sequestration amount: calculating the carbon sink C of the tree i by using the biomass variable quantity and the carbon-containing coefficient i ；
ΔM i ＝M Ag -M Ag ′ (3)
plant＝Ag-Ag′ (4)
Wherein, Δ M i Representing the amount of change in biomass over a period of time, M, of the tree i Ag Representation treeWood current biomass, M Ag ' represents biomass when the tree is transplanted, and plant represents planting age of the tree; ag represents the current age of the tree; ag' represents the tree age of the tree in transplanting and planting, and the unit of carbon sink is kilogram;
C i ＝ΔM i ×Tc i (5)
wherein Tc is i Is the carbon content coefficient of tree i;
step4.3: and (3) calculating the planting area of the tree i: calculating the tree crown width of the V1 data set in Step3, regarding the tree crown width as a circle, wherein the tree crown width is the diameter of the circle, and making an externally tangent square of the circle, wherein the diameter of the circle is equal to the side length of the square, and the area of the square is the planting area of a single tree in the corresponding planting year limit;
S i ＝W i 2 (6)
wherein S is i The planting area of the tree corresponding to the planting age limit is expressed in square meter, W i Is the crown of tree i;
step4.4: carbon sink density ρ calculation: the carbon sink density is defined as the amount of carbon dioxide absorbed by trees planted in a unit land area in a unit time and is used for measuring the condition that the trees in the area absorb the carbon dioxide within a specific time and range;
wherein the unit of the carbon sink density is kg/(m) 2 ·y)；
Step4.5: according to the data set V1, calculating the carbon sink density of all sample trees of different tree species under different planting years by using a breast height, tree height and crown width combined formula (7), and then respectively calculating corresponding average carbon sink densities;
where ρ is b Representing the carbon sink density of the sample tree;
constructing a new data set V2:
step4.6: and (3) calculating the total carbon sink density f (rho) under the target operational age by combining the data set V2 in Step4.5, wherein f (rho) represents the total value of the carbon sink density of the trees planted within the target operational age, and solving under the condition of the maximum total carbon sink density so that:
wherein the content of the first and second substances, tplan the time limit of the operation is shown,
denotes the G th j The planting age and the length of the plant>
Denotes the G th j Average carbon sink density corresponding to planting age of the seed tree;
calculating planning planting years corresponding to different tree species according to a formula (9), and obtaining a planting planning analysis set V3 of the forestry within a target operation year:
2. a method for planning a carbon sink in forestry based on carbon sink density as claimed in claim 1, wherein: the specific steps of Step1 are as follows:
step1.1: collecting regional information: collecting the planting area S of a tree planting area, unit: square meter; width, unit: rice; the length of the length,unit: rice; operational age t plan The unit: year;
step1.2: climate information acquisition: obtain planting regional climate information from the internet platform, include: annual average temperature, unit: c, centigrade degree; annual average precipitation, unit: millimeter; annual average sunshine time, unit: hours; a climate type;
step1.3: geographic information acquisition: using a positioning device, collecting geographical information of a planting area, the collected information of the planting area comprising: longitude, latitude, altitude, units of altitude: rice;
step1.4: collecting information of the Internet tree species; according to climate and geographic information of planting areas collected in Step1.2 and Step1.3, the tree species suitable for planting under the conditions of the same latitude, the same altitude and the same climate type are collected based on the Internet or an expert judgment method, and the tree species is set as a tree species set K1.
3. A method for forestry carbon sequestration planning based on carbon sequestration density, as claimed in claim 2, wherein: the specific steps of Step2 are as follows:
step2.1: collecting field tree species: collecting tree species information near a target planting area in a field by adopting an automatic or manual mode, and setting the tree species information as a tree species set K2;
step2.2: screening an alternative tree species set: combining the tree species sets K1 and K2 collected in Step1.4 and Step2.1, screening out alternative tree species sets K3 and G j Represents different tree species:
K3＝K1∩K2＝{G 1 ,G 2 ,...,G j } (1)
step2.3: measuring a sample tree: sampling the same number m of samples of different tree species and trees corresponding to the tree species with different ages on the spot near a planting area according to the alternative tree species set K3, and measuring detailed information of the samples, wherein the detailed information specifically comprises the following steps: tree species G j Age Ag, height H, diameter D, and crown width W.
4. A method for forestry carbon sequestration planning based on carbon sequestration density, as claimed in claim 3, wherein: the concrete steps of Step3 are as follows:
step3.1: data clustering: processing and analyzing the tree information measured in field in Step2.3 according to tree species by adopting a clustering analysis method and using a K-means or DBSCAN algorithm, removing abnormal data, and repeatedly measuring and processing the data until the number of the trees of the sample with each age corresponding to each required tree species reaches m;
step3.2: and (3) data set construction: constructing the data subjected to noise reduction processing into a data set V1:
v1= { tree number i, tree species G j Age of tree Ag i Height of tree H i Diameter at breast height D i Crown width W i }。
5. A method for forestry carbon sink planning based on carbon sink density as claimed in claim 4, wherein: the specific steps of Step5 are as follows:
step5.1: and (3) calculating planting intervals: according to the analysis set V3 of Step4.6, correspondingly carrying out planting planning based on G j The planting age limit of the tree species
Calculating the average crown width of all the sample trees j Corresponding planting interval of tree species>
Wherein, W b Representing the crown width of the sample tree;
step5.2: planting planning: according to the length and width of the target planting area collected in Step1.1 and the planting interval in Step5.1
Calculate G j Corresponding number of plantable trees of a tree species->
Wherein, length is the length of the planting area, width is the width of the planting area,
is G j The planting interval of the tree species is set, 2 [ 2 ]]Representing rounding;
step5.3: based on the calculation results of Step5.1 and Step5.2, correspondingly constructing a tree planting data set V4:
6. a method for forestry carbon sequestration planning based on carbon sequestration density, as claimed in claim 2, wherein: in Step1.1: the planting area in the area information acquisition refers to the area obtained directly or calculated through a defined planting range, the planting area is in any shape, and the planting area is converted into a rectangle through a pattern cutting mode.
7. A method for forestry carbon sequestration planning based on carbon sequestration density, as claimed in claim 2, wherein: in Step1.1: the operational age refers to the land use age or the maximum planting time, and in Step1.4: the method for collecting the trees suitable for planting in the same latitude, altitude and climate types by using the Internet refers to that a latitude range is defined according to the climate type of a local region and the climate similarity condition of the adjacent region of the latitude.
8. A method for planning a carbon sink in forestry based on carbon sink density as claimed in claim 3, wherein: in Step2.1: the range near the planting area is collected, the specific radius is not limited, and the radius range is defined according to the difference situation of the surrounding climate and the altitude.
9. A method for forestry carbon sink planning based on carbon sink density as claimed in claim 1, wherein: in Step4.1: biomass, referring to aboveground biomass of the tree, in step4.2: and (4) transplanting, namely selecting different tree age data of the same tree species in the area to carry out simulated transplanting, or transplanting the trees from the original growing area to a target planting area to carry out planting.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210419419.5A CN114881805B (en) | 2022-04-21 | 2022-04-21 | Forestry carbon sink planning method based on carbon sink density |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210419419.5A CN114881805B (en) | 2022-04-21 | 2022-04-21 | Forestry carbon sink planning method based on carbon sink density |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114881805A CN114881805A (en) | 2022-08-09 |
| CN114881805B true CN114881805B (en) | 2023-04-07 |
Family
ID=82671905
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210419419.5A Active CN114881805B (en) | 2022-04-21 | 2022-04-21 | Forestry carbon sink planning method based on carbon sink density |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114881805B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115796624B (en) * | 2023-02-13 | 2023-04-04 | 广东省木链网科技股份有限公司 | Environment-friendly tree planting management platform based on carbon neutralization target |
| CN117314020B (en) * | 2023-11-28 | 2024-02-27 | 生态环境部华南环境科学研究所(生态环境部生态环境应急研究所) | Wetland carbon sink data monitoring system of plankton |
| CN117455127B (en) * | 2023-12-26 | 2024-03-15 | 临沂市园林环卫保障服务中心 | Plant carbon sink dynamic data monitoring system based on wisdom gardens |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105335547A (en) * | 2015-09-18 | 2016-02-17 | 北京林业大学 | Optimal tree seed diameter grade density distribution method with biggest impact on carbon sink and biomass increment by stand density |
| CN106022652A (en) * | 2016-06-12 | 2016-10-12 | 刘兆刚 | Processing method of forest carbon sink operating plan and processing device of forest carbon sink operating plan |
| CN113705015A (en) * | 2021-09-03 | 2021-11-26 | 知晓（北京）通信科技有限公司 | Forestry carbon remittance method for forest land multi-data fusion |
-
2022
- 2022-04-21 CN CN202210419419.5A patent/CN114881805B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN114881805A (en) | 2022-08-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114881805B (en) | Forestry carbon sink planning method based on carbon sink density | |
| Gilmanov et al. | Productivity, respiration, and light-response parameters of world grassland and agroecosystems derived from flux-tower measurements | |
| CN110376167A (en) | Rice leaf nitrogen content monitoring method based on unmanned plane EO-1 hyperion | |
| Polglase et al. | Nitrogen and phosphorus cycling in relation to stand age of Ecucalyptus regnans F. Muell: I. Return from plant to soil in litterfall | |
| Hudiburg et al. | Evaluation and improvement of the Community Land Model (CLM4) in Oregon forests | |
| CN113268923B (en) | Summer corn yield estimation method based on simulated multispectral | |
| CN113268703A (en) | Nitrogen fertilizer deficiency rapid detection and precision topdressing method applied to wheat field management | |
| Shvidenko et al. | Carbon budget of the Russian boreal forests: a systems analysis approach to uncertainty | |
| Eastham et al. | Planting density effects on water use efficiency of trees and pasture in an agroforestry experiment | |
| CN111160799A (en) | Natural resource database construction method | |
| CN104820065B (en) | A kind of carbon remittance measuring method of city individual plant arbor | |
| CN111781200A (en) | Research method for carbon reserves of forest soil | |
| Ranasinghe et al. | Estimation of carbon sequestration potential in coconut plantations under different agro-ecological regions and land suitability classes | |
| CN115952702A (en) | Forest NEP calculation method based on FORCCHN model and remote sensing data | |
| Yang et al. | Modeling canopy photosynthesis and light interception partitioning among shoots in bi-axis and single-axis apple trees (Malus domestica Borkh.) | |
| CN113269450B (en) | County scale crop resource utilization and environmental effect evaluation method and system | |
| Käller | Vegetation-environment interactions in a boreonemoral forest in east central Sweden | |
| Sharpe | Methods of assessing the primary production of regions | |
| CN108154320B (en) | Clonal fir leaf surface photosynthetic benefit approximation calculation method | |
| CN112034111B (en) | Method for identifying high yield and low emission of rice varieties | |
| CN114792274A (en) | Layout method for artificial vegetation construction in arid region and mountain land | |
| Pate et al. | Uptake, partitioning and utilization of carbon and nitrogen in the phloem bleeding tree, Tasmanian blue gum (Eucalyptus globulus) | |
| CN109118382B (en) | Method for establishing relation model of soil moisture content and fertilizing amount and application | |
| CN109601219B (en) | Method for comprehensively evaluating drought resistance of Minjiang juniper in seedling stage | |
| Nonhebel | Harvesting the sun's energy using agro-ecosystems. |
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 |