CN114969665A - Method and system for estimating carbon emission and carbon reserve of mineral resource base - Google Patents
Method and system for estimating carbon emission and carbon reserve of mineral resource base Download PDFInfo
- Publication number
- CN114969665A CN114969665A CN202210620129.7A CN202210620129A CN114969665A CN 114969665 A CN114969665 A CN 114969665A CN 202210620129 A CN202210620129 A CN 202210620129A CN 114969665 A CN114969665 A CN 114969665A
- Authority
- CN
- China
- Prior art keywords
- carbon
- area
- sample
- vegetation
- content
- 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.)
- Granted
Links
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 437
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 427
- 229910052500 inorganic mineral Inorganic materials 0.000 title claims abstract description 159
- 239000011707 mineral Substances 0.000 title claims abstract description 159
- 238000000034 method Methods 0.000 title claims abstract description 68
- 239000011435 rock Substances 0.000 claims abstract description 108
- 238000005265 energy consumption Methods 0.000 claims abstract description 42
- UQSXHKLRYXJYBZ-UHFFFAOYSA-N Iron oxide Chemical compound [Fe]=O UQSXHKLRYXJYBZ-UHFFFAOYSA-N 0.000 claims description 108
- 239000002689 soil Substances 0.000 claims description 105
- 239000000395 magnesium oxide Substances 0.000 claims description 47
- CPLXHLVBOLITMK-UHFFFAOYSA-N magnesium oxide Inorganic materials [Mg]=O CPLXHLVBOLITMK-UHFFFAOYSA-N 0.000 claims description 47
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 46
- AXZKOIWUVFPNLO-UHFFFAOYSA-N magnesium;oxygen(2-) Chemical compound [O-2].[Mg+2] AXZKOIWUVFPNLO-UHFFFAOYSA-N 0.000 claims description 45
- 239000000292 calcium oxide Substances 0.000 claims description 42
- ODINCKMPIJJUCX-UHFFFAOYSA-N calcium oxide Inorganic materials [Ca]=O ODINCKMPIJJUCX-UHFFFAOYSA-N 0.000 claims description 42
- 238000003860 storage Methods 0.000 claims description 40
- 238000005070 sampling Methods 0.000 claims description 37
- BRPQOXSCLDDYGP-UHFFFAOYSA-N calcium oxide Chemical compound [O-2].[Ca+2] BRPQOXSCLDDYGP-UHFFFAOYSA-N 0.000 claims description 28
- 238000012216 screening Methods 0.000 claims description 26
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 23
- 239000001569 carbon dioxide Substances 0.000 claims description 22
- KMUONIBRACKNSN-UHFFFAOYSA-N potassium dichromate Chemical compound [K+].[K+].[O-][Cr](=O)(=O)O[Cr]([O-])(=O)=O KMUONIBRACKNSN-UHFFFAOYSA-N 0.000 claims description 22
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 21
- 239000002028 Biomass Substances 0.000 claims description 20
- 235000019402 calcium peroxide Nutrition 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 12
- 238000010438 heat treatment Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 claims description 8
- LHJQIRIGXXHNLA-UHFFFAOYSA-N calcium peroxide Chemical compound [Ca+2].[O-][O-] LHJQIRIGXXHNLA-UHFFFAOYSA-N 0.000 claims description 8
- 230000010354 integration Effects 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 6
- SPAGIJMPHSUYSE-UHFFFAOYSA-N Magnesium peroxide Chemical compound [Mg+2].[O-][O-] SPAGIJMPHSUYSE-UHFFFAOYSA-N 0.000 claims description 4
- 238000005303 weighing Methods 0.000 claims description 4
- 238000011161 development Methods 0.000 abstract description 25
- 230000000694 effects Effects 0.000 abstract description 3
- 238000005065 mining Methods 0.000 description 37
- 238000005553 drilling Methods 0.000 description 8
- 239000000843 powder Substances 0.000 description 4
- 238000000926 separation method Methods 0.000 description 4
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 3
- FYYHWMGAXLPEAU-UHFFFAOYSA-N Magnesium Chemical compound [Mg] FYYHWMGAXLPEAU-UHFFFAOYSA-N 0.000 description 3
- 230000037396 body weight Effects 0.000 description 3
- 239000011575 calcium Substances 0.000 description 3
- 229910001576 calcium mineral Inorganic materials 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 229910001608 iron mineral Inorganic materials 0.000 description 3
- 239000011777 magnesium Substances 0.000 description 3
- 229910001607 magnesium mineral Inorganic materials 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000009919 sequestration Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000009841 combustion method Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000003306 harvesting Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 229910052751 metal Inorganic materials 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000004876 x-ray fluorescence Methods 0.000 description 2
- BPQQTUXANYXVAA-UHFFFAOYSA-N Orthosilicate Chemical compound [O-][Si]([O-])([O-])[O-] BPQQTUXANYXVAA-UHFFFAOYSA-N 0.000 description 1
- 239000003082 abrasive agent Substances 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000007605 air drying Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- PTKRHFQQMJPPJN-UHFFFAOYSA-N dipotassium;oxido-(oxido(dioxo)chromio)oxy-dioxochromium;sulfuric acid Chemical compound [K+].[K+].OS(O)(=O)=O.[O-][Cr](=O)(=O)O[Cr]([O-])(=O)=O PTKRHFQQMJPPJN-UHFFFAOYSA-N 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000008213 purified water Substances 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 238000005527 soil sampling Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/188—Vegetation
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computational Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Business, Economics & Management (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Bioinformatics & Computational Biology (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Economics (AREA)
- Mining & Mineral Resources (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Probability & Statistics with Applications (AREA)
- Marine Sciences & Fisheries (AREA)
- Evolutionary Biology (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Multimedia (AREA)
- Algebra (AREA)
Abstract
The invention discloses a method and a system for estimating carbon emission and carbon reserve of a mineral resource base, wherein the carbon emission and the carbon reserve of comprehensive energy consumption are obtained by calculating the ore amount corresponding to each type of ore in an area to be estimated; simultaneously calculating the area and the carbon reserve of each region in the region to be estimated to obtain the loss amount of the ecological carbon; obtaining mineral carbonated carbon reserves according to the average content of each mineral in the surface rock sample and the drill core sample; and obtaining the carbon emission and the carbon reserve of the area to be estimated based on the ecological carbon loss, the comprehensive energy consumption carbon emission and the mineral carbonation carbon reserve. According to the invention, the carbon emission and the carbon reserve of the mineral resource base are integrally estimated, so that the estimation efficiency of the carbon emission and the carbon reserve before the development of the mineral resource base is improved, the blank of the estimation system of the carbon emission and the carbon reserve before the development of the mineral industry is filled, and important green support is provided for the subsequent activities of environment protection, resource utilization, geological exploration, mine planning and the like.
Description
Technical Field
The invention relates to the technical field of estimation of carbon emission and carbon reserve of mineral resource bases, in particular to a method and a system for estimating the carbon emission and the carbon reserve of the mineral resource bases.
Background
Along with the development of the country, the demand of various metal resources is increasing day by day, the external dependence degree of a large number of metal minerals in China is high for a long time, and the development of a mineral resource base is promoted in a plateau area in order to solve the self-supply capacity of the mineral resources in China. The problem of carbon dioxide emission and absorption is beginning to be emphasized by various industries, especially energy systems are used as the largest carbon emission source, and in mineral development, estimation of carbon emission and carbon reserve should be emphasized.
At present, the international research on the carbon dioxide emission in mineral development mainly researches the carbon dioxide absorbed by the respiration of soil and vegetation in the occupied land range, the indirect carbon emission of energy consumed in the production process of enterprises and the like. The studies on carbon dioxide sequestration are mostly carried out by using abandoned mine tunnels, oil and gas field rock strata sequestration, ocean sequestration and mineral carbonation. However, when the related data of carbon dioxide in mineral development is estimated, the comprehensive quantitative research is less, and the carbon dioxide emission and the carbon reserve need to be acquired by splitting, so that the acquisition efficiency of the carbon dioxide emission and the carbon reserve is reduced.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: the method and the system for estimating the carbon emission and the carbon reserve of the mineral resource base improve the estimation efficiency of the carbon emission and the carbon reserve before the development of the mineral resources by integrating and estimating the carbon emission and the carbon reserve of the mineral resource base, fill the blank of a carbon dioxide emission and carbon reserve estimation system before the development of the mineral industry at present, and provide important green support for the follow-up activities of environment protection, resource utilization, geological exploration, mine planning and the like.
In order to solve the technical problem, the invention provides a method for estimating carbon emission and carbon reserve of a mineral resource base, which comprises the following steps:
constructing an ore body model of a region to be estimated, dividing the ore body model into a plurality of sub-modules, and calculating and obtaining the comprehensive energy consumption carbon emission amount based on the ore amount corresponding to each type of ore according to the volume and small weight parameters of each sub-module;
acquiring a normalized vegetation index map corresponding to the area to be estimated according to the remote sensing image map of the area to be estimated, carrying out area division on the area to be estimated according to the normalized vegetation index map, and calculating the area of each area;
acquiring a sampling sample of each region, measuring the carbon reserves of the sampling sample to obtain the carbon reserves corresponding to each region, and calculating the ecological carbon loss of the region to be estimated according to the area of the region and the carbon reserves;
collecting the surface rock sample and the drill core sample of the area to be estimated, and calculating and obtaining the mineral carbonation carbon reserves of the area to be estimated according to the average content of each mineral in the surface rock sample and the drill core sample;
and combining the ecological carbon loss amount and the comprehensive energy consumption carbon emission amount to obtain the carbon emission amount of the area to be estimated, and obtaining the carbon storage amount of the area to be estimated according to the mineral carbonation carbon storage amount.
In a possible implementation manner, the obtaining a sampling sample of each region, and performing carbon reserve measurement on the sampling sample to obtain a carbon reserve corresponding to each region specifically includes:
setting and sampling vegetation and soil in each area according to a sampling line corresponding to each area to obtain a vegetation sample and a soil sample corresponding to each area, wherein the vegetation sample comprises an overground vegetation sample and an underground vegetation sample;
according to a dry burning method, determining the carbon content of the vegetation sample, and calculating the carbon reserve of the vegetation in each area according to the determined carbon content of the vegetation sample;
and measuring the carbon content of the soil sample according to a potassium dichromate external heating method, and calculating the carbon reserve of the soil in each area according to the measured carbon content of the soil sample.
In a possible implementation manner, the area division is performed on the area to be estimated according to the normalized vegetation index map, and the area of each area is calculated, which specifically includes:
dividing the area to be estimated into a desertification and salinization area, a low vegetation coverage area, a medium vegetation coverage area and a high vegetation coverage area according to the normalized vegetation index map;
setting and counting the number of pixel units corresponding to each area and the precision value of the remote sensing image according to the normalized vegetation index threshold;
and calculating the area of each region according to the counted number of the pixel units and the precision value of the remote sensing image.
In a possible implementation manner, the calculating and obtaining the mineral carbon carbonate carbon reserves of the area to be estimated according to the average content of each mineral in the surface rock sample and the borehole core sample specifically includes:
performing main quantity element analysis on the obtained surface rock sample and the obtained drill core sample respectively to obtain mineral contents corresponding to each surface rock sample and each drill core sample, wherein the mineral contents comprise calcium oxide content, magnesium oxide content and ferrous oxide content;
obtaining the average content of each mineral in each surface rock sample and each drill core sample according to the mineral content corresponding to each surface rock sample and each drill core sample;
and respectively calculating the relative molecular weight ratio of the carbon dioxide to the calcium oxide content, the magnesium oxide content and the ferrous oxide content, and obtaining the mineral carbonation carbon reserves of the area to be estimated according to the relative molecular weight ratio and the average content of each mineral.
In a possible implementation manner, the method for measuring the carbon content of the soil sample according to the potassium dichromate external heating method and calculating the carbon reserve of the soil per unit area in each area according to the measured carbon content of the soil sample specifically comprises the following steps:
randomly selecting three soil samples from all soil samples, carrying out water content measurement on the three soil samples to obtain water content coefficients corresponding to the three soil samples, and measuring the carbon content of the three soil samples according to a potassium dichromate external heating method to obtain the carbon content corresponding to the three soil samples;
and calculating the carbon reserve of the soil in unit area of each area based on the water content coefficient and the carbon content, wherein the calculation formula of the carbon reserve of the soil in unit area is as follows:
wherein, C y1 The carbon reserve per unit area of the soil in the desertification and salinization areas, C y2 Carbon reserve per unit area in areas of low vegetation coverage, C y3 The reserve of carbon per unit area in the area of medium vegetation coverage, C y4 The unit area carbon storage capacity in the area with high vegetation coverage is disclosed, and M represents the corresponding unit area soil quality in the collected area; k 1 Is the water content coefficient, K, of the first soil sample 2 Is the water content coefficient, K, of the second soil sample 3 Is the water content coefficient, C, of a third soil sample n1 Carbon content, C, of the first soil sample n2 Carbon content, C, of the second soil sample n3 Carbon content of the third soil sample.
In a possible implementation manner, the method includes the steps of measuring the carbon content of the vegetation sample according to a dry combustion method, and calculating the carbon reserve of vegetation in each area according to the measured carbon content of the vegetation sample, and specifically includes:
respectively weighing the overground vegetation sample and the underground vegetation sample to obtain overground vegetation biomass and underground vegetation biomass;
respectively measuring the carbon content of the overground vegetation sample and the underground vegetation sample according to a dry combustion method to obtain the carbon content of the overground vegetation and the carbon content of the underground vegetation;
and calculating the carbon reserve of the vegetation in each area according to the carbon content of the vegetation on the ground, the carbon content of the underground vegetation, the biomass of the vegetation on the ground and the biomass of the underground vegetation, wherein the carbon reserve of the vegetation in each area is calculated according to the following formula:
C xi(i＝1,2,3,4) ＝C 1 ×N 1 +C 2 ×N 2 ；
wherein, C x1 The carbon reserve of vegetation in unit area of desertification and salinization areas, C x2 For the carbon reserve of vegetation in the area of low vegetation coverage, C x3 Carbon reserve per unit area of vegetation in medium vegetation coverage area, C x4 For carbon reserves of vegetation per unit area in areas of high vegetation coverage, N 1 Carbon content of overground vegetation, N 2 Carbon content of underground vegetation, C 1 Is the biomass of overground vegetation C 2 Is the biomass of underground vegetation.
In a possible implementation manner, the obtaining of the average content of each mineral in each surface rock sample and each drill core sample according to the mineral content corresponding to each surface rock sample and each drill core sample specifically includes:
obtaining a first calcium oxide content, a first magnesium oxide content, a first ferrous oxide content and a first loss on ignition in each surface rock sample, screening all the surface rock samples according to the first loss on ignition, and calculating the average value of each mineral content of the screened surface rock samples according to the following formula:
wherein, T Cover layer CaO The average content of the first calcium oxide in the screened surface rock sample is obtained; t is Capping layer of MgO The average content of the first magnesium oxide in the screened surface rock sample is obtained; t is Capping layer of FeO The average content of the first ferrous oxide in the screened surface rock sample, T CaO1 For the first calcium oxide content, T, of each surface rock sample after screening MgO1 For the first magnesium oxide content, T, of each surface rock sample after screening FeO1 The first ferrous oxide content of each screened surface rock sample, and n is the number of the screened surface rock samples;
obtaining a second calcium dioxide content, a second magnesium dioxide content, a second ferrous oxide content and a second loss on ignition in each drill core sample, screening all drill core samples according to the second loss on ignition, and calculating the average value of each mineral content of the drill core samples subjected to screening according to the following formula:
wherein, T Core CaO The average content of the second calcium dioxide in the screened drill core sample is obtained; t is a unit of Rock core MgO The average content of the second magnesium oxide in the screened drill core sample; t is Core FeO The average content of the second ferrous oxide, T, in the screened core samples CaO2 For the second calcium oxide content, T, of each drill core sample after screening MgO2 For the second magnesium oxide content, T, of each drill core sample after screening FeO2 And m is the second ferrous oxide content of each screened drill core sample, and the number of the screened drill core samples.
The invention also provides a system for estimating carbon emission and carbon reserve of a mineral resource base, which comprises: the system comprises a comprehensive energy consumption carbon emission estimation module, a region division module, an ecological carbon loss estimation module, a mineral carbonation carbon storage estimation module and a data integration module;
the comprehensive energy consumption carbon emission estimation module is used for constructing an ore body model of an area to be estimated, dividing the ore body model into a plurality of sub-modules, and calculating and obtaining comprehensive energy consumption carbon emission based on the ore amount corresponding to each type of ore according to the volume and small weight parameters of each sub-module;
the area division module is used for acquiring a normalized vegetation index map corresponding to the area to be estimated according to the remote sensing image map of the area to be estimated, carrying out area division on the area to be estimated according to the normalized vegetation index map, and calculating the area of each area;
the ecological carbon loss estimation module is used for acquiring a sampling sample of each region, measuring the carbon reserves of the sampling sample to obtain the carbon reserves corresponding to each region, and calculating the ecological carbon loss of the region to be estimated according to the area of the region and the carbon reserves;
the mineral carbonation carbon reserve estimation module is used for collecting the surface rock sample and the drill core sample of the area to be estimated, calculating and obtaining the mineral carbonation carbon reserve of the area to be estimated according to the average content of each mineral in the surface rock sample and the drill core sample;
and the data integration module is used for combining the ecological carbon loss amount and the comprehensive energy consumption carbon emission amount to obtain the carbon emission amount of the area to be estimated, and obtaining the carbon storage amount of the area to be estimated according to the mineral carbonation carbon storage amount.
In a possible implementation manner, the ecological carbon loss estimation module is configured to obtain a sampling sample of each area, and perform carbon reserve measurement on the sampling sample to obtain a carbon reserve corresponding to each area, and specifically includes:
setting and sampling vegetation and soil in each area according to a sampling line corresponding to each area to obtain a vegetation sample and a soil sample corresponding to each area, wherein the vegetation sample comprises an overground vegetation sample and an underground vegetation sample;
according to a dry burning method, determining the carbon content of the vegetation sample, and calculating the carbon reserve of the vegetation in each area according to the determined carbon content of the vegetation sample;
and measuring the carbon content of the soil sample according to a potassium dichromate external heating method, and calculating the carbon reserve of the soil in each area according to the measured carbon content of the soil sample.
In a possible implementation manner, the area dividing module is configured to perform area division on the area to be estimated according to the normalized vegetation index map, and calculate an area of each area, and specifically includes:
dividing the area to be estimated into a desertification and salinization area, a low vegetation coverage area, a medium vegetation coverage area and a high vegetation coverage area according to the normalized vegetation index map;
setting and counting the number of pixel units corresponding to each area and the precision value of the remote sensing image according to the normalized vegetation index threshold;
and calculating the area of each region according to the counted number of the pixel units and the precision value of the remote sensing image.
In a possible implementation manner, the mineral carbon carbonate reserve estimation module is configured to calculate and obtain the mineral carbon carbonate reserve of the area to be estimated according to an average content of each mineral in the surface rock sample and the borehole core sample, and specifically includes:
performing main quantity element analysis on the obtained surface rock sample and the obtained drill core sample respectively to obtain mineral contents corresponding to each surface rock sample and each drill core sample, wherein the mineral contents comprise calcium oxide content, magnesium oxide content and ferrous oxide content;
obtaining the average content of each mineral in each surface rock sample and each drill core sample according to the mineral content corresponding to each surface rock sample and each drill core sample;
and respectively calculating the relative molecular weight ratio of the carbon dioxide to the calcium oxide content, the magnesium oxide content and the ferrous oxide content, and obtaining the mineral carbonation carbon reserves of the area to be estimated according to the relative molecular weight ratio and the average content of each mineral.
The invention also provides a terminal device comprising a processor, a memory and a computer program stored in the memory and configured to be executed by the processor, wherein the processor executes the computer program to realize the carbon emission and carbon reserve estimation method of the mineral resource base.
The invention also provides a computer readable storage medium comprising a stored computer program, wherein the computer program when executed controls a device in which the computer readable storage medium is located to perform the method for estimating carbon emission and carbon reserve in a mineral resource base as described in any one of the above.
Compared with the prior art, the invention has the following beneficial effects:
according to the method, the ore amount corresponding to each type of ore in the area to be estimated is calculated through modeling, and the comprehensive energy consumption carbon emission amount of the area to be estimated is calculated based on the ore amount; meanwhile, on the basis of ecological damage in the development process, the area to be estimated is divided by obtaining a normalized vegetation index map of the area to be estimated, the carbon reserves of different areas are calculated by sampling samples, and the ecological carbon reserve of the area to be estimated is used as the ecological loss; and simultaneously, respectively obtaining an earth surface rock sample and a drill core sample of the area to be estimated, and calculating the mineral carbonation carbon storage capacity of the area to be estimated based on the carbon fixation capacity of each mineral by calculating the average content of each mineral in the earth surface rock sample and the drill core sample. Compared with the prior art, the method comprehensively considers the problems of main carbon emission and carbon reserve in the current mining development process, comprehensively and quantitatively estimates the comprehensive energy consumption carbon emission, ecological loss and mineral carbonation carbon reserve of the area to be estimated respectively to obtain the carbon emission and carbon reserve of the area to be estimated, fills the blank of a carbon dioxide emission and carbon reserve estimation system before mining development, improves the estimation efficiency of the carbon emission and carbon reserve before mining resource base development, and can provide important green support for the follow-up activities of environmental protection, resource utilization, geological exploration, mine planning and the like.
Drawings
FIG. 1 is a schematic flow chart diagram illustrating one embodiment of a method for estimating carbon emissions and carbon reserves in a mineral resource base according to the present invention;
FIG. 2 is a schematic diagram of a carbon emission and carbon reserve estimation system for a mineral resource base according to an embodiment of the present invention;
FIG. 3 is a representation of the mining integrated energy consumption limits for one embodiment of the present invention;
fig. 4 is a view illustrating the ore dressing comprehensive energy consumption limit according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without inventive step based on the embodiments of the present invention, are within the scope of protection of the present invention.
Example 1
Referring to fig. 1, fig. 1 is a schematic flow chart of an embodiment of a method for estimating carbon emission and carbon reserve of a mineral resource base according to the present invention, as shown in fig. 1, the method includes steps 101-104, which are as follows:
step 101: and constructing an ore body model of the area to be estimated, dividing the ore body model into a plurality of sub-modules, and calculating and obtaining the comprehensive energy consumption carbon emission amount based on the ore amount corresponding to each type of ore according to the volume and small weight parameters of each sub-module.
In one embodiment, the mining area is defined according to known mines within the mining right range of the mining company having the mining right, and the defined mining area is used as the area to be estimated in the embodiment, wherein the area to be estimated includes: a main mining area of the mine that has been drilled for determination, a known boundary area of the mine defined by geological conditions, and an area of the earth that the mine construction needs to occupy.
In one embodiment, the estimation work of the ore amount of the area to be estimated is completed through a geostatistical method, a distance inverse ratio method and the like provided in Micomie software, specifically, an ore body model is constructed through an ore body to be estimated in the area to be estimated, the ore body model is divided into a plurality of sub-modules, the size of each divided sub-module is set to be consistent with the minimum mining block section designed in mining, and the module volume corresponding to each sub-module can be directly estimated and obtained based on the setting.
In one embodiment, mining the deposit can be divided into open-pit mining and underground mining, so in this embodiment, each sub-module is classified according to the open-pit mining mode and the underground mining mode, drilling operation is performed on each classified sub-module to obtain a drilling sample corresponding to each sub-module, small weight parameters of the drilling sample are determined according to a wax sealing drainage method, humidity determination is performed at the same time to obtain small ore weight data corresponding to the sub-module, and ore weight data corresponding to the sub-module are obtainedCarrying out humidity correction on the ore small body weight data, and respectively calculating the ore small body weight average values corresponding to various sub-modules
Taking the average ore small weight of the sub-module corresponding to the open-pit mining mode as the average small weight of the open-pit ore; and taking the average value of the ore small weight of the sub-module corresponding to the underground mining mode as the average small weight of the underground ore.
Obtaining the volume V of various ores by calculating the sub-module volume corresponding to various sub-modules Ore ore And according to the volume V of each type of ore Ore ore And average small body weight of various minerals
Calculating the ore amount M of various ores Amount of ore The ore amount calculation formula is as follows:
in one embodiment, mining industry mining is divided into two parts of mining energy consumption and mineral separation energy consumption, and mining is divided into two parts of surface mining and underground mining, and mining energy consumption and mineral separation energy consumption of an area to be estimated are calculated according to a mining comprehensive energy consumption quota table and a mineral separation comprehensive energy consumption quota table specified by related industry standards, wherein the mining comprehensive energy consumption quota table is shown in fig. 3, and the mineral separation comprehensive energy consumption quota table is shown in fig. 4; the carbon emission of the consumed energy is also calculated according to the relevant standard of the supplied power, such as 3.465kg of CO2 emitted by 1kg of standard coal. Based on this, calculating the comprehensive energy consumption carbon emission P of the area to be estimated requires calculating the open-cut mining energy consumption carbon emission P in the mining energy consumption carbon emission P 1 And underground mining energy consumption carbon emission P 2 And the beneficiation energy consumption carbon emission P 3 The calculation formula is as follows:
P＝P 1 +P 2 +P 3 。
step 102: and acquiring a normalized vegetation index map corresponding to the area to be estimated according to the remote sensing image map of the area to be estimated, carrying out area division on the area to be estimated according to the normalized vegetation index map, and calculating the area of each area.
In one embodiment, the remote sensing image map of the area to be estimated is obtained, and the normalized vegetation index data of the area to be estimated is analyzed through the multispectral image to generate a normalized vegetation index map corresponding to the area to be estimated.
In one embodiment, the area to be estimated is divided into a desertification and salinization area, a low vegetation coverage area, a medium vegetation coverage area and a high vegetation coverage area according to the normalized vegetation index map.
In one embodiment, the number of pixel units and the precision value of the remote sensing image corresponding to each region are set and counted according to a normalized vegetation index threshold, and the region area of each region is calculated according to the counted number of pixel units and the precision value of the remote sensing image, wherein a region area calculation formula is as follows:
L (1,2,3,4) ＝I number of pixel units ×E Accuracy of remote sensing image 2 ；
Wherein L is 1 Is the area of the desertification and the salinization area; l is 2 Is the area of the area with low vegetation coverage; l is 3 The area of the vegetation coverage area is shown; l is 4 Is the area of the area with high vegetation coverage.
Step 103: and acquiring a sampling sample of each region, measuring the carbon reserves of the sampling sample to obtain the carbon reserves corresponding to each region, and calculating the ecological carbon loss of the region to be estimated according to the area of the region and the carbon reserves.
In one embodiment, vegetation and soil sampling is carried out on each area according to a sampling line corresponding to each area, so as to obtain a vegetation sample and a soil sample corresponding to each area, wherein the vegetation sample comprises an overground vegetation sample and an underground vegetation sample.
In each area range, a sampling line is arranged according to Arcgis and DGSS software, and equidistant experimental sample points are established on the sampling line, and the sampling interval can be adjusted according to the size of the area range, in the embodiment, 50m is used as the point interval, the number of the experimental sample points in each area is at most not more than 10, and the design specification of the sample is 1m × 1 m; because the vegetation is divided into overground vegetation and underground vegetation, the overground vegetation is harvested by a mowing method, namely harvesting the vegetation on the ground at the same time, and the vegetation harvested on the ground is filled into a pollution-free sample bag to finish the harvesting of the overground vegetation sample; digging a soil sample which is 0-15cm downward from the ground in a sample square, screening the underground roots of the vegetation in the soil sample, filling the soil in a pollution-free sample bag, weighing the soil which is screened and obtained from the roots of the vegetation, and then filling the soil in the sample square again to finish the collection of the underground vegetation sample; randomly drilling 3 times in a sample square by using a soil drill with the diameter of 5cm, and respectively taking three layers of soil samples of 0-5cm, 5-10cm and 10-15cm to finish the collection of the soil samples; because vegetation does not grow in the desertification and salinization areas, vegetation samples are not collected in the areas, and only soil samples are collected.
In one embodiment, the carbon content of the vegetation sample is measured according to a dry burning method, and the carbon reserve per unit area of the vegetation in each area is calculated according to the measured carbon content of the vegetation sample.
After the overground vegetation samples and the underground vegetation samples are collected, the overground vegetation samples and the underground vegetation samples are respectively cleaned for 2-3 times by using purified water, the samples are crushed after water is dried and screened through a 100-mesh sieve, the screened samples are dried for 24 hours at 80 ℃ to constant weight, the samples are weighed and recorded and numbered, and the weights of the samples are recorded as overground vegetation biomass C 1 And underground plantingQuilt biomass C 2 。
The carbon content of the treated overground vegetation sample and the carbon content of the underground vegetation sample are measured by a dry burning method and are respectively recorded as the carbon content N of the overground vegetation 1 And carbon content N of underground vegetation 2 Then, the carbon reserve of vegetation in unit area in the corresponding region is calculated as follows:
C xi(i＝1,2,3,4) ＝C 1 ×N 1 +C 2 ×N 2 ；
wherein, C x1 The carbon reserve of vegetation in unit area of desertification and salinization areas; c x2 The carbon reserves of the vegetation in the unit area of the area with low vegetation coverage; c x3 The carbon reserve of vegetation in a unit area in a medium vegetation coverage area; c x4 The carbon reserve of the vegetation in the unit area of the high vegetation coverage area.
After the vegetation samples in each area are measured, integrating the carbon reserves of the vegetation in the area corresponding to each area, calculating the mean carbon reserve of the vegetation in the unit area corresponding to each area, and recording the mean carbon reserve as the mean carbon reserve
In one embodiment, the carbon content of the soil sample is measured according to a potassium dichromate external heating method, and the carbon reserve per unit area of the soil in each area is calculated according to the measured carbon content of the soil sample.
Numbering the collected soil samples according to the samples collected in layers, randomly selecting three-part soil samples from the three soil samples at each sampling point, respectively measuring the water content of the three soil samples at each sampling point, and recording as a water content coefficient K 1 、K 2 、K 3 Air-drying the soil sample, and measuring the carbon content of three randomly selected soil samples by using a concentrated sulfuric acid-potassium dichromate external heating method, and recording the carbon content as C n1 、C n2 、C n3 Then, the calculation formula of the carbon storage of the soil per unit area of each area is as follows:
wherein, C y1 The carbon reserve of the soil in unit area in the desertification and salinization areas; c y2 The carbon reserve per unit area in the area with low vegetation coverage; c y3 The unit area carbon reserve in the vegetation coverage area is obtained; c y4 The unit area carbon reserve in the area with high vegetation coverage is shown, and M represents the corresponding unit area soil quality in the collected area.
After the soil samples in each area are measured, integrating the data of the carbon reserves of the soil in the unit area of the area corresponding to each area, calculating the average value of the carbon reserves of the soil in the unit area corresponding to each area, and recording the average value as the average value
In an embodiment, the ecological carbon reserves of the corresponding regions are calculated according to the region areas of the regions, the average value of the carbon reserves of the soil in the unit area corresponding to each region, and the average value of the carbon reserves of the vegetation, which are calculated in step 102, and the calculation formula is as follows:
wherein, C Desertification and salinization The ecological carbon reserves of the desertification and salinization areas; c Low vegetation coverage For areas of low vegetation coverageEcological carbon reserves of the domain; c Coverage of medium vegetation The ecological carbon reserve of the vegetation coverage area is shown; c High vegetation coverage The ecological carbon reserve of the high vegetation coverage area.
In one embodiment, the obtained ecological carbon reserves corresponding to each region are counted and used as the ecological carbon reserves of the region to be estimated; because vegetation and soil in the development area are damaged when the mineral resource base is developed, and the ecological carbon reserve of the area to be estimated is lost, in this embodiment, the ecological carbon reserve of the area to be estimated is set as the ecological carbon loss of the area to be estimated, and the calculation formula is as follows:
C loss of ecological carbon ＝C Ecological carbon reserve ＝C Desertification and salinization +C Low vegetation coverage +C Coverage of medium vegetation +
C High vegetation coverage 。
Step 104: and acquiring the surface rock sample and the drill core sample of the area to be estimated, and calculating and obtaining the mineral carbonation carbon reserves of the area to be estimated according to the average content of each mineral in the surface rock sample and the drill core sample.
In one embodiment, 1:2000 map filling is carried out on a main mining area of a mineral deposit of an area to be estimated through DGSS software, geological map filling routes are designed, wherein the geological map filling routes are designed to be 100m in distance of each route, point distances on each route are designed to be about 50m, and the route length is determined according to the size of a mineral deposit range. And (3) acquiring a corresponding surface rock sample at each point position on the route, if the surface rock sample acquired at the point position has a relevant stratum outcrop, stripping a surface weathered layer from the acquired surface rock sample, and acquiring a surface scattered rock sample which has no obvious crack and is not easy to crack for the point position without the stratum outcrop, wherein the sample size is not smaller than 10cm multiplied by 10 cm.
In one embodiment, the collection of the drill core samples is completed by utilizing the drilling work of a mining company on a mineral deposit, the related logging work is carried out on the core before the collection, the sampling work is carried out according to different lithological characters, the lithological layering part is taken as an initial point, samples are taken at an interval of 50m, the sample size is about 8-12cm, the sampling quantity is adjusted according to the drilling depth, when the drill core samples are collected, the collection of the drill core samples with cracks or loose holes is avoided, the depths of the drill core samples are marked according to the drilling hole numbers, and the related lithological characters are marked according to the information of the logging work.
In one embodiment, the obtained surface rock samples and the obtained borehole core samples are subjected to main quantity element analysis to obtain mineral contents corresponding to each surface rock sample and each borehole core sample, wherein the mineral contents include calcium oxide content, magnesium oxide content and ferrous oxide content, and the average contents of minerals in the surface rock samples and the borehole core samples are obtained according to the mineral contents corresponding to each surface rock sample and each borehole core sample.
After obtaining the surface rock sample and the drill core sample, respectively cleaning the surface rock sample and the drill core sample, airing the surface moisture, after the moisture is aired, classifying the surface rock sample and the drill core sample according to different lithologies, and after the classification is finished, carrying out optical sheet manufacturing on the surface rock sample and the drill core sample, wherein the surface rock sample is reserved about 5cm multiplied by 5cm as a backup sample, the drill core sample is reserved about 3cm long as the backup sample, the rest parts are crushed, and the crushing is carried out by using quartz abrasive materials, so that the surface rock sample powder and the drill core sample powder are obtained.
Respectively obtaining 10g of surface rock sample powder, and performing principal component element analysis by using a wavelength dispersion X-ray fluorescence spectrometer to determine the first calcium oxide FeO content, the first magnesium oxide MgO content, the first ferrous oxide FeO content and the first loss on ignition LOI in the surface rock sample.
According to the first loss on ignition LOI, removing samples with the first loss on ignition LOI more than 10% in the surface rock sample data, and calculating the average value of the mineral contents of the surface rock samples subjected to removal treatment according to the following formula:
wherein, T Cover layer CaO The average content of CaO in the screened surface rock sample is the first calcium oxide; t is Capping layer of MgO The average content of first magnesium oxide (MgO) in the screened surface rock sample is shown; t is Capping layer of FeO The average content of the first ferrous oxide FeO in the screened surface rock sample, T CaO1 The first calcium oxide CaO content, T, of each ground rock sample after screening MgO1 The first magnesium oxide MgO content, T, of each surface rock sample after screening FeO1 And n is the quantity of the screened surface rock samples, wherein the content of the first ferrous oxide FeO of each screened surface rock sample is the quantity of the screened surface rock samples.
Respectively obtaining 10g of drill core sample powder, performing principal component element analysis by using a wavelength dispersion X-ray fluorescence spectrometer, and determining the content of second calcium dioxide FeO, the content of second magnesium dioxide MgO, the content of second ferrous oxide FeO and a second loss on ignition LOI in the surface rock sample.
According to the second loss on ignition LOI, removing samples with the second loss on ignition LOI more than 10% in the drill core sample data, and calculating the average value of the mineral content of the drill core samples after removal according to the following formula:
wherein, T Core CaO The average content of CaO in the screened drill core sample is the average content of CaO in the second calcium oxide; t is Rock core MgO The average content of second magnesium oxide MgO in the screened drill core sample is shown; t is Core FeO The average content of the second ferrous oxide FeO in the screened core sample, T CaO2 The CaO content of the second calcium dioxide of each drill core sample after screening, T MgO2 The second magnesium oxide MgO content, T, of each drill core sample after screening FeO2 And m is the number of the screened drill core samples, wherein the content of the second ferrous oxide FeO of each screened drill core sample is the content of the second ferrous oxide FeO of each screened drill core sample.
In one embodiment, the relative molecular weight ratios of carbon dioxide to the calcium oxide content, the magnesium oxide content and the ferrous oxide content are calculated respectively, and the mineral carbonation carbon reserves of the area to be estimated are obtained according to the relative molecular weight ratios and the average contents of the minerals.
Since the carbon energy storage capacity of the mineral is related to the chemical components of the mineral, wherein calcium, magnesium and iron minerals exist in the rock in the form of silicate and are good carbon fixing raw materials, and the corresponding oxide content is determined by performing principal element analysis on an earth surface rock sample and a drill core sample, the chemical formula of the mineral carbonation reaction is as follows:
CaO+CO 2 ＝CaCO 3 ；
MgO+CO 2 ＝MgCO 3 ；
FeO+CO 2 ＝FeCO 3 。
calculating relative molecular weight ratios of carbon dioxide to the calcium oxide content, the magnesium oxide content, and the ferrous oxide content based on the mineral carbonation reaction formula, wherein the relative molecular weight ratio of carbon dioxide to calcium oxide content is
The relative molecular weight ratio of the contents of carbon dioxide and magnesium oxide is
The relative molecular weight ratio of the carbon dioxide to the ferrous oxide content is
Estimating the carbon fixed quantity on the earth surface, namely the carbon reserve of the cover layer, according to the mineral quantity of the open-pit ore obtained in the step 101, the average content of each mineral in the screened earth surface rock sample and the relative molecular weight ratio of carbon dioxide to each mineral, wherein the estimation formula is as follows:
since calcium, magnesium and iron minerals are not the main economic minerals in the deposit, the calcium, magnesium and iron minerals are used as the tailings, and therefore the carbon fixed quantity of the core, namely the carbon reserve of the tailings, is estimated according to the mineral quantity of the ore under the well obtained in the step 101, the average content of each mineral in the screened drilling core sample and the relative molecular weight ratio of carbon dioxide to each mineral, and the estimation formula is as follows:
in one embodiment, the carbon reserves of the cap layer and the tailings are integrated to obtain the carbon reserves of the mineral carbonation in the area to be estimated, wherein the integrated formula is as follows:
C mineral carbon carbonate reserves ＝C Carbon reserve of cap layer +C Carbon reserve of tailings 。
Step 105: and combining the ecological carbon loss amount and the comprehensive energy consumption carbon emission amount to obtain the carbon emission amount of the area to be estimated, and obtaining the carbon storage amount of the area to be estimated according to the mineral carbonation carbon storage amount.
In one embodiment, the comprehensive energy consumption carbon emission of the area to be estimated obtained in step 101, the ecological carbon loss obtained in step 103, and the mineral carbonation carbon storage obtained in step 104 are correspondingly integrated and unitized to obtain the overall carbon emission and carbon storage numerical value of the deposit, so that the basic carbon emission and carbon storage data related to the carbon trading cost of the subsequent enterprise can be visually displayed in front of the ore enterprise, the control of the carbon emission cost in the development process of the ore enterprise is facilitated, and a decision basis can be provided for the related development process of the ore enterprise. The integration is as follows:
C discharge capacity ＝C Comprehensive energy consumption and carbon emission +C Loss of ecological carbon ；
C Carbon reserve ＝C Mineral carbon carbonate reserves 。
In the embodiment, the acquired data are integrated, so that the basic carbon emission and carbon reserve data related to the carbon transaction cost of the follow-up enterprise can be visually displayed in front of the mining enterprise, the control of the carbon emission cost in the development process of the mining enterprise is facilitated, and a decision basis can be provided for the related development process of the mining enterprise.
In summary, the method for estimating carbon emission and carbon reserve of a mineral resource base provided by this embodiment is based on the development characteristics of the mineral resource base, and can perform comprehensive quantitative estimation on the comprehensive energy consumption carbon emission, ecological loss and mineral carbonation carbon reserve of an area to be estimated respectively by comprehensively considering the problems of main carbon emission and carbon reserve faced in the current mineral development process, so as to obtain the carbon emission and carbon reserve of the area to be estimated, fill up the blank of a carbon dioxide emission and carbon reserve estimation system before mineral development, improve the estimation efficiency of the carbon emission and carbon reserve before mineral development, and greatly reduce the related estimation investment of a mineral enterprise; similarly, the method for estimating carbon emission and carbon reserve of the mineral resource base provided by this embodiment may also perform corresponding adjustment according to the characteristics of the mineral resource base, such as geographic conditions, ecological conditions, and dominant mineral species, and only perform single estimation on the main carbon emission or carbon reserve in the mining development process, thereby satisfying the estimation requirements of enterprises.
Example 2
Referring to fig. 2, fig. 2 is a schematic structural diagram of an embodiment of a system for estimating carbon emission and carbon reserve in a mineral resource base according to the present invention, as shown in fig. 2, the system includes: the comprehensive energy consumption carbon emission amount estimation module 201, the area division module 202, the ecological carbon loss amount estimation module 203, the mineral carbon carbonation amount estimation module 204 and the data integration module 205 are specifically as follows:
the comprehensive energy consumption carbon emission estimation module 201 is used for constructing an ore body model of an area to be estimated, dividing the ore body model into a plurality of sub-modules, and calculating and obtaining the comprehensive energy consumption carbon emission based on the ore amount corresponding to each type of ore according to the volume and small weight parameters of each sub-module.
The area division module 202 is configured to obtain a normalized vegetation index map corresponding to the area to be estimated according to the remote sensing image map of the area to be estimated, perform area division on the area to be estimated according to the normalized vegetation index map, and calculate an area of each area.
And the ecological carbon loss estimation module 203 is configured to obtain a sampling sample of each region, determine a carbon reserve of the sampling sample, obtain a carbon reserve corresponding to each region, and calculate an ecological carbon loss of the region to be estimated according to the region area and the carbon reserve.
And the mineral carbon carbonate storage estimation module 204 is configured to collect the surface rock sample and the drill core sample of the area to be estimated, calculate and obtain the mineral carbon carbonate storage of the area to be estimated according to the average content of each mineral in the surface rock sample and the drill core sample.
And the data integration module 205 is configured to combine the ecological carbon loss amount and the integrated energy consumption carbon emission amount to obtain the carbon emission amount of the area to be estimated, and obtain the carbon storage amount of the area to be estimated according to the mineral carbon carbonation carbon storage amount.
In an embodiment, the ecological carbon loss estimating module 203 is configured to obtain a sampling sample of each area, and perform carbon storage measurement on the sampling sample to obtain a carbon storage corresponding to each area, and specifically includes: setting and sampling vegetation and soil in each area according to a sampling line corresponding to each area to obtain a vegetation sample and a soil sample corresponding to each area, wherein the vegetation sample comprises an overground vegetation sample and an underground vegetation sample; according to a dry burning method, determining the carbon content of the vegetation sample, and calculating the carbon reserve of the vegetation in each area according to the determined carbon content of the vegetation sample; and measuring the carbon content of the soil sample according to a potassium dichromate external heating method, and calculating the carbon reserve of the soil in each area according to the measured carbon content of the soil sample.
In an embodiment, the area dividing module 202 is configured to perform area division on the area to be estimated according to the normalized vegetation index map, and calculate an area of each area, and specifically includes: dividing the region to be estimated into a desertification and salinization region, a low vegetation coverage region, a medium vegetation coverage region and a high vegetation coverage region according to the normalized vegetation index map; setting and counting the number of pixel units corresponding to each area and the precision value of the remote sensing image according to the normalized vegetation index threshold; and calculating the area of each region according to the counted number of the pixel units and the precision value of the remote sensing image.
In an embodiment, the mineral carbon carbonate storage estimation module 204 is configured to calculate and obtain the mineral carbon carbonate storage of the area to be estimated according to an average content of each mineral in the surface rock sample and the borehole core sample, and specifically includes: performing main quantity element analysis on the obtained surface rock sample and the obtained drill core sample respectively to obtain mineral contents corresponding to each surface rock sample and each drill core sample, wherein the mineral contents comprise calcium oxide content, magnesium oxide content and ferrous oxide content; obtaining the average content of each mineral in each surface rock sample and each drill core sample according to the mineral content corresponding to each surface rock sample and each drill core sample; and respectively calculating the relative molecular weight ratio of the carbon dioxide to the calcium oxide content, the magnesium oxide content and the ferrous oxide content, and obtaining the mineral carbonation carbon reserves of the area to be estimated according to the relative molecular weight ratio and the average content of each mineral.
In one embodiment, the ecological carbon loss estimation module 203 is configured to determine the carbon content of the soil sample according to a potassium dichromate external heating method, and calculate the ecological carbon loss estimation module for the soil carbon storage per unit area in each region according to the determined carbon content of the soil sample, and is specifically configured to: randomly selecting three soil samples from all soil samples, carrying out water content measurement on the three soil samples to obtain water content coefficients corresponding to the three soil samples, and measuring the carbon content of the three soil samples according to a potassium dichromate external heating method to obtain the carbon content corresponding to the three soil samples; and calculating the carbon reserve of the soil in unit area of each area based on the water content coefficient and the carbon content, wherein the calculation formula of the carbon reserve of the soil in unit area is as follows:
wherein, C y1 The carbon reserve per unit area of the soil in the desertification and salinization areas, C y2 Carbon reserve per unit area of soil in areas of low vegetation coverage, C y3 Is the unit area carbon reserve in the medium vegetation coverage area, C y4 The unit area carbon storage capacity in the area with high vegetation coverage is disclosed, and M represents the corresponding unit area soil quality in the collected area; k 1 Is the water content coefficient, K, of the first soil sample 2 Is the water content coefficient, K, of the second soil sample 3 Is the water content coefficient, C, of a third soil sample n1 Carbon content, C, of the first soil sample n2 Carbon content, C, of the second soil sample n3 Carbon content of the third soil sample.
In an embodiment, the ecological carbon loss estimation module 203 is configured to determine the carbon content of the vegetation sample according to a dry burning method, and calculate the carbon reserve of vegetation per unit area in each area according to the determined carbon content of the vegetation sample, and specifically includes: respectively weighing the overground vegetation sample and the underground vegetation sample to obtain overground vegetation biomass and underground vegetation biomass; respectively measuring the carbon content of the overground vegetation sample and the underground vegetation sample according to a dry burning method to obtain the carbon content of the overground vegetation and the carbon content of the underground vegetation; and calculating the carbon reserve of the vegetation in each area according to the carbon content of the vegetation on the ground, the carbon content of the underground vegetation, the biomass of the vegetation on the ground and the biomass of the underground vegetation, wherein the carbon reserve of the vegetation in each area is calculated according to the following formula:
C xi(i＝1,2,3,4) ＝C 1 ×N 1 +C 2 ×N 2 ；
wherein, C x1 The carbon reserve of vegetation in unit area of desertification and salinization areas, C x2 For the carbon reserve of vegetation in the area of low vegetation coverage, C x3 Carbon reserve per unit area of vegetation in medium vegetation coverage area, C x4 For carbon reserves of vegetation per unit area in areas of high vegetation coverage, N 1 Carbon content of overground vegetation, N 2 Carbon content of underground vegetation, C 1 Is the biomass of overground vegetation C 2 Is the biomass of underground vegetation.
In an embodiment, the mineral carbonated carbon reserve estimation module 204 is configured to obtain an average content of each mineral in each surface rock sample and each drill core sample according to a mineral content corresponding to each surface rock sample and each drill core sample, and specifically includes: obtaining a first calcium oxide content, a first magnesium oxide content, a first ferrous oxide content and a first loss on ignition in each surface rock sample, screening all the surface rock samples according to the first loss on ignition, and calculating the average value of each mineral content of the screened surface rock samples according to the following formula:
wherein, T Cover layer CaO The average content of the first calcium oxide in the screened surface rock sample is obtained; t is Capping layer of MgO The average content of the first magnesium oxide in the screened surface rock sample is obtained; t is Capping layer of FeO The average content of the first ferrous oxide in the screened surface rock sample, T CaO1 For the first calcium oxide content, T, of each surface rock sample after screening MgO1 The first magnesium oxide content, T, of each surface rock sample after screening FeO1 The first ferrous oxide content of each screened surface rock sample, and n is the number of the screened surface rock samples;
obtaining a second calcium dioxide content, a second magnesium dioxide content, a second ferrous oxide content and a second loss on ignition in each drill core sample, screening all drill core samples according to the second loss on ignition, and calculating the average value of each mineral content of the drill core samples subjected to screening according to the following formula:
wherein, T Core CaO The average content of the second calcium dioxide in the screened drill core sample is determined; t is a unit of Rock core MgO The average content of the second magnesium oxide in the screened drill core sample; t is FeO core The average content of the second ferrous oxide, T, in the screened core samples CaO2 For the second calcium oxide content, T, of each drill core sample after screening MgO2 For the second magnesium oxide content, T, of each drill core sample after screening FeO2 And m is the second ferrous oxide content of each screened drill core sample, and the number of the screened drill core samples.
It can be clearly understood by those skilled in the art that, for convenience and brevity of description, the specific working process of the system described above may refer to the corresponding process in the foregoing method embodiment, and is not described herein again.
It should be noted that the above-mentioned embodiments of the system for estimating carbon emissions and carbon reserves of a mineral resource base are merely illustrative, wherein the modules described as separate components may or may not be physically separate, and the components displayed as modules may or may not be physical units, i.e. may be located in one place, or may be distributed over a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment.
On the basis of the embodiment of the carbon emission and carbon storage estimation method for the mineral resource base, another embodiment of the invention provides a carbon emission and carbon storage estimation terminal device for the mineral resource base, which comprises a processor, a memory and a computer program stored in the memory and configured to be executed by the processor, wherein the processor executes the computer program to realize the carbon emission and carbon storage estimation method for the mineral resource base according to any one embodiment of the invention.
On the basis of the above embodiment of the method for estimating carbon emission and carbon reserve of a mineral resource base, another embodiment of the present invention provides a storage medium including a stored computer program, wherein when the computer program runs, a device in which the storage medium is located is controlled to execute the method for estimating carbon emission and carbon reserve of a mineral resource base according to any one of the embodiments of the present invention.
In summary, according to the method and the system for estimating carbon emission and carbon reserve of the mineral resource base, provided by the invention, the ore body model of the area to be estimated is constructed and divided into a plurality of sub-modules, and the corresponding ore amount of each type of ore is calculated and based on the volume and small weight parameters of each sub-module, so that the comprehensive energy consumption carbon emission is obtained; acquiring a normalized vegetation index map corresponding to the area to be estimated according to the remote sensing image map of the area to be estimated, carrying out area division on the area to be estimated according to the normalized vegetation index map, and calculating the area of each area; acquiring a sampling sample of each region, measuring the carbon reserves of the sampling sample to obtain the carbon reserves corresponding to each region, and calculating the ecological carbon loss of the region to be estimated according to the area of the region and the carbon reserves; collecting the surface rock sample and the drill core sample of the area to be estimated, and calculating and obtaining the mineral carbonation carbon reserves of the area to be estimated according to the average content of each mineral in the surface rock sample and the drill core sample; and combining the ecological carbon loss amount and the comprehensive energy consumption carbon emission amount to obtain the carbon emission amount of the area to be estimated, and obtaining the carbon storage amount of the area to be estimated according to the mineral carbonation carbon storage amount. Compared with the prior art, the method has the advantages that the carbon emission and the carbon reserve of the mineral resource base are integrally estimated, so that the estimation efficiency of the carbon emission and the carbon reserve before the development of the mineral resources is improved, and the blank of a carbon dioxide emission and carbon reserve estimation system before the development of the mineral industry is filled at present.
The above description is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, many modifications and substitutions can be made without departing from the technical principle of the present invention, and these modifications and substitutions should also be regarded as the protection scope of the present invention.
Claims (10)
1. A method for estimating carbon emissions and carbon reserves in a mineral resource base, comprising:
constructing an ore body model of a region to be estimated, dividing the ore body model into a plurality of sub-modules, and calculating and obtaining the comprehensive energy consumption carbon emission amount based on the ore amount corresponding to each type of ore according to the volume and small weight parameters of each sub-module;
acquiring a normalized vegetation index map corresponding to the area to be estimated according to the remote sensing image map of the area to be estimated, carrying out area division on the area to be estimated according to the normalized vegetation index map, and calculating the area of each area;
acquiring a sampling sample of each region, measuring the carbon reserves of the sampling sample to obtain the carbon reserves corresponding to each region, and calculating the ecological carbon loss of the region to be estimated according to the area of the region and the carbon reserves;
collecting the surface rock sample and the drill core sample of the area to be estimated, and calculating and obtaining the mineral carbonation carbon reserves of the area to be estimated according to the average content of each mineral in the surface rock sample and the drill core sample;
and combining the ecological carbon loss amount and the comprehensive energy consumption carbon emission amount to obtain the carbon emission amount of the area to be estimated, and obtaining the carbon storage amount of the area to be estimated according to the mineral carbonation carbon storage amount.
2. The method for estimating carbon emissions and carbon reserves in a mineral resource base as set forth in claim 1, wherein the step of obtaining a sample for each area and measuring the carbon reserves of the sample to obtain the carbon reserves corresponding to each area comprises:
setting and sampling vegetation and soil in each area according to a sampling line corresponding to each area to obtain a vegetation sample and a soil sample corresponding to each area, wherein the vegetation sample comprises an overground vegetation sample and an underground vegetation sample;
according to a dry burning method, determining the carbon content of the vegetation sample, and calculating the carbon reserve of the vegetation in each area according to the determined carbon content of the vegetation sample;
and measuring the carbon content of the soil sample according to a potassium dichromate external heating method, and calculating the carbon reserve of the soil in each area according to the measured carbon content of the soil sample.
3. The method for estimating carbon emission and carbon reserve in a mineral resource base according to claim 2, wherein the dividing the area to be estimated according to the normalized vegetation index map and calculating the area of each area comprises:
dividing the area to be estimated into a desertification and salinization area, a low vegetation coverage area, a medium vegetation coverage area and a high vegetation coverage area according to the normalized vegetation index map;
setting and counting the number of pixel units corresponding to each area and the precision value of the remote sensing image according to the normalized vegetation index threshold;
and calculating the area of each region according to the counted number of the pixel units and the precision value of the remote sensing image.
4. The method for estimating carbon emissions and carbon reserves in a mineral resource base of claim 1, wherein the calculating and obtaining the mineral carbonation carbon reserves in the area to be estimated based on the average mineral contents in the surface rock sample and the borehole core sample comprises:
performing main quantity element analysis on the obtained surface rock sample and the obtained drill core sample respectively to obtain mineral contents corresponding to each surface rock sample and each drill core sample, wherein the mineral contents comprise calcium oxide content, magnesium oxide content and ferrous oxide content;
obtaining the average content of each mineral in each earth surface rock sample and each drill hole core sample according to the mineral content corresponding to each earth surface rock sample and each drill hole core sample;
and respectively calculating the relative molecular weight ratio of the carbon dioxide to the calcium oxide content, the magnesium oxide content and the ferrous oxide content, and obtaining the mineral carbonated carbon reserve of the area to be estimated according to the relative molecular weight ratio and the average content of each mineral.
5. The method for estimating carbon emission and carbon reserve in a mineral resource base according to claim 3, wherein said measuring the carbon content of said soil sample by potassium dichromate external heating and calculating the carbon reserve per unit area of each area based on the measured carbon content of said soil sample comprises:
randomly selecting three soil samples from all soil samples, carrying out water content measurement on the three soil samples to obtain water content coefficients corresponding to the three soil samples, and measuring the carbon content of the three soil samples according to a potassium dichromate external heating method to obtain the carbon content corresponding to the three soil samples;
and calculating the carbon reserve of the soil in unit area of each area based on the water content coefficient and the carbon content, wherein the calculation formula of the carbon reserve of the soil in unit area is as follows:
wherein, C y1 The carbon reserve per unit area of the soil in the desertification and salinization areas, C y2 Carbon reserve per unit area of soil in areas of low vegetation coverage, C y3 The reserve of carbon per unit area in the area of medium vegetation coverage, C y4 The unit area carbon storage capacity in the area with high vegetation coverage is disclosed, and M represents the corresponding unit area soil quality in the collected area; k 1 Is the water content coefficient, K, of the first soil sample 2 Is the water content coefficient, K, of the second soil sample 3 Is the third soilWater content coefficient of sample, C n1 Carbon content, C, of the first soil sample n2 Carbon content, C, of the second soil sample n3 Carbon content of the third soil sample.
6. The method of estimating carbon emissions and carbon reserves of a mineral resource base of claim 3, wherein said determining the carbon content of the vegetation sample according to a dry-fire method and calculating the carbon reserve of vegetation per unit area of each area based on the determined carbon content of the vegetation sample comprises:
respectively weighing the overground vegetation sample and the underground vegetation sample to obtain overground vegetation biomass and underground vegetation biomass;
respectively measuring the carbon content of the overground vegetation sample and the underground vegetation sample according to a dry burning method to obtain the carbon content of the overground vegetation and the carbon content of the underground vegetation;
and calculating the carbon reserve of the vegetation in each area according to the carbon content of the vegetation on the ground, the carbon content of the underground vegetation, the biomass of the vegetation on the ground and the biomass of the underground vegetation, wherein the carbon reserve of the vegetation in each area is calculated according to the following formula:
C xi(i＝1,2,3,4) ＝C 1 ×N 1 +C 2 ×N 2 ；
wherein, C x1 The carbon reserve of vegetation in unit area of desertification and salinization areas, C x2 For the carbon reserve of vegetation in the area of low vegetation coverage, C x3 Carbon reserve per unit area of vegetation in medium vegetation coverage area, C x4 For carbon reserves of vegetation per unit area in areas of high vegetation coverage, N 1 Carbon content of overground vegetation, N 2 Carbon content of underground vegetation, C 1 Is the biomass of overground vegetation C 2 Is the biomass of underground vegetation.
7. The method for estimating carbon emission and carbon reserve in a mineral resource base according to claim 4, wherein the obtaining of the average content of each mineral in the surface rock sample and the borehole core sample according to the mineral content corresponding to each surface rock sample and each borehole core sample comprises:
obtaining a first calcium oxide content, a first magnesium oxide content, a first ferrous oxide content and a first loss on ignition in each surface rock sample, screening all the surface rock samples according to the first loss on ignition, and calculating the average value of each mineral content of the screened surface rock samples according to the following formula:
wherein, T Cover layer CaO The average content of the first calcium oxide in the screened surface rock sample is obtained; t is Capping layer of MgO The average content of the first magnesium oxide in the screened surface rock sample is obtained; t is Capping layer of FeO The average content of the first ferrous oxide in the screened surface rock sample, T CaO1 For the first calcium oxide content, T, of each surface rock sample after screening MgO1 The first magnesium oxide content, T, of each surface rock sample after screening FeO1 The first ferrous oxide content of each screened surface rock sample, and n is the number of the screened surface rock samples;
obtaining a second calcium dioxide content, a second magnesium dioxide content, a second ferrous oxide content and a second loss on ignition in each drill core sample, screening all drill core samples according to the second loss on ignition, and calculating the average value of each mineral content of the drill core samples subjected to screening according to the following formula:
wherein, T Core CaO The average content of the second calcium dioxide in the screened drill core sample is obtained; t is Rock core MgO The average content of the second magnesium oxide in the screened drill core sample; t is FeO core The average content of the second ferrous oxide, T, in the screened core samples CaO2 For the second calcium oxide content, T, of each drill core sample after screening MgO2 For the second magnesium oxide content, T, of each drill core sample after screening FeO2 And m is the second ferrous oxide content of each screened drill core sample, and the number of the screened drill core samples.
8. A system for estimating carbon emissions and carbon reserves in a mineral resource base, comprising: the system comprises a comprehensive energy consumption carbon emission estimation module, a region division module, an ecological carbon loss estimation module, a mineral carbonation carbon storage estimation module and a data integration module;
the comprehensive energy consumption carbon emission estimation module is used for constructing an ore body model of an area to be estimated, dividing the ore body model into a plurality of sub-modules, and calculating and obtaining the comprehensive energy consumption carbon emission based on the ore amount corresponding to each type of ore according to the volume and small weight parameters of each sub-module;
the area division module is used for acquiring a normalized vegetation index map corresponding to the area to be estimated according to the remote sensing image map of the area to be estimated, carrying out area division on the area to be estimated according to the normalized vegetation index map, and calculating the area of each area;
the ecological carbon loss estimation module is used for acquiring a sampling sample of each region, measuring the carbon reserves of the sampling sample to obtain the carbon reserves corresponding to each region, and calculating the ecological carbon loss of the region to be estimated according to the area of the region and the carbon reserves;
the mineral carbonation carbon reserve estimation module is used for collecting the surface rock sample and the drill core sample of the area to be estimated, calculating and obtaining the mineral carbonation carbon reserve of the area to be estimated according to the average content of each mineral in the surface rock sample and the drill core sample;
and the data integration module is used for combining the ecological carbon loss amount and the comprehensive energy consumption carbon emission amount to obtain the carbon emission amount of the area to be estimated, and obtaining the carbon storage amount of the area to be estimated according to the mineral carbonation carbon storage amount.
9. A terminal device comprising a processor, a memory and a computer program stored in the memory and configured to be executed by the processor, the processor when executing the computer program implementing the method of carbon emission and carbon reserve estimation of a mineral resource base according to any one of claims 1 to 7.
10. A computer-readable storage medium, comprising a stored computer program, wherein the computer program, when executed, controls an apparatus in which the computer-readable storage medium is located to perform the method of estimating carbon emissions and carbon reserves of a mineral resource base according to any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210620129.7A CN114969665B (en) | 2022-06-02 | 2022-06-02 | Method and system for estimating carbon emission and carbon reserve of mineral resource base |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210620129.7A CN114969665B (en) | 2022-06-02 | 2022-06-02 | Method and system for estimating carbon emission and carbon reserve of mineral resource base |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114969665A true CN114969665A (en) | 2022-08-30 |
| CN114969665B CN114969665B (en) | 2023-02-10 |
Family
ID=82960028
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210620129.7A Active CN114969665B (en) | 2022-06-02 | 2022-06-02 | Method and system for estimating carbon emission and carbon reserve of mineral resource base |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114969665B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010119962A (en) * | 2008-11-20 | 2010-06-03 | National Institute Of Advanced Industrial Science & Technology | Carbon dioxide storage in shallow aquifer |
| CN109544038A (en) * | 2018-12-20 | 2019-03-29 | 西北农林科技大学 | A kind of biomass carbon discharge estimating system and method based on life cycle analysis |
| CN111488902A (en) * | 2020-01-14 | 2020-08-04 | 沈阳农业大学 | Method and system for quantitatively estimating carbon reserves of ecosystem of primary coastal wetland |
| CN112836610A (en) * | 2021-01-26 | 2021-05-25 | 平衡机器科技(深圳)有限公司 | Land use change and carbon reserve quantitative estimation method based on remote sensing data |
-
2022
- 2022-06-02 CN CN202210620129.7A patent/CN114969665B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010119962A (en) * | 2008-11-20 | 2010-06-03 | National Institute Of Advanced Industrial Science & Technology | Carbon dioxide storage in shallow aquifer |
| CN109544038A (en) * | 2018-12-20 | 2019-03-29 | 西北农林科技大学 | A kind of biomass carbon discharge estimating system and method based on life cycle analysis |
| CN111488902A (en) * | 2020-01-14 | 2020-08-04 | 沈阳农业大学 | Method and system for quantitatively estimating carbon reserves of ecosystem of primary coastal wetland |
| CN112836610A (en) * | 2021-01-26 | 2021-05-25 | 平衡机器科技(深圳)有限公司 | Land use change and carbon reserve quantitative estimation method based on remote sensing data |
Non-Patent Citations (4)
| Title |
|---|
| KENING 等: "Land use change induced by the implementation of ecological restoration Programs increases future terrestrial ecosystem carbon sequestration in red soil hilly region of China", 《SCIENCEDIRECT》 * |
| 刘远林等: "内蒙古白音高勒有色金属矿区基于DIMINE的三维地质建模及矿产储量估算", 《矿产与地质》 * |
| 张敬宗等: "基于平行断面法的矿产资源储量可视化动态估算及应用", 《甘肃科技》 * |
| 李蓓蓓 等: "农业土地利用与碳排放研究综述", 《中国农机化学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114969665B (en) | 2023-02-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Shao et al. | Petrology and geochemistry of the high-sulphur coals from the Upper Permian carbonate coal measures in the Heshan Coalfield, southern China | |
| CN104612675B (en) | A kind of carbonate formation Lithology while drilling method for quickly identifying | |
| Donovan et al. | Sequence stratigraphy of the Eagle Ford (Boquillas) Formation in the subsurface of South Texas and outcrops of West Texas | |
| Morgans-Bell et al. | Integrated stratigraphy of the Kimmeridge Clay Formation (Upper Jurassic) based on exposures and boreholes in south Dorset, UK | |
| CN106257310B (en) | Sedimentary basin oil gas strong reducing action area Prospecting Sandstone-type Uranium Deposits method for establishing model | |
| Fuenlabrada et al. | Geochemistry of the Ediacaran–Early Cambrian transition in Central Iberia: Tectonic setting and isotopic sources | |
| Hazen et al. | Mineral evolution: Episodic metallogenesis, the supercontinent cycle, and the coevolving geosphere and biosphere | |
| Olander et al. | Surface soil changes following selective logging in an eastern Amazon forest | |
| CN106599377A (en) | Method for quantitative division of coal body structure based on logging data and coal body structure parameters | |
| Oerter et al. | Early to middle Miocene climate in the Atacama Desert of northern Chile | |
| Ji et al. | A novel method for estimating methane emissions from underground coal mines: The Yanma coal mine, China | |
| Marshall et al. | Reinterpreting the age of the uppermost ‘Old Red Sandstone'and Early Carboniferous in Scotland | |
| Pereira et al. | Provenance of upper Triassic sandstone, southwest Iberia (Alentejo and Algarve basins): tracing variability in the sources | |
| CN114969665B (en) | Method and system for estimating carbon emission and carbon reserve of mineral resource base | |
| Zhao et al. | Lithic technological responses to environmental change during the penultimate glacial cycle (MIS 7–6) at the Yangshang site, western Chinese Loess Plateau | |
| CN116794746A (en) | Geophysical quantitative prediction method for sea shale lithofacies | |
| CN111045105B (en) | Method for predicting geochemical connate zonation model for deep exploration of volcanic rock type uranium ores | |
| Asiedu et al. | Palaeoclimatic control on the composition of Palaeozoic shales from southern Ghana, West Africa | |
| Zhang et al. | Multivariate relationships and spatial distribution of geochemical features of soils in Jamaica | |
| Jiménez-Madrid et al. | Comparative analysis of intrinsic groundwater vulnerability assessment methods for carbonate aquifers | |
| CN110728074A (en) | Method for continuously explaining lithology of mixed-volume fine-grained rock and modeling method of model thereof | |
| Edwards et al. | The importance of unsaturated zone biogeochemical processes in determining groundwater composition, southeastern Australia | |
| Maithel | Characterization of Cross-bed Depositional Processes in the Coconino Sandstone | |
| Allen | Structurally controlled Cu-Zn-(Co-Pb) mineralization in the Neoproterozoic Ombombo Subgroup, Kaokoland, Namibia | |
| Lee | Metallogeny and characterization of late cretaceous superimposed porphyry Cu-Au-Mo and epithermal Au-Ag systems in the Dawson Range, Yukon, Canada: case study on the Klaza deposit |
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 |