CN112362556A - Method for obtaining permeability coefficient continuous function of coal mine mining stable area - Google Patents

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

Info

Publication number
CN112362556A
CN112362556A CN202011269283.1A CN202011269283A CN112362556A CN 112362556 A CN112362556 A CN 112362556A CN 202011269283 A CN202011269283 A CN 202011269283A CN 112362556 A CN112362556 A CN 112362556A
Authority
CN
China
Prior art keywords
crushed stone
pressure
seepage
darcy
permeability coefficient
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011269283.1A
Other languages
Chinese (zh)
Other versions
CN112362556B (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

Images

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 permeability coefficient continuous function of a coal mining stable area, which comprises the following steps: 1) collecting a surrounding rock sample of a target working surface, measuring the uniaxial compressive strength and density of the complete rock sample, and setting the pressure-bearing strength of the crushed stone according to the uniaxial compressive strength; 2) preparing not less than three same lithologic crushed stone samples with grain sizes, and drying; 3) carrying out pressure-bearing seepage tests 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 the pore pressure gradient and the 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 permeability coefficient and a non-Darcy flow factor; and then carrying out binary fitting 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.

Description

Method for obtaining permeability coefficient continuous function of coal mine mining stable area
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 particularly relates to a method for obtaining a permeability coefficient continuous function of a coal mine mining stable region, which can be directly applied to selection of a ground position of a gas extraction well of the coal mine mining stable region.
Background
In the research of the coal bed gas migration rule in the extraction process of the mining stable area, the distribution condition of the permeability coefficient is a key parameter for research. At present, the permeability coefficient is mostly selected to be an even distribution function or a sectional even distribution function, which is different from the actual situation to a certain extent, and the movement rule of the coal bed methane obtained by research is possibly inaccurate.
Disclosure of Invention
In view of the above, the present invention provides a method for obtaining a continuous function of permeability coefficients in a coal mining stable area, so as to solve the technical problem of inaccurate permeability coefficient selection.
The invention discloses a method for obtaining a permeability coefficient continuous function of a coal mining stable area, which comprises the following steps:
1) collecting a surrounding rock sample of a target working surface, measuring the uniaxial compressive strength and density of the complete rock sample, and setting the pressure-bearing strength of the crushed stone according to the uniaxial compressive strength;
2) preparing not less than three lithotripsy samples with the same particle size and according 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 compressing the crushed stone sample so as to simulate the porosity change of the caving rock in different places of the mining stable area due to the action of mine pressure, and gas for testing seepage conditions is introduced into the crushed stone sample in the process of compressing the crushed stone sample; the method comprises the steps that 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 are obtained by changing the displacement of a pressure head and the pressure of introduced air;
4) after the test is finished, selecting a plurality of groups of typical data, making a change curve of the pore pressure gradient and the seepage velocity, comparing the change curve of the pore pressure gradient and the seepage velocity obtained by fitting with a pore pressure gradient and seepage velocity relation curve expressed by a Forchheimer type non-Darcy equation, and judging whether the change curve accords with the Forchheimer type non-Darcy equation; the Forchheimer type NodeDarcy equation is as follows:
Figure BDA0002777119100000021
wherein J is a pore pressure gradient in Pa/m; rho is the fluid density in kg/m3(ii) a g is the gravity acceleration with the unit of N/kg; beta is a factor of Fidexi flow in m-1(ii) a v is the seepage velocity in m/s; k is the coefficient of permeability of the Darcy, the unit is m/s,
Figure BDA0002777119100000022
where k is the permeability in m-2(ii) a u is the dynamic viscosity coefficient 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 conforms to a 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 permeability coefficient K and a non-Darcy flow factor beta; then binary fitting is carried out to obtain a function relational expression f (n, d) of K, the porosity n of the crushed stone sample and the granularity d of the crushed stone sample in the pressure-bearing seepage test and a relational expression g (n, d) of beta, n and d,
Figure BDA0002777119100000023
Figure BDA0002777119100000024
in the formula: k-permeability coefficient; beta-noncardiac flow factor; n-the porosity of the crushed stone sample in the pressure bearing seepage test; d-crushed stone sample granularity; a is1、a2And a3Fitting the phase of the curve of the permeability coefficient K with respect to the variation of the porosity n and the particle size d to experimental dataA fitting constant is closed; b1、b2And b3The relevant fitting constants for the curves of the non-darcy flow factor beta with respect to the variation of porosity n and particle size d are fitted by experimental data.
The invention has the beneficial effects that:
the invention relates to a method for obtaining a continuous function of permeability coefficient of a coal mining stable area, which comprises the steps of carrying out a multi-particle-size crushed stone pressure-bearing seepage test to obtain experimental data such as seepage velocity, pore pressure gradient and particle size, fitting the data to obtain a function relation of the permeability coefficient K on porosity and particle size.
Drawings
FIG. 1 is a graph of pore pressure gradient versus seepage velocity.
Detailed Description
The method for obtaining the permeability coefficient continuous function of the coal mining stable area comprises the following steps:
1) and collecting a surrounding rock sample of the target working surface, measuring the uniaxial compressive strength and the density of the complete rock sample, and setting the pressure-bearing strength of the crushed stone according to the uniaxial compressive strength.
2) Preparing the same lithotripsy sample with not less than three particle sizes which is consistent with the test result of the step 1), and drying.
3) Carrying out pressure-bearing seepage tests on gravel samples with different particle sizes: the porosity of the crushed stone sample is reduced by compressing the crushed stone sample so as to simulate the porosity change of the caving rock in different places of the mining stable area due to the action of mine pressure, and gas for testing seepage conditions is introduced into the crushed stone sample in the process of compressing the crushed stone sample; the 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 are obtained by changing the pressure head displacement and the pressure of the introduced air.
In the embodiment, the pressure-bearing seepage test of the crushed stone samples with different particle sizes is to load the crushed stone samples into the device for testing the permeability parameters of the crushed rock disclosed by the application number 202010014728.5 and carry out the pressure-bearing seepage test 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 the pore pressure gradient and the seepage velocity, comparing the change curve of the pore pressure gradient and the seepage velocity obtained by fitting with a pore pressure gradient and seepage velocity relation curve expressed by a Forchheimer type non-Darcy equation, and judging whether the change curve accords with the Forchheimer type non-Darcy equation; the Forchheimer type NodeDarcy equation is as follows:
Figure BDA0002777119100000031
wherein J is a pore pressure gradient in Pa/m; rho is the fluid density in kg/m3(ii) a g is the gravity acceleration with the unit of N/kg; beta is a factor of Fidexi flow in m-1(ii) a v is the seepage velocity in m/s; k is the coefficient of permeability of the Darcy, the unit is m/s,
Figure BDA0002777119100000041
where k is the permeability in m-2(ii) a u is the dynamic viscosity coefficient of the fluid, in Pa · s.
5) After determining that the change curve of the pore pressure gradient and the seepage velocity conforms to a 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 permeability coefficient K and a non-Darcy flow factor beta; then binary fitting is carried out to obtain a function relational expression f (n, d) of K, the porosity n of the crushed stone sample and the granularity d of the crushed stone sample in the pressure-bearing seepage test and a relational expression g (n, d) of beta, n and d,
Figure BDA0002777119100000042
Figure BDA0002777119100000043
in the formula: k-permeability coefficient; beta-noncardiac flow factor; n-the porosity of the crushed stone sample in the pressure bearing seepage test; d-crushed stone sample granularity; a is1、a2And a3Fitting the relevant fitting constants of the permeability coefficient K on the variation curves of the porosity n and the particle size d through the test data; b1、b2And b3The relevant fitting constants for the curves of the non-darcy flow factor beta with respect to the variation of porosity n and particle size d are fitted by experimental data.
The permeability coefficient K obtained in the embodiment is more consistent with the actual situation, has a positive effect on the research of the coal bed gas migration rule in the extraction process of the mining stable area, and can help to select the position of the gas extraction well of the mining stable area of the coal mine.
Finally, the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting, 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 or equivalent substitutions may be made to the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention, and all of them should be covered in the claims of the present invention.

Claims (1)

1. The method for obtaining the permeability coefficient continuous function of the coal mining stable area is characterized by comprising the following steps of:
1) collecting a surrounding rock sample of a target working surface, measuring the uniaxial compressive strength and density of the complete rock sample, and setting the pressure-bearing strength of the crushed stone according to the uniaxial compressive strength;
2) preparing not less than three lithotripsy samples with the same particle size and according with the test result of the step 1), and drying;
3) carrying out pressure-bearing seepage tests on gravel samples with different particle sizes: the porosity of the crushed stone sample is reduced by compressing the crushed stone sample so as to simulate the porosity change of the caving rock in different places of the mining stable area due to the action of mine pressure, and gas for testing seepage conditions is introduced into the crushed stone sample in the process of compressing the crushed stone sample; the method comprises the steps that 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 are obtained by changing the displacement of a pressure head and the pressure of introduced air;
4) after the test is finished, selecting a plurality of groups of typical data, making a change curve of the pore pressure gradient and the seepage velocity, comparing the change curve of the pore pressure gradient and the seepage velocity obtained by fitting with a pore pressure gradient and seepage velocity relation curve expressed by a Forchheimer type non-Darcy equation, and judging whether the change curve accords with the Forchheimer type non-Darcy equation; the Forchheimer type NodeDarcy equation is as follows:
Figure FDA0002777119090000011
wherein J is a pore pressure gradient in Pa/m; rho is the fluid density in kg/m3(ii) a g is the gravity acceleration with the unit of N/kg; beta is a factor of Fidexi flow in m-1(ii) a v is the seepage velocity in m/s; k is the coefficient of permeability of the Darcy, the unit is m/s,
Figure FDA0002777119090000012
where k is the permeability in m-2(ii) a u is the dynamic viscosity coefficient 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 conforms to a 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 permeability coefficient K and a non-Darcy flow factor beta; then binary fitting is carried out to obtain a function relational expression f (n, d) of K, the porosity n of the crushed stone sample and the granularity d of the crushed stone sample in the pressure-bearing seepage test and a relational expression g (n, d) of beta, n and d,
Figure FDA0002777119090000021
Figure FDA0002777119090000022
in the formula: k-permeability coefficient; beta-noncardiac flow factor; n-the porosity of the crushed stone sample in the pressure bearing seepage test; d-crushed stone sample granularity; a is1、a2And a3Fitting the relevant fitting constants of the permeability coefficient K on the variation curves of the porosity n and the particle size d through the test data; b1、b2And b3The relevant fitting constants for the curves of the non-darcy flow factor beta with respect to the variation of porosity n and particle size d are fitted by 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 true CN112362556A (en) 2021-02-12
CN112362556B 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 (16)

* 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
US20130054157A1 (en) * 2010-01-22 2013-02-28 Centre National De La Recherche Scientifique- Cnrs Measurement of parameters linked to the flow of fluids in a porous material
US20140350906A1 (en) * 2011-12-16 2014-11-27 Landmark Graphics Corporation System and Method for Simulation of Gas Desorption in a Reservoir Using a Multi-Porosity Approach
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
CN111189755A (en) * 2020-01-07 2020-05-22 重庆大学 Numerical simulation obtaining method for effective storage space of coal mine mining stable area gas
CN111189754A (en) * 2020-01-07 2020-05-22 重庆大学 Device for testing permeability parameters of broken rock
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

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130054157A1 (en) * 2010-01-22 2013-02-28 Centre National De La Recherche Scientifique- Cnrs Measurement of parameters linked to the flow of fluids in a porous material
CN102323196A (en) * 2011-05-23 2012-01-18 申坤 Design manufacturing and application of triaxial penetration tester
US20140350906A1 (en) * 2011-12-16 2014-11-27 Landmark Graphics Corporation System and Method for Simulation of Gas Desorption in a Reservoir Using a Multi-Porosity Approach
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
CN111189755A (en) * 2020-01-07 2020-05-22 重庆大学 Numerical simulation obtaining method for effective storage space of coal mine mining stable area gas
CN111189754A (en) * 2020-01-07 2020-05-22 重庆大学 Device for testing permeability parameters of broken rock
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 (2)

* Cited by examiner, † Cited by third party
Title
王丽: "废弃煤矿采空区抽水蓄能水库初步可行性研究", 重庆大学学报, 30 April 2020 (2020-04-30), pages 47 - 54 *
贺香兰;周佳庆;魏凯;王敏;: "破碎花岗岩非达西渗流的试验研究", 中国农村水利水电, no. 09, 15 September 2017 (2017-09-15), pages 155 - 160 *

Also Published As

Publication number Publication date
CN112362556B (en) 2024-03-29

Similar Documents

Publication Publication Date Title
Jasinge et al. Effects of effective stress changes on permeability of latrobe valley brown coal
Davudov et al. Impact of pore compressibility and connectivity loss on shale permeability
Ma et al. An experimental investigation of permeability measurement of water flow in crushed rocks
Ma et al. Compaction and seepage properties of crushed limestone particle mixture: an experimental investigation for Ordovician karst collapse pillar groundwater inrush
CN108590601B (en) Experimental method for optimizing water injection expansion construction parameters
Bohnsack et al. Stress sensitivity of porosity and permeability under varying hydrostatic stress conditions for different carbonate rock types of the geothermal Malm reservoir in Southern Germany
Feng et al. Experimental investigation on the seepage property of saturated broken red sandstone of continuous gradation
Ma et al. Effect of particle mixture on seepage properties of crushed mudstones
CN107462936B (en) Utilize the method for pressure monitoring Data Inversion low permeability reservoir non-Darcy percolation law
CN112727534B (en) Gas drilling hole arrangement method based on true three-dimensional stress and permeability dynamic change
CN110501758B (en) Glutenite reservoir longitudinal continuous brittleness index prediction method
Klimov et al. Experimental study of the influence of a triaxial stress state with unequal components on rock permeability
Cao et al. Reservoir petrofacies—A tool for characterization of reservoir quality and pore structures in a tight sandstone reservoir: A study from the sixth member of Upper Triassic Yanchang Formation, Ordos Basin, China
M'Jahad et al. Characterization of transport and water retention properties of damaged Callovo-Oxfordian claystone
Davudov* et al. Impact of pore connectivity and topology on gas productivity in Barnett and Haynesville Shale plays
Toll et al. Triaxial testing of unsaturated samples of undisturbed residual soil from Singapore
CN112362556A (en) Method for obtaining permeability coefficient continuous function of coal mine mining stable area
Papamichos et al. A volumetric sand production experiment
Popov et al. Experimental and Numerical Assessment of the Influence of Bottomhole Pressure Drawdown on Terrigenous Reservoir Permeability and Well Productivity.
Shilova et al. Estimation of propped fracture compressibility in coal seam
Cheshomi et al. An empirical approach to estimate EPMT and PL of silty clays based on SPT
CN112417663B (en) Method for representing coal bed gas migration path in extraction process of coal mining stable area
CN114819314B (en) Multi-factor weakly cemented overburden water-guiding fracture height prediction method
Zhou et al. Quantitative Prediction of Sand Production in Weakly Consolidated Sandstone Through Numerical Simulation
Rybalkin et al. Effect of compressive stresses on permeability of coal: Experimental and modeling estimates

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