WO2020006818A1 - Method for area-based coal rock water injection seepage-damage-stress coupling value simulation - Google Patents

Method for area-based coal rock water injection seepage-damage-stress coupling value simulation Download PDF

Info

Publication number
WO2020006818A1
WO2020006818A1 PCT/CN2018/100462 CN2018100462W WO2020006818A1 WO 2020006818 A1 WO2020006818 A1 WO 2020006818A1 CN 2018100462 W CN2018100462 W CN 2018100462W WO 2020006818 A1 WO2020006818 A1 WO 2020006818A1
Authority
WO
WIPO (PCT)
Prior art keywords
coal
simulation
calculation
determine whether
seepage
Prior art date
Application number
PCT/CN2018/100462
Other languages
French (fr)
Chinese (zh)
Inventor
周刚
尹文婧
张文政
魏星
张国宝
Original Assignee
山东科技大学
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 山东科技大学 filed Critical 山东科技大学
Priority to RU2020110457A priority Critical patent/RU2743121C1/en
Publication of WO2020006818A1 publication Critical patent/WO2020006818A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • E21B49/006Measuring wall stresses in the borehole
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21CMINING OR QUARRYING
    • E21C41/00Methods of underground or surface mining; Layouts therefor
    • E21C41/16Methods of underground mining; Layouts therefor
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Definitions

  • the present invention relates to the field of mine rock mechanics, and in particular, to a numerical simulation method for coupling water seepage-damage-stress coupling of coal rock masses containing geological structures.
  • Coal seam water injection is the process of injecting pressurized water and aqueous solution into the coal body through drilling.
  • the wetting effect of water makes the coal body plastic and the brittleness weakens.
  • many brittle fractures become plastic deformation, which greatly reduces it.
  • the possibility that the coal body is broken into dust particles is reduced, and the amount of coal dust generated is reduced.
  • the whole process involves many disciplines such as computational fluid mechanics, fracture rock mechanics and fluid-structure coupling. The process of numerical simulation is more complicated.
  • the object of the present invention is to provide a numerical simulation method for seepage-damage-stress coupled numerical simulation of water injection in coal and rock masses containing geological structures, which can more accurately simulate coal body damage during the water injection process And water transport laws, providing data for coal seam water injection.
  • the solution of the present invention includes:
  • a coupled numerical simulation method for seepage-damage-stress coupling of water injection in coal and rock masses comprising the following steps:
  • Step 1 Establish a coal and rock mass model including faults based on geological survey results
  • Step 2 Scan the corresponding coal body in the coal rock body obtained from the geological survey, and combine the FDK 3D reconstruction algorithm to construct a 3D digital core model with real pore structure characteristics;
  • Step 3 In the boundary element environment, execute the "mining fissure" generation algorithm to determine whether it is a mining affected area, that is, determine whether the distance from the mining face is less than 80m. If it is a mining affected area, go to step 4; if it is not mining Move the affected area, go to step five;
  • Step 4 Determine whether the cross-sectional area of the watershed is greater than 30 ⁇ m 2. If the cross-sectional area of the watershed is greater than 30 ⁇ m 2 , perform NS calculation processing, store and output the result; if less than 30 ⁇ m 2 , perform Darcy calculation processing, store and output the result;
  • Step 5 Determine whether the cross-sectional area of the watershed is greater than 30 ⁇ m 2. If the cross-sectional area of the watershed is greater than 30 ⁇ m 2 , go to step 6; if less than 30 ⁇ m 2 , go to step 12;
  • Step 6 Based on the C # language, the programming generates initial stress concentration points based on burial depth, geological conditions, lithology of overlying strata, and coal seam inclination.
  • Step 7 Mesh the 3D digital core model processed in Step 6;
  • Step 8 Calculate the tensile and shearing moments of the coal body grid, compare the tensile and shearing moments of the grid with the tensile and shearing moments, and find the most vulnerable grid set; The grid of moments is marked as "Native Crack", and then step 9 is performed; if there is no grid with a tensile-shear moment greater than the tensile-shear moment, the calculation is terminated and step 11 is performed;
  • Step 9 Replace the material in the "primary fissure" grid with the gas in the coal seam;
  • Step 10 Determine whether the total area of the "primary fissures" after the processing in step 9 exceeds twice the fault area. If not, generate a new model after deleting the "primary fissure" grid, and return to repeat step 8. If it exceeds , The calculation is terminated and step 11 is executed;
  • Step 11 Construct a model of two types of coal seam fractures that ultimately include primary fractures in the geological structure and fractures generated during the calculation and re-mesh the mesh, derive the STL general geometry, perform N-S calculation processing, and store and output the results;
  • Step 12 In the boundary element environment, after coupling the custom equation to calculate the seepage field and stress parameters, go to Step 13.
  • Step 13 In the meshless simulation environment, determine whether there are points where the tensile-shear moment is greater than the relevant parameters of the coal and rock mass. If so, connect such points in turn, mark the closed area as "invalid coal and rock mass", and execute Step 14; if not, go to Step 15;
  • Step 14 Add the water inlet conditions to the marked "invalid coal and rock mass" boundary in step 13.
  • Step 15 Perform turbulence simulation separately, calculate the stress distribution, and store the simulation results of the node;
  • Step 16 Determine whether the fracture separation from the simulation result in step 15 reaches the surface of the coal body. If not, go to step 17; if it reaches the surface of the coal body, go to step 18.
  • Step 17 Check whether the accumulated storage time reaches the preset simulation time. If not, go back to Step 12. If the calculation time is reached, go to Step 18.
  • Step 18 Stop the calculation in the meshless simulation environment. After calculating turbulence and stress only in the boundary element environment, perform Darcy calculation processing, store and output the results;
  • Step 19 The output results of Step 4, Step 11 and Step 18 are integrated and output as an independent file to obtain a quantitative statistical result.
  • the numerical simulation method for the coupled seepage-damage-stress numerical simulation of coal-rock mass injection includes:
  • Step A N-S initialization
  • Step B Calculate water pressure and gas pressure
  • Step C determine whether the water pressure is greater than the gas pressure, and if it is greater, change the next grid material to water, and execute step D; if it is less, directly perform step D;
  • Step D Determine whether the calculation time setting is reached. If the calculation time is not reached, return to step A. If the calculation time is reached, directly store and output the result.
  • the coupled numerical simulation method for seepage-damage-stress coupled injection of coal and rock mass, wherein the Darcy calculation process includes:
  • Step E Darcy initialization
  • Step F Calculate water pressure and gas pressure
  • Step J Determine whether the water pressure is greater than the gas pressure, and if it is greater, change the next grid material to water, and perform step H; if it is less, directly perform step H;
  • Step H Determine whether the calculation time setting is reached. If the calculation time is not reached, return to step E. If the calculation time is reached, directly store and output the result.
  • step 1 specifically further includes: reserve "fault start point” and "fault end point” in the process of establishing a coal rock mass model including fault ",” Fault turning point “,” Azimuth difference ",” Drop “and” Dip angle "are provided for user input.
  • the method for numerically coupling seepage-damage-stress coupled numerical simulation of coal rock mass injection wherein the above-mentioned step dimer further includes: using an image processing algorithm to perform a filtering operation on a three-dimensional digital core model, smoothing the model edges, and then implementing based on the threshold. Data segmentation was used to obtain the micropore structure of the coal body. The island pores with poor connectivity were eliminated, and the final numerical model of coal body pores in STL format was derived.
  • the numerical simulation method for coupled seepage-damage-stress numerical simulation of coal-rock mass injection wherein the above-mentioned quantitative statistical results include simulation results of hydraulic pressure field, simulation results of gas pressure, simulation results of coal pressure, simulation results of water flow velocity, and simulation of gas flow velocity result.
  • the present invention provides a numerical simulation method for coupling water seepage-damage-stress coupled numerical simulation of coal rock masses containing geological structures, using two mathematical models of Navier-Stokes equation and Darcy's law, respectively, based on two methods: boundary element method and meshless method.
  • the boundary element method is used to simulate the seepage process, and the meshless method is used to simulate the coal body fracture process.
  • the advantages of the simulation method can more accurately simulate the coal body damage and water movement during the regional injection of coal and rock masses, providing a solid data foundation for coal seam water injection, and greatly reducing the possibility of coal bodies being broken into dust particles. It reduces the amount of coal dust and ensures the safety of coal seam mining.
  • FIG. 1 is a schematic diagram of a numerical simulation method for seepage-damage-stress coupling of water injection in a coal rock mass according to the present invention
  • FIG. 3 is a schematic diagram of a Darcy calculation process in the present invention.
  • FIG. 5 is a schematic view of a coal body model after removing islands in the present invention.
  • FIG. 6 is a simulation diagram of water pressure results in the present invention.
  • FIG. 8 is a simulation diagram of coal pressure results of the present invention.
  • FIG. 9 is a simulation diagram of the water flow velocity result of the present invention.
  • FIG. 10 is a simulation diagram of gas flow velocity results of the present invention.
  • the present invention provides a numerical simulation method for coupling water seepage-damage-stress coupled numerical simulation of coal rock mass containing geological structures.
  • the present invention is further described in detail below. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.
  • the invention provides a numerical simulation method for seepage-damage-stress coupling numerical simulation of water injection in coal and rock bodies, which includes the following steps:
  • Step 1 Establish a coal and rock mass model including faults based on geological survey results
  • Step 2 Scan the corresponding coal body in the coal rock body obtained from the geological survey, and combine the FDK 3D reconstruction algorithm to construct a 3D digital core model with real pore structure characteristics;
  • Step 3 In the boundary element environment, execute the "mining fissure" generation algorithm to determine whether it is a mining affected area, that is, determine whether the distance from the mining face is less than 80m. If it is a mining affected area, go to step 4; if it is not mining Move the affected area, go to step five;
  • Step 4 Determine whether the cross-sectional area of the watershed is greater than 30 ⁇ m 2. If the cross-sectional area of the watershed is greater than 30 ⁇ m 2 , as shown in FIG. 2, perform NS calculation processing, store and output the result; if less than 30 ⁇ m 2 , as shown in FIG. 3 , Then perform Darcy calculation processing, store and output the results;
  • Step 5 Determine whether the cross-sectional area of the watershed is greater than 30 ⁇ m 2. If the cross-sectional area of the watershed is greater than 30 ⁇ m 2 , go to step 6; if less than 30 ⁇ m 2 , go to step 12;
  • Step 6 Based on the C # language, the programming generates initial stress concentration points based on burial depth, geological conditions, lithology of overlying strata, and coal seam inclination.
  • Step 7 Mesh the 3D digital core model processed in Step 6;
  • Step 8 Calculate the tensile and shearing moments of the coal body grid, compare the tensile and shearing moments of the grid with the tensile and shearing moments, and find the most vulnerable grid set; if there is a tensile and shearing moment greater than the tensile and shearing, The grid of moments is marked as "Native Crack", and then step 9 is performed; if there is no grid with a tensile-shear moment greater than the tensile-shear moment, the calculation is terminated and step 11 is performed;
  • Step 9 Replace the material in the "primary fissure" grid with the gas in the coal seam;
  • Step 10 Determine whether the total area of the "primary fissures" after the processing in step 9 exceeds twice the fault area. If not, generate a new model after deleting the "primary fissure" grid, and return to repeat step 8. If it exceeds , The calculation is terminated and step 11 is executed;
  • Step 11 Construct a model of two types of coal seam fractures that ultimately include primary fractures in the geological structure and fractures generated during calculation, and re-mesh the mesh to derive the STL general geometry for N-S calculation processing, and store and output the results;
  • Step 12 In the boundary element environment, after coupling the custom equation to calculate the seepage field and stress parameters, go to Step 13.
  • Step 13 In the meshless simulation environment, determine whether there are points where the tensile-shear moment is greater than the relevant parameters of the coal and rock mass. If so, connect such points in turn, mark the closed area as "invalid coal and rock mass", and execute Step 14; if not, go to Step 15;
  • Step 14 Add the water inlet conditions to the marked "invalid coal and rock mass" boundary in step 13.
  • Step 15 Perform turbulence simulation separately, calculate the stress distribution, and store the simulation results of the node;
  • Step 16 Determine whether the fracture separation from the simulation result in step 15 reaches the surface of the coal body. If not, go to step 17; if it reaches the surface of the coal body, go to step 18.
  • Step 17 Check whether the accumulated storage time reaches the preset simulation time. If not, go back to Step 12. If the calculation time is reached, go to Step 18.
  • Step 18 Stop the calculation in the meshless simulation environment. After calculating turbulence and stress only in the boundary element environment, perform Darcy calculation processing, store and output the results;
  • Step 19 The output results of Step 4, Step 11 and Step 18 are integrated and output as an independent file to obtain a quantitative statistical result.
  • the N-S calculation process includes:
  • Step A N-S initialization
  • Step B Calculate water pressure and gas pressure
  • Step C determine whether the water pressure is greater than the gas pressure, and if it is greater, change the next grid material to water, and execute step D; if it is less, directly perform step D;
  • Step D Determine whether the calculation time setting is reached. If the calculation time is not reached, return to step A. If the calculation time is reached, directly store and output the result.
  • Darcy calculation processing includes:
  • Step E Darcy initialization
  • Step F Calculate water pressure and gas pressure
  • Step J determine whether the water pressure is greater than the gas pressure, and if it is greater, change the next grid material to water, and execute step H; if it is less, directly execute step H;
  • Step H Determine whether the calculation time setting is reached. If the calculation time is not reached, return to step E. If the calculation time is reached, directly store and output the result.
  • step 1 specifically includes: reserve "fault start point”, “fault end point”, “fault turning point”, “azimuth difference”, “fall difference” and Six parameters of "tilt angle" are provided for user input.
  • step dimer also includes: using an image processing algorithm to perform a filtering operation on the three-dimensional digital core model, smoothing the edges of the model, and then performing data segmentation based on the threshold to obtain the microscopic pore structure of the coal body, and removing island pores with poor connectivity,
  • the final numerical model of coal pores in STL format was derived.
  • the above-mentioned quantitative statistical results include simulation results of water pressure field, simulation results of gas pressure, simulation results of coal pressure, simulation results of water flow velocity and simulation results of gas flow velocity.
  • a coupled numerical simulation method for seepage-damage-stress coupling of water injection in a coal rock mass containing a geological structure includes the following steps:
  • Step 1 Based on the geological survey results, establish a coal and rock mass model including faults;
  • coal seam model is programmed and programmed, in which six parameters of "starting point of the fault", “end point of the fault”, “turning point of the fault”, “azimuth difference”, “fall”, and “dip angle” are reserved for user input. ;
  • Step 2 CT array model
  • Step 3 In the boundary element environment, execute the "mining fissure" generation algorithm to determine whether it is a mining affected area, that is, determine whether the distance from the mining surface is less than 80m. If it is a mining-affected area, go to step 4; if it is not a mining-affected area, go to step 5.
  • Step 4 Determine whether the cross-sectional area of the watershed is greater than 30 ⁇ m 2. If the cross-sectional area of the watershed is greater than 30 ⁇ m 2 , go to step 6; if it is less than 30 ⁇ m 2 , go to step 7.
  • Step 5 Determine whether the cross-sectional area of the watershed is greater than 30 ⁇ m 2. If the cross-sectional area of the watershed is greater than 30 ⁇ m 2 , go to step 8; if it is less than 30 ⁇ m 2 , go to step 14.
  • Step 6 Perform N-S calculations, store and output the results.
  • Step 7. Perform Darcy calculation, save and output the result.
  • Step 8 Based on the C # language, program to generate initial stress concentration points based on burial depth, geological conditions, lithology of overlying strata, and inclination of coal seams.
  • Step 10 Calculate the tensile and shearing moments of the coal body grid, compare the tensile and shearing moments of the grid with the tensile and shearing moments, and find the most vulnerable grid set. If there is a mesh with a tensile-shear moment greater than the tensile-shear moment, mark it as "primary crack", and then perform step 11. If there is no mesh with a tensile-shear moment greater than the tensile-shear moment, the calculation is terminated and step 13 is performed.
  • Step 11 Replace the material in the "primary fissure" grid with the gas in the coal seam.
  • Step 12 Determine whether the total area of "primary fissures" exceeds twice the fault area. If not, generate a new model that deletes the original fissure mesh and return to repeat step 10. If it exceeds, terminate the calculation and perform step 10. three.
  • Step 13 Construct a model of the two types of coal seam fractures that ultimately include the primary fractures of the geological structure and the fractures generated during the calculation, and then re-mesh the grid to perform the sixth step after exporting the STL general geometry.
  • Step 14 In the boundary element environment, after coupling the custom equation to calculate the seepage field and stress parameters, go to step 15.
  • Step 15 In the meshless simulation environment, write the corresponding algorithm to determine whether there are points where the tensile-shear moment is greater than the relevant parameters of the coal rock mass. If so, connect such points and mark the closed area as "invalid coal rock mass.” ", Go to step sixteen; if not, go to step seventeen.
  • Step 16 Add the water inlet conditions to the marked "invalid coal rock body" boundary in step 15.
  • Step 17 Perform turbulence simulation separately, calculate the stress distribution, and store the simulation results of this node.
  • Step 18 Determine whether the fracture separation has reached the surface of the coal body. If not, go to step 19; if it has reached the surface of the coal body, go to step 20.
  • Step 19 Check whether the accumulated storage time reaches the preset simulation time. If not, go back to step 14. If the calculation time is reached, go to step 20.
  • Step 20 Stop the calculation in the meshless environment and perform step 7 only after calculating turbulence and stress in the boundary element environment.
  • Step 21 The integrated result is output and stored as an independent file, and the statistical results are quantified, such as the hydraulic pressure field simulation result diagram of FIG. 6, the gas pressure simulation result diagram of FIG. 7, and the coal body pressure simulation diagram of FIG. 8. The simulation results of the water flow velocity in FIG. 9 and the simulation results of the gas flow velocity in FIG. 10.
  • Step 61 N-S initialization.
  • Step 62 Calculate water pressure and gas pressure.
  • Step 63 Determine whether the water pressure is greater than the gas pressure. If it is greater, change the next grid material to water and perform step 64. If it is less, proceed directly to step 64.
  • Step 64 Determine whether the calculation time setting is reached. If the calculation time is not reached, return to step 61 again. If the calculation time is reached, directly store and output the result.
  • Step 71 Initialize Darcy.
  • Step 72 Calculate water pressure and gas pressure.
  • Step 73 Determine whether the water pressure is greater than the gas pressure. If it is greater, change the next grid material to water, and perform step 74. If it is less, directly perform step 74.
  • Step 74 Determine whether the calculation time setting is reached. If it cannot be reached, return to step 71 again. If the calculation time is reached, directly store and output the result.
  • Step 1 Based on the results of the geological survey, a coal and rock mass model containing faults is established, which specifically includes the following steps:
  • Step 2CT array model which specifically includes the following steps:
  • a nanoVoxel-2000 series-X-ray 3D microscope was used to scan the long-flame coal with an accuracy of 0.5 ⁇ m, and a 3D digital core model with real pore structure characteristics was constructed in combination with the FDK 3D reconstruction algorithm;
  • the image processing algorithm is used to filter the three-dimensional digital core model to reduce noise, smooth the model edges, and then perform volume rendering based on the gray value, and then implement data segmentation based on the threshold to obtain the micro-pore structure of the coal body;
  • Step 3 In the boundary element environment, execute the "mining fissure" generation algorithm to determine whether it is a mining affected area, that is, whether the distance from the mining face is less than 80m. If it is a mining-affected area, go to step 4; if it is not a mining-affected area, go to step 5. In order to complete the theory itself, the following assumptions are introduced:
  • Step 4 determines whether the cross-sectional area of the watershed is greater than 30 ⁇ m 2. If the cross-sectional area of the watershed is greater than 30 ⁇ m 2 , go to step 6; if it is smaller than 30 ⁇ m 2 , go to step 7.
  • Step 5 determines whether the cross-sectional area of the watershed is greater than 30 ⁇ m 2. If the cross-sectional area of the watershed is greater than 30 ⁇ m 2 , go to step 8; if it is smaller than 30 ⁇ m 2 , go to step 14.
  • Step 6 Perform N-S calculation, store and output the result.
  • N-S is used to describe the free flow of water in water injection holes and cracks. It is expressed in Cartesian coordinate system as:
  • t time
  • water pressure in coal body, MPa
  • v kinetic viscosity coefficient
  • u x , u y , u z axial unit mass forces of x, y, and z axes, respectively, mg ⁇ s -2
  • F x , F y , F z are components of external force, N
  • u, v, w are the velocity components of the fluid at point (x, y, z) at time t.
  • Step 7 performs Darcy calculation, stores and outputs the result.
  • P is the water pressure MPa in the coal body
  • V x and V y are the velocity components in the x and y axis directions
  • K is the permeability
  • is the dynamic viscosity of water Pa ⁇ s
  • g is the acceleration constant of gravity.
  • is the density of the liquid.
  • V is the velocity vector m / s of the coal bed gas movement;
  • gradp is the coal bed methane pore pressure gradient Pa / m;
  • is the coal bed methane pore pressure function;
  • k is the coal bed methane permeability coefficient is the absolute viscosity of the coal bed methane Pa ⁇ s;
  • Step 8 is based on the C # language and is programmed to generate initial stress concentration points based on burial depth, geological conditions, lithology of overlying strata and inclination of coal seams.
  • Step 9 Divides the grid, which includes the following steps:
  • Tetrahedral meshing is performed on the coal body geometric model, and the meshing accuracy is increased as much as possible within the range allowed by the computing power of the computer to ensure the accuracy and reliability of the fracture calculation results.
  • Step 10 In the boundary element environment, calculate and compare the tensile and shear moment values of each mesh, which specifically includes the following steps:
  • the Coulomb criterion states that the tensile shear failure value of the rock mass during the tensile-shear failure process in a certain plane exceeds the tensile-shear moment value, that is, the cohesion of the material and the constant plane method stress.
  • C + ⁇ tan ⁇ , where ⁇ is the tensile-shear moment value; C is the cohesive force or cohesive force, which is the tensile-shear strength without positive pressure; ⁇ is the internal friction angle; ⁇ is a fixed constant.
  • Step 11 If there is no mesh with a tensile-shear moment greater than the tensile-shear moment, the calculation is terminated and Step 13 is executed.
  • Step 11 replaces the material in the "primary fissure" grid with gas in the coal seam.
  • Step 12 Determine whether the total area of the "primary fissure" exceeds twice the fault area, which specifically includes the following steps:
  • Step 13 If not, generate a new model that deletes the original fissure mesh, and return to repeat Step 10. If it exceeds, terminate the calculation and execute Step 13.
  • Step 13 Build a model that finally includes two types of coal seam fractures, including primary fractures in the geological structure and fractures generated during the calculation.
  • the specific steps include the following steps:
  • Step 14 In the boundary element environment, coupled with custom equations to calculate the seepage field and stress parameters.
  • Step 15 In the meshless simulation environment, write an algorithm to determine whether there are points where the tensile-shear moment is greater than the relevant parameters of the coal and rock mass, which specifically include the following steps:
  • Step 17 If there is no point where the tensile-shear moment is greater than the relevant parameters of the coal rock mass, then Step 17 is continued.
  • Step 16 adds the water inlet conditions to the marked "invalid coal rock body" boundary in Step 15.
  • Step 17 Perform turbulence simulation separately, calculate the stress distribution, and store the simulation results of this node.
  • Step 18 Determine whether the fracture separation has reached the surface of the coal body, which specifically includes the following steps:
  • Step 19 is performed;
  • Step 20 If the fracture breaks off and reaches the surface of the coal body, perform Step 20.
  • Step 19 Whether the accumulated storage time reaches the preset simulation time, which specifically includes the following steps:
  • Step 20 stops the calculation in the meshless method environment, calculates turbulence and stress only in the boundary element environment, and then executes Step 7.
  • Step 21 outputs the integrated result and stores it as an independent file to quantify the statistical results, such as the hydraulic pressure field simulation result chart in Fig. 6, the gas pressure simulation result chart in Fig. 7, the coal pressure result simulation chart in Fig. 8, and Fig. 9

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Computer Hardware Design (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • Remote Sensing (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

Disclosed is a method for area-based coal rock water injection seepage-damage-stress coupling value simulation, comprising: establishing a coal seam model according to a geological survey result; performing area division and respectively simulating coal rock water injection seepage of a non-mining induced disturbance area and a mining induced disturbance area; in the non-mining induced disturbance area, performing programming to compare tension-shear torque values of meshes to determine whether deformation or even fracture occurs to coal rock under the action of water pressure stress; simultaneously simulating a coal fracture process using a meshfree method and simulating a seepage process using a boundary element method, wherein the advantages of the two simulation methods are combined on a micro scale; and in the mining induced disturbance area, respectively performing simulation using two mathematic models, i.e., the N-S equation and the Darcy's law, to accurately simulate coal damage and a water transfer law during an area-based coal rock water injection process. The method provides a solid data basis for coal seam water injection, greatly reduces the possibility that coal is crushed into dust particles, decreases the production amount of coal dust, and ensures the security of coal seam mining.

Description

一种煤岩体分区注水渗流-损伤-应力耦合数值模拟方法Numerical simulation method for seepage-damage-stress coupling of water injection in coal and rock mass 技术领域Technical field
本发明涉及矿山岩石力学领域,尤其涉及一种含地质构造的煤岩体分区注水渗流-损伤-应力耦合数值模拟方法。The present invention relates to the field of mine rock mechanics, and in particular, to a numerical simulation method for coupling water seepage-damage-stress coupling of coal rock masses containing geological structures.
背景技术Background technique
煤层注水是通过钻孔,将压力水和水溶液注入煤体的过程,水的湿润作用使煤体塑性增强,脆性减弱,当煤体受外力作用时,许多脆性破碎变为塑性形变,因而大量减少了煤体被破碎为尘粒的可能性,降低了煤尘的产生量。整个过程涉及计算流体力学、断裂岩石力学与流固耦合等多门学科,对其进行数值模拟的过程较为复杂。Coal seam water injection is the process of injecting pressurized water and aqueous solution into the coal body through drilling. The wetting effect of water makes the coal body plastic and the brittleness weakens. When the coal body is exposed to external forces, many brittle fractures become plastic deformation, which greatly reduces it. The possibility that the coal body is broken into dust particles is reduced, and the amount of coal dust generated is reduced. The whole process involves many disciplines such as computational fluid mechanics, fracture rock mechanics and fluid-structure coupling. The process of numerical simulation is more complicated.
目前针对煤层注水的数值模拟研究大都以宏观角度进行,通过模拟宏观层面的渗流速度场及渗流压力场等研究水分在煤体内的运移规律。但作为一种典型的多孔介质材料,煤体润湿的本质是水分进入煤体内的众多细微孔隙,宏观角度的模拟无法复现这一过程以供深入研究。且现有技术中针对水分渗流的研究大都基于有限元或离散元之一单独进行分析,但是有限元难以实现煤体断裂,离散元无法准确描述渗流过程与水分增量数据。At present, most of the numerical simulation studies on coal seam water injection are carried out from a macro perspective. The simulation of the seepage velocity field and seepage pressure field at the macro level is used to study the water transport law in the coal body. However, as a typical porous media material, the essence of coal body wetting is that water enters many fine pores in the coal body, and the simulation of macroscopic angle cannot reproduce this process for further research. In addition, most of the researches on water seepage in the prior art are based on separate analysis of finite element or discrete element, but it is difficult for finite element to achieve coal fracture, and discrete element cannot accurately describe the seepage process and moisture increment data.
因此,现有技术有待于更进一步的改进和发展。Therefore, the existing technology needs to be further improved and developed.
发明内容Summary of the invention
鉴于上述现有技术的不足,本发明的目的在于提供一种含地质构造的煤岩体分区注水渗流-损伤-应力耦合数值模拟方法,更准确的模拟出煤岩体分区注水过程中煤体损伤及水分的运移规律,为煤层注水提供数据。In view of the above-mentioned shortcomings of the prior art, the object of the present invention is to provide a numerical simulation method for seepage-damage-stress coupled numerical simulation of water injection in coal and rock masses containing geological structures, which can more accurately simulate coal body damage during the water injection process And water transport laws, providing data for coal seam water injection.
为解决上述技术问题,本发明方案包括:To solve the above technical problems, the solution of the present invention includes:
一种煤岩体分区注水渗流-损伤-应力耦合数值模拟方法,其包括以下步骤:A coupled numerical simulation method for seepage-damage-stress coupling of water injection in coal and rock masses, comprising the following steps:
步骤一,依据地质勘测结果,建立包含断层的煤岩体模型;Step 1: Establish a coal and rock mass model including faults based on geological survey results;
步骤二,对地质勘测得到之煤岩体中的相应煤体进行扫描,并结合FDK三维重建算法构建具有真实孔隙结构特征的三维数字岩心模型;Step 2: Scan the corresponding coal body in the coal rock body obtained from the geological survey, and combine the FDK 3D reconstruction algorithm to construct a 3D digital core model with real pore structure characteristics;
步骤三、在边界元环境下,执行“采动裂隙”生成算法,判断是否为采动影响区,即判断距离采掘面是否小于80m;如果是采动影响区,则执行步骤四;如果不是采动影响区,则执行步骤五;Step 3: In the boundary element environment, execute the "mining fissure" generation algorithm to determine whether it is a mining affected area, that is, determine whether the distance from the mining face is less than 80m. If it is a mining affected area, go to step 4; if it is not mining Move the affected area, go to step five;
步骤四、判断流域的截面积是否大于30μm 2,若流域截面积大于30μm 2,则进行N-S计算处理,储存并输出结果;若小于30μm 2,则进行Darcy计算处理,储存并输出结果; Step 4. Determine whether the cross-sectional area of the watershed is greater than 30 μm 2. If the cross-sectional area of the watershed is greater than 30 μm 2 , perform NS calculation processing, store and output the result; if less than 30 μm 2 , perform Darcy calculation processing, store and output the result;
步骤五、判断流域的截面积是否大于30μm 2,若流域截面积大于30μm 2,则执行步骤六,若小于30μm 2,则执行步骤十二; Step 5. Determine whether the cross-sectional area of the watershed is greater than 30 μm 2. If the cross-sectional area of the watershed is greater than 30 μm 2 , go to step 6; if less than 30 μm 2 , go to step 12;
步骤六、以C#语言为基础,编程以埋深、地质条件、上覆岩层岩性及煤层倾角生成初始应力集中点。 Step 6. Based on the C # language, the programming generates initial stress concentration points based on burial depth, geological conditions, lithology of overlying strata, and coal seam inclination.
步骤七、对经过步骤六处理的三维数字岩心模型进行剖分网格;Step 7. Mesh the 3D digital core model processed in Step 6;
步骤八、通过计算煤体网格所受拉剪力矩,将网格所受拉剪力矩与抗拉剪力矩比对,寻找最易被破坏的网格集合;若存在拉剪力矩大于抗拉剪力矩的网格,则标记为“原生裂隙”,然后执行步骤九;若不存在拉剪力矩大于抗拉剪力矩的网格,则终止计算,执行步骤十一; Step 8. Calculate the tensile and shearing moments of the coal body grid, compare the tensile and shearing moments of the grid with the tensile and shearing moments, and find the most vulnerable grid set; The grid of moments is marked as "Native Crack", and then step 9 is performed; if there is no grid with a tensile-shear moment greater than the tensile-shear moment, the calculation is terminated and step 11 is performed;
步骤九、替换“原生裂隙”网格中的材料为煤层内气体;Step 9. Replace the material in the "primary fissure" grid with the gas in the coal seam;
步骤十、判断经过步骤九处理后之“原生裂隙”的总面积是否超过断层面积的两倍,如果没有超过,则生成删除“原生裂隙”网格后的新模型,返回重复步骤八;如果超过,则终止计算,执行步骤十一;Step 10. Determine whether the total area of the "primary fissures" after the processing in step 9 exceeds twice the fault area. If not, generate a new model after deleting the "primary fissure" grid, and return to repeat step 8. If it exceeds , The calculation is terminated and step 11 is executed;
步骤十一、构建最终包含地质构造原生裂隙、计算过程中生成裂隙两类煤层裂隙的模型并重新剖分网格,导出STL通用几何后进行N-S计算处理,储存并输出结果;Step 11: Construct a model of two types of coal seam fractures that ultimately include primary fractures in the geological structure and fractures generated during the calculation and re-mesh the mesh, derive the STL general geometry, perform N-S calculation processing, and store and output the results;
步骤十二、在边界元环境下,耦合自定义方程计算渗流场与应力参数后,执行步骤十三; Step 12. In the boundary element environment, after coupling the custom equation to calculate the seepage field and stress parameters, go to Step 13.
步骤十三、在无网格仿真环境中,判断是否存在拉剪力矩大于煤岩体相关参数的点,如果有,则将此类点依次连接,闭合区域标记为“无效煤岩体”,执行步骤十四;如果没有,则执行步骤十五;Step 13. In the meshless simulation environment, determine whether there are points where the tensile-shear moment is greater than the relevant parameters of the coal and rock mass. If so, connect such points in turn, mark the closed area as "invalid coal and rock mass", and execute Step 14; if not, go to Step 15;
步骤十四、将步骤十三中,标记的“无效煤岩体”边界重新添加水分入口条件; Step 14. Add the water inlet conditions to the marked "invalid coal and rock mass" boundary in step 13.
步骤十五、单独进行湍流模拟,并计算应力分布,并储存该节点的模拟结果;Step 15: Perform turbulence simulation separately, calculate the stress distribution, and store the simulation results of the node;
步骤十六、判断步骤十五中模拟结果的断裂脱离是否达到煤体表面,如果没有,则执行步骤十七;如果到达煤体表面,则执行步骤十八; Step 16. Determine whether the fracture separation from the simulation result in step 15 reaches the surface of the coal body. If not, go to step 17; if it reaches the surface of the coal body, go to step 18.
步骤十七、累计存储时间是否达到预设模拟时间,如果没有达到,则返回执行步骤十二;如果达到计算时间,则执行步骤十八;Step 17. Check whether the accumulated storage time reaches the preset simulation time. If not, go back to Step 12. If the calculation time is reached, go to Step 18.
步骤十八、停止无网格仿真环境中的运算,仅在边界元环境中计算湍流与应力后则进行Darcy计算处理,储存并输出结果;Step 18. Stop the calculation in the meshless simulation environment. After calculating turbulence and stress only in the boundary element environment, perform Darcy calculation processing, store and output the results;
步骤十九、将步骤四、步骤十一与步骤十八的输出结果整合后输出并存储为独立文件,得到量化统计结果。Step 19: The output results of Step 4, Step 11 and Step 18 are integrated and output as an independent file to obtain a quantitative statistical result.
所述的煤岩体分区注水渗流-损伤-应力耦合数值模拟方法,其中,N-S计算处理包括:The numerical simulation method for the coupled seepage-damage-stress numerical simulation of coal-rock mass injection, wherein the N-S calculation processing includes:
步骤A、N-S初始化;Step A, N-S initialization;
步骤B、计算水压和瓦斯压;Step B: Calculate water pressure and gas pressure;
步骤C、判断水压是否大于瓦斯压,若大于,则将下一网格材料变更为水,执行步骤D,若小于,则直接执行步骤D;Step C, determine whether the water pressure is greater than the gas pressure, and if it is greater, change the next grid material to water, and execute step D; if it is less, directly perform step D;
步骤D、判断是否达到计算时间设定,若不能达到,则返回重新执行步骤A,若达到计算时间,则直接储存并输出结果。Step D: Determine whether the calculation time setting is reached. If the calculation time is not reached, return to step A. If the calculation time is reached, directly store and output the result.
所述的煤岩体分区注水渗流-损伤-应力耦合数值模拟方法,其中,Darcy计算处理包括:The coupled numerical simulation method for seepage-damage-stress coupled injection of coal and rock mass, wherein the Darcy calculation process includes:
步骤E、Darcy初始化;Step E, Darcy initialization;
步骤F、计算水压和瓦斯压;Step F: Calculate water pressure and gas pressure;
步骤J、判断水压是否大于瓦斯压,若大于,则将下一网格材料变更为水,执行步骤H,若小于,则直接执行步骤H;Step J: Determine whether the water pressure is greater than the gas pressure, and if it is greater, change the next grid material to water, and perform step H; if it is less, directly perform step H;
步骤H、判断是否达到计算时间设定,若不能达到,则返回重新执行步骤E,若达到计算时间,则直接储存并输出结果。Step H: Determine whether the calculation time setting is reached. If the calculation time is not reached, return to step E. If the calculation time is reached, directly store and output the result.
所述的煤岩体分区注水渗流-损伤-应力耦合数值模拟方法,其中,上述步骤一具体的还包括:建立包含断层的煤岩体模型的过程中预留出“断层起点”、“断层终点”、“断层转折点”、“方位差”、“落差”与“倾角”六个参数供用户输入。The numerical simulation method for coupled injection and seepage-damage-stress simulation of coal and rock masses, wherein the above-mentioned step 1 specifically further includes: reserve "fault start point" and "fault end point" in the process of establishing a coal rock mass model including fault "," Fault turning point "," Azimuth difference "," Drop "and" Dip angle "are provided for user input.
所述的煤岩体分区注水渗流-损伤-应力耦合数值模拟方法,其中,上述步骤二聚体的还包括:利用图像处理算法对三维数字岩心模型进行滤波操作,平滑模型边缘,进而基于阈值实现数据分割得到煤体的微观孔隙结构,剔除连通性较差的孤岛孔隙,导出得到最终的STL格式煤体孔隙数字模型。The method for numerically coupling seepage-damage-stress coupled numerical simulation of coal rock mass injection, wherein the above-mentioned step dimer further includes: using an image processing algorithm to perform a filtering operation on a three-dimensional digital core model, smoothing the model edges, and then implementing based on the threshold. Data segmentation was used to obtain the micropore structure of the coal body. The island pores with poor connectivity were eliminated, and the final numerical model of coal body pores in STL format was derived.
所述的煤岩体分区注水渗流-损伤-应力耦合数值模拟方法,其中,上述量化统计结果包括水压场模拟结果、瓦斯压力模拟结果、煤体压力模拟结果、水分流速模拟结果与瓦斯流速模拟结果。The numerical simulation method for coupled seepage-damage-stress numerical simulation of coal-rock mass injection, wherein the above-mentioned quantitative statistical results include simulation results of hydraulic pressure field, simulation results of gas pressure, simulation results of coal pressure, simulation results of water flow velocity, and simulation of gas flow velocity result.
本发明提供了一种含地质构造的煤岩体分区注水渗流-损伤-应力耦合数值模拟方法,分别采用Navier-Stokes方程和Darcy定律两种数学模型,基于边界元和无网格法两种方法进行模拟;边界元法与有限元相比,具有单元个数少,数据准备简单等优点,用边界元法模拟渗流过程,采用无网格法模拟煤体断裂过程,在微观尺度上结合两种模拟方法的优点,可以更准确的模拟出煤岩体分区注水过程中煤体损伤及水分的运移规律,为煤层注水提供坚实的数据基础,在大量减少了煤体被破碎为尘粒的可能性,降低了煤尘的产生量,保证了煤层开采的安全性。The present invention provides a numerical simulation method for coupling water seepage-damage-stress coupled numerical simulation of coal rock masses containing geological structures, using two mathematical models of Navier-Stokes equation and Darcy's law, respectively, based on two methods: boundary element method and meshless method. Perform simulation; compared with finite element, the boundary element method has the advantages of fewer elements and simple data preparation. The boundary element method is used to simulate the seepage process, and the meshless method is used to simulate the coal body fracture process. The advantages of the simulation method can more accurately simulate the coal body damage and water movement during the regional injection of coal and rock masses, providing a solid data foundation for coal seam water injection, and greatly reducing the possibility of coal bodies being broken into dust particles. It reduces the amount of coal dust and ensures the safety of coal seam mining.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
图1为本发明中煤岩体分区注水渗流-损伤-应力耦合数值模拟方法的示意图;FIG. 1 is a schematic diagram of a numerical simulation method for seepage-damage-stress coupling of water injection in a coal rock mass according to the present invention;
图2为本发明中N-S计算处理的示意图;2 is a schematic diagram of N-S calculation processing in the present invention;
图3为本发明中Darcy计算处理的示意图;3 is a schematic diagram of a Darcy calculation process in the present invention;
图4为本发明中去除孤岛煤块示意图;4 is a schematic diagram of removing island coal blocks in the present invention;
图5为本发明中去除孤岛后的煤体模型示意图;5 is a schematic view of a coal body model after removing islands in the present invention;
图6为本发明中为水压结果模拟图;FIG. 6 is a simulation diagram of water pressure results in the present invention;
图7为本发明瓦斯压力结果模拟图;7 is a simulation diagram of gas pressure results of the present invention;
图8为本发明煤体压力结果模拟图;FIG. 8 is a simulation diagram of coal pressure results of the present invention; FIG.
图9为本发明水分流速结果模拟图;FIG. 9 is a simulation diagram of the water flow velocity result of the present invention; FIG.
图10为本发明瓦斯流速结果模拟图。FIG. 10 is a simulation diagram of gas flow velocity results of the present invention.
具体实施方式detailed description
本发明提供了一种含地质构造的煤岩体分区注水渗流-损伤-应力耦合数值模拟方法,为使本发明的目的、技术方案及效果更加清楚、明确,以下对本发明进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。The present invention provides a numerical simulation method for coupling water seepage-damage-stress coupled numerical simulation of coal rock mass containing geological structures. In order to make the object, technical scheme and effect of the present invention clearer and more specific, the present invention is further described in detail below. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.
本发明提供了一种煤岩体分区注水渗流-损伤-应力耦合数值模拟方法,如图1所示的,其包括以下步骤:The invention provides a numerical simulation method for seepage-damage-stress coupling numerical simulation of water injection in coal and rock bodies, which includes the following steps:
步骤一,依据地质勘测结果,建立包含断层的煤岩体模型;Step 1: Establish a coal and rock mass model including faults based on geological survey results;
步骤二,对地质勘测得到之煤岩体中的相应煤体进行扫描,并结合FDK三维重建算法构建具有真实孔隙结构特征的三维数字岩心模型;Step 2: Scan the corresponding coal body in the coal rock body obtained from the geological survey, and combine the FDK 3D reconstruction algorithm to construct a 3D digital core model with real pore structure characteristics;
步骤三、在边界元环境下,执行“采动裂隙”生成算法,判断是否为采动影响区,即判断距离采掘面是否小于80m;如果是采动影响区,则执行步骤四;如果不是采动影响区,则执行步骤五;Step 3: In the boundary element environment, execute the "mining fissure" generation algorithm to determine whether it is a mining affected area, that is, determine whether the distance from the mining face is less than 80m. If it is a mining affected area, go to step 4; if it is not mining Move the affected area, go to step five;
步骤四、判断流域的截面积是否大于30μm 2,若流域截面积大于30μm 2,如图2所示的,则进行N-S计算处理,储存并输出结果;若小于30μm 2,如图3所示的,则进行Darcy计算处理,储存并输出结果; Step 4. Determine whether the cross-sectional area of the watershed is greater than 30 μm 2. If the cross-sectional area of the watershed is greater than 30 μm 2 , as shown in FIG. 2, perform NS calculation processing, store and output the result; if less than 30 μm 2 , as shown in FIG. 3 , Then perform Darcy calculation processing, store and output the results;
步骤五、判断流域的截面积是否大于30μm 2,若流域截面积大于30μm 2,则执行步骤六,若小于30μm 2,则执行步骤十二; Step 5. Determine whether the cross-sectional area of the watershed is greater than 30 μm 2. If the cross-sectional area of the watershed is greater than 30 μm 2 , go to step 6; if less than 30 μm 2 , go to step 12;
步骤六、以C#语言为基础,编程以埋深、地质条件、上覆岩层岩性及煤层倾角生成初始应力集中点。 Step 6. Based on the C # language, the programming generates initial stress concentration points based on burial depth, geological conditions, lithology of overlying strata, and coal seam inclination.
步骤七、对经过步骤六处理的三维数字岩心模型进行剖分网格;Step 7. Mesh the 3D digital core model processed in Step 6;
步骤八、通过计算煤体网格所受拉剪力矩,将网格所受拉剪力矩与抗拉剪力矩比对,寻找最易被破坏的网格集合;若存在拉剪力矩大于抗拉剪力矩的网格,则标记为“原生裂隙”,然后执行步骤九;若不存在拉剪力矩大于抗拉剪力矩的网格,则终止计算,执行步骤十一; Step 8. Calculate the tensile and shearing moments of the coal body grid, compare the tensile and shearing moments of the grid with the tensile and shearing moments, and find the most vulnerable grid set; if there is a tensile and shearing moment greater than the tensile and shearing, The grid of moments is marked as "Native Crack", and then step 9 is performed; if there is no grid with a tensile-shear moment greater than the tensile-shear moment, the calculation is terminated and step 11 is performed;
步骤九、替换“原生裂隙”网格中的材料为煤层内气体;Step 9. Replace the material in the "primary fissure" grid with the gas in the coal seam;
步骤十、判断经过步骤九处理后之“原生裂隙”的总面积是否超过断层面积的两倍,如果没有超过,则生成删除“原生裂隙”网格后的新模型,返回重复步骤八;如果超过,则终止计算,执行步骤十一;Step 10. Determine whether the total area of the "primary fissures" after the processing in step 9 exceeds twice the fault area. If not, generate a new model after deleting the "primary fissure" grid, and return to repeat step 8. If it exceeds , The calculation is terminated and step 11 is executed;
步骤十一、构建最终包含地质构造原生裂隙、计算过程中生成裂隙两类煤层裂隙的模型,并重新剖分网格,导出STL通用几何后进行N-S计算处理,储存并输出结果;Step 11: Construct a model of two types of coal seam fractures that ultimately include primary fractures in the geological structure and fractures generated during calculation, and re-mesh the mesh to derive the STL general geometry for N-S calculation processing, and store and output the results;
步骤十二、在边界元环境下,耦合自定义方程计算渗流场与应力参数后,执行步骤十三; Step 12. In the boundary element environment, after coupling the custom equation to calculate the seepage field and stress parameters, go to Step 13.
步骤十三、在无网格仿真环境中,判断是否存在拉剪力矩大于煤岩体相关参数的点,如果有,则将此类点依次连接,闭合区域标记为“无效煤岩体”,执行步骤十四;如果没有,则执行步骤十五;Step 13. In the meshless simulation environment, determine whether there are points where the tensile-shear moment is greater than the relevant parameters of the coal and rock mass. If so, connect such points in turn, mark the closed area as "invalid coal and rock mass", and execute Step 14; if not, go to Step 15;
步骤十四、将步骤十三中,标记的“无效煤岩体”边界重新添加水分入口条件; Step 14. Add the water inlet conditions to the marked "invalid coal and rock mass" boundary in step 13.
步骤十五、单独进行湍流模拟,并计算应力分布,并储存该节点的模拟结果;Step 15: Perform turbulence simulation separately, calculate the stress distribution, and store the simulation results of the node;
步骤十六、判断步骤十五中模拟结果的断裂脱离是否达到煤体表面,如果没有,则执行步骤十七;如果到达煤体表面,则执行步骤十八; Step 16. Determine whether the fracture separation from the simulation result in step 15 reaches the surface of the coal body. If not, go to step 17; if it reaches the surface of the coal body, go to step 18.
步骤十七、累计存储时间是否达到预设模拟时间,如果没有达到,则返回执行步骤十二;如果达到计算时间,则执行步骤十八;Step 17. Check whether the accumulated storage time reaches the preset simulation time. If not, go back to Step 12. If the calculation time is reached, go to Step 18.
步骤十八、停止无网格仿真环境中的运算,仅在边界元环境中计算湍流与应力后则进行Darcy计算处理,储存并输出结果;Step 18. Stop the calculation in the meshless simulation environment. After calculating turbulence and stress only in the boundary element environment, perform Darcy calculation processing, store and output the results;
步骤十九、将步骤四、步骤十一与步骤十八的输出结果整合后输出并存储为独立文件,得到量化统计结果。Step 19: The output results of Step 4, Step 11 and Step 18 are integrated and output as an independent file to obtain a quantitative statistical result.
在本发明的另一较佳实施例中,N-S计算处理包括:In another preferred embodiment of the present invention, the N-S calculation process includes:
步骤A、N-S初始化;Step A, N-S initialization;
步骤B、计算水压和瓦斯压;Step B: Calculate water pressure and gas pressure;
步骤C、判断水压是否大于瓦斯压,若大于,则将下一网格材料变更为水,执行步骤D,若小于,则直接执行步骤D;Step C, determine whether the water pressure is greater than the gas pressure, and if it is greater, change the next grid material to water, and execute step D; if it is less, directly perform step D;
步骤D、判断是否达到计算时间设定,若不能达到,则返回重新执行步骤A,若达到计算时间,则直接储存并输出结果。Step D: Determine whether the calculation time setting is reached. If the calculation time is not reached, return to step A. If the calculation time is reached, directly store and output the result.
还有Darcy计算处理包括:Also Darcy calculation processing includes:
步骤E、Darcy初始化;Step E, Darcy initialization;
步骤F、计算水压和瓦斯压;Step F: Calculate water pressure and gas pressure;
步骤J、判断水压是否大于瓦斯压,若大于,则将下一网格材料变更为水,执行步骤H, 若小于,则直接执行步骤H;Step J: determine whether the water pressure is greater than the gas pressure, and if it is greater, change the next grid material to water, and execute step H; if it is less, directly execute step H;
步骤H、判断是否达到计算时间设定,若不能达到,则返回重新执行步骤E,若达到计算时间,则直接储存并输出结果。Step H: Determine whether the calculation time setting is reached. If the calculation time is not reached, return to step E. If the calculation time is reached, directly store and output the result.
更进一步的,上述步骤一具体的还包括:建立包含断层的煤岩体模型的过程中预留出“断层起点”、“断层终点”、“断层转折点”、“方位差”、“落差”与“倾角”六个参数供用户输入。Furthermore, the above-mentioned step 1 specifically includes: reserve "fault start point", "fault end point", "fault turning point", "azimuth difference", "fall difference" and Six parameters of "tilt angle" are provided for user input.
并且上述步骤二聚体的还包括:利用图像处理算法对三维数字岩心模型进行滤波操作,平滑模型边缘,进而基于阈值实现数据分割得到煤体的微观孔隙结构,剔除连通性较差的孤岛孔隙,导出得到最终的STL格式煤体孔隙数字模型。In addition, the above-mentioned step dimer also includes: using an image processing algorithm to perform a filtering operation on the three-dimensional digital core model, smoothing the edges of the model, and then performing data segmentation based on the threshold to obtain the microscopic pore structure of the coal body, and removing island pores with poor connectivity, The final numerical model of coal pores in STL format was derived.
上述量化统计结果包括水压场模拟结果、瓦斯压力模拟结果、煤体压力模拟结果、水分流速模拟结果与瓦斯流速模拟结果。The above-mentioned quantitative statistical results include simulation results of water pressure field, simulation results of gas pressure, simulation results of coal pressure, simulation results of water flow velocity and simulation results of gas flow velocity.
为了更进一步的描述本发明,以下列举更为详尽的实施例。In order to further describe the present invention, more detailed examples are listed below.
一种含地质构造的煤岩体分区注水渗流-损伤-应力耦合数值模拟方法,其包括以下步骤:A coupled numerical simulation method for seepage-damage-stress coupling of water injection in a coal rock mass containing a geological structure includes the following steps:
步骤一、依据地质勘测结果,建立包含断层的煤岩体模型;Step 1. Based on the geological survey results, establish a coal and rock mass model including faults;
依据实际地质勘测的结果,编程建立煤层模型,其中,预留出“断层起点”、“断层终点”、“断层转折点”、“方位差”、“落差”、“倾角”六个参数供用户输入;According to the results of actual geological surveys, a coal seam model is programmed and programmed, in which six parameters of "starting point of the fault", "end point of the fault", "turning point of the fault", "azimuth difference", "fall", and "dip angle" are reserved for user input. ;
步骤二、CT阵列模型;Step 2: CT array model;
对地质勘测得到之煤岩体中的相应煤体进行扫描实验,结合FDK三维重建算法构建具有真实孔隙结构特征的三维数字岩心模型。利用图像处理算法对三维数字岩心模型进行滤波操作,平滑模型边缘,进而基于阈值实现数据分割得到煤体的微观孔隙结构,剔除连通性较差的孤岛孔隙,如图4中非浅色区域,处理后的煤体孔隙结构如图5所示,导出即为最终的STL格式煤体孔隙数字模型。Scanning experiments were performed on the corresponding coal bodies in the coal and rock bodies obtained from the geological survey, and a 3D digital core model with real pore structure characteristics was constructed in combination with the FDK 3D reconstruction algorithm. The image processing algorithm is used to filter the 3D digital core model to smooth the edges of the model, and then the data is segmented based on the threshold to obtain the micro-pore structure of the coal body. The island pores with poor connectivity are eliminated, as shown in Figure 4. The resulting coal body pore structure is shown in Figure 5, which is derived as the final digital model of coal body pores in STL format.
步骤三、在边界元环境下,执行“采动裂隙”生成算法,判断是否为采动影响区,即判断距离采掘面是否小于80m。如果是采动影响区,则执行步骤四;如果不是采动影响区,则执行步骤五。Step 3: In the boundary element environment, execute the "mining fissure" generation algorithm to determine whether it is a mining affected area, that is, determine whether the distance from the mining surface is less than 80m. If it is a mining-affected area, go to step 4; if it is not a mining-affected area, go to step 5.
步骤四、判断流域的截面积是否大于30μm 2,若流域截面积大于30μm 2,则执行步骤六,若小于30μm 2,则执行步骤七。 Step 4. Determine whether the cross-sectional area of the watershed is greater than 30 μm 2. If the cross-sectional area of the watershed is greater than 30 μm 2 , go to step 6; if it is less than 30 μm 2 , go to step 7.
步骤五、判断流域的截面积是否大于30μm 2,若流域截面积大于30μm 2,则执行步骤八,若小于30μm 2,则执行步骤十四。 Step 5. Determine whether the cross-sectional area of the watershed is greater than 30 μm 2. If the cross-sectional area of the watershed is greater than 30 μm 2 , go to step 8; if it is less than 30 μm 2 , go to step 14.
步骤六、进行N-S计算,储存并输出结果。 Step 6. Perform N-S calculations, store and output the results.
步骤七、进行Darcy计算,储存并输出结果。Step 7. Perform Darcy calculation, save and output the result.
步骤八、以C#语言为基础,编程以埋深、地质条件、上覆岩层岩性及煤层倾角生成初始 应力集中点。 Step 8. Based on the C # language, program to generate initial stress concentration points based on burial depth, geological conditions, lithology of overlying strata, and inclination of coal seams.
步骤九、剖分网格。Step 9. Meshing.
步骤十、通过计算煤体网格所受拉剪力矩,将网格所受拉剪力矩与抗拉剪力矩比对,寻找最易被破坏的网格集合。若存在拉剪力矩大于抗拉剪力矩的网格,则标记为“原生裂隙”,然后执行步骤十一。若不存在拉剪力矩大于抗拉剪力矩的网格,则终止计算,执行步骤十三。Step 10. Calculate the tensile and shearing moments of the coal body grid, compare the tensile and shearing moments of the grid with the tensile and shearing moments, and find the most vulnerable grid set. If there is a mesh with a tensile-shear moment greater than the tensile-shear moment, mark it as "primary crack", and then perform step 11. If there is no mesh with a tensile-shear moment greater than the tensile-shear moment, the calculation is terminated and step 13 is performed.
步骤十一、替换“原生裂隙”网格中的材料为煤层内气体。Step 11: Replace the material in the "primary fissure" grid with the gas in the coal seam.
步骤十二、判断“原生裂隙”的总面积是否超过断层面积的两倍,如果没有超过,则生成删除原生裂隙网格的新模型,返回重复步骤十;如果超过,则终止计算,执行步骤十三。 Step 12. Determine whether the total area of "primary fissures" exceeds twice the fault area. If not, generate a new model that deletes the original fissure mesh and return to repeat step 10. If it exceeds, terminate the calculation and perform step 10. three.
步骤十三、构建最终包含地质构造原生裂隙、计算过程中生成裂隙两类煤层裂隙的模型,并重新剖分网格,导出STL通用几何后执行步骤六。Step 13: Construct a model of the two types of coal seam fractures that ultimately include the primary fractures of the geological structure and the fractures generated during the calculation, and then re-mesh the grid to perform the sixth step after exporting the STL general geometry.
步骤十四、在边界元环境下,耦合自定义方程计算渗流场与应力参数后,执行步骤十五。 Step 14. In the boundary element environment, after coupling the custom equation to calculate the seepage field and stress parameters, go to step 15.
步骤十五、在无网格仿真环境中,编写相应算法,判断是否存在拉剪力矩大于煤岩体相关参数的点,如果有,则将此类点连接,闭合区域标记为“无效煤岩体”,执行步骤十六;如果没有,则执行步骤十七。Step 15: In the meshless simulation environment, write the corresponding algorithm to determine whether there are points where the tensile-shear moment is greater than the relevant parameters of the coal rock mass. If so, connect such points and mark the closed area as "invalid coal rock mass." ", Go to step sixteen; if not, go to step seventeen.
步骤十六、将步骤十五中,标记的“无效煤岩体”边界重新添加水分入口条件。 Step 16. Add the water inlet conditions to the marked "invalid coal rock body" boundary in step 15.
步骤十七、单独进行湍流模拟,并计算应力分布,并储存该节点的模拟结果。Step 17: Perform turbulence simulation separately, calculate the stress distribution, and store the simulation results of this node.
步骤十八、判断断裂脱离是否达到煤体表面,如果没有,则执行步骤十九;如果到达煤体表面,则执行步骤二十。Step 18. Determine whether the fracture separation has reached the surface of the coal body. If not, go to step 19; if it has reached the surface of the coal body, go to step 20.
步骤十九、累计存储时间是否达到预设模拟时间,如果没有达到,则返回执行步骤十四;如果达到计算时间,则执行步骤二十。Step 19: Check whether the accumulated storage time reaches the preset simulation time. If not, go back to step 14. If the calculation time is reached, go to step 20.
步骤二十、停止无网格法环境中的运算,仅在边界元环境中计算湍流与应力后执行步骤七。 Step 20. Stop the calculation in the meshless environment and perform step 7 only after calculating turbulence and stress in the boundary element environment.
步骤二十一、将整合后的结果输出并存储为独立文件,量化统计结果,如图6的水压场模拟结果图、图7的瓦斯压力模拟结果图、图8的煤体压力结果模拟图、图9的水分流速模拟结果图与图10的瓦斯流速模拟结果图。Step 21: The integrated result is output and stored as an independent file, and the statistical results are quantified, such as the hydraulic pressure field simulation result diagram of FIG. 6, the gas pressure simulation result diagram of FIG. 7, and the coal body pressure simulation diagram of FIG. 8. The simulation results of the water flow velocity in FIG. 9 and the simulation results of the gas flow velocity in FIG. 10.
还包括N-S计算处理和Darcy计算处理的步骤:It also includes the steps of N-S calculation processing and Darcy calculation processing:
N-S计算处理步骤;N-S calculation processing steps;
步骤六一、N-S初始化。Step 61: N-S initialization.
步骤六二、计算水压和瓦斯压。Step 62: Calculate water pressure and gas pressure.
步骤六三、判断水压是否大于瓦斯压,若大于,则将下一网格材料变更为水,执行步骤六四,若小于,则直接执行步骤六四。Step 63: Determine whether the water pressure is greater than the gas pressure. If it is greater, change the next grid material to water and perform step 64. If it is less, proceed directly to step 64.
步骤六四、判断是否达到计算时间设定,若不能达到,则返回重新执行步骤六一,若达到计算时间,则直接储存并输出结果。Step 64: Determine whether the calculation time setting is reached. If the calculation time is not reached, return to step 61 again. If the calculation time is reached, directly store and output the result.
Darcy计算处理步骤;Darcy calculation processing steps;
步骤七一、Darcy初始化。Step 71: Initialize Darcy.
步骤七二、计算水压和瓦斯压。Step 72: Calculate water pressure and gas pressure.
步骤七三、判断水压是否大于瓦斯压,若大于,则将下一网格材料变更为水,执行步骤七四,若小于,则直接执行步骤七四。Step 73: Determine whether the water pressure is greater than the gas pressure. If it is greater, change the next grid material to water, and perform step 74. If it is less, directly perform step 74.
步骤七四、判断是否达到计算时间设定,若不能达到,则返回重新执行步骤七一,若达到计算时间,则直接储存并输出结果。Step 74: Determine whether the calculation time setting is reached. If it cannot be reached, return to step 71 again. If the calculation time is reached, directly store and output the result.
实施例1Example 1
Step 1依据地质勘测结果,建立包含断层的煤岩体模型,其具体包括如下步骤:Step 1 Based on the results of the geological survey, a coal and rock mass model containing faults is established, which specifically includes the following steps:
(1)在笛卡尔坐标系下,依据注水煤层高度、走向长度及倾向长度建立煤体几何模型,模型维度为三维;(1) In the Cartesian coordinate system, a geometric model of the coal body is established based on the height, strike length, and inclination length of the water injection coal seam, and the model dimension is three-dimensional;
(2)其中,预留出“断层起点”、“断层终点”、“断层转折点”、“方位差”、“落差”、“倾角”六个参数供用户输入。(2) Among them, six parameters of "starting point of the fault", "end point of the fault", "turning point of the fault", "azimuth difference", "fall" and "inclination angle" are reserved for user input.
(3)然后导出STL通用几何。(3) Then export the STL general geometry.
Step 2CT阵列模型,其具体包括如下步骤:Step 2CT array model, which specifically includes the following steps:
(1)采用精度0.5μm的nanoVoxel-2000系列-X射线三维显微镜对长焰煤进行扫描实验,结合FDK三维重建算法构建了具有真实孔隙结构特征的三维数字岩心模型;(1) A nanoVoxel-2000 series-X-ray 3D microscope was used to scan the long-flame coal with an accuracy of 0.5 μm, and a 3D digital core model with real pore structure characteristics was constructed in combination with the FDK 3D reconstruction algorithm;
(2)利用图像处理算法对三维数字岩心模型进行滤波操作,减少噪声,平滑模型边缘,然后基于灰度值,进行体渲染,进而基于阈值实现数据分割得到煤体的微观孔隙结构;(2) The image processing algorithm is used to filter the three-dimensional digital core model to reduce noise, smooth the model edges, and then perform volume rendering based on the gray value, and then implement data segmentation based on the threshold to obtain the micro-pore structure of the coal body;
(3)剔除连通性较差的孤岛孔隙,如图4中非浅蓝色区域(具体参见其他证明材料中的图4),处理后的煤体孔隙结构如图5所示(具体参见其他证明材料中的图5),导出即为最终的STL格式煤体孔隙数字模型(3) Remove the island pores with poor connectivity, as shown in the non-light blue area in Figure 4 (see Figure 4 in other certification materials for details), and the pore structure of the coal body after processing is shown in Figure 5 (see other certifications for details) Figure 5) in the material, which is derived as the final digital model of coal pores in STL format
Step 3在边界元环境下,执行“采动裂隙”生成算法,判断是否为采动影响区,即距离采掘面是否小于80m。如果是采动影响区,则执行步骤四;如果不是采动影响区,则执行步骤五。为了理论本身的完备性,引入以下假设:Step 3 In the boundary element environment, execute the "mining fissure" generation algorithm to determine whether it is a mining affected area, that is, whether the distance from the mining face is less than 80m. If it is a mining-affected area, go to step 4; if it is not a mining-affected area, go to step 5. In order to complete the theory itself, the following assumptions are introduced:
(1)忽略或不考虑基质岩体的透水性和储水性(因为相比较于“采动裂隙”,基质岩体的透水性和储水性能均特别弱)。(1) Ignore or disregard the water permeability and water storage capacity of the matrix rock mass (because the water permeability and water storage performance of the matrix rock mass are particularly weak compared to "mining fractures").
(2)“采动裂隙”渗流服从达西定律
Figure PCTCN2018100462-appb-000001
其中k=-b 2/12
(2) Seepage of "mining fracture" obeys Darcy's law
Figure PCTCN2018100462-appb-000001
Where k = -b 2/12
(3)“采动裂隙”变形规律服从Goodman节理模型。(3) The deformation law of "mining fracture" obeys the Goodman joint model.
Step 4判断流域的截面积是否大于30μm 2,若流域截面积大于30μm 2,则执行步骤六,若小于30μm 2,则执行步骤七。 Step 4 determines whether the cross-sectional area of the watershed is greater than 30 μm 2. If the cross-sectional area of the watershed is greater than 30 μm 2 , go to step 6; if it is smaller than 30 μm 2 , go to step 7.
Step 5判断流域的截面积是否大于30μm 2,若流域截面积大于30μm 2,则执行步骤八,若小于30μm 2,则执行步骤十四。 Step 5 determines whether the cross-sectional area of the watershed is greater than 30 μm 2. If the cross-sectional area of the watershed is greater than 30 μm 2 , go to step 8; if it is smaller than 30 μm 2 , go to step 14.
Step 6进行N-S计算,储存并输出结果。Step 6: Perform N-S calculation, store and output the result.
(1)N-S初始化。(1) N-S initialization.
以N-S来描述水分注水孔及裂隙中水分的自由流动,在笛卡尔坐标系下表示为:N-S is used to describe the free flow of water in water injection holes and cracks. It is expressed in Cartesian coordinate system as:
Figure PCTCN2018100462-appb-000002
Figure PCTCN2018100462-appb-000002
Figure PCTCN2018100462-appb-000003
Figure PCTCN2018100462-appb-000003
Figure PCTCN2018100462-appb-000004
Figure PCTCN2018100462-appb-000004
式中,t—时间;ρ—煤体中的水压,MPa;v—运动黏滞系数;u x、u y、u z—分别为x、y、z轴轴向单位质量力,mg·s -2;F x、F y、F z是外力的分量,N;
Figure PCTCN2018100462-appb-000005
—那勃勒算子;u、v、w是流体在t时刻,在点(x、y、z)处的速度分量。
In the formula, t—time; ρ—water pressure in coal body, MPa; v—kinetic viscosity coefficient; u x , u y , u z — axial unit mass forces of x, y, and z axes, respectively, mg · s -2 ; F x , F y , F z are components of external force, N;
Figure PCTCN2018100462-appb-000005
-Naples operator; u, v, w are the velocity components of the fluid at point (x, y, z) at time t.
(2)计算水压和瓦斯压。(2) Calculate water pressure and gas pressure.
(3)判断水压是否大于瓦斯压,若大于,则将下一网格材料变更为水,执行(4),若小于,则直接执行(4)。(3) Determine whether the water pressure is greater than the gas pressure. If it is greater, change the next grid material to water and perform (4). If it is less, directly perform (4).
(4)判断是否达到计算时间设定,若不能达到,则返回重新执行(1),若达到计算时间,则直接储存并输出结果。(4) Determine whether the calculation time setting is reached. If it cannot be reached, return to re-execute (1). If the calculation time is reached, directly store and output the result.
Step 7进行Darcy计算,储存并输出结果。Step 7 performs Darcy calculation, stores and outputs the result.
(1)Darcy初始化。(1) Darcy initialization.
(2)计算水压和瓦斯压。(2) Calculate water pressure and gas pressure.
Darcy定律描述水在煤体内的渗流运动,其微分方程在笛卡尔坐标系中表示为:Darcy's law describes the percolation motion of water in coal. Its differential equation is expressed in Cartesian coordinate system as:
Figure PCTCN2018100462-appb-000006
Figure PCTCN2018100462-appb-000006
Figure PCTCN2018100462-appb-000007
Figure PCTCN2018100462-appb-000007
式中,P为煤体中的水压MPa,V x与V y分别为x、y轴方向上的速度分量;K为渗透率,μ为水的动力粘度Pa·s,g为重力加速度常数;ρ为液体密度。 In the formula, P is the water pressure MPa in the coal body, V x and V y are the velocity components in the x and y axis directions; K is the permeability, μ is the dynamic viscosity of water Pa · s, and g is the acceleration constant of gravity. ; Ρ is the density of the liquid.
Darcy定律认为煤层内煤层气运动基本符合线性渗透规律,即:Darcy's law states that the movement of coalbed methane in coal seams basically conforms to the law of linear permeability, namely:
Figure PCTCN2018100462-appb-000008
Figure PCTCN2018100462-appb-000008
其中,V为煤层气运动的速度矢量m/s;gradp为煤层内煤层气孔隙压力梯度Pa/m;μ为煤层内煤层气孔隙压力函数;k为煤层气渗透性系数为煤层气的绝对黏度,Pa·s;Among them, V is the velocity vector m / s of the coal bed gas movement; gradp is the coal bed methane pore pressure gradient Pa / m; μ is the coal bed methane pore pressure function; k is the coal bed methane permeability coefficient is the absolute viscosity of the coal bed methane Pa · s;
瓦斯沿煤层的流动一般服从Darcy定律,即瓦斯的渗流速度与瓦斯压力梯度成正比:The flow of gas along the coal seam generally obeys Darcy's law, that is, the seepage velocity of gas is proportional to the pressure gradient of gas:
Figure PCTCN2018100462-appb-000009
Figure PCTCN2018100462-appb-000009
其中,q为瓦斯渗流速度,cm/s;λ为煤层透气系数,m 2/(MPa 2·d);μ为流体的动力黏度系数,Pa·s;k为渗透率,m 2
Figure PCTCN2018100462-appb-000010
为瓦斯压力梯度,P/cm。
Among them, q is gas seepage velocity, cm / s; λ is coal seam permeability coefficient, m 2 / (MPa 2 · d); μ is fluid dynamic viscosity coefficient, Pa · s; k is permeability, m 2 ;
Figure PCTCN2018100462-appb-000010
Is the gas pressure gradient, P / cm.
(3)判断水压是否大于瓦斯压,若大于,则将下一网格材料变更为水,执行(4),若小于,则直接执行(4)。(3) Determine whether the water pressure is greater than the gas pressure. If it is greater, change the next grid material to water and perform (4). If it is less, directly perform (4).
(4)判断是否达到计算时间设定,若不能达到,则返回重新执行(1),若达到计算时间,则直接储存并输出结果。(4) Determine whether the calculation time setting is reached. If it cannot be reached, return to re-execute (1). If the calculation time is reached, directly store and output the result.
Step 8以C#语言为基础,编程以埋深、地质条件、上覆岩层岩性及煤层倾角生成初始应力集中点。 Step 8 is based on the C # language and is programmed to generate initial stress concentration points based on burial depth, geological conditions, lithology of overlying strata and inclination of coal seams.
Step 9剖分网格,其具体包括如下步骤:Step 9 Divides the grid, which includes the following steps:
对煤体几何模型进行四面体网格剖分,在计算机运算力允许的范围内尽可能提高网格剖分精度,以保证裂隙计算结果的准确性与可靠性。Tetrahedral meshing is performed on the coal body geometric model, and the meshing accuracy is increased as much as possible within the range allowed by the computing power of the computer to ensure the accuracy and reliability of the fracture calculation results.
Step 10边界元环境下,计算并比较各网格的拉剪力矩值,其具体包括如下步骤:Step 10 In the boundary element environment, calculate and compare the tensile and shear moment values of each mesh, which specifically includes the following steps:
(1)依据但不限于Coulomb-Mohr强度准则描述岩体的拉剪破裂过程。Coulomb准则认为岩体的拉剪破裂过程发生在某一平面产生的破坏拉剪力矩值超过了抗拉剪力矩值即材料的内聚力和乘以常数的平面法应力。其数学表达式为|τ|=C+γtanμ其中τ是拉剪力矩值;C是内聚力或者黏结力,是无正压力时的抗拉剪强度;μ是内摩擦角;γ是一个固定常数。(1) Describe the tensile-shear failure process of rock masses based on but not limited to the Coulomb-Mohr strength criterion. The Coulomb criterion states that the tensile shear failure value of the rock mass during the tensile-shear failure process in a certain plane exceeds the tensile-shear moment value, that is, the cohesion of the material and the constant plane method stress. Its mathematical expression is | τ | = C + γtanμ, where τ is the tensile-shear moment value; C is the cohesive force or cohesive force, which is the tensile-shear strength without positive pressure; μ is the internal friction angle; γ is a fixed constant.
(2)通过计算煤体网格所受拉剪力矩,将网格所受拉剪力矩与抗拉剪力矩比对,寻找最易被破坏的网格集合。若存在拉剪力矩大于抗拉剪力矩的网格,则标记为“采动裂隙”,然后执行Step 11。若不存在拉剪力矩大于抗拉剪力矩的网格,则终止计算,执行Step 13。(2) By calculating the tensile and shearing moments of the coal body grid, comparing the tensile and shearing moments of the grid with the tensile and shearing moments, find the most vulnerable grid set. If there is a mesh with a tensile-shear moment greater than the tensile-shear moment, it is marked as "moving crack", and then Step 11 is performed. If there is no mesh with a tensile-shear moment greater than the tensile-shear moment, the calculation is terminated and Step 13 is executed.
Step 11替换“原生裂隙”网格中的材料为煤层内气体。Step 11 replaces the material in the "primary fissure" grid with gas in the coal seam.
(1)在煤体几何模型中剔除Step 10中被标记为“采动裂隙”的网格;(1) In the coal body geometric model, remove the grid marked as "mining fracture" in Step 10;
(2)将被剔除的几何单元材料(网格)变更为煤层内气体,孔隙率变更为1;(2) Change the removed geometric element material (grid) to gas in the coal seam, and change the porosity to 1;
(3)对因该单元被剔除而新生成的边界设定与相邻边界相同的边界条件;(3) Set the same boundary conditions as the adjacent boundary to the newly generated boundary due to the unit being eliminated;
Step 12判断“原生裂隙”的总面积是否超过断层面积的两倍,其具体包括如下步骤:Step 12: Determine whether the total area of the "primary fissure" exceeds twice the fault area, which specifically includes the following steps:
(1)如果没有超过,则生成删除原生裂隙网格的新模型,返回重复Step 10;(2)如果 超过,则终止计算,执行Step 13。(1) If not, generate a new model that deletes the original fissure mesh, and return to repeat Step 10. If it exceeds, terminate the calculation and execute Step 13.
Step 13构建最终包含地质构造原生裂隙、计算过程中生成裂隙两类煤层裂隙的模型,其具体包括如下步骤:Step 13 Build a model that finally includes two types of coal seam fractures, including primary fractures in the geological structure and fractures generated during the calculation.The specific steps include the following steps:
(1)生成包含断层的煤岩体模型;(1) Generate a coal rock model containing faults;
(2)重新剖分网格;(2) Re-mesh the mesh;
(3)导出STL通用几何;(3) Export STL general geometry;
Step 14在边界元环境下,耦合自定义方程计算渗流场与应力参数。 Step 14 In the boundary element environment, coupled with custom equations to calculate the seepage field and stress parameters.
(1)依据但不限于达西定律描述水分在煤体中的运移过程,即以
Figure PCTCN2018100462-appb-000011
作为渗流模拟的数学模型,其中,t为时间;φk为裂隙比空隙度;s为沿裂隙长度的坐标;k f为裂隙渗透系数;W为源汇项;p为煤体中水压(MPa)。若研究人员有需要添加的其它流动控制方程,可手动添加;
(1) Describe the transport process of water in coal body according to but not limited to Darcy's law.
Figure PCTCN2018100462-appb-000011
As a mathematical model for seepage simulation, where t is time; φk is the fracture-to-void ratio; s is the coordinate along the length of the fracture; k f is the fracture permeability coefficient; W is the source-sink term; p is the water pressure in the coal body (MPa ). If researchers have other flow control equations that they need to add, they can add them manually;
(2)岩体在孔隙流体的作用下,遵循修正的应力规律。依据但不限于基于流体渗流压力的应力方程σ ij=σ' ij+αpσ ij计算煤体应力分布,其中α称为有效应力系数。实践证明,α是孔隙压p和体积应力θ的函数,α=f(p,θ)。 (2) Under the action of pore fluid, the rock mass follows the modified stress law. The coal body stress distribution is calculated according to, but not limited to, the stress equation σ ij = σ ' ij + αpσ ij based on the fluid seepage pressure, where α is called the effective stress coefficient. Practice has proved that α is a function of pore pressure p and volume stress θ, α = f (p, θ).
在包含其它应力分布规律或特殊情况的煤层,可编程加入其它本构方程,对更加贴近实际的煤体应力分布进行模拟分析。In coal seams that contain other stress distribution laws or special conditions, other constitutive equations can be programmed to simulate and analyze the stress distribution of the coal body closer to the actual.
Step 15在无网格仿真环境中,编写算法,判断是否存在拉剪力矩大于煤岩体相关参数的点,其具体包括如下步骤:Step 15 In the meshless simulation environment, write an algorithm to determine whether there are points where the tensile-shear moment is greater than the relevant parameters of the coal and rock mass, which specifically include the following steps:
(1)如果存在,将拉剪力矩大于煤岩体相关参数的点,连接起来,形成闭合区域标记为“无效煤岩体”,继续执行Step 16;(1) If it exists, connect the points where the tensile-shear moment is greater than the relevant parameters of the coal rock mass to form a closed area marked as "invalid coal rock mass" and continue with Step 16;
(2)如果没有拉剪力矩大于煤岩体相关参数的点,则继续执行Step 17。(2) If there is no point where the tensile-shear moment is greater than the relevant parameters of the coal rock mass, then Step 17 is continued.
Step 16将Step 15中,标记的“无效煤岩体”边界重新添加水分入口条件。 Step 16 adds the water inlet conditions to the marked "invalid coal rock body" boundary in Step 15.
Step 17单独进行湍流模拟,并计算应力分布,并储存该节点的模拟结果。Step 17: Perform turbulence simulation separately, calculate the stress distribution, and store the simulation results of this node.
Step 18判断断裂脱离是否达到煤体表面,其具体包括如下步骤:Step 18: Determine whether the fracture separation has reached the surface of the coal body, which specifically includes the following steps:
(1)如果断裂脱离没有到达煤体表面,则执行Step 19;(1) If the fracture does not reach the surface of the coal body, then Step 19 is performed;
(2)如果断裂脱离到达煤体表面,则执行Step 20。(2) If the fracture breaks off and reaches the surface of the coal body, perform Step 20.
Step 19累计存储时间是否达到预设模拟时间,其具体包括如下步骤:Step 19 Whether the accumulated storage time reaches the preset simulation time, which specifically includes the following steps:
(1)如果存储时间没有达到预设模拟时间,则返回执行Step 14;(1) If the storage time does not reach the preset simulation time, return to Step 14;
(2)如果存储时间达到预设模拟时间,则执行Step 20。(2) If the storage time reaches the preset simulation time, execute Step 20.
Step 20停止无网格法环境中的运算,仅在边界元环境中计算湍流与应力,后执行Step 7。 Step 20 stops the calculation in the meshless method environment, calculates turbulence and stress only in the boundary element environment, and then executes Step 7.
Step 21将整合后的结果输出并存储为独立文件,量化统计结果,如图6的水压场模拟结果图、图7的瓦斯压力模拟结果图、图8的煤体压力结果模拟图、图9的水分流速模拟结果图与图10的瓦斯流速模拟结果图。Step 21 outputs the integrated result and stores it as an independent file to quantify the statistical results, such as the hydraulic pressure field simulation result chart in Fig. 6, the gas pressure simulation result chart in Fig. 7, the coal pressure result simulation chart in Fig. 8, and Fig. 9 The simulation results of water flow velocity and the simulation results of gas flow velocity in FIG.
依本数值模拟方法结果运算得到的结果可以看出,随着更大范围的煤体被水分润湿,水分运移的主要动力来源已变为毛细作用力而非注水压力,因此水分与瓦斯的平均压力均逐渐下降;但水分压力的下降梯度较为恒定,而瓦斯压力收到不规则孔、裂隙结构的影响,下降梯度差异明显。煤体随注水作业进行而逐渐被软化,因此综合作用压力总体上呈现出下降趋势。同时,水分与瓦斯的流速变化与各自的压力变化基本一致,但瓦斯流速高于水分流速。According to the results calculated by the results of this numerical simulation method, it can be seen that as a larger range of coal bodies are wetted by water, the main source of power for water transport has become capillary force rather than water injection pressure. The average pressure gradually decreases; but the gradient of the water pressure is relatively constant, while the gas pressure is affected by the structure of irregular pores and fissures, and the gradient of the gradient is obvious. The coal body is gradually softened with the water injection operation, so the overall action pressure shows a downward trend. At the same time, the changes in the flow rate of water and gas are basically consistent with the changes in their respective pressures, but the flow rate of gas is higher than the flow rate of water.
上述现象与工程现场应用结果相一致,说明本发明的模拟方法所得结果是可靠性的。The above phenomena are consistent with the application results in the engineering field, which shows that the results obtained by the simulation method of the present invention are reliable.
当然,以上说明仅仅为本发明的较佳实施例,本发明并不限于列举上述实施例,应当说明的是,任何熟悉本领域的技术人员在本说明书的教导下,所做出的所有等同替代、明显变形形式,均落在本说明书的实质范围之内,理应受到本发明的保护。Of course, the above description is only a preferred embodiment of the present invention, and the present invention is not limited to the above embodiments. It should be noted that all equivalent substitutions made by those skilled in the art under the teaching of this specification The obvious deformation forms all fall within the substantive scope of this specification and should be protected by the present invention.

Claims (6)

  1. 一种煤岩体分区注水渗流-损伤-应力耦合数值模拟方法,其包括以下步骤:A coupled numerical simulation method for seepage-damage-stress coupling of water injection in coal and rock masses, comprising the following steps:
    步骤一,依据地质勘测结果,建立包含断层的煤岩体模型;Step 1: Establish a coal and rock mass model including faults based on geological survey results;
    步骤二,对地质勘测得到之煤岩体中的相应煤体进行扫描,并结合FDK三维重建算法构建具有真实孔隙结构特征的三维数字岩心模型;Step 2: Scan the corresponding coal body in the coal rock body obtained from the geological survey, and combine the FDK 3D reconstruction algorithm to construct a 3D digital core model with real pore structure characteristics;
    步骤三、在边界元环境下,执行“采动裂隙”生成算法,判断是否为采动影响区,即判断距离采掘面是否小于80m;如果是采动影响区,则执行步骤四;如果不是采动影响区,则执行步骤五;Step 3: In the boundary element environment, execute the "mining fissure" generation algorithm to determine whether it is a mining affected area, that is, determine whether the distance from the mining face is less than 80m. If it is a mining affected area, go to step 4; if it is not mining Move the affected area, go to step five;
    步骤四、判断流域的截面积是否大于30μm 2,若流域截面积大于30μm 2,则进行N-S计算处理,储存并输出结果;若小于30μm 2,则进行Darcy计算处理,储存并输出结果; Step 4. Determine whether the cross-sectional area of the watershed is greater than 30 μm 2. If the cross-sectional area of the watershed is greater than 30 μm 2 , perform NS calculation processing, store and output the result; if less than 30 μm 2 , perform Darcy calculation processing, store and output the result;
    步骤五、判断流域的截面积是否大于30μm 2,若流域截面积大于30μm 2,则执行步骤六,若小于30μm 2,则执行步骤十二; Step 5. Determine whether the cross-sectional area of the watershed is greater than 30 μm 2. If the cross-sectional area of the watershed is greater than 30 μm 2 , go to step 6; if less than 30 μm 2 , go to step 12;
    步骤六、以C#语言为基础,编程以埋深、地质条件、上覆岩层岩性及煤层倾角生成初始应力集中点。Step 6. Based on the C # language, the programming generates initial stress concentration points based on burial depth, geological conditions, lithology of overlying strata, and coal seam inclination.
    步骤七、对经过步骤六处理的三维数字岩心模型进行剖分网格;Step 7. Mesh the 3D digital core model processed in Step 6;
    步骤八、通过计算煤体网格所受拉剪力矩,将网格所受拉剪力矩与抗拉剪力矩比对,寻找最易被破坏的网格集合;若存在拉剪力矩大于抗拉剪力矩的网格,则标记为“原生裂隙”,然后执行步骤九;若不存在拉剪力矩大于抗拉剪力矩的网格,则终止计算,执行步骤十一;Step 8. Calculate the tensile and shearing moments of the coal body grid, compare the tensile and shearing moments of the grid with the tensile and shearing moments, and find the most vulnerable grid set; if there is a tensile and shearing moment greater than the tensile and shearing, The grid of moments is marked as "Native Crack", and then step 9 is performed; if there is no grid with a tensile-shear moment greater than the tensile-shear moment, the calculation is terminated and step 11 is performed;
    步骤九、替换“原生裂隙”网格中的材料为煤层内气体;Step 9. Replace the material in the "primary fissure" grid with the gas in the coal seam;
    步骤十、判断经过步骤九处理后之“原生裂隙”的总面积是否超过断层面积的两倍,如果没有超过,则生成删除“原生裂隙”网格后的新模型,返回重复步骤八;如果超过,则终止计算,执行步骤十一;Step 10. Determine whether the total area of the "primary fissures" after the processing in step 9 exceeds twice the fault area. If not, generate a new model after deleting the "primary fissure" grid, and return to repeat step 8. If it exceeds , The calculation is terminated and step 11 is executed;
    步骤十一、构建最终包含地质构造原生裂隙、计算过程中生成裂隙两类煤层裂隙的模型,并重新剖分网格,导出STL通用几何后进行N-S计算处理,储存并输出结果;Step 11: Construct a model of two types of coal seam fractures that ultimately include primary fractures in the geological structure and fractures generated during calculation, and re-mesh the mesh to derive the STL general geometry for N-S calculation processing, and store and output the results;
    步骤十二、在边界元环境下,耦合自定义方程计算渗流场与应力参数后,执行步骤十三;Step 12. In the boundary element environment, after coupling the custom equation to calculate the seepage field and stress parameters, go to Step 13.
    步骤十三、在无网格仿真环境中,判断是否存在拉剪力矩大于煤岩体相关参数的点,如果有,则将此类点依次连接,闭合区域标记为“无效煤岩体”,执行步骤十四;如果没有,则执行步骤十五;Step 13. In the meshless simulation environment, determine whether there are points where the tensile-shear moment is greater than the relevant parameters of the coal and rock mass. If so, connect such points in turn, mark the closed area as "invalid coal and rock mass", and execute Step 14; if not, go to Step 15;
    步骤十四、将步骤十三中,标记的“无效煤岩体”边界重新添加水分入口条件;Step 14. Add the water inlet conditions to the marked "invalid coal and rock mass" boundary in step 13.
    步骤十五、单独进行湍流模拟,并计算应力分布,并储存该节点的模拟结果;Step 15: Perform turbulence simulation separately, calculate the stress distribution, and store the simulation results of the node;
    步骤十六、判断步骤十五中模拟结果的断裂脱离是否达到煤体表面,如果没有,则执行步骤十七;如果到达煤体表面,则执行步骤十八;Step 16. Determine whether the fracture separation from the simulation result in step 15 reaches the surface of the coal body. If not, go to step 17; if it reaches the surface of the coal body, go to step 18.
    步骤十七、累计存储时间是否达到预设模拟时间,如果没有达到,则返回执行步骤十二;如果达到计算时间,则执行步骤十八;Step 17. Check whether the accumulated storage time reaches the preset simulation time. If not, go back to Step 12. If the calculation time is reached, go to Step 18.
    步骤十八、停止无网格仿真环境中的运算,仅在边界元环境中计算湍流与应力后则进行Darcy计算处理,储存并输出结果;Step 18. Stop the calculation in the meshless simulation environment. After calculating turbulence and stress only in the boundary element environment, perform Darcy calculation processing, store and output the results;
    步骤十九、将步骤四、步骤十一与步骤十八的输出结果整合后输出并存储为独立文件,得到量化统计结果。Step 19: The output results of Step 4, Step 11 and Step 18 are integrated and output as an independent file to obtain a quantitative statistical result.
  2. 根据权利要求1所述的煤岩体分区注水渗流-损伤-应力耦合数值模拟方法,其特征在于,N-S计算处理包括:The method for numerically coupling seepage-damage-stress coupled numerical simulation of water injection in a coal rock mass according to claim 1, wherein the N-S calculation process includes:
    步骤A、N-S初始化;Step A, N-S initialization;
    步骤B、计算水压和瓦斯压;Step B: Calculate water pressure and gas pressure;
    步骤C、判断水压是否大于瓦斯压,若大于,则将下一网格材料变更为水,执行步骤D,若小于,则直接执行步骤D;Step C, determine whether the water pressure is greater than the gas pressure, and if it is greater, change the next grid material to water, and execute step D; if it is less, directly perform step D;
    步骤D、判断是否达到计算时间设定,若不能达到,则返回重新执行步骤A,若达到计算时间,则直接储存并输出结果。Step D: Determine whether the calculation time setting is reached. If the calculation time is not reached, return to step A. If the calculation time is reached, directly store and output the result.
  3. 根据权利要求1所述的煤岩体分区注水渗流-损伤-应力耦合数值模拟方法,其特征在于,Darcy计算处理包括:The numerical simulation method for seepage-damage-stress coupled numerical simulation of water injection in a coal rock mass according to claim 1, wherein the Darcy calculation process includes:
    步骤E、Darcy初始化;Step E, Darcy initialization;
    步骤F、计算水压和瓦斯压;Step F: Calculate water pressure and gas pressure;
    步骤J、判断水压是否大于瓦斯压,若大于,则将下一网格材料变更为水,执行步骤H,若小于,则直接执行步骤H;Step J: Determine whether the water pressure is greater than the gas pressure, and if it is greater, change the next grid material to water, and perform step H; if it is less, directly perform step H;
    步骤H、判断是否达到计算时间设定,若不能达到,则返回重新执行步骤E,若达到计算时间,则直接储存并输出结果。Step H: Determine whether the calculation time setting is reached. If the calculation time is not reached, return to step E. If the calculation time is reached, directly store and output the result.
  4. 根据权利要求1所述的煤岩体分区注水渗流-损伤-应力耦合数值模拟方法,其特征在于,上述步骤一具体的还包括:建立包含断层的煤岩体模型的过程中预留出“断层起点”、“断层终点”、“断层转折点”、“方位差”、“落差”与“倾角”六个参数供用户输入。The method for numerical simulation of seepage-damage-stress coupling of zoned water injection in a coal rock mass according to claim 1, wherein a specific step of the above step further comprises: "faults are reserved in the process of establishing a coal rock mass model including faults" The starting point, the end point of the fault, the turning point of the fault, the azimuth difference, the drop, and the inclination are provided for user input.
  5. 根据权利要求1所述的煤岩体分区注水渗流-损伤-应力耦合数值模拟方法,其特征在于,上述步骤二聚体的还包括:利用图像处理算法对三维数字岩心模型进行滤波操作,平滑模型边缘,进而基于阈值实现数据分割得到煤体的微观孔隙结构,剔除连通性较差的孤岛孔隙,导出得到最终的STL格式煤体孔隙数字模型。The method for numerically coupling water injection seepage-damage-stress coupling of a coal rock mass according to claim 1, wherein the step dimer further comprises: using an image processing algorithm to perform a filtering operation on a three-dimensional digital core model to smooth the model. The edge, and then based on the threshold to achieve data segmentation to obtain the micro-pore structure of the coal body, remove the island pores with poor connectivity, and derive the final digital model of coal body pores in STL format.
  6. 根据权利要求1所述的煤岩体分区注水渗流-损伤-应力耦合数值模拟方法,其特征在于,上述量化统计结果包括水压场模拟结果、瓦斯压力模拟结果、煤体压力模拟结果、水分 流速模拟结果与瓦斯流速模拟结果。The method for numerical simulation of seepage-damage-stress coupling of water injection in a coal rock mass according to claim 1, wherein the quantitative statistical results include simulation results of hydraulic pressure field, simulation results of gas pressure, simulation results of coal pressure, and water flow velocity Simulation results and gas velocity simulation results.
PCT/CN2018/100462 2018-07-02 2018-08-14 Method for area-based coal rock water injection seepage-damage-stress coupling value simulation WO2020006818A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2020110457A RU2743121C1 (en) 2018-07-02 2018-08-14 Method of numerical modeling of filtration/damage/stress communication during water injection into coal-bearing mass during regionalization

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201810705645.3A CN109063257B (en) 2018-07-02 2018-07-02 A kind of coal and rock subregion water filling seepage flow-damage-stress coupling method for numerical simulation
CN201810705645.3 2018-07-02

Publications (1)

Publication Number Publication Date
WO2020006818A1 true WO2020006818A1 (en) 2020-01-09

Family

ID=64818148

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2018/100462 WO2020006818A1 (en) 2018-07-02 2018-08-14 Method for area-based coal rock water injection seepage-damage-stress coupling value simulation

Country Status (3)

Country Link
CN (1) CN109063257B (en)
RU (1) RU2743121C1 (en)
WO (1) WO2020006818A1 (en)

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111307687A (en) * 2020-03-05 2020-06-19 中煤科工集团重庆研究院有限公司 Method for evaluating permeability and strength of bonding of high polymer material and coal rock mass
CN111666699A (en) * 2020-04-30 2020-09-15 山东大学 Rock mass engineering cross-scale simulation calculation method based on REV full-area coverage
CN111860952A (en) * 2020-06-16 2020-10-30 重庆大学 Method for quickly optimizing key mining parameters of outburst coal seam
CN111898187A (en) * 2020-07-23 2020-11-06 武汉大学 Method for determining longitudinal model range value of tunnel excavation seepage simulation analysis
CN112131709A (en) * 2020-08-25 2020-12-25 山东大学 Jointed rock mechanics simulation method and system based on near-field dynamics constitutive model
CN112507575A (en) * 2020-11-03 2021-03-16 辽宁工程技术大学 Fault grouting reinforcement effect evaluation method based on FLAC3D numerical simulation
CN112560226A (en) * 2020-11-25 2021-03-26 国家能源投资集团有限责任公司 Coal mine underground reservoir environment risk assessment method
CN112800597A (en) * 2021-01-13 2021-05-14 安徽马钢张庄矿业有限责任公司 Intelligent fine efficient ecological mining analysis method for high and middle segments of mine resources
CN112926270A (en) * 2021-03-18 2021-06-08 西安科技大学 Gas multi-factor coupling relation analysis and early warning model construction method
CN112924331A (en) * 2021-01-12 2021-06-08 江苏师范大学 Method for establishing water-rock coupling model of compressive strength of coal rock after water solution soaking
CN113533042A (en) * 2021-07-07 2021-10-22 北京科技大学 Comprehensive index calculation method for representing rock stress and fracture and application
CN113654477A (en) * 2021-08-16 2021-11-16 中国矿业大学 Coal body deformation testing device, testing system and testing method
CN113742940A (en) * 2021-09-16 2021-12-03 重庆大学 Mining pressure relief boundary space-time curve cluster determination method
CN114002129A (en) * 2021-12-02 2022-02-01 河北省交通规划设计研究院有限公司 High-water-pressure-crack rock mass seepage test platform
CN114217048A (en) * 2021-12-10 2022-03-22 国家能源投资集团有限责任公司 Mining three-dimensional simulation experiment model experiment method
CN114495679A (en) * 2022-01-25 2022-05-13 中国矿业大学 Method for manufacturing real coal two-dimensional microfluidic model
CN114544437A (en) * 2022-02-28 2022-05-27 中国矿业大学 Coal rock flow field in-situ fluorescent bacteria microscopic tracing method
CN114638137A (en) * 2022-03-31 2022-06-17 王永亮 Hot-dry rock heat production prediction method based on heat-flow-solid-damage coupling
CN114660268A (en) * 2022-03-22 2022-06-24 中铁水利水电规划设计集团有限公司 Water retaining layer seepage monitoring system for farmland of farmland lifting area of reservoir inundation area
CN114809992A (en) * 2022-04-20 2022-07-29 太原理工大学 Low-permeability reservoir coal-based gas full-life-cycle efficient extraction method
CN115014982A (en) * 2022-08-09 2022-09-06 中国矿业大学(北京) Method for evaluating coal rock damage inferiority based on gas migration pressure fluctuation characteristics
CN115201898A (en) * 2022-05-20 2022-10-18 中国地震局地球物理研究所 Numerical simulation method and system for three-dimensional injection-production induced seismic fault fracture sliding
CN115524275A (en) * 2022-10-18 2022-12-27 中国水利水电科学研究院 Fractured rock mass permeability tensor determination method considering rock mass seepage and fracture seepage
CN117408191A (en) * 2023-12-15 2024-01-16 山东大学 Fracture rock mass permeation-flow grouting simulation method and system

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109117589B (en) * 2018-09-11 2022-10-04 中煤科工集团重庆研究院有限公司 Quantitative description method for coal seam roof fracture field
CN109883920A (en) * 2019-03-08 2019-06-14 西南石油大学 Characterize the method and device of rock thermal damage
CN110043259B (en) * 2019-04-18 2020-05-26 中国矿业大学 Regional stoping process for composite top working face with large mining depth and high mining height
CN110118991B (en) * 2019-05-16 2020-06-23 中国矿业大学 Mining induced stress assessment method based on microseismic damage reconstruction
CN110489832B (en) * 2019-07-31 2023-05-23 中国航发沈阳发动机研究所 Simulation test method for pneumatic performance of turbulence control screen unit body
CN110390176B (en) * 2019-07-31 2020-05-19 西南交通大学 Ballastless track freezing and damage behavior calculation method
CN110872942B (en) * 2019-10-24 2021-09-28 中国石油化工股份有限公司 Research method for contribution of different acting forces to recovery ratio in oil reservoir injection-production coupling mode
CN110866337B (en) * 2019-11-12 2021-06-01 中南大学 Differential stress-based mining fault activation tendency judgment method
CN111695285B (en) * 2020-06-17 2023-12-22 大连海事大学 Anisotropic rock stress-damage-seepage coupling numerical simulation method
CN112329312B (en) * 2020-11-10 2022-07-26 河海大学 Rapid generation method of three-dimensional seepage stress coupling cohesion unit
CN112903966A (en) * 2021-01-20 2021-06-04 中国矿业大学(北京) Coal mining damage range determining method based on energy transfer conservation
CN113722960A (en) * 2021-08-31 2021-11-30 中国地质大学(武汉) Three-dimensional numerical simulation method for humidity migration and drying shrinkage cracking
CN116756985B (en) * 2022-11-29 2024-01-30 华东师范大学 COMSOL Multiphysics-based field multi-medium environment organic pollutant migration simulation method
CN117627669A (en) * 2024-01-26 2024-03-01 中交一航局第三工程有限公司 Rectangular jacking pipe construction method based on blind area advanced pilot tunnel expansion and excavation

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU796420A1 (en) * 1977-04-04 1981-01-15 Всесоюзный Научно-Исследовательскийинститут Горной Геомеханики Имаркшейдерского Дела Method of simulating strained state of rock body with models made of similar materials
CN102252957A (en) * 2011-04-19 2011-11-23 河南理工大学 Experimental measurement device and method of solid-liquid conversion parameters of decompressed coal
CN103034765A (en) * 2012-12-14 2013-04-10 天津大学 Goaf grouting dynamic overall process simulation method based on numerical simulation
CN105401939A (en) * 2015-11-30 2016-03-16 中国石油大学(北京) Coal seam well wall stability analysis method under multi-factor coupled effect
CN105787220A (en) * 2016-04-22 2016-07-20 山东科技大学 Coal bed high-pressure water injection fracturing-flow seeping value simulation method
CN106960070A (en) * 2016-12-28 2017-07-18 山东科技大学 A kind of seepage simulation method that coal body is reconstructed based on finite element discretization member CT

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012079009A2 (en) * 2010-12-10 2012-06-14 Baker Hughes Incorporated A method to improve reservoir simulation and recovery from fractured reservoirs
CN104653226B (en) * 2014-12-26 2017-07-14 中国矿业大学 A kind of division methods of the coal mine rock burst danger zone based on stress gradient
CN108038282A (en) * 2017-11-30 2018-05-15 安徽理工大学 Load transmission method for numerical simulation is exploited under a kind of unconsolidated confined aquifer

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU796420A1 (en) * 1977-04-04 1981-01-15 Всесоюзный Научно-Исследовательскийинститут Горной Геомеханики Имаркшейдерского Дела Method of simulating strained state of rock body with models made of similar materials
CN102252957A (en) * 2011-04-19 2011-11-23 河南理工大学 Experimental measurement device and method of solid-liquid conversion parameters of decompressed coal
CN103034765A (en) * 2012-12-14 2013-04-10 天津大学 Goaf grouting dynamic overall process simulation method based on numerical simulation
CN105401939A (en) * 2015-11-30 2016-03-16 中国石油大学(北京) Coal seam well wall stability analysis method under multi-factor coupled effect
CN105787220A (en) * 2016-04-22 2016-07-20 山东科技大学 Coal bed high-pressure water injection fracturing-flow seeping value simulation method
CN106960070A (en) * 2016-12-28 2017-07-18 山东科技大学 A kind of seepage simulation method that coal body is reconstructed based on finite element discretization member CT

Cited By (43)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111307687A (en) * 2020-03-05 2020-06-19 中煤科工集团重庆研究院有限公司 Method for evaluating permeability and strength of bonding of high polymer material and coal rock mass
CN111307687B (en) * 2020-03-05 2022-07-19 中煤科工集团重庆研究院有限公司 Method for evaluating permeability and strength of bonding of high polymer material and coal rock mass
CN111666699A (en) * 2020-04-30 2020-09-15 山东大学 Rock mass engineering cross-scale simulation calculation method based on REV full-area coverage
CN111666699B (en) * 2020-04-30 2023-06-02 山东大学 Rock mass engineering trans-scale simulation calculation method based on REV full-area coverage
CN111860952A (en) * 2020-06-16 2020-10-30 重庆大学 Method for quickly optimizing key mining parameters of outburst coal seam
CN111860952B (en) * 2020-06-16 2024-02-27 重庆大学 Method for rapidly optimizing key mining parameters of outburst coal seam
CN111898187A (en) * 2020-07-23 2020-11-06 武汉大学 Method for determining longitudinal model range value of tunnel excavation seepage simulation analysis
CN111898187B (en) * 2020-07-23 2022-10-11 武汉大学 Method for determining longitudinal model range value of tunnel excavation seepage simulation analysis
CN112131709A (en) * 2020-08-25 2020-12-25 山东大学 Jointed rock mechanics simulation method and system based on near-field dynamics constitutive model
CN112131709B (en) * 2020-08-25 2024-04-19 山东大学 Joint rock mass mechanical simulation method and system based on near field dynamics constitutive model
CN112507575A (en) * 2020-11-03 2021-03-16 辽宁工程技术大学 Fault grouting reinforcement effect evaluation method based on FLAC3D numerical simulation
CN112560226B (en) * 2020-11-25 2023-12-26 国家能源投资集团有限责任公司 Coal mine underground reservoir environment risk assessment method
CN112560226A (en) * 2020-11-25 2021-03-26 国家能源投资集团有限责任公司 Coal mine underground reservoir environment risk assessment method
CN112924331A (en) * 2021-01-12 2021-06-08 江苏师范大学 Method for establishing water-rock coupling model of compressive strength of coal rock after water solution soaking
CN112800597A (en) * 2021-01-13 2021-05-14 安徽马钢张庄矿业有限责任公司 Intelligent fine efficient ecological mining analysis method for high and middle segments of mine resources
CN112800597B (en) * 2021-01-13 2022-09-13 安徽马钢张庄矿业有限责任公司 Intelligent fine efficient ecological mining analysis method for high and middle segments of mine resources
CN112926270A (en) * 2021-03-18 2021-06-08 西安科技大学 Gas multi-factor coupling relation analysis and early warning model construction method
CN112926270B (en) * 2021-03-18 2023-05-05 西安科技大学 Gas multi-factor coupling relation analysis and early warning model construction method
CN113533042A (en) * 2021-07-07 2021-10-22 北京科技大学 Comprehensive index calculation method for representing rock stress and fracture and application
CN113654477A (en) * 2021-08-16 2021-11-16 中国矿业大学 Coal body deformation testing device, testing system and testing method
CN113742940B (en) * 2021-09-16 2023-09-08 重庆大学 Mining pressure relief boundary space-time curve cluster determination method
CN113742940A (en) * 2021-09-16 2021-12-03 重庆大学 Mining pressure relief boundary space-time curve cluster determination method
CN114002129B (en) * 2021-12-02 2022-08-09 河北省交通规划设计研究院有限公司 High-water-pressure-crack rock mass seepage test platform
CN114002129A (en) * 2021-12-02 2022-02-01 河北省交通规划设计研究院有限公司 High-water-pressure-crack rock mass seepage test platform
CN114217048A (en) * 2021-12-10 2022-03-22 国家能源投资集团有限责任公司 Mining three-dimensional simulation experiment model experiment method
CN114217048B (en) * 2021-12-10 2024-05-24 国家能源投资集团有限责任公司 Mining three-dimensional simulation experiment model experiment method
CN114495679B (en) * 2022-01-25 2022-10-28 中国矿业大学 Method for manufacturing real coal two-dimensional microfluidic model
CN114495679A (en) * 2022-01-25 2022-05-13 中国矿业大学 Method for manufacturing real coal two-dimensional microfluidic model
CN114544437B (en) * 2022-02-28 2023-12-19 中国矿业大学 Microscopic tracing method for in-situ fluorescent bacteria of coal-rock flow field
CN114544437A (en) * 2022-02-28 2022-05-27 中国矿业大学 Coal rock flow field in-situ fluorescent bacteria microscopic tracing method
CN114660268A (en) * 2022-03-22 2022-06-24 中铁水利水电规划设计集团有限公司 Water retaining layer seepage monitoring system for farmland of farmland lifting area of reservoir inundation area
CN114638137A (en) * 2022-03-31 2022-06-17 王永亮 Hot-dry rock heat production prediction method based on heat-flow-solid-damage coupling
CN114638137B (en) * 2022-03-31 2023-03-28 王永亮 Hot-dry rock heat production prediction method based on heat-flow-solid-damage coupling
CN114809992B (en) * 2022-04-20 2023-08-08 太原理工大学 Low permeability reservoir coal gas full life cycle efficient extraction method
CN114809992A (en) * 2022-04-20 2022-07-29 太原理工大学 Low-permeability reservoir coal-based gas full-life-cycle efficient extraction method
CN115201898B (en) * 2022-05-20 2023-06-20 中国地震局地球物理研究所 Numerical simulation method and system for three-dimensional injection and production induced seismic fault fracture sliding
CN115201898A (en) * 2022-05-20 2022-10-18 中国地震局地球物理研究所 Numerical simulation method and system for three-dimensional injection-production induced seismic fault fracture sliding
CN115014982B (en) * 2022-08-09 2022-10-11 中国矿业大学(北京) Method for evaluating coal rock damage inferiority based on gas migration pressure fluctuation characteristics
CN115014982A (en) * 2022-08-09 2022-09-06 中国矿业大学(北京) Method for evaluating coal rock damage inferiority based on gas migration pressure fluctuation characteristics
CN115524275B (en) * 2022-10-18 2023-08-22 中国水利水电科学研究院 Fractured rock mass permeability tensor determination method considering rock mass seepage and fracture seepage
CN115524275A (en) * 2022-10-18 2022-12-27 中国水利水电科学研究院 Fractured rock mass permeability tensor determination method considering rock mass seepage and fracture seepage
CN117408191A (en) * 2023-12-15 2024-01-16 山东大学 Fracture rock mass permeation-flow grouting simulation method and system
CN117408191B (en) * 2023-12-15 2024-04-02 山东大学 Fracture rock mass permeation-flow grouting simulation method and system

Also Published As

Publication number Publication date
CN109063257B (en) 2019-04-26
CN109063257A (en) 2018-12-21
RU2743121C1 (en) 2021-02-15

Similar Documents

Publication Publication Date Title
WO2020006818A1 (en) Method for area-based coal rock water injection seepage-damage-stress coupling value simulation
CN106960070B (en) Seepage simulation method for reconstructing coal body based on finite element-discrete element CT
CN107060746B (en) A kind of method of complex fracture oil deposit flow simulation
WO2016192077A1 (en) Method for establishing and solving numerical well-testing model of horizontal well for tight gas hydraulic fracturing
CN105019877A (en) Calculating method for size of fracturing and broken region of horizontal well for shale
Segura et al. Coupled HM analysis using zero‐thickness interface elements with double nodes—Part II: Verification and application
CA2739590A1 (en) Sand and fluid production and injection modeling methods
Jin et al. Two-phase PFEM with stable nodal integration for large deformation hydromechanical coupled geotechnical problems
CN110579400B (en) Measuring and calculating method for micro-scale strength and residual strength of brittle rock
CN110984973A (en) Determination method for single-well controlled reserve of fracture-cave carbonate gas reservoir
CN109162701B (en) Coal seam open hole well fracture pressure prediction method
CN112012712A (en) Numerical simulation method and device for water injection growth seam of embedded discrete seam
CN107423466A (en) A kind of proppant embedment and the method for numerical simulation of fracture condudtiviy quantitative forecast
Xiong et al. Influence of natural fractures on hydraulic fracture propagation behaviour
CN105114065A (en) Method for simulating hydraulic fracturing fracture extension process
CN110738001A (en) Calculation method for unconventional reservoir fracturing yield-increasing transformation area
Cheng et al. 3-D dynamic evolution analysis of coal-rock damaged field and gas seepage field during the gas extraction process
CN115510778A (en) Continental facies shale reservoir infinite stage fracturing process optimization method and system
CN113591420B (en) Simulation method and processor for true triaxial hydraulic fracturing experiment
Damirchi et al. Coupled hydro-mechanical modelling of saturated fractured porous media with unified embedded finite element discretisations
Bertrand et al. Anisotropic modelling of Opalinus Clay behaviour: From triaxial tests to gallery excavation application
CN109558614B (en) Simulation method and system for gas flow in shale gas reservoir multi-scale fracture
Um et al. Geochemical and geomechanical effects on wellbore cement fractures
CN107066744B (en) Earthquake-resistant analysis method for subway station in water-rich stratum
Yang et al. Coupled hydro-mechanical analysis of gas production in fractured shale reservoir by random fracture network modeling

Legal Events

Date Code Title Description
NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 18925208

Country of ref document: EP

Kind code of ref document: A1