CN109284848A - The long-range casting planing method of interior row's opencut and computing system - Google Patents

The long-range casting planing method of interior row's opencut and computing system Download PDF

Info

Publication number
CN109284848A
CN109284848A CN201710600403.3A CN201710600403A CN109284848A CN 109284848 A CN109284848 A CN 109284848A CN 201710600403 A CN201710600403 A CN 201710600403A CN 109284848 A CN109284848 A CN 109284848A
Authority
CN
China
Prior art keywords
stage
dump
dumps
advance rate
slope surface
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.)
Pending
Application number
CN201710600403.3A
Other languages
Chinese (zh)
Inventor
黄甫
肖双双
马力
丁小华
刘武
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN201710600403.3A priority Critical patent/CN109284848A/en
Publication of CN109284848A publication Critical patent/CN109284848A/en
Pending legal-status Critical Current

Links

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
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Theoretical Computer Science (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Quality & Reliability (AREA)
  • Animal Husbandry (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Agronomy & Crop Science (AREA)
  • Development Economics (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Mining & Mineral Resources (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

The invention discloses the long-range casting planing method of interior row's opencut and computing systems, include basic data acquisition, calculation processing, result output module.By constantly adjusting mined out sector width of each stage, tentative calculation respectively dumps the stage and dumps height this method, and further determines that and respectively dump stage mass center.According to stripping stage mass center is respectively adopted with the corresponding stage mass center that dumps the sum of horizontal transport distance and vertical-lift height is solved, and further determines that each stage overburden freight.During constantly adjusting mined out sector width of each stage, if the sum of each stage overburden freight calculated is less than last calculated result, it is iterated until reaching the calculation times of setting or utilizing all combined situations of the exhaustive complete each stage mined out sector width of the method for exhaustion.Finally calculate the minimum cost that overburden in calculating cycle transports, at the same can provide respectively dump the stage dump parameter.The present invention can instruct the design of opencut internal dumping, save the transportation cost of overburden.

Description

The long-range casting planing method of interior row's opencut and computing system
Technical field
The present invention relates to the long-term plan of opencut internal dumping designs, are suitable for opencut internal dumping under the conditions of interior row Long-range dumps engineered locations design, can be realized the purpose that overburden transportation cost is minimum in calculating cycle.
Background technique
The transportation cost of opencut overburden accounts for the significant proportion of cost of winning.Casting engineered locations design (including Refuse dump dumps height, mined out sector width) transportation range of the overburden from stope through end side transport to refuse dump is directly influenced, Thus rationally economically plan that refuse dump can reduce the transportation cost of overburden.
Currently, the mode that opencut casting design generallys use manual tentative calculation carries out the design of refuse dump, with casting capacity Meeting production requirement is target, does not consider casting transportation cost in this process.
Time this for " monthly plan " shorter " short-term casting design ", if by stope, refuse dump several digging points, Dry dump a little is reduced to fixed point, and weight is arranged according to distance, gradient etc. in the existing route of mine, utilizes the side such as network Method is calculated overburden and transports the specific optimal flux flow direction transported to refuse dump from stope.When this for " yearly plan " Between longer " mid-term casting design ", then be to divide overburden according to step, obtain the volume and mass center of each step, together Each step mass center in the refuse dump that sample is designed and rock amount determine the flow of overburden using the methods of network on this basis And flow direction.
Currently, being still for " long-range casting design " (time span was in 3 years, 5 years or even entire mine life cycle) It is time-consuming and laborious using the design method of similar " year design ", and there is obvious drawback.First, " short-term/mid-term casting design " Refuse dump design when only consider that capacity does not consider cost, transport optimizing is carried out after the completion of refuse dump is designed.It is " short The casting engineered locations of phase/mid-term casting design " can not determine whether casting design economic benefit is optimal, and it is each that it is even more impossible to determinations Can stage casting engineered locations make entire mine life cycle In-pit dumping optimal economic benefit.Second, " short-term/mid-term casting Design " frequently adjusted due to only considering that the shorter time cycle may cause open-pit mine blast height, be unfavorable for it is subsequent reclaim, Third, casting space can not in due course reasonable utilization, transport Dumping cost is high.
Summary of the invention
The present invention is exactly in order to solve the above-mentioned technical problem, to provide the long-range row of open-pit mine blast under the conditions of arranging in one kind The computing system of soil planning, can make the overburden transportation cost in entire calculating cycle minimum, and provides and each dump the stage Casting engineered locations (including refuse dump dump height, dump stage slope surface position).
Computing system is planned in the long-range casting of open-pit mine blast under the conditions of a kind of interior row provided by the invention, includes basic number Three modules are exported according to acquisition, calculation processing, result.It is characterized in that basic data, which obtains module, obtains three kinds of data packets of stope It includes and adopts stripping stage bottom of slope position, adopt stripping stage mass center and rock amount volume, three kinds of refuse dump data include that fixed advance rate dumps rank Section cutting solid volume and mass center, fixed advance rate dump stage bottom of slope position, specifically include the following steps:
Step 1: establishing customized one-dimensional coordinate system for characterizing the position of stope slope surface Yu refuse dump slope surface.Choose one For the straight line vertical with stope direction of propulsion as reference line, the distance in plane a little to the straight line is the point customized one The coordinate for tieing up coordinate system, does not consider elevation.Regulation straight line is negative along the coordinate side of stope direction of propulsion, the other side be positive ( It can make opposite regulations)；
Step 2: the stope Block Model comprising ore Rock information is established, it is " medium-term and long-term by 2013 editions 3DMine softwares Planning "-" cutting of knife amount " function is respectively adopted the rock amount volume in stripping stage and generates each stage entity file, and 3DMine is called VBA interface reads these barycenter datas for adopting stripping stage entity；
Step 3: being established in CAD or 3DMine and respectively adopt stripping stage slope surface DTM, call 3DMine VBA interface to obtain each It adopts stripping stage slope surface DTM and (point is most for the maximum point of coordinate in customized one-dimensional coordinate system on (DTM) intersection of final mining area hole Close to goaf and refuse dump side, when the positive negative direction of customized one-dimensional coordinate system makees opposite regulations, take coordinate minimum Point), as adopting stripping stage bottom of slope position；
Step 4: establishing internal dumping Block Model, the boundary of internal dumping exists using earth's surface boundary as boundary along side Connect when outer dump with along side outer dump." differential " is carried out after establishing internal dumping, i.e., is pressed along refuse dump developing direction Setting dump slope angle and fixed advance rate carries out " cutting " and obtains dumping the stage under the conditions of fixed advance rate, in order to mould What type to be solved dump, and the stage is mutually distinguished, this is dumped stage definitions and dumps the stage for fixed advance rate.Simultaneously in vertical direction On carry out vertical " cutting " according to certain intervals and obtain the reality dumped after the cutting of stage vertical direction under the conditions of fixed advance rate Body, these cutting substantial definitions are that fixed advance rate dumps stage cutting entity.It calls 3DMine VBA interface to read these to cut Cut the volume and barycenter data of entity；
Step 5: establishing each fixed advance rate in CAD or 3DMine and dump the stage and dump slope surface DTM, call 3DMine VBA interface obtains each fixed advance rate and dumps that the stage dumps slope surface DTM and final mining area is cheated on (DTM) intersection customized one-dimensional The smallest point of coordinate in coordinate system is (point near goaf and stope side, when the positive negative direction of customized one-dimensional coordinate system When making opposite regulations, coordinate maximum point is taken), stage bottom of slope position is dumped as each fixed advance rate.
Computing system is planned in the long-range casting of open-pit mine blast under the conditions of a kind of interior row provided by the invention, includes basic number According to acquisition, calculation processing, result export three modules, can be calculated in entire calculating cycle overburden transport it is minimum at Sheet and each casting engineered locations for dumping the stage.It is characterized by:
Calculation processing module is realized using following steps:
Step 1: reading the data in basic data acquisition module；
Step 2: initiation parameter dumps stage mined out sector width L (i)=Lmax(i=1to N0), Flag (i) (mark I-th dumps whether the stage meets constraint condition), Eflag mark the case where having traversed all mined out sector widths, in calculating cycle Freight Cf；
L (i)=Lmax(i=1to N0), Flag (i)=False (i=1to N0), Eflag=Flase, Cf=1e100, J=1 (calculates first and dumps the stage)；
Step 3: judging Eflag=False, be, enter step 4, otherwise enter step 17；
Step 4: judging Flag (j)=False, be to enter step 5, be to enter step 10；
Step 5: according to jth stage mined out sector width LjWith adopt stripping stage bottom of slope positionIt determines and dumps stage slope surface solid Determine advance rate and dump position in the stage, that is, dumps stage slope surface and be in fixed advance rate and dump which in stage and apart from this in stage Stage casting slope surface distance further obtains fixed advance rate and dumps stage cutting entity and adopted the stripping stage dumping slope surface and cutting It is segmented into the mass center and volume of each section after two parts；
Step 6: " integrating " in vertical direction --- it successively adds up upwards from lowest part and dumps stage slope surface in j-th Position dumps stage slope surface position with j-1 and (dumps stage slope surface position and model calculating cycle starting point according to the 1st when j=1 When initial refuse dump position or initial land form) between fixation advance rate dump the volume of stage cutting entity, until big In current required casting space or it is added to one layer entity of topmost, it is empty according to cumulative solid volume and current required casting Between difference, preceding once accumulative solid volume, the height of a most upper layer entity carries out difference and determines that dumping the stage for j-th dumps Height Hj.To in dumping stage slope surface position for j-th and j-1 is a dumps stage slope surface position and (dump when j=1 according to the 1st Initial refuse dump position or initial land form when stage slope surface position and model calculating cycle starting point) between, dump rank j-th The mass center that section dumps the entity between height and boundary hole bottom plate is weighted to obtain j-th according to the volume of each entity dumps rank Section dumps mass center；
Step 7: calculating is adjacent to dump stage height difference Δ H=Hj-Hj-1If height difference Δ H, which is less than, gives high difference h, Flag (j)=True, otherwise Flag (j)=False；
Step 8: judging Flag (j)=True, 10 are entered step when being true, enter step 9 for fictitious time；
Step 9: judging L (j) > LminEven the mined out sector width in jth stage reduces step pitch Lstep L (j)=L when being true (j) 3 are entered step-Lstep and, enters step 10 for fictitious time；
Step 10: judging Flag (j)=True, 11 are entered when being true, enter 12 for fictitious time；
Step 11: judging j=N0, that is, judge whether that having calculated the last one dumps the stage.When being true according to respectively adopt stripping rank Duan Zhixin and the horizontal haul distance of overburden and vertical-lift height in stage centroid calculation each stage are dumped accordingly, and further counted Target function value C is calculated, enters step 13；J=j+1 (carry out next dump stage calculate) is enabled for fictitious time, enters step 3；
Step 12: judging j=1, that is, judge currently to calculate dumps whether the stage is first to dump the stage in stage, when being true Eflag=True is enabled, enters step 3, enters step 14 for fictitious time；
Step 13: judging C < Cf, it is very season Cf=C, iteration respectively dump stage parameter, are directly entered 14 for fictitious time；
Step 14: judging L (j)≤Lmin, 15 are entered step when being true, are included in 16 for fictitious time；
Step 15: judging j=1, i.e., whether be first to dump the stage, be that very season Eflag=True enters step 16, be Fictitious time enables j=j -1, carries out the previous calculating for dumping the stage, enters step 14；
Step 16: judging Eflag=False, the mined out sector width in jth stage reduces step pitch Lstep when being true, enables L (j) =L (j)-Lstep, to partial parameters assignment Flag (i)=False (i=j to N0), L (i)=Lmax(i=j+1to N0), 3 are entered step, is directly entered step 3 for fictitious time；
Step 17: export following data (including but not limited to following data): respectively dump stage overburden freight it With respectively dump the absolute altitude that dumps in stage, respectively dump fixation advance rate locating for the stage and dump the stage, respectively dump the stage and dump slope surface It dumps the distance for dumping slope surface in stage apart from locating fixation advance rate, each stage overburden horizontal transport distance and vertically mentions Rise.
The detailed process of the step 5 are as follows:
5-1, stripping stage bottom of slope position is adopted according to i-thI-th of mined out sector width LiIt calculates i-th and dumps stage bottom of slope Position(adopting stripping stage label i to reduce along customized one-dimensional coordinate system positive direction)；
5-2, by traverse fixed advance rate dump stage bottom of slope position can determine there are goaf from Adopt stripping stage bottom of slope positionNearest fixation advance rate dumps the stage, labelPosition in array is Mi(regulation fixation pushes away The label j that progress dumps the stage increases along customized one-dimensional coordinate system positive direction),Wherein:For Fixed advance rate jth dumps the bottom of slope position in stage；
5-3, fixation advance rate locating for stage slope surface is dumped it is possible to further obtaining i-th dump stage mi
Wherein MiMeetFor fixed advance rate MiA stage bottom of slope position that dumps is one Tie up the coordinate in coordinate system, m；ncStage bottom of slope position is dumped for fixed advance rate be in dump slope surface and MiA fixed propulsion Degree dumps the number between slope surface, dimensionless；
5-4 is according to miIt calculates i-th and dumps stage slope surface and miA fixed advance rate dumps stage slope surface and adopts in horizontal plane Distance on the direction of propulsion of field
Wherein: lDThe fixation advance rate in stage, m are dumped for fixed advance rate；For fixed advance rate miIt is a to dump rank Coordinate of the section bottom of slope position in one-dimensional coordinate system, m；For fixed advance rate mi+ 1 dumps stage bottom of slope position one Tie up the coordinate in coordinate system, m；
5-5 calculates miA fixed advance rate dumps in stage vertical direction j-th of entity and dumps stage slope surface by i-th The mass center seat of k-th of the entity (k=1 is to be in currently to dump the stage, and k=2 be to be in next to dump the stage) formed after cutting Mark
WhereinFor miA fixed advance rate dumps in stage vertical direction j-th of entity mass center Centroid vector form coordinate, m；Be customized one-dimensional coordinate system positive direction unit vector in three-dimensional system of coordinate to Measure expression-form, m；
5-6, dump stage slope surface by the fixation advance rate where it dump stage splitting object be two-part each section Volume then is weighted to obtain by cross-sectional area:
Wherein Vi,j,kFor miA fixed advance rate dumps in stage vertical direction j-th of entity and dumps rank by i-th K-th of the entity (k=1 is to be in currently to dump the stage, and k=2 is to be in next to dump the stage) formed after section slope surface cutting Volume, m3；Vmi,jFor miA fixed advance rate dumps j-th of solid volume in stage vertical direction, m3。
Detailed process in the step 6 are as follows:
6-1, it is dumped in the case that the stage slope surface position i determines in steps of 5, from i stage slope surface is dumped to dumping i-1 rank Fixation advance rate between section slope surface dumps stage cutting entity and is successively accumulated to n-th layer entity from top to bottom in vertical direction When volume Vn:
Wherein: n0The number of stages in stage, dimensionless are dumped for fixed advance rate；Vi,j,1V when being k=1i,j,k, i.e., miA fixed advance rate dumps j-th of entity in stage vertical direction and is dumped being in of being formed after the cutting of stage slope surface for i-th I-th of volume for dumping the entity in the stage, m3；Vi-1,j,2V when being k=2i-1,j,k, i.e. mi-1A fixed advance rate dumps In stage vertical direction j-th of entity dumped by (i-1)-th formed after the cutting of stage slope surface be in i-th and dump in the stage The volume of entity, m3；
Work as VnMeetOr n=nhWhen (being accumulated to a most upper layer entity), calculates i-th and dump stage row Abandon the height difference of n-th layer entity Yu n-th layer entity starting altitude
In formula: HeRefuse dump Block Model highest elevation, m；HbFor cutting when starting elevation, i.e., from the elevation the following are First entity, m；HsThe fixed cutting spacing of stage vertical direction, m are dumped for fixed advance rate；nhIt is dumped for fixed advance rate Stage vertically cut after number of entities, dimensionless.(note: n=nhWhen formula in " % " be take the remainder operation)；
6-2, it calculates i-th and dumps the stage and dump slope surface height Hi:
6-3, determining that dumping the stage dumps height HiAfterwards, it is also necessary to it determines and dumps stage center-of-mass coordinate, thus needs pair Dump in the stage that n-th of entity exists in vertical directionIt is cut at height.Stage casting mark is dumped after cutting at i-th The x of high its mass center of part below, y-coordinate is constant, and z coordinate is that the coordinate of entity mass center before cutting subtractsIts Middle HnFor the fixed advance rate of the n-th (n >=2) layer dump the stage cut entity thickness (" % " be complementation in formula, " and " to divide exactly):
Dumping at i-th the stage after 6-4, cutting dumps the solid volume of absolute altitude or less part and multiplies for solid volume before cutting With
6-5, the volume for dumping i-th each entity in phase range can obtain the weighting of its center-of-mass coordinate as weight Stage center-of-mass coordinate is dumped to i-th
Computing system is planned in the long-range casting of open-pit mine blast under the conditions of a kind of interior row provided by the invention, includes basic number According to acquisition, calculation processing, result export three modules, can be calculated in entire calculating cycle overburden transport it is minimum at Sheet and each casting engineered locations for dumping the stage.It is characterized by:
As a result the absolute altitude that dumps for respectively dumping the stage that output module is exported according to calculation processing module, respectively dumped locating for the stage Fixation advance rate dump the stage, respectively dump the stage and dump slope surface and dump slope surface apart from what locating fixation advance rate dumped the stage Distance, that internal dumping is drawn out in mapping software (including but not limited to 3DMine software) plot space respectively dumps the stage Step bench crest, step bench toe.
Compared with the prior art, the advantages of the present invention are as follows:
1) traditional refuse dump design does not consider cost, in the present invention design of refuse dump consider overburden transport at This, several is divided into entire calculating cycle adopts the stripping stage and with corresponding dump the stage, and the present invention can make in calculating cycle respectively The sum of the overburden transportation cost in stage is minimum；
2) design of traditional refuse dump needs by hand repeatedly tentative calculation until casting space meets error requirements, and the present invention can Automatically the engineered locations that dump in each stage are provided, including each stage dumps absolute altitude, casting slope surface position (step bench crest/step Bench toe)；
3) the adjacent height difference maximum value for dumping the stage can be set in the present invention, make respectively to dump the stage dump absolute altitude fluctuation compared with It is small, conducive to casting operation and subsequent reclaim.
Detailed description of the invention
Straight line a is the reference line of customized one-dimensional coordinate system in Fig. 1, C be stope slope surface and final mining area cheat on intersection c away from The farthest point of reference line, distance are LC.D is that refuse dump slope surface and final mining area cheat point nearest away from reference line on intersection d, distance For LD。
Fig. 2 is that fixed advance rate dumps the stage and fixed advance rate dumps stage cutting entity, lDFor fixed advance rate, Hb The top absolute altitude of elevation, H are cut for the 1st layereFor internal dumping highest absolute altitude, nhFor the total number of plies for cutting entity, HsIt is arrived for 2 nhThe height of -1 entity.
Fig. 3 is the calculating schematic diagram for dumping the stage for i-th, and overstriking solid line is to dump the stage for i-th to dump range in figure.mi Fixation advance rate locating for the stage, which is dumped, for i-th dumps the stage.
Entity is that the fixation advance rate of shade filling in Fig. 3 dumps stage cutting entity in Fig. 4.H point is that fixation pushes away in Fig. 4 Progress dumps the mass center of stage cutting entity A BCD, and G is the mass center of AEFD, and I is the mass center of BCFE,The stage is dumped for i-th Dump the distance for dumping slope surface that slope surface dumps the stage apart from locating advance rate.
Fig. 5 is the algorithm flow chart of calculation processing module.
Fig. 6 is the structure and data flow diagram of interior row's open-pit mine blast Long-term planning computing system.
Specific embodiment
The long-range casting planing method of interior row's opencut of the invention and computing system, by taking the opencut of Inner Mongol as an example, specifically Steps are as follows:
Basic data obtains module:
Step 1: establishing customized one-dimensional coordinate system, straight line is chosen (with stope propulsion side in stope side outside boundary To vertically) it is used as reference line, it is some coordinate of this in customized one-dimensional coordinate system to the distance of the straight line in plane, no Consider elevation.Regulation straight line is negative along the coordinate side of stope direction of propulsion, and the other side is positive.
Step 2: establishing stope Block Model, obtain each stage rock amount using 3DMine software, call VBA interface to read each Stage mass center.
Step 3: according to each stage stope slope surface in 956 horizontal (for positioning each stage, other levels also may be selected) The position line makes each stage stope slope surface DTM with CAD or 3DMine, and 3DMine VBA interface is called to obtain each stage stope slope On face and final mining area DTM intersection in a customized coordinate system the maximum point of coordinate.
Step 4: internal dumping Block Model is established, substantially using earth's surface boundary as boundary, when existing along side outer dump Connect with outer dump, refuse dump highest absolute altitude is set as higher numerical value.Fixed advance rate is obtained using 3DMine software to arrange Abandoning stage each stage entity, and (cutting is complete in 3DMine for the volume and mass center for calling VBA interface to read each stage cutting entity At reading is completed in VBA environment).
Step 5: dumping stage slope surface according to each fixed advance rate (for positioning each stage, also may be selected it in 960 levels It is horizontal) the position line make each fixed advance rate with CAD or 3DMine and dump stage slope surface DTM, call 3DMine VBA Interface obtains each stage and dumps on slope surface and final mining area DTM intersection the maximum point of coordinate in a customized coordinate system.
Calculation processing module:
Step 1: read basic data obtain module data, wherein stope data include the rock amount volume for adopting the stripping stage, It adopts stripping stage center-of-mass coordinate, adopt stripping stage bottom of slope position, refuse dump data include that fixed advance rate dumps stage cutting entity Volume, fixed advance rate dump the volume of stage cutting entity, fixed advance rate dumps the center-of-mass coordinate, solid of stage cutting entity Determine the bottom of slope position that advance rate dumps the stage.
Step 2: parameter assignment: the coefficient of volumetric expansion takes 1.15, goaf maximum width Lmax=200m, goaf minimum widith LminThe goaf=20m adjusts step pitch Lstep=20m, adjacent to dump stage maximum height difference Δ H=48m.Initiation parameter: it dumps Stage mined out sector width L (i), (mark has traversed by Flag (i) (mark the i-th dumps whether the stage meets constraint condition), Eflag The case where all mined out sector widths), the freight C in calculating cyclef。
L (i)=Lmax(i=1to N0), Flag (i)=False (i=1to N0), Eflag=Flase, Cf=1e100, J=1 (calculates first and dumps the stage)
Step 3: judging Eflag=False, be, enter step 4, otherwise enter step 17
Step 4: judging Flag (j)=False, be to enter step 5, be to enter step 10
Step 5: according to mined out sector width LjWith adopt stripping stage bottom of slope positionDetermine that dump stage slope surface promotes fixed Degree dumps position in the stage, that is, dumps stage slope surface and be in fixed advance rate and dump which in stage in stage and arrange apart from the stage Slight slope identity distance is from further, obtaining fixed advance rate and dump stage cutting entity and adopted the stripping stage dumping slope surface and being cut into two The mass center and volume of each section behind part.
Step 6: " integrating " in vertical direction --- it successively adds up upwards from lowest part and dumps stage slope surface in j-th Position dumps stage slope surface position with j-1 and (dumps stage slope surface position and model calculating cycle starting point according to the 1st when j=1 When initial refuse dump position or initial land form) between fixation advance rate dump the volume of stage cutting entity, until big In current required casting space or it is added to one layer entity of topmost, it is empty according to cumulative solid volume and current required casting Between difference, preceding once accumulative solid volume, the height of a most upper layer entity carries out difference and determines that dumping the stage for j-th dumps Height Hj.To in dumping stage slope surface position for j-th and j-1 is a dumps stage slope surface position and (dump when j=1 according to the 1st Initial refuse dump position or initial land form when stage slope surface position and model calculating cycle starting point) between, dump rank j-th The mass center that section dumps the entity between height and boundary hole bottom plate is weighted to obtain j-th according to the volume of each entity dumps rank Section dumps mass center.
Step 7: calculating is adjacent to dump stage height difference Δ H=Hj-Hj-1If height difference Δ H, which is less than, gives high difference h, Flag (j)=True, otherwise Flag (j)=False.
Step 8: judging Flag (j)=True, 10 are entered step when being true, enter step 9 for fictitious time
Step 9: judging L (j) > LminEven the mined out sector width in jth stage reduces step pitch Lstep L (j)=L when being true (j) 3 are entered step-Lstep and, enters step 10 for fictitious time
Step 10: judging Flag (j)=True, 11 are entered step when being true, enter step 12 for fictitious time
Step 11: judging j=N0, that is, judge whether that having calculated the last one dumps the stage.Evidence respectively adopts the stripping stage when being true Mass center and the horizontal haul distance of overburden and vertical-lift height in stage centroid calculation each stage are dumped accordingly, and further calculated The sum of stage overburden freight C is respectively dumped, enters step 13；J=j+1 is enabled (to carry out next dumping the stage for fictitious time Calculate), enter step 3.
C calculation formula are as follows:
In formula: λ is the overburden coefficient of volumetric expansion, dimensionless；cNFor level road freight (round-trip), cN=5 yuan/kmm3；cLu For heavily loaded upward slope freight, cLu=3 yuan/kmm3；cUdFor empty wagons descending freight, cUd=1.6 yuan/kmm3；cLdFor Heavily loaded descending freight, cLd=2.4 yuan/kmm3；cUuFor empty wagons upward slope freight, cUu=1.8 yuan/kmm3；G is oblique Ramp limits slope, takes g=0.06.
Step 12: judging j=1, that is, judge currently to calculate dumps whether the stage is first to dump the stage in stage, when being true Eflag=True is enabled, enters step 3, enters step 14 for fictitious time.
Step 13: judging C < Cf, it is very season Cf=C, iteration respectively dump stage parameter, are directly entered step 14 for fictitious time.
Step 14: judging L (j)≤Lmin, 15 are entered step when being true, enter step 16 for fictitious time.
Step 15: judging j=1, i.e., whether be first to dump the stage, be that very season Eflag=True enters step 16, be Fictitious time enables j=j -1, carries out the previous calculating for dumping the stage, enters step 14.
Step 16: judging Eflag=False, the mined out sector width in jth stage reduces step pitch Lstep when being true, enables L (j) =L (j)-Lstep, to partial parameters assignment Flag (i)=False (i=j to N0), L (i)=Lmax(i=j+1to N0), 3 are entered step, is directly entered step 3 for fictitious time.
Step 17: output data --- the minimum freight finally determined is 8.20426E+10 member.The part in each stage Dump the following (note: N of parameteriIndicating to dump for i-th the stage dumps height and dumps stage cutting entity N in fixed advance ratei Layer):
As a result output module:
Engineered locations data are dumped according to each stage exported in calculation processing module, call 3DMine VBA interface will be defeated The engineered locations data that dump out are drawn out in 3DMine plot space and respectively dump stage step bench crest/step bench toe.

Claims (5)

  1. The long-range casting planing method of row's opencut and computing system in 1., include basic data acquisition, calculation processing, result output Three modules, it is characterised in that it includes following operation that basic data, which obtains module:
    Customized one-dimensional coordinate system is established for characterizing the position of stope slope surface Yu refuse dump slope surface, one is chosen and is promoted with stope For the vertical straight line in direction as reference line, the distance in plane a little to the straight line is this in customized one-dimensional coordinate system Coordinate does not consider elevation, it is specified that straight line is negative along the coordinate side of stope direction of propulsion, and the other side is positive (can also advise on the contrary It is fixed)；
    Stripping stage slope surface DTM is respectively adopted in foundation, respectively to adopt on stripping stage slope surface DTM and final mining area hole DTM intersection customized one It is (positive and negative when customized one-dimensional coordinate system the point near goaf and refuse dump side to tie up the maximum point of coordinate in coordinate system When opposite regulations are made in direction, coordinate smallest point is taken), as adopting stripping stage bottom of slope position；
    Internal dumping Block Model is established, the boundary of internal dumping is using earth's surface boundary as boundary, when existing along side outer dump Connect with along side outer dump, carries out " differential " after establishing internal dumping, i.e., in the dumping by setting along refuse dump developing direction Slope angle and fixed advance rate carry out " cutting " and obtain dumping the stage under the conditions of fixed advance rate, in order to model needs solution Dumping the stage mutually distinguishes, this is dumped stage definitions and dumps the stage for fixed advance rate, while in vertical direction according to certain Interval carries out vertical " cutting " and obtains the entity dumped after the cutting of stage vertical direction under the conditions of fixed advance rate, these cuttings Substantial definition is that fixed advance rate dumps stage cutting entity, and 3DMine VBA interface is called to read the volume of these cutting entities And barycenter data；
    Establish each fixed advance rate and dump the stage and dump slope surface DTM, with each fixed advance rate dump the stage dump slope surface DTM with most On the DTM intersection of whole boundary hole in customized one-dimensional coordinate system the smallest point of coordinate (point near goaf and stope side, When the positive negative direction of customized one-dimensional coordinate system makees opposite regulations, coordinate maximum point is taken), rank is dumped as each fixed advance rate Section bottom of slope position.
  2. The long-range casting planing method of row's opencut and computing system in 2., include basic data acquisition, calculation processing, result output Three modules, it is characterised in that calculation processing module is realized using following steps:
    Step 1: reading the data in basic data acquisition module；
    Step 2: initiation parameter dumps stage mined out sector width L (i)=Lmax(i=1to N0), Flag (i) (mark the i-th row Whether the abandoning stage meets constraint condition), Eflag indicates the case where having traversed all mined out sector widths, the transport in calculating cycle Expense Cf,
    L (i)=Lmax(i=1to N0), Flag (i)=False (i=1to N0), Eflag=Flase, Cf=1e100, j=1 (calculate first and dump the stage)；
    Step 3: judging Eflag=False, be, enter step 4, otherwise enter step 17；
    Step 4: judging Flag (j)=False, be to enter step 5, be to enter step 10；
    Step 5: according to jth stage mined out sector width LjWith adopt stripping stage bottom of slope positionDetermine that dump stage slope surface pushes away in fixation Progress dumps the position in the stage, that is, dumps stage slope surface and be in fixed advance rate and dump which in stage and apart from the stage in stage Casting slope surface distance further obtains fixed advance rate and dumps stage cutting entity and adopted the stripping stage dumping slope surface and being cut into The mass center and volume of each section after two parts；
    Step 6: " integrating " in vertical direction --- it successively adds up upwards from lowest part and dumps stage slope surface position in j-th Stage slope surface position is dumped (when dumping stage slope surface position and model calculating cycle starting point according to the 1st when j=1 with j-1 Initial refuse dump position or initial land form) between fixation advance rate dump the volume of stage cutting entity, work as until being greater than Casting space needed for phase is added to one layer entity of topmost, according to cumulative solid volume and current required casting space Difference, preceding once accumulative solid volume, the height of a most upper layer entity carry out difference and determine that dumping the stage for j-th dumps height Hj, in dumping stage slope surface position for j-th and j-1 dump stage slope surface position and (dump the stage according to the 1st when j=1 Initial refuse dump position or initial land form when slope surface position and model calculating cycle starting point) between, dump stage row for j-th The mass center for the entity abandoned between height and boundary hole bottom plate is weighted to obtain j-th according to the volume of each entity dumps the stage Dump mass center；
    Step 7: calculating is adjacent to dump stage height difference Δ H=Hj-Hj-1If height difference Δ H, which is less than, gives high difference h, Flag (j) =True, otherwise Flag (j)=False；
    Step 8: judging Flag (j)=True, 10 are entered step when being true, enter step 9 for fictitious time；
    Step 9: judging L (j) > LminEven the mined out sector width in jth stage reduces step pitch Lstep L (j)=L (j)-when being true Lstep simultaneously enters step 3, enters step 10 for fictitious time；
    Step 10: judging Flag (j)=True, 11 are entered step when being true, enter step 12 for fictitious time；
    Step 11: judging j=N0, that is, judge whether that having calculated the last one dumps the stage, according to respectively adopting stripping stage mass center when being true And the horizontal haul distance of overburden and vertical-lift height in stage centroid calculation each stage are dumped accordingly, and further calculate target Functional value C, enters step 13；J=j+1 (carry out next dump stage calculate) is enabled for fictitious time, enters step 3；
    Step 12: judging j=1, that is, judge currently to calculate dumps whether the stage is first to dump the stage in stage, is very seasonal Eflag=True enters step 3, enters step 14 for fictitious time；
    Step 13: judging C < Cf, it is very season Cf=C, iteration respectively dump stage parameter, are directly entered step 14 for fictitious time；
    Step 14: judging L (j)≤Lmin, 15 are entered step when being true, enter step 16 for fictitious time；
    Step 15: judging j=1, i.e., whether be first to dump the stage, be that very season Eflag=True enters step 16, be fictitious time J=j -1 is enabled, the previous calculating for dumping the stage is carried out, enters step 14；
    Step 16: judging Eflag=False, the mined out sector width in jth stage reduces step pitch Lstep when being true, enables L (j)=L (j)-Lstep, to partial parameters assignment Flag (i)=False (i=j to N0), L (i)=Lmax(i=j+1to N0), enter Step 3, step 3 is directly entered for fictitious time；
    Step 17: exporting following data (including but not limited to following data): respectively dumping the sum of stage overburden freight, respectively The absolute altitude that dumps for dumping the stage respectively dumps fixation advance rate locating for the stage and dumps the stage, respectively dumps the stage and dump slope surface distance Locating fixation advance rate dumps the distance for dumping slope surface in stage, and each stage overburden horizontal transport distance and vertical-lift are high Degree.
  3. 3. the step 5 of the calculation processing module as described in right 2, it is characterised in that the detailed process of the step 5:
    5-1, stripping stage bottom of slope position is adopted according to i-thI-th of mined out sector width LiIt calculates i-th and dumps stage bottom of slope position (adopting stripping stage label i to reduce along customized one-dimensional coordinate system positive direction)；
    5-2, stage bottom of slope position is dumped by traversing fixed advance rate can determine there are goaf from adopting stripping Stage bottom of slope positionNearest fixation advance rate dumps the stage, labelPosition in array is Mi(regulation is fixed to be promoted Degree dumps the label j in stage and increases along customized one-dimensional coordinate system positive direction), Mi=j,Wherein:It is solid Determine the bottom of slope position that advance rate jth dumps the stage；
    5-3, fixation advance rate locating for stage slope surface is dumped it is possible to further obtaining i-th dump stage mi
    Wherein MiMeet For fixed advance rate MiA stage bottom of slope position that dumps is in one-dimensional seat Coordinate in mark system, m；ncStage bottom of slope position is dumped for fixed advance rate be in dump slope surface and MiA fixed advance rate row Abandon the number between slope surface, dimensionless；
    5-4 is according to miIt calculates i-th and dumps stage slope surface and miA fixed advance rate dumps stage slope surface and pushes away in horizontal plane stope Into the distance on direction
    Wherein: lDThe fixation advance rate in stage, m are dumped for fixed advance rate；For fixed advance rate miIt is a to dump stage bottom of slope Coordinate of the position in one-dimensional coordinate system, m；For fixed advance rate mi+ 1 dumps stage bottom of slope position in one-dimensional coordinate Coordinate in system, m；
    5-5 calculates miA fixed advance rate dumps in stage vertical direction j-th of entity and dumps the cutting of stage slope surface by i-th The center-of-mass coordinate of k-th of the entity (k=1 is to be in currently to dump the stage, and k=2 is to be in next to dump the stage) formed afterwards
    WhereinFor miA fixed advance rate dump the mass center of j-th of entity mass center in stage vertical direction to Amount form coordinate, m；The vector for being customized one-dimensional coordinate system positive direction unit vector in three-dimensional system of coordinate expression Form, m；
    5-6, dump stage slope surface by the fixation advance rate where it dump stage splitting object be two-part each section volume Then weight to obtain by cross-sectional area:
    Wherein Vi,j,kFor miA fixed advance rate dumps in stage vertical direction j-th of entity and dumps stage slope by i-th The volume of k-th of the entity (k=1 is to be in currently to dump the stage, and k=2 is to be in next to dump the stage) formed after the cutting of face, m3；For miA fixed advance rate dumps j-th of solid volume in stage vertical direction, m3。
  4. 4. the step 6 of the calculation processing module as described in right 2, it is characterised in that the detailed process of the step 6:
    6-1, it is dumped in the case that the stage slope surface position i determines in steps of 5, from i stage slope surface is dumped to dumping the stage slope i-1 Fixation advance rate between face dumps when stage cutting entity is successively accumulated to n-th layer entity from top to bottom in vertical direction Volume Vn:
    Wherein: n0The number of stages in stage, dimensionless are dumped for fixed advance rate；Vi,j,1V when being k=1i,j,k, i.e. miIt is a Fixed advance rate dumps j-th of entity in stage vertical direction and is dumped being in i-th of being formed after the cutting of stage slope surface for i-th Dump the volume of the entity in the stage, m3；Vi-1,j,2V when being k=2i-1,j,k, i.e. mi-1It is vertical that a fixed advance rate dumps the stage Upward j-th of the entity of histogram is dumped by (i-1)-th the entity that i-th dumped in the stage of being in formed after the cutting of stage slope surface Volume, m3；
    Work as VnMeet Vn-1< λ Vi s< VnOr n=nhWhen (being accumulated to a most upper layer entity), calculates i-th and dump the stage and dump The height difference of n-th layer entity and n-th layer entity starting altitude
    In formula: HeRefuse dump Block Model highest elevation, m；HbStarting elevation when to cut, i.e., from the elevation, the following are first A entity, m；HsThe fixed cutting spacing of stage vertical direction, m are dumped for fixed advance rate；nhThe stage is dumped for fixed advance rate Number of entities after vertical cutting, dimensionless.(note: n=nhWhen formula in " % " be take the remainder operation)；
    6-2, it calculates i-th and dumps the stage and dump slope surface height Hi:
    6-3, determining that dumping the stage dumps height HiAfterwards, it is also necessary to determine and dump stage center-of-mass coordinate, thus need to dumping rank N-th of entity exists in vertical direction in sectionIt is cut at height.Stage casting absolute altitude or less is dumped after cutting at i-th Its mass center of part x, y-coordinate is constant, and z coordinate is that the coordinate of entity mass center before cutting subtractsWherein HnFor The fixed advance rate of n-th (n >=2) layer dump stage cutting entity thickness (" % " is complementation in formula, " and " be to divide exactly):
    Dump after 6-4, cutting at i-th the stage dump the solid volume of absolute altitude or less part for solid volume before cutting multiplied by/Hn；
    6-5, i-th is dumped the volume of each entity in phase range as weight the can be obtained to the weighting of its center-of-mass coordinate I dump stage center-of-mass coordinate
  5. The long-range casting planing method of row's opencut and computing system in 5., include basic data acquisition, calculation processing, result output Three modules.It is characterized in that being:
    As a result the absolute altitude that dumps for respectively dumping the stage that output module is exported according to calculation processing module, respectively dumps and consolidates locating for the stage Determine advance rate dump the stage, respectively dump the stage dump slope surface apart from locating fixation advance rate dump the stage dump slope surface away from From that draws out internal dumping in mapping software (including but not limited to 3DMine software) plot space respectively dumps stage step Bench crest, step bench toe.
CN201710600403.3A 2017-07-21 2017-07-21 The long-range casting planing method of interior row's opencut and computing system Pending CN109284848A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710600403.3A CN109284848A (en) 2017-07-21 2017-07-21 The long-range casting planing method of interior row's opencut and computing system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710600403.3A CN109284848A (en) 2017-07-21 2017-07-21 The long-range casting planing method of interior row's opencut and computing system

Publications (1)

Publication Number Publication Date
CN109284848A true CN109284848A (en) 2019-01-29

Family

ID=65184743

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710600403.3A Pending CN109284848A (en) 2017-07-21 2017-07-21 The long-range casting planing method of interior row's opencut and computing system

Country Status (1)

Country Link
CN (1) CN109284848A (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110751315A (en) * 2019-09-23 2020-02-04 中南大学 Man-machine interactive line selection method, system and controller for strip mine road system
CN112214867A (en) * 2020-08-28 2021-01-12 辽宁工程技术大学 Strip mine mining boundary and mining program collaborative optimization method under complex coal seam condition
CN112302657A (en) * 2020-09-18 2021-02-02 中煤科工集团沈阳设计研究院有限公司 Truck transportation discharge method for strip mine
CN113344328A (en) * 2021-05-07 2021-09-03 北方魏家峁煤电有限责任公司 Method and device for delimiting outer soil discharge field of opencast coal mine
CN114266379A (en) * 2021-11-22 2022-04-01 中煤科工集团沈阳设计研究院有限公司 Short-term internal discharge method for truck transportation of strip mine

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102966355A (en) * 2012-12-06 2013-03-13 马钢集团设计研究院有限责任公司 Method for creating internal dumping conditions of strip mine
CN104082056A (en) * 2014-07-01 2014-10-08 中国矿业大学(北京) Screening method for replacing surface soil of strip mine reclamation
US20150206433A1 (en) * 2014-01-21 2015-07-23 Hitachi Construction Machinery Co., Ltd. Vehicle control system
CN105023189A (en) * 2015-07-29 2015-11-04 中国神华能源股份有限公司 Land reclamation monitoring method and device for opencast mine area

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102966355A (en) * 2012-12-06 2013-03-13 马钢集团设计研究院有限责任公司 Method for creating internal dumping conditions of strip mine
US20150206433A1 (en) * 2014-01-21 2015-07-23 Hitachi Construction Machinery Co., Ltd. Vehicle control system
CN104082056A (en) * 2014-07-01 2014-10-08 中国矿业大学(北京) Screening method for replacing surface soil of strip mine reclamation
CN105023189A (en) * 2015-07-29 2015-11-04 中国神华能源股份有限公司 Land reclamation monitoring method and device for opencast mine area

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
万蒙蒙: "基于露天矿斜坡道的排土规划研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110751315A (en) * 2019-09-23 2020-02-04 中南大学 Man-machine interactive line selection method, system and controller for strip mine road system
CN110751315B (en) * 2019-09-23 2023-09-19 中南大学 Man-machine interaction type route selection method, system and controller for strip mine road system
CN112214867A (en) * 2020-08-28 2021-01-12 辽宁工程技术大学 Strip mine mining boundary and mining program collaborative optimization method under complex coal seam condition
CN112214867B (en) * 2020-08-28 2024-04-02 辽宁工程技术大学 Collaborative optimization method for mining boundaries and mining procedures of strip mine under complex coal seam condition
CN112302657A (en) * 2020-09-18 2021-02-02 中煤科工集团沈阳设计研究院有限公司 Truck transportation discharge method for strip mine
CN112302657B (en) * 2020-09-18 2024-05-17 中煤科工集团沈阳设计研究院有限公司 Method for discharging truck transportation of strip mine
CN113344328A (en) * 2021-05-07 2021-09-03 北方魏家峁煤电有限责任公司 Method and device for delimiting outer soil discharge field of opencast coal mine
CN114266379A (en) * 2021-11-22 2022-04-01 中煤科工集团沈阳设计研究院有限公司 Short-term internal discharge method for truck transportation of strip mine

Similar Documents

Publication Publication Date Title
CN109284848A (en) The long-range casting planing method of interior row&#39;s opencut and computing system
CN111814369B (en) Strip division method capable of accurately calculating force inclination angle between strips and stability safety coefficient of side slope
Tung et al. Optimal risk‐based design of flood levee systems
CN112214867A (en) Strip mine mining boundary and mining program collaborative optimization method under complex coal seam condition
CN112149314B (en) Multi-sand reservoir capacity silt flushing simulation method based on virtual reservoir capacity correction
CN110068573A (en) Based on photogrammetric and particle contact theory blast fragmentation size measurement method
CN112035932B (en) Arch dam intelligent progress simulation method
CN111767602B (en) High arch dam progress simulation method based on Internet of things
CN105804739A (en) Prediction method and prediction system for gas emission quantity
CN108978573A (en) A kind of method of terrain data quick visualization auxiliary river bed change research
CN112364474A (en) Mining scheme optimization method based on strip mine partition mining process
CN109740183A (en) Tunnel tunnel face estimation of stability model and benching tunnelling method excavation height design method
Stevanović et al. Application of MCDA in selection of different mining methods and solutions
CN109190294B (en) Design and optimization method for slope cutting and reinforcing scheme of huge thick landslide under reservoir water level change condition
Adamchuk et al. Substantiation of the method of determination the open-cast mine final con-tours taking into account the transport parameters
Wen et al. Efficiency analysis of the marine economy in the Guangdong–Hong Kong–Macao greater bay area based on a DEA model
CN109978256B (en) Multi-station earth and stone side loading and transporting machine optimal configuration method
CN110751315A (en) Man-machine interactive line selection method, system and controller for strip mine road system
CN115130169A (en) Comprehensive site selection method and system for bridge abutment in mountainous area
CN110728438B (en) Stock ground reserve analysis and evaluation method
JP2003020648A (en) Sedimentary construction method for water bottom ground
CN110852551B (en) Stock yard reserve calculation method
CN116399211B (en) Method for calculating slope surface sliding based on two-point Beidou monitoring and positioning
CN117145464A (en) Method and device for predicting sweet spots of tight reservoir based on sediment transport system theory
CN117057633B (en) Method and system for determining grade index of low-grade multi-metal ore

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20190129