CN112362556B - Method for obtaining continuous function of permeability coefficient of coal mine mining stable region - Google Patents

Method for obtaining continuous function of permeability coefficient of coal mine mining stable region Download PDF

Info

Publication number
CN112362556B
CN112362556B CN202011269283.1A CN202011269283A CN112362556B CN 112362556 B CN112362556 B CN 112362556B CN 202011269283 A CN202011269283 A CN 202011269283A CN 112362556 B CN112362556 B CN 112362556B
Authority
CN
China
Prior art keywords
pressure
seepage
darcy
test
porosity
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011269283.1A
Other languages
Chinese (zh)
Other versions
CN112362556A (en
Inventor
任奕玮
姜德义
陈结
刘文浩
贾栋
杜俊生
陈波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University
Original Assignee
Chongqing University
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 Chongqing University filed Critical Chongqing University
Priority to CN202011269283.1A priority Critical patent/CN112362556B/en
Publication of CN112362556A publication Critical patent/CN112362556A/en
Application granted granted Critical
Publication of CN112362556B publication Critical patent/CN112362556B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/082Investigating permeability by forcing a fluid through a sample
    • G01N15/0826Investigating permeability by forcing a fluid through a sample and measuring fluid flow rate, i.e. permeation rate or pressure change
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088Investigating volume, surface area, size or distribution of pores; Porosimetry

Landscapes

  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Dispersion Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Fluid Mechanics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The invention discloses a method for obtaining a continuous function of permeability coefficient of a mining stable region of a coal mine, which comprises the following steps: 1) Collecting a target working face surrounding rock sample, measuring the uniaxial compressive strength and the uniaxial compressive strength of the complete rock sample, and setting the crushing pressure-bearing strength by the uniaxial compressive strength; 2) Preparing the same lithotripsy stone sample with the particle size of not less than three, and drying; 3) Carrying out pressure-bearing seepage test on crushed stone samples with different particle sizes; 4) After the test is finished, selecting a plurality of groups of typical data, and making a change curve of pore pressure gradient and seepage velocity; 5) Processing the change curve of the pore pressure gradient and the seepage velocity obtained by fitting according to a Forchheimer type non-Darcy equation to obtain a non-Darcy osmotic coefficient and a non-Darcy flow factor; and then binary fitting is carried out to obtain a functional relation between the permeability coefficient and the porosity and granularity of the crushed stone sample in the pressure-bearing seepage test, wherein the permeability coefficient obtained by the method is more in line with the actual situation.

Description

Method for obtaining continuous function of permeability coefficient of coal mine mining stable region
Technical Field
The invention relates to the technical field of coal mine gas extraction, in particular to the technical field of coal mine tunneling, and in particular relates to a method for obtaining a continuous function of permeability coefficient of a coal mine mining stability zone, which can be directly applied to ground position selection of a coal mine mining stability zone gas extraction well.
Background
In the research of the coal bed gas migration rule in the extraction process of the mining stable region, the permeability coefficient distribution condition is a key parameter of the research. At present, most of permeability coefficient selection is a uniform distribution function or a segmented uniform distribution function, which is different from the actual situation to a certain extent, and may cause inaccurate coal bed gas migration rules obtained through research.
Disclosure of Invention
In view of the above, the invention aims to provide a method for obtaining a continuous function of permeability coefficient of a mining stable region of a coal mine so as to solve the technical problem of inaccurate permeability coefficient selection.
The invention discloses a method for obtaining a continuous function of permeability coefficient of a mining stable region of a coal mine, which comprises the following steps:
1) Collecting a target working face surrounding rock sample, measuring the uniaxial compressive strength and the uniaxial compressive strength of the complete rock sample, and setting the crushing pressure-bearing strength by the uniaxial compressive strength;
2) Preparing identical lithotripsy stone samples with the particle sizes not less than three and consistent with the test results of the step 1), and drying;
3) Carrying out pressure-bearing seepage test on crushed stone samples with different particle sizes: the porosity of the crushed stone sample is reduced by compression, so that the porosity change of the rock falling from different places in a mining stable region due to the action of mine pressure is simulated, and gas for testing seepage conditions is introduced into the crushed stone sample in the process of compressing the crushed stone sample; obtaining experimental data of a plurality of groups of axial loads, pressure head displacement, air inlet pressure, air outlet pressure, air inlet flow and air outlet flow by changing the pressure head displacement and the pressure of the introduced gas;
4) After the test is finished, selecting a plurality of groups of typical data, making a change curve of pore pressure gradient and seepage velocity, comparing the change curve of pore pressure gradient and seepage velocity obtained by fitting with a relation curve of pore pressure gradient and seepage velocity expressed by a Forchheimer type non-Darcy equation, and judging whether the relation curve accords with the Forchheimer type non-Darcy equation; the Forchheimer's non-Darcy equation is as follows:
wherein J is pore pressure gradient, and the unit is Pa/m; ρ is the fluid density in kg/m 3 The method comprises the steps of carrying out a first treatment on the surface of the g is gravity acceleration, and the unit is N/kg; beta is the fidaxflow factor in m -1 The method comprises the steps of carrying out a first treatment on the surface of the v is the seepage velocity, the unit is m/s; k is the non-Darcy osmotic coefficient, the unit is m/s,wherein k is permeability in m -2 The method comprises the steps of carrying out a first treatment on the surface of the u is the dynamic coefficient of viscosity of the fluid, and the unit is Pa.s;
5) After determining that the change curve of the pore pressure gradient and the seepage velocity accords with the Forchheimer type non-Darcy equation, processing the change curve of the pore pressure gradient and the seepage velocity obtained by fitting according to a formula (1) to obtain a non-Darcy osmotic coefficient K and a non-Darcy flow factor beta; then binary fitting is carried out to obtain a functional relation f (n, d) of K and the broken stone sample porosity n and the broken stone sample granularity d in the pressure-bearing seepage test and a relation g (n, d) of beta and n and d,
wherein: k-permeability coefficient; beta-fidaxb factor; porosity of crushed stone sample in n-pressure-bearing seepage test; d-particle size of crushed stone sample; a, a 1 、a 2 And a 3 Fitting a relative fitting constant of the permeability coefficient K to the porosity n and particle size d curves for the experimental data; b 1 、b 2 And b 3 The relevant fitting constants of the fidaxflow factor β with respect to the porosity n and particle size d change curve were fitted to the experimental data.
The invention has the beneficial effects that:
the method for obtaining the continuous function of the permeability coefficient of the mining stable region of the coal mine obtains experimental data such as seepage velocity, pore pressure gradient, particle size and the like by carrying out a multi-particle-size broken stone pressure-bearing seepage test, and fits the data to obtain a functional relation of the permeability coefficient K with respect to the porosity and the particle size.
Drawings
FIG. 1 is a graph of pore pressure gradient versus seepage velocity.
Detailed Description
The method for obtaining the continuous function of the permeability coefficient of the mining stable region of the coal mine comprises the following steps:
1) And collecting a target working face surrounding rock sample, measuring the uniaxial compressive strength and the density of the complete rock sample, and setting the crushing pressure-bearing strength by the uniaxial compressive strength.
2) Preparing identical lithotripsy stone samples with not less than three particle sizes and consistent with the test result of the step 1), and drying.
3) Carrying out pressure-bearing seepage test on crushed stone samples with different particle sizes: the porosity of the crushed stone sample is reduced by compression, so that the porosity change of the rock falling from different places in a mining stable region due to the action of mine pressure is simulated, and gas for testing seepage conditions is introduced into the crushed stone sample in the process of compressing the crushed stone sample; and obtaining experimental data of a plurality of groups of axial loads, pressure head displacement, air inlet pressure, air outlet pressure, air inlet flow and air outlet flow by changing the pressure head displacement and the pressure of the introduced air.
In the embodiment, the pressure-bearing seepage test is carried out on the crushed stone samples with different particle sizes, wherein the crushed stone samples are put into a device for testing the permeability parameters of crushed rock disclosed by the application number 202010014728.5, and the pressure-bearing seepage test is carried out on a GCTS test system platform.
4) After the test is finished, selecting a plurality of groups of typical data, making a change curve of pore pressure gradient and seepage velocity, comparing the change curve of pore pressure gradient and seepage velocity obtained by fitting with a relation curve of pore pressure gradient and seepage velocity expressed by a Forchheimer type non-Darcy equation, and judging whether the relation curve accords with the Forchheimer type non-Darcy equation; the Forchheimer's non-Darcy equation is as follows:
wherein J is pore pressure gradient, and the unit is Pa/m; ρ is the fluid density in kg/m 3 The method comprises the steps of carrying out a first treatment on the surface of the g is gravity acceleration, and the unit is N/kg; beta is the fidaxflow factor in m -1 The method comprises the steps of carrying out a first treatment on the surface of the v is the seepage velocity, the unit is m/s; k is the non-Darcy osmotic coefficient, the unit is m/s,wherein k is permeability in m -2 The method comprises the steps of carrying out a first treatment on the surface of the u is the dynamic coefficient of viscosity of the fluid in Pa.s.
5) After determining that the change curve of the pore pressure gradient and the seepage velocity accords with the Forchheimer type non-Darcy equation, processing the change curve of the pore pressure gradient and the seepage velocity obtained by fitting according to a formula (1) to obtain a non-Darcy osmotic coefficient K and a non-Darcy flow factor beta; then binary fitting is carried out to obtain a functional relation f (n, d) of K and the broken stone sample porosity n and the broken stone sample granularity d in the pressure-bearing seepage test and a relation g (n, d) of beta and n and d,
wherein: k-permeability coefficient; beta-fidaxb factor; porosity of crushed stone sample in n-pressure-bearing seepage test; d-particle size of crushed stone sample; a, a 1 、a 2 And a 3 Fitting a relative fitting constant of the permeability coefficient K to the porosity n and particle size d curves for the experimental data; b 1 、b 2 And b 3 The relevant fitting constants of the fidaxflow factor β with respect to the porosity n and particle size d change curve were fitted to the experimental data.
The permeability coefficient K obtained in the embodiment is more in line with the actual situation, has positive effects on researching the coal bed gas migration rule in the extraction process of the mining stabilization zone, and can help to select the position of the gas extraction well of the mining stabilization zone of the coal mine.
Finally, it is noted that the above embodiments are only for illustrating the technical solution of the present invention and not for limiting the same, and although the present invention has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that modifications and equivalents may be made thereto without departing from the spirit and scope of the technical solution of the present invention, which is intended to be covered by the scope of the claims of the present invention.

Claims (1)

1. The method for obtaining the continuous function of the permeability coefficient of the mining stable region of the coal mine is characterized by comprising the following steps of:
1) Collecting a target working face surrounding rock sample, measuring the uniaxial compressive strength and the uniaxial compressive strength of the complete rock sample, and setting the crushing pressure-bearing strength by the uniaxial compressive strength;
2) Preparing identical lithotripsy stone samples with the particle sizes not less than three and consistent with the test results of the step 1), and drying;
3) Carrying out pressure-bearing seepage test on crushed stone samples with different particle sizes: the porosity of the crushed stone sample is reduced by compression, so that the porosity change of the rock falling from different places in a mining stable region due to the action of mine pressure is simulated, and gas for testing seepage conditions is introduced into the crushed stone sample in the process of compressing the crushed stone sample; obtaining experimental data of a plurality of groups of axial loads, pressure head displacement, air inlet pressure, air outlet pressure, air inlet flow and air outlet flow by changing the pressure head displacement and the pressure of the introduced gas;
4) After the test is finished, selecting a plurality of groups of typical data, making a change curve of pore pressure gradient and seepage velocity, comparing the change curve of pore pressure gradient and seepage velocity obtained by fitting with a relation curve of pore pressure gradient and seepage velocity expressed by a Forchheimer type non-Darcy equation, and judging whether the relation curve accords with the Forchheimer type non-Darcy equation; the Forchheimer's non-Darcy equation is as follows:
wherein J is pore pressure gradient, and the unit is Pa/m; ρ is the fluid density in kg/m 3 The method comprises the steps of carrying out a first treatment on the surface of the g is gravity acceleration, and the unit is N/kg; beta is the fidaxflow factor in m -1 The method comprises the steps of carrying out a first treatment on the surface of the v is the seepage velocity, the unit is m/s; k is the non-Darcy osmotic coefficient, the unit is m/s,wherein k is permeability in m -2 The method comprises the steps of carrying out a first treatment on the surface of the u is the dynamic coefficient of viscosity of the fluid, and the unit is Pa.s;
5) After determining that the change curve of the pore pressure gradient and the seepage velocity accords with the Forchheimer type non-Darcy equation, processing the change curve of the pore pressure gradient and the seepage velocity obtained by fitting according to a formula (1) to obtain a non-Darcy osmotic coefficient K and a non-Darcy flow factor beta; then binary fitting is carried out to obtain a functional relation f (n, d) of K and the broken stone sample porosity n and the broken stone sample granularity d in the pressure-bearing seepage test and a relation g (n, d) of beta and n and d,
wherein: k-permeability coefficient; beta-fidaxb factor; porosity of crushed stone sample in n-pressure-bearing seepage test; d-particle size of crushed stone sample; a, a 1 、a 2 And a 3 Fitting the permeability coefficient K to the variation of the porosity n and the particle size d by means of the test dataA correlation fit constant for the curve; b 1 、b 2 And b 3 The relevant fitting constants of the fidaxflow factor β with respect to the porosity n and particle size d change curve were fitted to the experimental data.
CN202011269283.1A 2020-11-13 2020-11-13 Method for obtaining continuous function of permeability coefficient of coal mine mining stable region Active CN112362556B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011269283.1A CN112362556B (en) 2020-11-13 2020-11-13 Method for obtaining continuous function of permeability coefficient of coal mine mining stable region

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011269283.1A CN112362556B (en) 2020-11-13 2020-11-13 Method for obtaining continuous function of permeability coefficient of coal mine mining stable region

Publications (2)

Publication Number Publication Date
CN112362556A CN112362556A (en) 2021-02-12
CN112362556B true CN112362556B (en) 2024-03-29

Family

ID=74514955

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011269283.1A Active CN112362556B (en) 2020-11-13 2020-11-13 Method for obtaining continuous function of permeability coefficient of coal mine mining stable region

Country Status (1)

Country Link
CN (1) CN112362556B (en)

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102323196A (en) * 2011-05-23 2012-01-18 申坤 Design manufacturing and application of triaxial penetration tester
CN102620996A (en) * 2012-04-11 2012-08-01 江苏师范大学 Operational method for measuring creep parameters and infiltration parameters of smashed rocks simultaneously
CN106706492A (en) * 2016-11-21 2017-05-24 大连理工大学 Penetration device for researching porous medium whole-zone penetration mechanism
CN107462936A (en) * 2017-08-28 2017-12-12 中国石油大学(北京) Utilize the method for pressure monitoring Data Inversion low permeability reservoir non-Darcy percolation law
WO2017217565A1 (en) * 2016-06-14 2017-12-21 한국가스공사 Method for calculating permeability of porous medium by means of tortuous hydraulic diameter
CN107976395A (en) * 2017-12-21 2018-05-01 西安科技大学 CO in loose coal2Osmosis system measurement device and its method
CN108614910A (en) * 2018-03-06 2018-10-02 江西理工大学 The computational methods of ion type rareearth mine in_situ leaching critical groundwater table
CN109612907A (en) * 2018-12-27 2019-04-12 太原理工大学 Fractured coal and rock testing permeability experimental rig and method
CN110119574A (en) * 2019-05-14 2019-08-13 中国矿业大学 A kind of non linear fluid flow through porous medium system stability distinguishing method of filling coal mining water-resisting key strata
CN110940610A (en) * 2019-11-27 2020-03-31 山东科技大学 Broken rock nonlinear seepage test system and method
CN111189754A (en) * 2020-01-07 2020-05-22 重庆大学 Device for testing permeability parameters of broken rock
CN111189755A (en) * 2020-01-07 2020-05-22 重庆大学 Numerical simulation obtaining method for effective storage space of coal mine mining stable area gas
KR20200098420A (en) * 2019-02-11 2020-08-20 한국가스공사 Method of calculating permeability variation of porous material using pore-scale simulation
CN111753414A (en) * 2020-06-16 2020-10-09 河海大学 Seepage characteristic simulation method for rock progressive destruction process based on non-Darcy law

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2955662B1 (en) * 2010-01-22 2014-08-22 Total Sa MEASUREMENT OF PARAMETERS RELATING TO THE FLOW OF FLUIDS IN POROUS MATERIAL
EP2791858A4 (en) * 2011-12-16 2016-11-16 Landmark Graphics Corp System and method for simulation of gas desorption in a reservoir using a multi-porosity approach

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102323196A (en) * 2011-05-23 2012-01-18 申坤 Design manufacturing and application of triaxial penetration tester
CN102620996A (en) * 2012-04-11 2012-08-01 江苏师范大学 Operational method for measuring creep parameters and infiltration parameters of smashed rocks simultaneously
WO2017217565A1 (en) * 2016-06-14 2017-12-21 한국가스공사 Method for calculating permeability of porous medium by means of tortuous hydraulic diameter
CN106706492A (en) * 2016-11-21 2017-05-24 大连理工大学 Penetration device for researching porous medium whole-zone penetration mechanism
CN107462936A (en) * 2017-08-28 2017-12-12 中国石油大学(北京) Utilize the method for pressure monitoring Data Inversion low permeability reservoir non-Darcy percolation law
CN107976395A (en) * 2017-12-21 2018-05-01 西安科技大学 CO in loose coal2Osmosis system measurement device and its method
CN108614910A (en) * 2018-03-06 2018-10-02 江西理工大学 The computational methods of ion type rareearth mine in_situ leaching critical groundwater table
CN109612907A (en) * 2018-12-27 2019-04-12 太原理工大学 Fractured coal and rock testing permeability experimental rig and method
KR20200098420A (en) * 2019-02-11 2020-08-20 한국가스공사 Method of calculating permeability variation of porous material using pore-scale simulation
CN110119574A (en) * 2019-05-14 2019-08-13 中国矿业大学 A kind of non linear fluid flow through porous medium system stability distinguishing method of filling coal mining water-resisting key strata
CN110940610A (en) * 2019-11-27 2020-03-31 山东科技大学 Broken rock nonlinear seepage test system and method
CN111189754A (en) * 2020-01-07 2020-05-22 重庆大学 Device for testing permeability parameters of broken rock
CN111189755A (en) * 2020-01-07 2020-05-22 重庆大学 Numerical simulation obtaining method for effective storage space of coal mine mining stable area gas
CN111753414A (en) * 2020-06-16 2020-10-09 河海大学 Seepage characteristic simulation method for rock progressive destruction process based on non-Darcy law

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
废弃煤矿采空区抽水蓄能水库初步可行性研究;王丽;重庆大学学报;20200430;47-54 *
破碎花岗岩非达西渗流的试验研究;贺香兰;周佳庆;魏凯;王敏;;中国农村水利水电;20170915(09);155-160 *
范鹏宏 ; 聂百胜 ; .含瓦斯煤Forchheimer型非达西渗流实验研究.矿业安全与环保.(06),23-26. *

Also Published As

Publication number Publication date
CN112362556A (en) 2021-02-12

Similar Documents

Publication Publication Date Title
US6725930B2 (en) Conductive proppant and method of hydraulic fracturing using the same
AU2015345707B2 (en) Method for determining earth surface interpenetrated crack distribution and air leakage characteristics in shallow burial coal mining
Xiao et al. Experimental investigation of the effects of suffusion on physical and geomechanic characteristics of sandy soils
Dabbous et al. The permeability of coal to gas and water
US7562583B2 (en) Apparatus and method for measuring crush-resistance of granular materials
Feng et al. Experimental investigation on the seepage property of saturated broken red sandstone of continuous gradation
Kim et al. Soil conditioning of weathered granite soil used for EPB shield TBM: A laboratory scale study
Ma et al. Seepage properties of crushed coal particles
Ma et al. Effect of particle mixture on seepage properties of crushed mudstones
CN112727534B (en) Gas drilling hole arrangement method based on true three-dimensional stress and permeability dynamic change
Ma et al. Impact of particle transfer on flow properties of crushed mudstones
Klimov et al. Experimental study of the influence of a triaxial stress state with unequal components on rock permeability
CN108153947B (en) Soft rock field state-based disintegration characteristic description method
CN112461657B (en) Rapid prediction method for critical damage stress of roadbed soil
Pang et al. Measurement of the coefficient of seepage characteristics in pore-crushed coal bodies around gas extraction boreholes
Vo et al. Cohesion and suction induced hang-up in ore passes
CN112362556B (en) Method for obtaining continuous function of permeability coefficient of coal mine mining stable region
Davudov* et al. Impact of pore connectivity and topology on gas productivity in Barnett and Haynesville Shale plays
Arasteh et al. A new model for calculation of the plastic compression index and porosity and permeability of gob materials in longwall mining
Zhang et al. Experimental study on seepage properties of postpeak fractured rocks under cyclic loading-unloading confining stress and axial stress
CN111946321B (en) Proppant parameter design method for sand-filling temporary plugging fracturing
CN114075974B (en) Temporary plugging steering fracturing method, device, control equipment and storage medium
CN112417663B (en) Method for representing coal bed gas migration path in extraction process of coal mining stable area
CN114611259A (en) Method for evaluating effective reconstruction volume of multi-scale artificial seam
Shilova et al. Estimation of propped fracture compressibility in coal seam

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant