CN114881298B - Optimization method and device for final boundary of opencut coal mine, medium and electronic equipment - Google Patents
Optimization method and device for final boundary of opencut coal mine, medium and electronic equipment Download PDFInfo
- Publication number
- CN114881298B CN114881298B CN202210414436.XA CN202210414436A CN114881298B CN 114881298 B CN114881298 B CN 114881298B CN 202210414436 A CN202210414436 A CN 202210414436A CN 114881298 B CN114881298 B CN 114881298B
- Authority
- CN
- China
- Prior art keywords
- boundary
- area
- land
- cost
- iteration number
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000003245 coal Substances 0.000 title claims abstract description 183
- 238000000034 method Methods 0.000 title claims abstract description 170
- 238000005457 optimization Methods 0.000 title claims abstract description 86
- 230000007613 environmental effect Effects 0.000 claims abstract description 126
- 238000004519 manufacturing process Methods 0.000 claims abstract description 85
- 230000008030 elimination Effects 0.000 claims abstract description 20
- 238000003379 elimination reaction Methods 0.000 claims abstract description 20
- 239000002689 soil Substances 0.000 claims description 87
- 239000011435 rock Substances 0.000 claims description 73
- 238000005065 mining Methods 0.000 claims description 51
- 238000005406 washing Methods 0.000 claims description 48
- 229910052799 carbon Inorganic materials 0.000 claims description 22
- 238000012423 maintenance Methods 0.000 claims description 18
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 17
- 238000013139 quantization Methods 0.000 claims description 14
- 238000003860 storage Methods 0.000 claims description 13
- 238000009412 basement excavation Methods 0.000 claims description 12
- 238000003892 spreading Methods 0.000 claims description 10
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 9
- 235000015097 nutrients Nutrition 0.000 claims description 9
- 229910052760 oxygen Inorganic materials 0.000 claims description 9
- 239000001301 oxygen Substances 0.000 claims description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 239000002028 Biomass Substances 0.000 claims description 8
- 238000004887 air purification Methods 0.000 claims description 8
- 239000005431 greenhouse gas Substances 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000011002 quantification Methods 0.000 claims description 5
- UIIMBOGNXHQVGW-UHFFFAOYSA-M Sodium bicarbonate Chemical compound [Na+].OC([O-])=O UIIMBOGNXHQVGW-UHFFFAOYSA-M 0.000 claims description 4
- 239000000758 substrate Substances 0.000 claims 1
- 230000008901 benefit Effects 0.000 description 13
- 238000004422 calculation algorithm Methods 0.000 description 11
- 238000004891 communication Methods 0.000 description 4
- 241000668854 Howardia biclavis Species 0.000 description 2
- 230000008094 contradictory effect Effects 0.000 description 2
- 238000005265 energy consumption Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007717 exclusion Effects 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002459 sustained effect Effects 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
-
- 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2219/00—Indexing scheme relating to application aspects of data processing equipment or methods
- G06F2219/10—Environmental application, e.g. waste reduction, pollution control, compliance with environmental legislation
-
- 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/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Marketing (AREA)
- Mathematical Physics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Data Mining & Analysis (AREA)
- Agronomy & Crop Science (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Mathematical Optimization (AREA)
- Health & Medical Sciences (AREA)
- Mathematical Analysis (AREA)
- General Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Primary Health Care (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The application provides an optimization method and device for the final boundary of an opencast coal mine, a medium and electronic equipment, wherein the method comprises the following steps: taking the iteration number to be 0, and optimizing the boundary by using a cone elimination method according to the unit production cost so as to determine the current boundary of the opencut coal mine; acquiring an iteration number of an optimization process; if the iteration number is not 0, judging whether the current boundary and the boundary obtained by the last iteration optimization meet the preset convergence condition, and if the current boundary and the boundary parameter are stored and output, and the optimization is finished; if the iteration number is 0, adding the environmental cost of each process unit production quantity to the production cost of each process unit production quantity to obtain the comprehensive cost, adding 1 to the iteration number, and re-optimizing the boundary by using a cone elimination method according to the comprehensive cost until the current boundary and the boundary obtained by the previous iteration optimization meet the preset convergence condition. The method of the application realizes that the environmental cost is included in the final boundary optimization like the production cost.
Description
Technical Field
The application relates to the technical field of open pit mining, in particular to an optimization method and device for the final boundary of an open pit coal mine, a medium and electronic equipment.
Background
In the related art, before the final boundary is determined, the land damage area is unknown, and the environmental cost cannot be calculated according to the land damage area in advance, so the environmental cost cannot be included in the final boundary optimization like the production cost.
Disclosure of Invention
In view of the above, the application provides an optimization method and device, medium and electronic equipment for the final boundary of an opencast coal mine, which realize the quantification of the environmental cost of ecological environment impact and bring the environmental cost into the final boundary optimization like the production cost.
According to one aspect of the application, there is provided a method for optimizing the final boundary of an opencast coal mine, comprising:
Determining unit production cost of unit production of each process in a plurality of processes according to the mining technical conditions of the opencast coal mine, wherein the plurality of processes comprise: raw coal mining, rock stripping, quaternary layer stripping and coal washing;
Determining the environmental cost of the unit area of the open pit coal mine according to the type of the ecological system;
taking the iteration number to be 0, and optimizing the boundary by using a cone elimination method according to the unit production cost so as to determine the current boundary of the opencut coal mine;
calculating the land damage area of each procedure in the current boundary by using the environmental impact quantization model and each land damage area allocation model;
calculating the environmental cost of unit production of each process according to the land damage area of each process;
Acquiring an iteration number of the optimization process, and judging whether the iteration number is 0;
if the iteration number is not 0, judging whether the current boundary and the boundary obtained by the last iteration optimization meet the preset convergence condition, and if the current boundary and the boundary parameter are stored and output, and the optimization is finished;
If the iteration number is 0, adding the environmental cost of each process unit production amount to the production cost of each process unit production amount to obtain the comprehensive cost, adding 1 to the iteration number, and re-optimizing the boundary by using a cone elimination method according to the comprehensive cost until the current boundary and the boundary obtained by the previous iteration optimization meet the preset convergence condition;
if the preset convergence condition is not met, judging whether the iteration number is greater than or equal to a preset iteration number threshold, and if the iteration number is greater than or equal to the preset iteration number threshold, storing and outputting current boundaries and boundary parameters, and finishing optimization;
if the iteration number is smaller than the preset iteration number threshold, adding the environmental cost of each process unit production quantity to the production cost of each process unit production quantity to obtain the comprehensive cost, adding 1 to the iteration number, and re-optimizing the boundary by using the cone elimination method according to the comprehensive cost until the current boundary and the boundary obtained by the previous iteration optimization meet the preset convergence condition.
Optionally, the step of calculating the land damage area of each process in the current boundary by using the environmental impact quantization model and each land damage area allocation model specifically comprises the following steps:
Determining a plurality of total land damage areas of the current boundary by using an environmental impact quantification model, wherein the plurality of total land damage areas comprise: the total area of the soil excavation loss, the total area of the soil discharge site, the total area of ground facilities, the total area of the gangue pile and the total area of the surface soil site;
spreading the total area of each land damage to each procedure, and respectively calculating each land damage area corresponding to each procedure after spreading by using each land damage area spreading model;
And calculating the land damage area of each procedure according to each land damage area.
Optionally, each total area of the land damage is allocated to each procedure, and each land damage area corresponding to each procedure after allocation is calculated by using each land damage area allocation model, which specifically comprises the following steps:
The ratio of the excavated area of the land at the boundary of the raw coal exploitation is Wherein V m is the volume of raw coal, V c is the first boundary volume, A c is the excavation area of the boundary land;
The ratio of the rock stripping to the excavated area of the land at the boundary is Wherein V y is rock volume;
The ratio of the digging area of the land at the boundary of the four-layer stripping is Wherein V 4 is the quaternary layer volume;
the occupation ratio of rock stripping in the land occupation area of the dumping site is Wherein A p is the occupied area of the soil discharge site;
the occupation ratio of the land occupation area of the dumping site by the stripping of the quaternary layer is
The ratio of the raw coal exploitation to the ground facility area isWherein A s is the ground facility area;
The ratio of rock stripping to the ground facility area is
The ratio of the area of the four-layer stripping to the ground facility area is
The ratio of the area of the coal washing on the ground facilities is
The occupation ratio of the surface soil field occupied area of raw coal exploitation isWherein h c is the average surface soil stripping thickness of the stope, h s is the average surface soil stripping thickness of the ground facility, V b is the surface soil field volume, and A b is the surface soil field occupied area;
the occupation ratio of rock stripping in the surface soil field is as follows:
wherein hp is the average stripping thickness of the surface soil of the dump;
the occupation ratio of the land occupation area of the surface soil field of the quaternary layer stripping is as follows:
The occupation ratio of the occupation area of the coal selective washing surface soil field is Wherein A g is the occupied area of the gangue pile, and h g is the surface soil average stripping thickness of the gangue pile.
Optionally, the step of calculating the land damage area of each procedure according to each land damage area specifically includes:
the land damage area of the raw coal exploitation is A m＝Acm+Asm+Abm;
The damage area of the rock stripped land is A y＝Apy+Acy+Asy+Aby;
the land damage area of the quaternary layer stripping is A 4＝Ap4+Ac4+As4+Ab4;
The damage area of the soil for coal washing is A x＝Ag+Asx+Abx.
Optionally, the step of calculating the environmental cost per unit throughput of each process according to the land damage area of each process specifically includes:
The first environmental cost of raw coal mining is:
Wherein Z yield is the biomass production value of the destroyed land, Z csoil is the soil carbon fixation value of the land ecosystem, C rec is the reclamation and maintenance cost, Z cplant is the plant carbon fixation value of the destroyed land ecosystem, Z O2 is the oxygen release value of the destroyed land ecosystem, Z air is the air purification value of the destroyed land ecosystem, Z soil is the soil maintenance value of the destroyed land ecosystem, Z H2O is the water conservation value of the destroyed land ecosystem, Z neut is the nutrient circulation value of the destroyed land ecosystem, F is the time length, F t is the time coefficient, C gm is the emission environment cost of raw coal mining, and Q m is the raw coal weight;
the first environmental cost of rock stripping is:
wherein, C gy is the environmental cost of the emission of rock stripping;
The first environmental cost of the quaternary layer stripping is:
wherein, C g4 is the environmental cost of the emission of the quaternary layer stripping;
The first environmental cost of coal washing is:
Wherein, C gx is the emission environmental cost of coal washing.
Optionally, the boundary parameters include at least one of: raw coal amount, rock amount, quaternary layer amount, multiple land damage areas, greenhouse gas emission amount and environmental cost of each process in the current world.
Optionally, the environmental cost per unit area includes at least one of: biomass production value, soil carbon fixation value, plant carbon fixation value, oxygen release value, air purification value, soil maintenance value, water source conservation value, nutrient circulation value and reclamation maintenance cost.
According to another aspect of the present application, there is provided an optimization apparatus for an open pit coal mine final boundary, comprising:
a first determining module, configured to determine a unit production cost of a unit production amount of each process in a plurality of processes according to a mining technical condition of an opencast coal mine, where the plurality of processes include: raw coal mining, rock stripping, quaternary layer stripping and coal washing;
The second determining module is used for determining the environmental cost of the unit area of the open pit coal mine according to the type of the ecological system;
The third determining module is used for taking the iteration number to be 0, optimizing the boundary by using a cone elimination method according to the unit production cost, and determining the current boundary of the opencast coal mine;
the first calculation module is used for calculating the land damage area of each procedure in the current boundary by using the environmental impact quantization model and each land damage area allocation model;
the second calculation module is used for calculating the environmental cost of unit production of each process according to the land damage area of each process;
the acquisition module is used for acquiring the iteration number of the optimization process;
the judging module is used for judging whether the iteration number is 0;
The judging module is further used for judging whether the current boundary and the boundary obtained by the last iterative optimization meet the preset convergence condition or not if the iterative number is not 0;
The output module is used for outputting the current boundary and boundary parameters if the preset convergence condition is met, and the optimization is finished;
the optimization module is used for adding the environmental cost of each process unit production amount to the production cost of each process unit production amount to obtain the comprehensive cost if the iteration number is 0, adding 1 to the iteration number, and re-optimizing the boundary by using the cone elimination method according to the comprehensive cost until the current boundary and the boundary obtained by the previous iteration optimization meet the preset convergence condition;
The judging module is further used for judging whether the iteration number is larger than or equal to a preset iteration number threshold value or not if the preset convergence condition is not met;
the output module is further used for outputting the current boundary and boundary parameters if the iteration number is greater than or equal to a preset iteration number threshold value, and optimizing is finished;
and the optimization module is also used for adding the environmental cost of each process unit production quantity to the production cost of each process unit production quantity to obtain the comprehensive cost if the iteration number is smaller than the preset iteration number threshold value, adding 1 to the iteration number, and re-optimizing the boundary by utilizing the cone elimination method according to the comprehensive cost until the current boundary and the boundary obtained by the previous iteration optimization meet the preset convergence condition.
According to a further aspect of the present application there is provided a storage medium having stored thereon a computer program which when executed by a processor implements the method of optimizing the final boundaries of an opencut mine as described above.
According to a further aspect of the application, there is provided an electronic device comprising a storage medium, a processor and a computer program stored on the storage medium and executable on the processor, the processor implementing the method for optimizing the final boundaries of an opencast mine as described above when executing the computer program.
By means of the technical scheme, the impact of the opencast coal mine on the ecological environment is quantified through the biochemical boundary optimization iterative algorithm in the environment cost, the environment cost of the ecological environment impact is quantified, the environment cost is included into the final boundary optimization like the production cost, and therefore comprehensive benefits are obtained, namely, after the environment cost is subtracted from the pure economic benefits, the highest boundary is used as the final boundary. The optimized final boundary achieves the best balance between improving the economic benefit of open pit mining as much as possible and reducing the impact on the ecological environment as much as possible, and the low-carbon ecological optimization of the final boundary of the open pit mine is realized.
The foregoing description is only an overview of the present application, and is intended to be implemented in accordance with the teachings of the present application in order that the same may be more clearly understood and to make the same and other objects, features and advantages of the present application more readily apparent.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this specification, illustrate embodiments of the application and together with the description serve to explain the application and do not constitute a limitation on the application. In the drawings:
fig. 1 shows a schematic flow chart of an optimization method for an open pit coal mine final boundary provided by an embodiment of the application;
Fig. 2 shows a block diagram of an optimization device for the final boundary of an opencast coal mine according to an embodiment of the present application.
Detailed Description
The application will be described in detail hereinafter with reference to the drawings in conjunction with embodiments. It should be noted that, without conflict, the embodiments of the present application and features of the embodiments may be combined with each other.
In this embodiment, an optimization method for the final boundary of an opencast coal mine is provided, as shown in fig. 1, and the method includes:
step 101, determining unit production cost of unit production of each process in a plurality of processes according to the exploitation technical conditions of an opencast coal mine;
Step 102, determining the environmental cost of the unit area of the open pit coal mine according to the type of the ecological system;
Step 103, taking the iteration number to be 0, and optimizing the boundary by using a cone elimination method according to the unit production cost so as to determine the current boundary of the opencast coal mine;
104, calculating the land damage area of each procedure in the current boundary by using the environmental impact quantization model and each land damage area allocation model;
step 105, calculating the environmental cost of unit production of each process according to the land damage area of each process;
step 106, obtaining the iteration number of the optimization process;
Step 107, judging whether the iteration number is 0, if so, proceeding to step 111, otherwise proceeding to step 108;
Step 108, judging whether the current boundary and the boundary obtained by the last iterative optimization meet the preset convergence condition, if yes, entering step 109, and if not, entering step 110;
Step 109, storing and outputting the current boundary and boundary parameters, and ending the optimization;
Step 110, judging whether the iteration number is greater than or equal to a preset iteration number threshold, if so, entering step 109, and if not, entering step 111;
And step 111, adding the environmental cost of each process unit production to the production cost of each process unit production to obtain the comprehensive cost, adding 1 to the iteration number, and re-optimizing the boundary by using a cone elimination method according to the comprehensive cost until the current boundary and the boundary obtained by the previous iteration optimization meet the preset convergence condition.
In this embodiment, the unit production cost per unit production amount of each process in a plurality of processes, namely four major processes of raw coal mining, rock stripping, quaternary layer stripping and coal washout in the mining process of the open pit coal mine, is determined according to the mining technical conditions of the newly built or extended open pit coal mine. Further, according to the type of the ecological system, the environmental cost per unit area of the open pit coal mine is determined. And taking the iteration number to be 0, optimizing the boundary by using a cone elimination method according to the unit production cost of each procedure, and further determining the current boundary of the opencast coal mine.
Further, an environmental impact quantization model and various land damage area allocation models are applied, and the land damage area of each procedure in the current boundary is calculated. And further, according to the land damage area of each process, calculating the environmental cost of each process unit production.
Further, in the iterative optimization process, judging whether the iterative ordinal number is 0, if the iterative ordinal number is not 0, judging whether the current boundary and the boundary obtained by the last iterative optimization meet the preset convergence condition, if the preset convergence condition is met, describing that the algorithm converges, and determining that the current boundary is a final boundary, and storing and outputting current boundary and boundary parameters, wherein the related parameters are parameters such as raw coal quantity, rock quantity, four-layer quantity, land damage area of each process, greenhouse gas emission quantity, environmental cost of each process and the like, and the iteration is ended at the moment. If the preset convergence condition is not met, judging whether the iteration number is larger than or equal to a preset ordinal threshold, if the iteration number is larger than or equal to the ordinal threshold, indicating that the algorithm reaches the maximum iteration number, storing and outputting the current boundary and relevant parameters thereof, and ending the iteration. If the iteration number is smaller than the preset ordinal threshold, the iteration number is adjusted, namely the iteration number i=i+1 is set, and the method returns to calculating a plurality of unit costs corresponding to the iteration number at the moment to perform the next boundary optimization.
When determining whether the obtained iteration number is 0, if the iteration number is 0, it is described that the iteration has not yet started at this time, the iteration number is adjusted, the iteration number i=i+1 is set, and a plurality of unit costs corresponding to the iteration number are calculated to perform boundary optimization. It will be appreciated that when the iteration number is 0, only the production cost need be considered, and no environmental cost need be considered. Therefore, at an iteration number of 0, the unit cost of each of the plurality of processes is equal to the production cost of each of the plurality of processes.
It should be noted that, the preset convergence condition is comparing the current boundary volume with the boundary volume of the boundary obtained by the previous iterative optimization, and judging whether the percentage of the two phase differences is smaller than dV (%), if so, the algorithm converges, wherein dV (%) is the preset convergence standard.
According to the embodiment of the invention, the impact of the opencast coal mine on the ecological environment is quantized by using a biochemical boundary optimization iterative algorithm in the environment cost, so that the environment cost of the ecological environment impact is quantized, the environment cost is included into the final boundary optimization like the production cost, and the comprehensive benefit is obtained, namely, the highest boundary is taken as the final boundary after the environment cost is subtracted from the pure economic benefit. The optimized final boundary achieves the best balance between improving the economic benefit of open pit mining as much as possible and reducing the impact on the ecological environment as much as possible, and the low-carbon ecological optimization of the final boundary of the open pit mine is realized.
In the embodiment of the application, the step of calculating the land damage area of each procedure in the current boundary by using the environmental impact quantization model and each land damage area allocation model specifically comprises the following steps: determining a plurality of total land damage areas of the current boundary by using an environmental impact quantification model, wherein the plurality of total land damage areas comprise: the total area of the soil excavation loss, the total area of the soil discharge site, the total area of ground facilities, the total area of the gangue pile and the total area of the surface soil site; spreading the total area of each land damage to each procedure, and respectively calculating each land damage area corresponding to each procedure after spreading by using each land damage area spreading model; and calculating the land damage area of each procedure according to each land damage area.
In the technical scheme, an environmental impact quantization model is applied, and a plurality of land damage areas of the current boundary are calculated, specifically, the total area of the plurality of land damage areas comprises: the total area of the soil excavation, the total area of the soil discharge site, the total area of ground facilities, the total area of the gangue pile and the total area of the surface soil site. And then apportioning the above-mentioned all land damage area parameters to raw coal mining, rock stripping, quaternary layer stripping and coal selective washing, applying all land damage area apportionment models, apportioning all land damage areas to each working procedure, and further, synthesizing the apportionment of all land damage areas, so as to obtain the land damage areas respectively caused by raw coal mining, rock stripping, quaternary layer stripping and coal selective washing.
According to the method, based on the boundary optimization iterative algorithm of the internal biochemistry of the environmental cost, the land damage area can be calculated before the final boundary is determined, the environmental cost is calculated, the internal biochemistry of the ecological cost is further included in the overall evaluation of the mine ecology, and the external generation cost is internalized. Compared with the prior art, no matter the final boundary is designed according to the economic reasonable stripping ratio or the whole optimization is carried out, only various investment and operation cost in mine construction and production are always considered, and most of ecological function loss caused by damage to an ecological system in mine production is 'external generation cost', and the technical problem is not considered.
In the embodiment of the application, further, each total area of the land damage is allocated to each procedure, each land damage area corresponding to each procedure after allocation is calculated by using each land damage area allocation model, and the method specifically comprises the following steps:
The ratio of the excavated area of the land at the boundary of the raw coal exploitation is Wherein V m is the volume of raw coal, V c is the first boundary volume, A c is the excavation area of the boundary land;
The ratio of the rock stripping to the excavated area of the land at the boundary is Wherein V y is rock volume;
The ratio of the digging area of the land at the boundary of the four-layer stripping is Wherein V 4 is the quaternary layer volume;
the occupation ratio of rock stripping in the land occupation area of the dumping site is Wherein A p is the occupied area of the soil discharge site;
the occupation ratio of the land occupation area of the dumping site by the stripping of the quaternary layer is
The ratio of the raw coal exploitation to the ground facility area isWherein A s is the ground facility area;
The ratio of rock stripping to the ground facility area is
The ratio of the area of the four-layer stripping to the ground facility area is
The ratio of the area of the coal washing on the ground facilities is
The occupation ratio of the surface soil field occupied area of raw coal exploitation isWherein h c is the average surface soil stripping thickness of the stope, h s is the average surface soil stripping thickness of the ground facility, V b is the surface soil field volume, and A b is the surface soil field occupied area;
the occupation ratio of rock stripping in the surface soil field is as follows:
wherein hp is the average stripping thickness of the surface soil of the dump;
the occupation ratio of the land occupation area of the surface soil field of the quaternary layer stripping is as follows:
The occupation ratio of the occupation area of the coal selective washing surface soil field is Wherein A g is the occupied area of the gangue pile, and h g is the surface soil average stripping thickness of the gangue pile.
In the technical scheme, the land excavation area of the boundary of the opencast coal mine is jointly caused by raw coal mining, rock stripping and quaternary layer stripping, so that the land excavation area is apportioned to the raw coal mining, rock stripping and quaternary layer stripping according to the proportion of the raw coal volume, the rock volume and the quaternary layer volume in the boundary to the total volume of the boundary.
Specifically, the ratio of the excavation area of the raw coal exploitation at the boundary land is as follows:
The ratio of the rock stripping to the excavating area of the land at the boundary is as follows:
the ratio of the digging area of the land at the boundary of the four-layer stripping is as follows:
Wherein V m is the volume of raw coal, V c is the first boundary volume, A c is the excavation area of the boundary land; v y is rock volume; v 4 is the quaternary layer volume.
Further, the damage to the land by the dump is due to rock and quaternary stripping, so the total area of the dump is apportioned to rock stripping and quaternary stripping in proportion to the total stripped volume of rock volume, quaternary volume in the boundary.
Specifically, the occupation ratio of rock stripping in the soil discharge site is:
the occupation ratio of the land occupation area of the dumping site by the stripping of the quaternary layer is as follows:
wherein A p is the occupied area of the dumping site.
Further, in ground facilities such as mine-dedicated access, factory building, warehouse, office, power supply, and drainage, some serve for stripping, for example, dedicated transportation routes from stopes to dumping grounds; some serve raw coal exploitation and coal washing at the same time, for example, a special transportation road from a stope to a coal washery; some service coal washings, e.g., coal washery footprints; some serve whole mines, such as office facilities. It is difficult to determine the area of each facility and its service objects before the final mining plan and the total map layout are determined. Thus, the total area of the ground facility is apportioned to raw coal mining, rock stripping, quaternary stripping and coal washout, only roughly in terms of the volume of raw coal, rock volume and quaternary layer volume to total volume of the boundary.
Specifically, the ratio of the raw coal exploitation to the ground facility area is as follows:
The ratio of rock stripping to the ground facility area is:
The ratio of the area of the quaternary layer stripping on the ground facility is as follows:
The ratio of the area of the coal washing on the ground facilities is as follows:
wherein A s is the ground facility area.
Further, the damage of the gangue pile to the land is completely caused by coal washing, so that the occupied area of the gangue pile is completely classified into the land damage area of the coal washing.
According to the mode, according to different working procedures corresponding to different land damage areas, each land damage area is allocated to raw coal exploitation, rock stripping, quaternary layer stripping and coal washing, so that the land damage areas caused by the four working procedures are obtained, the calculation accuracy of the land damage areas is ensured, and the accuracy of a final boundary optimization result is improved.
In the embodiment of the present application, further, according to each land damage area, the step of calculating the land damage area of each procedure specifically includes:
the land damage area of the raw coal exploitation is A m＝Acm+Asm+Abm;
The damage area of the rock stripped land is A y＝Apy+Acy+Asy+Aby;
the land damage area of the quaternary layer stripping is A 4＝Ap4+Ac4+As4+Ab4;
The damage area of the soil for coal washing is A x＝Ag+Asx+Abx.
In the technical scheme, the land damage areas respectively caused by raw coal exploitation, rock stripping, four-layer stripping and coal washing can be obtained by combining the above-mentioned apportionment of each land damage area.
In the embodiment of the present application, further, the step of calculating the environmental cost per unit throughput of each process according to the land damage area of each process specifically includes:
The first environmental cost of raw coal mining is:
Wherein Z yield is the biomass production value of the destroyed land, Z csoil is the soil carbon fixation value of the land ecosystem, C rec is the reclamation and maintenance cost, Z cplant is the plant carbon fixation value of the destroyed land ecosystem, Z O2 is the oxygen release value of the destroyed land ecosystem, Z air is the air purification value of the destroyed land ecosystem, Z soil is the soil maintenance value of the destroyed land ecosystem, Z H2O is the water conservation value of the destroyed land ecosystem, Z neut is the nutrient circulation value of the destroyed land ecosystem, F is the time length, F t is the time coefficient, C gm is the emission environment cost of raw coal mining, and Q m is the raw coal weight;
the first environmental cost of rock stripping is:
wherein, C gy is the environmental cost of the emission of rock stripping;
The first environmental cost of the quaternary layer stripping is:
wherein, C g4 is the environmental cost of the emission of the quaternary layer stripping;
The first environmental cost of coal washing is:
Wherein, C gx is the emission environmental cost of coal washing.
In this technical scheme, the environmental costs caused by the raw coal mining include the land damage environmental cost and the greenhouse gas emission environmental cost of the raw coal mining, and therefore, the unit environmental cost of the raw coal mining, that is, the first environmental cost is:
Similarly, the unit environmental costs of rock stripping, quaternary layer stripping and coal selective washing are respectively as follows:
Wherein, C em is the unit environmental cost of raw coal exploitation, and the unit is the unit per t; c ey is the unit environmental cost of rock stripping, the unit is the unit environmental cost of quaternary layer stripping, the unit is the unit environmental cost of coal selective washing, and the unit is the unit environmental cost of quaternary layer stripping; c gm is the emission environmental cost of raw coal exploitation, and the unit is this/t; c gy is the unit discharge environmental cost of rock stripping, the unit is the unit discharge environmental cost of quaternary layer stripping, the unit is the unit discharge environmental cost of coal selective washing, and the unit is the unit discharge environmental cost of coal selective washing; z yield is the biomass production value of the destroyed land, namely the characteristic land price, the unit is the soil carbon fixing value of the destroyed land ecological system, the unit is the soil/hm 2;Crec is the reclamation and maintenance cost, the unit is the plant carbon fixing value of the destroyed land ecological system, the unit is the soil/(hm 2·a);ZO2) is the oxygen releasing value of the destroyed land ecological system, the unit is the soil purifying value of the destroyed land ecological system/(hm 2·a);Zair), the unit is the soil holding value of the destroyed land ecological system/(hm 2·a);Zsoil), the unit is the water conservation value of the destroyed land ecological system/(hm 2·a);ZH2O), the unit is the nutrient circulating value of the destroyed land ecological system/(hm 2. A); F is the length of time from beginning to end of mining and restoring the ecological function of the land, which is estimated according to reasonable annual production capacity of raw coal, F t is a time coefficient, and Q m is the weight of raw coal, and the unit is 10 4 t.
It should be noted that, in the environmental cost, the losses of plant carbon fixation, oxygen release, air purification, soil maintenance, water conservation and nutrient circulation values of the land ecosystem are sustained, and the losses of these values occur each year over a time span from the start of land damage to the completion of reclamation and restoration of its ecological function, and the total amount of these values is not only related to the area of land damage but also to the length of time for which the damaged state is maintained, i.e., depending on the damaged time and the reclamation time. Since there is no specific mining plan and reclamation plan before the final design of the boundary is determined, the duration of each land reclamation unit in the reclamation state cannot be estimated, and therefore, in the optimization process, the average duration of the land in the reclamation state, namely Ff t years, is estimated using the mining life F of the boundary, including reclamation after the end of the mining and the maintenance time, and a time coefficient F t. For the life F, a reasonable raw coal annual production capacity is estimated in terms of shearable quantities in the boundary without production planning, and the life F of the boundary is calculated therefrom to be equal to the life plus reclamation and maintenance time, the latter typically being 3 years to 5 years. For f t, considering that a large area of land is damaged in a mine infrastructure, and it is possible to reclaim some land damage units after a long time of production, so that most of the land damage area is in a damage state for a long time, the value range of the time coefficient f t is set to 0.5-0.8.
Through the method, the unit environmental cost of raw coal exploitation, rock stripping, quaternary layer stripping and coal washing is calculated respectively, so that the ecological cost is internally biochemically integrated and is incorporated into the overall evaluation of mine ecology, and the economic and ecological total benefits of the final boundaries after optimization are maximized.
In the practical application scene, the coal mining inevitably causes impact on the ecological environment, and the larger the mining scale is, the larger the impact amount is. When adopting the open pit mining, the final boundary represents the overall mining scale of the open pit coal mine. Reducing the final boundary can correspondingly reduce the impact on the ecological environment, but the economic loss is brought by reducing the final boundary to be lower than the optimal economic boundary, namely the final boundary with the largest economic benefit. Thus, when considering the ecological problem in the final boundary optimization, the optimization objective is twofold, that is, it is necessary to improve the economic efficiency of the surface mining as much as possible and reduce the impact on the ecological environment as much as possible. However, in the prior art, two targets are contradictory or contradictory, and the application provides a border optimization method of an open pit coal mine in biochemical environment cost, so that the problem is successfully solved.
Specifically, the basic principle of the boundary optimization method of the biochemical open-pit coal mine in the environmental cost is to quantize the impact of the open-pit coal mine on the ecological environment, further quantize the environmental cost of the ecological environment impact, and bring the environmental cost into the boundary optimization like the production cost, thereby obtaining the comprehensive benefit, namely the highest boundary after subtracting the environmental cost from the pure economic benefit, and taking the highest boundary as the optimal boundary. This boundary achieves an optimal balance between maximizing the economic benefits of surface mining and minimizing impact on the ecological environment.
However, in the process of internally generating the environmental cost in the final boundary optimization, the environmental cost of the land ecosystem damage is related to the land damage area, which is unknown before the final boundary is determined, and the environmental cost thereof cannot be calculated. Therefore, the application provides a biochemical boundary optimization iterative algorithm in environmental cost, and the damage to a land ecological system and the emission of greenhouse gases can be greatly reduced by applying the algorithm, so that the low-carbon ecological optimization of the final boundary of the opencast coal mine is realized.
Specifically, let C m、Cy、C4 and C x denote unit production costs of raw coal mining, rock stripping, quaternary layer stripping, and coal washing, respectively; c m,i、Cy,i、C4,i and C x,i respectively represent unit costs of raw coal mining, rock stripping, quaternary stripping and coal washing after adding environmental cost in the ith iteration, wherein i represents the iteration number. The final boundary optimization algorithm for biochemistry within the environmental cost is as follows:
Step 1: according to the type of the ecological system of the mining area, an environmental cost quantification model is applied, and the environmental cost of the unit land ecological system damaged area is calculated, namely, a biomass production value Z yield, a soil carbon fixation value Z csoil, a plant carbon fixation value Z cplant, an oxygen release value Z O2, an air purification value Z air, a soil maintenance value Z soil, a water source conservation value Z H2O, a nutrient circulation value Z neut and a reclamation maintenance cost C rec. Further, according to the unit energy consumption and the explosive unit energy consumption of raw coal exploitation, rock stripping, quaternary strata stripping and coal washing, the environmental impact quantization model is applied to calculate the unit greenhouse gas emission amounts e m、ey、e4 and e x of the four working procedures respectively, and then the environmental cost quantization model is applied to calculate the unit emission environmental costs C gm、Cgy、Cg4 and C gx of the four working procedures respectively. A convergence criterion dV (%) and a maximum number of iterations n max are set.
Step 2: let i=0; c m,i＝Cm、Cy,i＝Cy、C4,i＝C4、Cx,i＝Cx, namely only the production cost is considered, and the environmental cost is not considered.
Step 3: with C m,i、Cy,i、C4,i and C x,i as cost parameters, the "final-boundary-optimized cone exclusion algorithm" is applied to optimize the boundary, which is the current boundary. The total volume, raw coal volume, rock volume, quaternary layer volume and raw coal weight of this boundary were V c、Vm、Vy、V4 and Q m, respectively, with volume units of 10 4m3 and weight units of 10 4 t.
Step 4: and calculating the total land excavation area A c, the total land occupation area A p, the total ground facility area A s, the total waste pile occupation area A g and the total surface soil field occupation area A b of the boundary by using an environment impact quantization model.
Step 5: and (3) spreading all the land damage areas obtained in the step (4) to raw coal mining, rock stripping, quaternary layer stripping and coal washing to obtain the land damage areas respectively caused by the four working procedures.
Specifically, the apportionment method is as follows:
The land excavated area A c of the boundary is commonly caused by raw coal mining, rock stripping and quaternary stripping, so that A c is apportioned to raw coal mining, rock stripping and quaternary stripping, and is marked as A cm、Acy and A c4, according to the proportion of the raw coal volume V m, the rock volume V y and the quaternary layer volume V 4 in the boundary to the total boundary volume V c:
Further, the damage to the land by the dump is due to rock and quaternary stripping, so the total area of the dump, a p, is apportioned to rock and quaternary stripping, noted as a py and a p4, in the proportion of the rock volume V y, quaternary volume V 4, to the total stripping volume in the boundary:
Further, ground facilities such as facilities of mine-dedicated roads and plants, warehouse, office, power supply, drainage, etc.; some are dedicated haul roads that serve for stripping, such as stopes to dumping grounds; some are simultaneously used for raw coal exploitation and coal washing, such as a special transportation road from a stope to a coal washery; some are used for coal washing, such as occupation of land in a coal washery; some serve whole mines, such as office facilities. It is difficult to determine the area of each facility and its service objects before the final mining plan and the total map layout are determined. Thus, the total surface facility area A s is apportioned to raw coal mining, rock stripping, quaternary stripping and coal washout, represented by A sm、Asy、As4 and A sx, respectively, only roughly in the proportions of raw coal volume V m, rock volume V y quaternary layer volume V 4 to total boundary volume V c:
Further, the damage of the gangue pile to the land is completely caused by coal washing, so that the occupied area A g of the gangue pile is completely classified as the land damage area of the coal washing.
Further, the topsoil originates from the quarry, the dump, the pile and the topsoil stripping on the ground for ground facilities. Assuming that the surface soil field occupation area is proportional to the soil amount, the surface soil field occupation area A b is apportioned to the raw coal exploitation, rock stripping, quaternary stripping and coal washing according to the proportion of the surface soil stripping volume to the total surface soil volume V b of the areas of raw coal exploitation, rock stripping, quaternary stripping and coal washing, and the surface soil field occupation area A b is respectively represented by A bm、Aby、Ab4 and A bx:
/>
Wherein h c、hp、hs、hg is the average stripping thickness of surface soil of a stope, a dumping site, a ground facility area and a gangue storage yard, and the unit is m.
By integrating the above-mentioned apportionment of various land damage areas, the land damage areas respectively caused by raw coal exploitation, rock stripping, four-layer stripping and coal washing and selecting can be obtained and respectively marked as A m、Ay、A4 and A x:
Am＝Acm+Asm+Abm；
Ay＝Apy+Acy+Asy+Aby；
A4＝Ap4+Ac4+As4+Ab4；
Ax＝Ag+Asx+Abx。
step 6: the unit environmental cost of raw coal mining, rock stripping, quaternary stripping and coal washout is calculated.
The environmental costs caused by the raw coal mining include the land damage environmental cost and the greenhouse gas emission environmental cost of the raw coal mining, so the unit environmental cost of the raw coal mining is:
Cem is the unit environmental cost of raw coal exploitation, and is/t; z yield is the biomass production value of the destroyed land, namely the land taking price, the unit is that of the soil carbon fixation value of the land ecological system is that of the soil/hm 2;Zcsoil, the unit is that of the reclamation and maintenance cost is that of the soil/hm 2;Crec, the unit is that of the soil/hm 2;Zcplant、ZO2、Zair、Zsoil、ZH2O、Zneut, the unit is that of the plant carbon fixation, oxygen release, air purification, soil maintenance, water conservation and nutrient circulation value of the destroyed land ecological system, and the unit is that of the soil/(hm 2. A); f is the estimated time length from the beginning of exploitation to the end of exploitation and restoring the ecological function of the land according to reasonable annual production capacity of raw coal; f t is a time coefficient; c gm is the unit emission environmental cost of raw coal mining, and the unit is the unit per t.
Similarly, the unit environmental costs of rock stripping, quaternary layer stripping and coal washing are:
Wherein, C ey and C e4 are the unit environmental cost of rock and quaternary layer stripping respectively, the unit is the unit environmental cost of coal washing and selecting per m 3;Cex, and the unit is the unit environmental cost of coal washing and selecting per t; c gy and C g4 are the unit discharge environmental costs of rock and quaternary stripping, respectively, in units of this/m 3;Cgx is the unit discharge environmental cost of coal washout, in units of this/t.
Step 7: if i=0, go to step 10; otherwise, i >0, step 8 is entered.
Step 8: comparing the total volume of the current boundary and the boundary obtained by the previous iteration, if the percentage of the two phase differences is smaller than dV, converging the algorithm, wherein the current boundary is the optimal boundary, and outputting the boundary and related parameters thereof, wherein the related parameters comprise raw coal quantity, rock quantity, quaternary layer quantity, various land damage areas, greenhouse gas emission quantity, various environmental costs and the like, and the iteration is ended; otherwise, step 9 is entered.
Step 9: if i is not less than n max, outputting the current boundary and related parameters thereof, and ending the iteration; if i < n max, go to step 10.
Step 10: setting i=i+1, returning Cm,i＝Cm+Cem,Cy,i＝Cy+Cey,C4,i＝C4+Ce4,Cx,i＝Cx+Cex, to the step 3, adding the environmental cost into the unit production cost, and performing the next boundary optimization.
Further, as a specific implementation of the optimization method of the final boundary of the opencast coal mine, the embodiment of the application provides an optimization device of the final boundary of the opencast coal mine, as shown in fig. 2, the optimization device 200 of the final boundary of the opencast coal mine includes:
A first determining module 201, configured to determine a unit production cost per unit production amount of each process in a plurality of processes according to a mining technical condition of an opencast coal mine, where the plurality of processes includes: raw coal mining, rock stripping, quaternary layer stripping and coal washing;
A second determining module 202, configured to determine an environmental cost per unit area of the opencast coal mine according to the type of the ecosystem;
A third determining module 203, configured to take the iteration number to 0, and optimize the boundary according to the unit production cost by using a cone elimination method, so as to determine the current boundary of the opencut mine;
a first calculation module 204, configured to calculate a land damage area of each process in the current boundary using the environmental impact quantization model and each land damage area allocation model;
A second calculation module 205 for calculating an environmental cost per unit throughput of each process based on the land damage area of each process;
an obtaining module 206, configured to obtain an iteration number of the optimization process;
a judging module 207, configured to judge whether the iteration number is 0;
the judging module 207 is further configured to judge whether the current boundary and the boundary obtained by the last iterative optimization meet the preset convergence condition if the iterative number is not 0;
The output module 208 is configured to output the current boundary and the boundary parameter if the preset convergence condition is satisfied, and the optimization is ended;
The optimizing module 209 is configured to, if the iteration number is 0, add the environmental cost of each process unit throughput to the production cost of each process unit throughput to obtain a comprehensive cost, add 1 to the iteration number, and re-optimize the boundary by using the cone elimination method according to the comprehensive cost until the current boundary and the boundary obtained by the previous iteration optimization meet the preset convergence condition;
The judging module 207 is further configured to judge whether the iteration number is greater than or equal to a preset iteration number threshold if the preset convergence condition is not satisfied;
The output module 208 is further configured to output the current boundary and the boundary parameter if the iteration number is greater than or equal to the preset iteration number threshold, and the optimization is ended;
The optimization module 209 is further configured to add the environmental cost of each process unit throughput to the production cost of each process unit throughput if the iteration number is smaller than the preset iteration number threshold, obtain a comprehensive cost, add 1 to the iteration number, and re-optimize the boundary by using the cone elimination method according to the comprehensive cost until the current boundary and the boundary obtained by the previous iteration optimization meet the preset convergence condition.
It should be noted that, other corresponding descriptions of each functional module related to the optimization device for the final boundary of the opencast coal mine provided by the embodiment of the present application may refer to the corresponding description of fig. 1, and are not repeated herein.
Based on the method shown in fig. 1, correspondingly, the embodiment of the application also provides a storage medium, on which a computer program is stored, and the program is executed by a processor to realize the optimization method of the final boundary of the opencast coal mine shown in fig. 1.
Based on such understanding, the technical solution of the present application may be embodied in a software product or a hardware product, or a combination of hardware and software, where the software product may be stored in a nonvolatile storage medium (may be a CD-ROM, a U-disk, a mobile hard disk, etc.), and includes several instructions for causing an electronic device (may be a personal computer, a server, or a network device, etc.) to execute the method of each implementation scenario of the present application.
Based on the method shown in fig. 1 and the embodiment of the optimizing device of the final boundary of the opencast mine shown in fig. 2, in order to achieve the above purpose, the embodiment of the application further provides an electronic device, which may be a personal computer, a server, a network device, etc., specifically, the electronic device includes a storage medium and a processor; a storage medium storing a computer program; a processor for executing a computer program for implementing the optimization method of the final boundary of the opencast mine as shown in fig. 1.
Optionally, the electronic device may also include a user interface, a network interface, a camera, radio Frequency (RF) circuitry, sensors, audio circuitry, WI-FI modules, and the like. The user interface may include a Display screen (Display), an input unit such as a Keyboard (Keyboard), etc., and the optional user interface may also include a USB interface, a card reader interface, etc. The network interface may optionally include a standard wired interface, a wireless interface (e.g., bluetooth interface, WI-FI interface), etc.
It will be appreciated by those skilled in the art that the structure of the electronic device provided in this embodiment is not limited to the electronic device, and may include more or fewer components, or may be combined with certain components, or may be arranged with different components.
The storage medium may also include an operating system, a network communication module. An operating system is a program that manages and saves electronic device hardware and software resources, supporting the execution of information handling programs, as well as other software and/or programs. The network communication module is used for realizing communication among all the controls in the storage medium and communication with other hardware and software in the entity equipment.
From the above description of the embodiments, it will be apparent to those skilled in the art that the present application may be implemented by means of software plus necessary general hardware platforms, or may be implemented by hardware.
Those skilled in the art will appreciate that the drawing is merely a schematic illustration of one preferred implementation scenario and that elements or processes in the drawing are not necessarily required to practice the application. Those skilled in the art will appreciate that elements of an apparatus in an implementation may be distributed throughout the apparatus in an implementation as described in the implementation, or that corresponding variations may be located in one or more apparatuses other than the present implementation. The units of the implementation scenario may be combined into one unit, or may be further split into a plurality of sub-units.
The above-mentioned inventive sequence numbers are merely for description and do not represent advantages or disadvantages of the implementation scenario. The foregoing disclosure is merely illustrative of some embodiments of the application, and the application is not limited thereto, as modifications may be made by those skilled in the art without departing from the scope of the application.
Claims (6)
1. A method for optimizing a final boundary of an opencast coal mine, the method comprising:
Determining unit production cost of unit production of each process in a plurality of processes according to mining technical conditions of an opencast coal mine, wherein the plurality of processes comprise: raw coal mining, rock stripping, quaternary layer stripping and coal washing;
determining the environmental cost of the open pit coal mine in unit area according to the type of the ecological system;
taking the iteration number to be 0, and optimizing the boundary by using a cone elimination method according to the unit production cost so as to determine the current boundary of the opencast coal mine;
Calculating the land damage area of each procedure in the current boundary by using an environmental impact quantization model and each land damage area allocation model;
Calculating the environmental cost of unit production of each process according to the land damage area of each process;
Acquiring an iteration number of an optimization process, and judging whether the iteration number is 0;
if the iteration number is not 0, judging whether the current boundary and the boundary obtained by the last iteration optimization meet a preset convergence condition, and if the preset convergence condition is met, storing and outputting the current boundary and boundary parameters, and finishing optimization;
If the iteration number is 0, adding the environmental cost of each process unit production amount to the production cost of each process unit production amount to obtain the comprehensive cost, adding 1 to the iteration number, and re-optimizing the boundary by using a cone elimination method according to the comprehensive cost until the current boundary and the boundary obtained by the previous iteration optimization meet the preset convergence condition;
If the preset convergence condition is not met, judging whether the iteration number is larger than or equal to a preset iteration number threshold, and if the iteration number is larger than or equal to the preset iteration number threshold, storing and outputting the current boundary and boundary parameters, and finishing optimization;
if the iteration number is smaller than the preset iteration number threshold, adding the environmental cost of each process unit production amount into the production cost of each process unit production amount to obtain the comprehensive cost, adding 1 to the iteration number, and optimizing the boundary again by using the cone elimination method according to the comprehensive cost until the current boundary and the boundary obtained by the previous iteration optimization meet the preset convergence condition;
the step of calculating the land damage area of each working procedure in the current boundary by using an environmental impact quantization model and each land damage area allocation model specifically comprises the following steps:
determining a plurality of total land damage areas of the current boundary by using the environmental impact quantification model, wherein the plurality of total land damage areas comprise: the total area of the soil excavation loss, the total area of the soil discharge site, the total area of ground facilities, the total area of the gangue pile and the total area of the surface soil site;
Spreading the total area of each land damage to each working procedure, and respectively calculating each land damage area corresponding to each working procedure after spreading by utilizing the spreading model of each land damage area;
according to the land damage areas, calculating the land damage area of each working procedure;
the step of allocating each land damage total area to each working procedure and calculating each land damage area corresponding to each working procedure after allocation by using each land damage area allocation model comprises the following steps:
The corresponding land excavated area of raw coal exploitation is Wherein V m is the volume of raw coal, V c is the current boundary volume, and A c is the total area of the land excavation;
the excavated area of the rock-stripped land is Wherein V y is the rock volume;
the land digging area corresponding to the quaternary layer stripping is Wherein V 4 is the quaternary layer volume;
The occupied area of the dumping site corresponding to the rock stripping is Wherein A p is the occupied area of the soil discharge site;
the land occupation area of the dumping site corresponding to the quaternary layer stripping is ；
The area of the ground facility corresponding to the raw coal exploitation isWherein a s is the ground facility area;
the area of the ground facility corresponding to the rock stripping is ；
The area of the ground facility corresponding to the quaternary layer stripping is；
The area of the ground facility corresponding to the coal washing is；
The surface soil field occupation area corresponding to the raw coal exploitation isWherein h c is the average surface soil stripping thickness of a stope, h s is the average surface soil stripping thickness of a ground facility, V b is the surface soil field volume, and A b is the surface soil field occupied area;
The surface soil field occupation area corresponding to the rock stripping is
Wherein hp is the average stripping thickness of the surface soil of the dump;
The surface soil field occupied area corresponding to the quaternary layer stripping is as follows:
The surface soil field occupation area corresponding to the coal selection and washing is Wherein A g is the occupied area of the gangue pile, and h g is the average stripping thickness of the surface soil of the gangue pile;
the step of calculating the land damage area of each procedure according to each land damage area specifically comprises the following steps:
the land damage area of the raw coal exploitation is A m=Acm+Asm+Abm;
the damage area of the rock stripped land is A y= Apy+Acy+Asy+Aby;
The land damage area of the quaternary layer stripping is A 4= Ap4+Ac4+As4+Ab4;
the damage area of the soil for coal washing is A x= Ag+Asx+Abx;
The step of calculating the environmental cost per unit production of each process according to the land damage area of each process specifically comprises the following steps:
the environmental cost of raw coal exploitation is:
，
Wherein Z yield is the biomass production value of the destroyed land, Z csoil is the soil carbon fixation value of the land ecosystem, C rec is the reclamation and maintenance cost, Z cplant is the plant carbon fixation value of the destroyed land ecosystem, Z O2 is the oxygen release value of the destroyed land ecosystem, Z air is the air purification value of the destroyed land ecosystem, Z soil is the soil maintenance value of the destroyed land ecosystem, Z H2O is the water conservation value of the destroyed land ecosystem, Z neut is the nutrient circulation value of the destroyed land ecosystem, F is the time length, F t is the time coefficient, C gm is the emission environment cost of raw coal mining, and Q m is the raw coal weight;
the environmental cost of the rock stripping is:
，
wherein, C gy is the environmental cost of the emission of rock stripping;
The environmental cost of the quaternary layer stripping is as follows:
，
wherein, C g4 is the environmental cost of the emission of the quaternary layer stripping;
The environmental cost of coal washing is as follows:
，
Wherein, C gx is the emission environmental cost of coal washing.
2. The method of claim 1, wherein the step of determining the position of the substrate comprises,
The boundary parameters include at least one of: the current boundary raw coal amount, rock amount, quaternary layer amount, multiple land damage areas, greenhouse gas emission amount and the environmental cost of each working procedure.
3. A method according to claim 1 or 2, characterized in that,
The environmental cost per unit area includes at least one of: biomass production value, soil carbon fixation value, plant carbon fixation value, oxygen release value, air purification value, soil maintenance value, water source conservation value, nutrient circulation value and reclamation maintenance cost.
4. An apparatus for implementing a method of optimizing the final boundaries of an opencast coal mine as claimed in claim 1, said apparatus comprising:
a first determining module, configured to determine a unit production cost of a unit production amount of each process in a plurality of processes according to a mining technical condition of an opencast coal mine, where the plurality of processes include: raw coal mining, rock stripping, quaternary layer stripping and coal washing;
The second determining module is used for determining the environmental cost per unit area of the open pit coal mine according to the type of the ecological system;
the third determining module is used for taking the iteration number to be 0, and optimizing the boundary by utilizing a cone elimination method according to the unit production cost so as to determine the current boundary of the opencast coal mine;
The first calculation module is used for calculating the land damage area of each procedure in the current boundary by using an environmental impact quantization model and each land damage area allocation model;
the second calculation module is used for calculating the environmental cost of the unit production of each process according to the land damage area of each process;
the acquisition module is used for acquiring the iteration number of the optimization process;
the judging module is used for judging whether the iteration number is 0;
The judging module is further configured to judge whether the current boundary and the boundary obtained by previous iterative optimization meet a preset convergence condition if the iterative ordinal is not 0;
the output module is used for outputting the current boundary and the boundary parameters if the preset convergence condition is met, and the optimization is finished;
The optimization module is used for adding the environmental cost of each process unit production amount to the production cost of each process unit production amount to obtain the comprehensive cost, adding 1 to the iteration number, and re-optimizing the boundary by using a cone elimination method according to the comprehensive cost until the current boundary and the boundary obtained by the previous iteration optimization meet the preset convergence condition;
The judging module is further configured to judge whether the iteration number is greater than or equal to a preset iteration number threshold if the preset convergence condition is not satisfied;
The output module is further configured to output the current boundary and the boundary parameter if the iteration number is greater than or equal to the preset iteration number threshold, and the optimization is ended;
and the optimization module is further configured to add the environmental cost of each process unit throughput to the production cost of each process unit throughput if the iteration number is smaller than the preset iteration number threshold, obtain the integrated cost, add 1 to the iteration number, and optimize the boundary again by using the cone elimination method according to the integrated cost until the current boundary and the boundary obtained by the last iteration optimization meet the preset convergence condition.
5. A storage medium having stored thereon a computer program, which when executed by a processor, implements the method of any of claims 1 to 3.
6. An electronic device comprising a storage medium, a processor and a computer program stored on the storage medium and executable on the processor, characterized in that the processor implements the method of any one of claims 1 to 3 when executing the computer program.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210414436.XA CN114881298B (en) | 2022-04-20 | 2022-04-20 | Optimization method and device for final boundary of opencut coal mine, medium and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210414436.XA CN114881298B (en) | 2022-04-20 | 2022-04-20 | Optimization method and device for final boundary of opencut coal mine, medium and electronic equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114881298A CN114881298A (en) | 2022-08-09 |
| CN114881298B true CN114881298B (en) | 2024-05-03 |
Family
ID=82672287
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210414436.XA Active CN114881298B (en) | 2022-04-20 | 2022-04-20 | Optimization method and device for final boundary of opencut coal mine, medium and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114881298B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112214867A (en) * | 2020-08-28 | 2021-01-12 | 辽宁工程技术大学 | Strip mine mining boundary and mining program collaborative optimization method under complex coal seam condition |
| CN113344328A (en) * | 2021-05-07 | 2021-09-03 | 北方魏家峁煤电有限责任公司 | Method and device for delimiting outer soil discharge field of opencast coal mine |
| CN114169629A (en) * | 2021-12-14 | 2022-03-11 | 东北大学 | Mining plan optimization method and system considering ecological environment |
-
2022
- 2022-04-20 CN CN202210414436.XA patent/CN114881298B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112214867A (en) * | 2020-08-28 | 2021-01-12 | 辽宁工程技术大学 | Strip mine mining boundary and mining program collaborative optimization method under complex coal seam condition |
| CN113344328A (en) * | 2021-05-07 | 2021-09-03 | 北方魏家峁煤电有限责任公司 | Method and device for delimiting outer soil discharge field of opencast coal mine |
| CN114169629A (en) * | 2021-12-14 | 2022-03-11 | 东北大学 | Mining plan optimization method and system considering ecological environment |
Non-Patent Citations (2)
| Title |
|---|
| 《Ultimate pit optimazation with ecological cost for open pit metal mines》;Xu Xiaochuan etal;《Science Direct》;20140227;全文 * |
| 《生成成本内生化的金属露天矿境界-开采计划整体优化研究》;胥孝川;《万方数据》;20150915;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114881298A (en) | 2022-08-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhou et al. | Assessment and planning of underground space use in Singapore | |
| CN107885960B (en) | Earthwork volume estimation system and method based on automatic line selection of construction roads in wind power plant | |
| Moselhi et al. | Optimization of earthmoving operations in heavy civil engineering projects | |
| US9589076B2 (en) | Area-based open pit mine designer | |
| Haiko et al. | DEVELOPMENT OF A MORPHOLOGICAL MODEL FOR TERRITORIAL DEVELOPMENT OF UNDERGROUND CITY SPACE. | |
| CN116402642A (en) | Carbon evaluation method, system and storage medium for land development project | |
| CN114881298B (en) | Optimization method and device for final boundary of opencut coal mine, medium and electronic equipment | |
| CN114491755A (en) | Mountain land photovoltaic power station road planning method, device, equipment and storage medium | |
| Devanand et al. | Assessment of embodied energy of masonry building materials in Hassan, India | |
| US9552445B2 (en) | Automatic change propagation in an area-based open pit mine designer | |
| CN102877846A (en) | Unit method based mined and reclaimed topsoil stripping time determining method | |
| Bao et al. | Optimization production scheduling of underground backfilling mining based on NSGA-II | |
| CN112053031A (en) | Method for determining slope-expanding mining boundary of deep open pit mine | |
| Kaliampakos et al. | Underground space development: Setting modern strategies | |
| CN114936679A (en) | Ecological impact quantitative accounting method and system for opencast coal mining | |
| Sari et al. | A landfill based approach to surface mine design | |
| Hou et al. | Optimization of mining methods for deep orebody of large phosphate mines | |
| CN111597696A (en) | Method for evaluating water delivery quantity of oasis in arid region based on ecological hydrological simulation and optimization | |
| CN112711916A (en) | Power transmission line engineering construction equipment model selection method based on improved genetic algorithm | |
| CN106326613A (en) | Economic floor area ratio calculation method and system based on land economic mathematical model | |
| CN103226738A (en) | Reclamation opportunity optimization method based on expected reclamation cultivation rates and reclamation cost | |
| CA2914118A1 (en) | Method of determining networks for resource recovery operations | |
| NL2032888B1 (en) | Method and device for optimizing final boundary of open-pit coal mine, medium and electronic apparatus | |
| CN117892560B (en) | High-geographic-precision intelligent city safety drainage pipe network simulation and prediction method | |
| Amar et al. | RETRACTED: Construction Schedule Management for Densely Populated Areas Using CPM-Crashing |
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 |