CN110533349A - A kind of coal seam gas-bearing capacity calculates and error analysis method - Google Patents
A kind of coal seam gas-bearing capacity calculates and error analysis method Download PDFInfo
- Publication number
- CN110533349A CN110533349A CN201910956488.8A CN201910956488A CN110533349A CN 110533349 A CN110533349 A CN 110533349A CN 201910956488 A CN201910956488 A CN 201910956488A CN 110533349 A CN110533349 A CN 110533349A
- Authority
- CN
- China
- Prior art keywords
- desorption
- gas
- solution
- time
- equation
- 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
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 71
- 239000003245 coal Substances 0.000 title claims abstract description 71
- 238000003795 desorption Methods 0.000 claims abstract description 196
- 238000000034 method Methods 0.000 claims abstract description 64
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 24
- 238000004364 calculation method Methods 0.000 claims abstract description 12
- 238000004088 simulation Methods 0.000 claims description 35
- 238000012937 correction Methods 0.000 claims description 10
- 241001074085 Scophthalmus aquosus Species 0.000 claims description 9
- 238000005553 drilling Methods 0.000 claims description 8
- 230000001186 cumulative effect Effects 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- 238000001179 sorption measurement Methods 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 238000007599 discharging Methods 0.000 claims description 3
- 238000010438 heat treatment Methods 0.000 claims description 3
- ONUFESLQCSAYKA-UHFFFAOYSA-N iprodione Chemical compound O=C1N(C(=O)NC(C)C)CC(=O)N1C1=CC(Cl)=CC(Cl)=C1 ONUFESLQCSAYKA-UHFFFAOYSA-N 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 7
- 230000000875 corresponding effect Effects 0.000 description 18
- 238000007872 degassing Methods 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 4
- 238000009825 accumulation Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N7/00—Analysing materials by measuring the pressure or volume of a gas or vapour
- G01N7/02—Analysing materials by measuring the pressure or volume of a gas or vapour by absorption, adsorption, or combustion of components and measurement of the change in pressure or volume of the remainder
- G01N7/04—Analysing materials by measuring the pressure or volume of a gas or vapour by absorption, adsorption, or combustion of components and measurement of the change in pressure or volume of the remainder by absorption or adsorption alone
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Operations Research (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Mathematical Optimization (AREA)
- Development Economics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Marketing (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Data Mining & Analysis (AREA)
- Primary Health Care (AREA)
- Mining & Mineral Resources (AREA)
- Marine Sciences & Fisheries (AREA)
- Game Theory and Decision Science (AREA)
- Animal Husbandry (AREA)
- Quality & Reliability (AREA)
- Agronomy & Crop Science (AREA)
- Pathology (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Algebra (AREA)
- Combustion & Propulsion (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
Abstract
The invention discloses a kind of Coal Bed Methane Content Estimation Method and error analysis method, calculation method includes: that S1 determines site desorption equation V2=ABT0 1/2/(1+B·T0 1/2)‑V1 Solution, the S2 coal heart desorption, S3 equation model;Error analysis method includes the error analysis of remaining desorption quantity and the error analysis of gas total content.The present invention propose site desorption equation combine specific live desorption of mash gas process can gas total content, gas loss amount and remaining desorption quantity to deep fractures calculate;Reduce gas bearing capacity relative error;Corresponding error analysis method is provided simultaneously to assist site desorption equation to carry out more accurate gas bearing capacity and calculate, to improve the accuracy and reliability of deep fractures gas bearing capacity calculating, be further conducive to make objective appraisal to deep gas bearing capacity and hosting geological condition.
Description
Technical field
The invention belongs to coal seam gas-bearing capacity determination techniques field, in particular to a kind of deep fractures gas bearing capacity calculate with
Error analysis method.
Background technique
Coal seam gas-bearing capacity is one of coal-bed gas major parameter, is important base necessary to coal and CBM exploration and development
Plinth parameter.
Ground is surveyed period coal seam gas-bearing capacity and is made of four parts: gas loss amount (V1), desorption of mash gas amount (V2), crush before
Heating in vacuum amount of gas evolved (V3) and crush after amount of gas evolved (V4).Deaerate the sum of gas amount after the preceding degassing gas amount of crushing and fine ore
It can be described as gas remaining quantity.Wherein desorption of mash gas amount is measured by site desorption tank, and gas remaining quantity passes through experimental determination, and
Gas loss amount mainly passes through coal heart site desorption rule and calculates to obtain.
Surveying the calculation method of period coal seam gas-bearing capacity domesticly at present is mainlyMethod (graphical method and analytic method): it is
According to drilling the coal heart start expose a period of time t in add up desorption quantity V andIt is in linear relation to be determined.It will start to solve
It inhales each point line in linear relation in a period of time and extends and intersect with ordinate, the intercept of straight line on longitudinal axis is required
Gas loss amount;Or the coordinate value put according to these, gas loss amount can also be found out by least square method.For graphical method,
Be withFor abscissa, whole measuring points are marked and drawed on graph paper, will start to desorb each point in linear relation in a period of time
Line, and extend related to axis of ordinates;The intercept of straight line on longitudinal axis is required gas loss amount.For analytic method,
Because coal sample start expose a period of time in, V withIt is in linear relation, i.e.,In formula a, b be etc. permanent number, whenWhen, V=a, a value is required gas loss amount.
Deep fractures are influenced by high geostress field, temperature field and Fluid Pressure Field, and hosting geological condition is complicated.It is shallow with routine
Portion compares in coal seam, when acquiring deep fractures product, gas loss time advance in drilling process, and conventional gas bearing capacity test
Method does not account for the loss amount of this part free gas when calculating loss tolerance, relatively low so as to cause test result.According to text
It offers report: when coal seam buried depth is less than 500m, there is 70% measured value relatively low 15~25%, 20% measured value higher 10~
15%;Coal seam buried depth be greater than 500m (when especially close to 800m), generally have the relatively low degree of measured value increase with hole depth and
The trend of increase has 85% measured value relatively low 30~40%, up to 50% or more, only 8% measured value higher 5~10%.
As shown in Fig. 2, correlative study show this method for its calculated value of deep fractures can less than normal, due to acquisition gas bearing capacity value and
For deep gas real content there are certain error, constrain deep gas bearing capacity and hosting geological condition objectively evaluates effect.
Summary of the invention
The present invention the above-mentioned prior art there are aiming at the problem that, a kind of coal seam gas-bearing capacity is provided and is calculated and error analysis
Method.
The present invention is realized by following technological means solves above-mentioned technical problem:
A kind of Coal Bed Methane Content Estimation Method, comprising the following steps:
S1 determines site desorption equation:
V2=ABT0 1/2/(1+B·T0 1/2)-V1 Solution
Wherein, A is the upper limit desorption of mash gas total amount of coal heart desorption, and B is the coefficient in relation to gas adsorption speed and heat of adsorption,
V1 SolutionFor the gas loss amount of solution, V2Add up desorption quantity, T for gas0=t0+ t, t0=1/2t1+t2,
V3 Solution=A-V1 Solution-V2, V3 SolutionFor remaining desorption quantity calculated value;
The desorption of the S2 coal heart: drilling as mud medium or water medium, when the coal heart that drills is mentioned to aperture half depth, as coal
The initial desorption time of the heart.It will mention and bore coal heart progress site desorption, during which, record, which mentions, bores time t1, tinning time t2, remember simultaneously
Under desorption time t and the accumulative desorption of mash gas amount V of t at any time in solution cucurbitula2;
S3 equation model: the gas of coal core sample at any time is added up into desorption quantity V2Throwing is drawn in a coordinate system, to data point
The fitting for carrying out site desorption equation, determines A and V1 SolutionValue, then A be coal-bed gas total content calculated value, V1 SolutionFor coal seam watt
The calculated value of this loss amount.
Further, be realize the gas bearing capacity based on field measurement data accurate calculating, by laboratory through crushing
The preceding remaining desorption quantity V to deaerate and acquisition is deaerated after crushing2 It is realIt is painted on graph paper as a part and throwing for accumulating desorption quantity,
Fitting is corrected to equation.
Further, the site desorption equation is suitable for different desorption rates, different accumulative desorption duration and not simultaneously
With the gas site desorption data of the modes such as observation interval.
Further, Data correction is carried out using single point correction or two-point calibration;
The method of the single point correction are as follows:
(1) establish coordinate (V2+V3+V4), throwing is drawn on graph paper, wherein V3It is de- to crush initial vacuum heating
Tolerance, V4For amount of gas evolved after crushing, V3+V4It is considered as in solution cucurbitula what uninterrupted desorption at any time continuously obtained, corresponding parsing
Time is ti;
(2) formula V is used2=ABT0 1/2/(1+B·T0 1/2)-V1 SolutionTo coordinate points (V2+V3+V4) carry out instead
Multiple fitting optimization, make coordinate points (V2+V3+V4) meet formula V as accurately as possible2=ABT0 1/2/(1+B·T0 1 /2)-V1 SolutionThe extension trend of curve;
(3) the desorption time coordinate for finally obtaining optimal fitting result, obtains A, B parameter value and the condition of fit equation
Under V1 Solution, V1 SolutionThe then gas loss amount to solve;
The method of the two-point calibration are as follows:
(1) establish coordinate (V2+V3) and (V2+V3+V4), throwing is drawn on graph paper, wherein V3For
Crush initial vacuum heat de-airing amount, V4For crush after amount of gas evolved, coordinate (V2+V3) correspond to V3It is considered as accumulative desorption quantity
Component part, the corresponding parsing time be ti-1, coordinate (V2+V3+V4) correspond to V3+V4It is considered as accumulative desorption quantity
Composition portion, corresponding parsing time are ti;
(2) formula V is used2=ABT0 1/2/(1+B·T0 1/2)-V1 SolutionTo coordinate points (V2+V3) and (V2+V3+V4) carry out repeatedly fitting optimization, make coordinate points (V2+V3) and (V2+V3+V4) to the greatest extent
Meet to possible accuracy formula V2=ABT0 1/2/(1+B·T0 1/2)-V1 SolutionThe extension trend of curve;
(3) the desorption time coordinate for finally obtaining optimal fitting result, obtains A, B parameter value and the condition of fit equation
Under V1 Solution, V1 SolutionThe then gas loss amount to solve;
Further, desorption of mash gas can be calculated to V2With V3+V4When corresponding theoretical observation interval be tΔ=ti-
t.At this point, t is cumulative maximum desorption quantity V2Corresponding accumulative desorption observation time.
The present invention also provides the error analysis method that a kind of coal seam gas-bearing capacity calculates, the error including remaining desorption quantity
The error analysis of analysis and gas total content;
Wherein, the error analysis equation of the remaining desorption quantity are as follows:
RE=100% × (V3 Solution-V3 It is real)/V3 It is real;
The error analysis equation of the gas total content are as follows:
R=100% × (A-VAlways Solution)/VAlways Solution, wherein VAlways Solution=V1 Solution+V2+V3 It is real。
It further, further include the different gas loss amount error analyses for mentioning and boring the time;
Wherein, the difference mentions the gas loss amount error analysis equation for boring the time are as follows:
RE=100% × [V0(t1’)-V0(t1)]/V0(t1);
V0(t1) it is that the smallest simulation mentions brill time t1' coal-bed gas loss amount, i.e. actual field desorption true exposure when
Between t1Under the conditions of according to site desorption equation V2=ABT0 1/2/(1+B·T0 1/2)-V1 SolutionThe V of acquisition1 SolutionValue;
V0(t1') it is that simulation borehole simulation mentions brill time point t1' when, can in the hope of exposure duration t1Coal-bed gas damage
Vector;
Wherein, V0(t1')=V1’(t1’)-V0’(t1');
In formula, V0’(t1') be simulation mention bore the time be t1' when reality coal seams gas discharging quantity;
V1’(t1') be simulation mention bore the time be t1' when gas loss amount, by site desorption equation V2=ABT0 1/2/(1
+B·T0 1/2)-V1 Solution, mentioned according to simulation and bore time t1' solve;
T0'=t0'+t, T0' it is total desorption time;
t0'=1/2t1’+t2 t0' be coal core sample board design the calculating loss amount time, t2For the tinning time;
t1'=t1+2×T0, t1' be simulation the coal heart mention be drilled into reach the earth's surface time, T0For in t1' under the conditions of the coal heart solving
The time desorbed in cucurbitula.
It further, further include loss amount caused by different number desorption time point and upper limit desorption quantity relative error minute
Analysis;
Wherein, loss amount relative error analysis equation caused by the different number desorption time point are as follows:
REV=100% × (V1 k-V1 n)/V1 n;
Wherein, upper limit desorption quantity relative error analysis equation caused by the different number desorption time point are as follows:
REA=100% × (Ak-An)/An;
In formula, V1 kAnd AkRespectively correspond desorption time point tnGas loss amount and upper limit desorption quantity;All in accordance with equation V2=
AB·T0 1/2/(1+B·T0 1/2)-V1 SolutionFitting obtains, n be observation time point, n=1,2,3,4 ... ...;Corresponding cumulative observations
Time is tn;K is to calculate time point, k=1,2,3,4 ... ... n, k≤n;
V1 nAnd AnFor the gas loss amount and upper limit desorption quantity of maximum duration observation point, i.e. tnValue when being maximized.
It further, further include different desorption time point tnWhen accumulative desorption of mash gas amount relative error analysis;
Wherein, the different desorption time point tnWhen accumulative desorption of mash gas amount relative error analysis equation are as follows:
RE=100% × (V2 n’-V2 n)/V2 n;
In formula, V2 nFor different number desorption time point tnThe accumulative desorption of mash gas amount of corresponding authentic testing, V2 n’For the time
Point tnIt is corresponding to be fitted accumulative desorption of mash gas amount, V2 n’By site desorption equation V2=ABT0 1/2/(1+B·T0 1/2)-V1 SolutionFitting
It finds out, V2 n’That is V2, n=1,2,3,4 ....
Further, the error analysis equation of the remaining desorption quantity, the error analysis equation of gas total content, difference mention
Bore loss amount relative error analysis side caused by the gas loss amount error analysis equation of time, different number desorption time point
Upper limit desorption quantity relative error analysis equation caused by journey, different number desorption time point and different desorption time point tnWhen it is tired
Meter desorption of mash gas amount relative error analysis equation is used equally for the live desorption of mash gas data and straight-line method, logarithm of different modes
The error analysis of the gas loss amount calculation method such as method, polynomial method, power function method, Amoco method.
The invention has the benefit that the present invention proposes that site desorption equation combines specific live desorption of mash gas process can
It is calculated with gas total content, gas loss amount and the remaining desorption quantity to deep fractures;It is opposite accidentally to reduce gas bearing capacity
Difference;There is provided corresponding error analysis method simultaneously to assist site desorption equation to carry out more accurate gas bearing capacity and calculate, from
And the accuracy and reliability of deep fractures gas bearing capacity calculating is improved, be further conducive to deep gas bearing capacity and preservation
Matter condition makes objective appraisal.
Detailed description of the invention
Fig. 1 is the corresponding desorption data of Examples 1 and 2;
Fig. 2 is using traditionThe Fitting Calculation result of method;
Fig. 3 is using equation V2=AB T0 1/2/(1+B T0 1/2)-V1 SolutionThe Fitting Calculation result;
Fig. 4 is using equation V2=AB T0 1/2/(1+B T0 1/2)-V1 SolutionIt is fitted and uses the calculated result of single point correction;
Fig. 5 is using equation V2=AB T0 1/2/(1+B T0 1/2)-V1 SolutionIt is fitted and uses the calculated result of two-point calibration;
Fig. 6 is the corresponding desorption data of embodiment 3~6;
Fig. 7 is that simulation difference proposes the relative error for boring the gas loss amount that the time obtains in embodiment 4;
Fig. 8 is the gas loss amount V of different desorption time observation point n in embodiment 51 kWith upper limit desorption quantity AkIt is opposite accidentally
Difference;
Fig. 9 is different desorption time point t in embodiment 6nWhen accumulative desorption of mash gas amount relative error;
Figure 10 be in embodiment 6 according to the calculation method that uses herein and other frequently with method acquire, simulation not
It is compared with the relative error for boring the gas loss amount that the time obtains is proposed.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art
Every other embodiment obtained without creative efforts, shall fall within the protection scope of the present invention.
Embodiment 1
The present embodiment combination gas site desorption experimental data is to equation V2=AB T0 1/2/(1+B T0 1/2)-V1 SolutionIt is said
Bright and explanation.
Gas site desorption experimental data, general desorption time are 120min.The acquirement of time interval point conforms generally to GB/
T 23249-2009 " survey period coal seam gas-bearing capacity measuring method ".Its method particularly includes: the continuous observation time, 120min stopped.
Reading time interval time rule is as follows: before starting to observe in 60min, first point to 1min-2min between thirdly, later every
2min-5min reading is primary;It is primary at interval of being read in 10min-20min in second hour.In addition, AQ 1046-2007 "
Survey period coal seam gas-bearing capacity measuring method " the reading time interval time is provided as follows: the continuous observation time, 120min stopped;Start
Before observing in 60min, the first point interval 2min is later primary every 3-5min reading;At interval of 10- in second hour
Reading is primary in 20min.
The volume data of desorption of mash gas amount used in the present embodiment and aftermentioned embodiment is corrected standard bar
Part lower body volume data.
With V2For ordinate, withFor abscissa, the throwing of whole measuring points is painted on graph paper, using equation V2=
AB·T0 1/2/(1+B·T0 1/2)-V1 SolutionIt is corrected, and reads V from Fitting curve equation1 SolutionValue, corresponding experimental data are detailed
See that Fig. 1, corresponding fitting result are shown in Fig. 3;In addition, being used using identical experiment dataThe corresponding fitting result of method is shown in Fig. 2.
It is gas loss amount site desorption data in Fig. 1, measured value and match value unit are mL in Fig. 1.
In conjunction with Fig. 2 and Fig. 3, using traditionalMethod and equation V2=ABT0 1/2/(1+B·T0 1/2)-V1 SolutionApproximating method
Obtained gas loss amount is respectively 516mL and 1213mL, traditionalMethod will lead to Lower result, and equation V2=AB T0 1 /2/(1+B T0 1/2)-V1 SolutionFitting is distorted result.
Embodiment 2
For in embodiment use equation V2=AB T0 1/2/(1+B T0 1/2)-V1 SolutionCarry out the correlometer of gasresolution amount
It calculates, although this method can accurately be fitted live gas bearing capacity, according to country and professional standard, when ground is surveyed
Scene desorption of mash gas time phase is generally 120 minutes.When using this method using 120 minutes desorption of mash gas data on site, meeting
Lead to the distortion of matched curve and the distortion of calculated result, calculating gas loss amount is made not conform to the actual conditions.Therefore, this method needs
Improvement is advanced optimized, the computational accuracy of gas bearing capacity is improved, and is convenient to field application.
In conjunction with the embodiments in 1 Fig. 1 experimental data, with V2For ordinate, withFor abscissa, whole measuring points are thrown
It is painted on graph paper.Laboratory is crushed into preceding degassing and the remaining desorption quantity V for the acquisition that deaerates after crushing3+V4It is considered as accumulation desorption
A part of amount, which is thrown, to be painted on graph paper, using equation V2=AB T0 1/2/(1+B T0 1/2)-V1 SolutionIt is corrected, and quasi- from curve
It closes and reads V on equation1 SolutionValue, corresponding fitting result are shown in Fig. 4.The gas loss amount for using single point correction approximating method to obtain for
794mL。
In conjunction with the embodiments in 1 Fig. 1 experimental data, with V2For ordinate, withFor abscissa, whole measuring points are thrown
It is painted on graph paper.Laboratory is crushed into preceding degassing and the remaining desorption quantity V for the acquisition that deaerates after crushing3And V3+V4It is respectively seen as
A part of accumulation desorption quantity, which is thrown, to be painted on graph paper, using equation V2=AB T0 1/2/(1+B T0 1/2)-V1 SolutionIt is corrected, and
V is read from Fitting curve equation1 SolutionValue, corresponding fitting result are shown in Fig. 5.The gas obtained using two-point calibration approximating method
Loss amount is 790mL.
In conjunction with the embodiments 1 and embodiment 2, knot is calculated using the gasresolution amount that single point correction and two-point calibration method obtain
Fruit is more bonded the practical desorption process of gas, eliminates using equation V2=AB T0 1/2/(1+B T0 1/2)-V1 SolutionIt is corrected fitting
The distortion of curve and the distortion of calculated result, while live 120 minutes desorption of mash gas data are also more suitable for, it is more advantageous to gas
Desorption quantity calculates.
Embodiment 3
For convenient for relative error analysis, the present embodiment uses the more another data sample of data point as the present embodiment
The relative error of remaining desorption quantity solves and the relative error of upper limit desorption quantity solves.Detailed data is shown in Fig. 6, corresponding fitting
As a result it is
V2=2524.1 × 0.04 × T0 1/2/(1+0.04×T0 1/2)-315.5
The relative error of remaining desorption quantity is solved: RE=100% × (V3 Solution-V3 It is real)/V3 It is real;
For gas site desorption data, A=2524.1mL, V2=489.7mL (it is read by practical desorption data, it is practical
It desorbs observation time and adds up 2280min, it is assumed that desorption observation time to 120min terminates, then this V2It is 120min for the corresponding time
When accumulative live desorption of mash gas amount), V1 Solution=315.5mL (is read) by fit equation;
Then V3 Solution=A-V1 Solution-V2=2524.1-315.5-489.7=1718.9mL;
After the 120min time total desorption of mash gas amount and laboratory degassing the sum of desorption quantity be considered as practical desorption quantity V3 Depending on real, then have:
V3 Depending on real=VAfter 120min Solution+V3 It is real=980.45+1072.2=2052.6mL;
Then RE=100% × (V3 Solution-V3 Depending on real)/V3 Depending on real=100% × (1718.9-2052.6)/2052.6=16.3%.
The relative error of upper limit desorption quantity is solved: R=100% × (A-VAlways Solution)/VAlways Solution, wherein VAlways Solution=V1 Solution+V2+V3 It is real,
V1 SolutionBy equation V2=AB T0 1/2/(1+B T0 1/2)-V1 SolutionIt reads, V2It is obtained by the experiment of solution cucurbitula site desorption, laboratory degassing solution
Suction obtains V3 It is real;
For gas site desorption experimental data, have:
V3 Depending on real=VAfter 120min Solution+V3 It is real=980.45+1072.15=2052.6mL
VAlways Solution=V1 Solution+V2+V3 Depending on real=315.5+489.7+2052.6=2857.8mL
R=100% × (A-VAlways Solution)/VAlways Solution=100% × (2524.1-2857.8)/2857.8=-11.68%
Embodiment 4
For convenient for relative error analysis, the present embodiment uses the more another data sample of data point as the present embodiment
Desorption experiment data, to equation V2=AB T0 1/2/(1+B T0 1/2)-V1 SolutionIt carries out gas amount error analysis and verifies this method, in detail
It counts accurately according to seeing Fig. 6;Simulation difference proposes the gas loss amount error analysis method for boring the time comprising following steps:
1) identical as graphical method and analytic method, it is specified that under operating condition drilling be mud medium or water medium, drill coal
It is the initial desorption time of the coal heart when heart is mentioned to aperture half depth.
2) t is set1' be simulation the coal heart mention be drilled into reach the earth's surface time (i.e. drilling in simulate exposure total time, t1’≥t1), T0
For in t1' under the conditions of the coal heart in solution cucurbitula time for having desorbed, then have t1' it is really to mention boring time t1It is being solved with the coal heart
The sum of 2 times of the time desorbed in cucurbitula, it may be assumed that
t1'=t1+2×T0
3) then the hypothesis of simulation different exposed total times in drilling is analyzed, the calculating loss of coal core sample board design
Measure time t0' it is equal to the drill interior half for simulating exposure total time and tinning time t2(the coal heart reach earth's surface to be packed into desorption
The sum of tank):
t0'=1/2t1’+t2
4) then total desorption time T0' it is simulation loss amount time t0' with solution cucurbitula in coal core sample desorption observation time
The sum of t, are as follows:
T0'=t0’+t
5) time is bored as t when borehole simulation mentions1' when, coal-bed gas simulates loss amount V1’(t1') are as follows:
V1’(t1')=V0(t1’)+V0’(t1’)
In formula: V0(t1') it is that simulation borehole simulation mentions brill time point t1' when, can in the hope of exposure duration t1Coal seam watt
This loss amount, V0’(t1') be simulation mention bore the time be t1' when reality coal seams gas discharging quantity
V1’(t1') be simulation mention bore the time be t1' when gas loss amount, by site desorption equation V2=ABT0 1/2/(1
+B·T0 1/2)-V1 Solution, mentioned according to simulation and bore time t1' solve.
6) time t is bored when simulation mentions1' it is really to mention boring time t1When, V0(t1')=V1’(t1');
7) time t is bored when simulation mentions1’>t1When, V0(t1')=V1’(t1’)-V0’(t1');
8) different simulations is mentioned and bores time t1', equation V can be used2=AB T0 1/2/(1+B T0 1/2)-V1 SolutionIt acquires
Difference simulation, which mentions, bores time point t1' when gas loss amount V1’(t1'), and then acquire exposure duration t1Coal-bed gas loss amount
V0(t1');
9) it defines the smallest simulation and mentions brill time t1' it (is at this time really to mention boring time t1) coal-bed gas loss amount V0
(t1) it is true value, other of board design mention the gas loss amount V for boring the time0(t1') it is observation.It then can calculating simulation drilling
True exposure duration t1When gas loss amount relative error and mapping analysis, relative error (RE) calculation formula are as follows:
RE=100% × [V0(t1’)-V0(t1)]/V0(t1);
10) Fig. 7 is using equation V2=AB T0 1/2/(1+B T0 1/2)-V1 SolutionThe different simulations acquired, which mention, bores what the time obtained
Gas loss amount and relative error;Fig. 3 is using equation V2=AB T0 1/2/(1+B T0 1/2)-V1 SolutionThe simulation difference acquired mentions brill
The relative error for the gas loss amount that time obtains, Fig. 3 content are detailed in attached drawing.
Simulation, which mentions, bores time t1' it is a dynamic value, it can change;As its value according to circumstances can be set in table 2
22min (is actually mentioned and is bored time t1)、30min、62min、112min、162min、222min。
Embodiment 5
The remaining solution that the present embodiment is obtained for degassing after deaerating and crush before in embodiment 2, being crushed using laboratory
Pipette V3+V4To solution cucurbitula desorption of mash gas amount V2It is modified, due to the uncertainty of gas site desorption length of time, uses
Equation V2=AB T0 1/2/(1+B T0 1/2)-V1 SolutionIt, can be due to the difference of desorption point quantity when being fitted solution to desorption data
Keep fitting acquired results inconsistent, there are errors to form the loss amount that calculating obtains.For convenient for relative error analysis, this reality
It applies example and is directed in embodiment 2 as embodiment offer because of different number solution using the data sample in the more Fig. 6 of data point
It is as follows to inhale loss amount relative error analysis method caused by time point:
(1) assume there be n observation time point, then adding up observation time is tn, n=1,2,3,4 ... ... n.
(2) the total desorption time T of coal sample0For the calculating exposure duration t before tinning0With desorption observation time t after tinningnIt
With i.e. T0=t0+tn。
(3) using mathematical program software (or Excel fitting, programming etc.) to observation time point tnAccumulative gas solution
Pipette adds up desorption quantity V2 nUsing equation V2=AB T0 1/2/(1+B T0 1/2)-V1 SolutionFitting, acquires different desorption time point tnIt is right
The gas loss amount V answered1 nWith upper limit desorption quantity An。
(4) the gas loss amount and upper limit desorption quantity for setting maximum duration observation point respectively are respectively V1 nAnd AnFor true value,
The gas loss amount V of different desorption time observation point n1 kWith upper limit desorption quantity Ak(k=1,2,3,4 ... ... relative error n)
(RE) it can be found out respectively by following formula:
REV=100% × (V1 k-V1 n)/V1 n
REA=100% × (Ak-An)/An
(5) Fig. 8 is for desorption of mash gas experimental data, to different number desorption time point tnAccumulation desorption of mash gas amount V2 n
Using Formula V2=AB T0 1/2/(1+B T0 1/2)-V1 SolutionFitting, the gas loss amount V of the different desorption time observation point n acquired1 kWith
Upper limit desorption quantity Ak(k=1,2,3,4 ... ... relative error (RE) n).
Embodiment 6
The present embodiment provides the remaining solutions obtained for degassing after deaerating and crush before being crushed in embodiment 1 by laboratory
Pipette V3+V4To solution cucurbitula desorption of mash gas amount V2It is modified, V2=AB T0 1/2/(1+B T0 1/2)-V1 SolutionTo true gas solution
The subduplicate changing rule of pipette at any time is fitted, and can be obtained the desorption of mash gas amount of fitting and be carried out error analysis.
For convenient for relative error analysis, the present embodiment uses the data sample in the more Fig. 6 of data point as analysis object.
For different number desorption time point tnAuthentic testing accumulative desorption of mash gas amount V2 nIf it corresponds to time point tn
Fitting add up desorption of mash gas amount be V2 n’, then different desorption time point t can be acquirednWhen accumulative desorption of mash gas amount relative error.
Expression formula are as follows:
RE=100% × (V2 n’-V2 n)/V2 n
For the fine desorption experiment data of gas, in each different desorption time point tnWhen accumulative desorption of mash gas amount it is opposite
Error is as shown in Figure 9.
It should be noted that the gas bearing capacity calculation method being mentioned herein, can not only be applied in coal geology exploratory bore-hole
The measurement of coal heart gas bearing capacity can also be applied to the coal bed gas content measurement of coal core sample in coal and cbm exploration.
It should be noted that the correlated error analysis method being mentioned herein can not only be to gas loss meter proposed in this paper
Calculation method carries out error analysis, can also be applied in coal and cbm exploration the desorption of mash gas data of coal core sample and
Method, logarithm method, polynomial method, Amoco and other common gas loss amount calculation method progress error analyses both at home and abroad.
As Figure 10 be according to the calculation method used herein with other frequently with method acquire, simulation is different mention brill when
Between obtain gas loss amount relative error comparison.
Particularly, desorption of mash gas point due to the time it is long, put it is relatively more, respective intervals time point therein can be extended,
It shortens, lack or deletes, but is little to influential effect is calculated.
In conclusion in conjunction with the embodiments 1 and embodiment 2, the present invention proposes site desorption equation V2=ABT0 1/2/(1+
B·T0 1/2)-V1 Solution, can gas total content, gas loss amount and remnants to deep fractures in conjunction with specific live resolving
Desorption quantity is calculated;The calculated result of 3~6 pairs of Examples 1 and 2 carries out error analysis in conjunction with the embodiments, passes through error analysis
Data the result shows that, site desorption equation V proposed by the present invention2=ABT0 1/2/(1+B·T0 1/2)-V1 SolutionCalculated result
Error it is smaller, and as the increase error of desorption time point quantity and observation time point can be smaller and smaller and tends towards stability, have
Accurate calculating, especially deep fractures gas are made conducive to coal seam gas-bearing capacity, is further conducive to deep fractures watt
This content and hosting geological condition make objective appraisal.
It should be noted that, in this document, such as first and second or the like relational terms are used merely to one if it exists
A entity or operation with another entity or operate distinguish, without necessarily requiring or implying these entities or operation it
Between there are any actual relationship or orders.Moreover, the terms "include", "comprise" or its any other variant are intended to
Cover non-exclusive inclusion, so that the process, method, article or equipment for including a series of elements not only includes those
Element, but also including other elements that are not explicitly listed, or further include for this process, method, article or setting
Standby intrinsic element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that
There is also other identical elements in the process, method, article or apparatus that includes the element.
The above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations;Although with reference to the foregoing embodiments
Invention is explained in detail, those skilled in the art should understand that: it still can be to aforementioned each implementation
Technical solution documented by example is modified or equivalent replacement of some of the technical features;And these modification or
Replacement, the spirit and scope for technical solution of various embodiments of the present invention that it does not separate the essence of the corresponding technical solution.
Claims (10)
1. a kind of Coal Bed Methane Content Estimation Method, which comprises the following steps:
S1 determines site desorption equation:
V2=ABT0 1/2/(1+B·T0 1/2)-V1 Solution
Wherein, A is the upper limit desorption of mash gas total amount of coal heart desorption, and B is the coefficient in relation to gas adsorption speed and heat of adsorption, V1 SolutionFor
Gas loss amount, V2Add up desorption quantity, T for gas0=t0+ t, t0=1/2t1+t2,
V3 Solution=A-V1 Solution-V2, V3 SolutionFor remaining desorption quantity calculated value;
The desorption of the S2 coal heart: drilling as mud medium or water medium, when the coal heart that drills is mentioned to aperture half depth, as the coal heart
Initial desorption time;It will mention and bore coal heart progress site desorption, during which, record, which mentions, bores time t1, tinning time t2, write down simultaneously
Desorption time t and the accumulative desorption quantity V of t at any time in solution cucurbitula2;And remaining desorption quantity V is recorded by experimental determination3 It is real;
S3 equation model: by the accumulative desorption quantity V of coal core sample at any time2Throwing is drawn in a coordinate system, carries out scene to data point
The fitting of Desorption equation, determines A and V1 SolutionValue, then A be coal-bed gas total content calculated value, V1 SolutionFor coal-bed gas loss amount
Calculated value.
2. a kind of Coal Bed Methane Content Estimation Method according to claim 1, which is characterized in that remaining desorption quantity V2 It is realMake
It is painted on graph paper to accumulate a part and the throwing of desorption quantity, fitting is corrected to equation.
3. a kind of Coal Bed Methane Content Estimation Method according to claim 2, which is characterized in that the site desorption equation
Calculating suitable for the desorption of mash gas data that coal seam gas-bearing capacity measuring method obtains.
4. a kind of Coal Bed Methane Content Estimation Method according to claim 3, which is characterized in that use single point correction or double
Point correction carries out Data correction;
The method of the single point correction are as follows:
(1) coordinate is establishedThrowing is drawn on graph paper, wherein V3To crush initial vacuum heating
Amount of gas evolved, V4For amount of gas evolved after crushing, V3+V4It is considered as in solution cucurbitula what uninterrupted desorption at any time continuously obtained, corresponding solution
The analysis time is ti;
(2) formula V is used2=ABT0 1/2/(1+B·T0 1/2)-V1 SolutionTo coordinate points It carries out anti-
Multiple fitting optimization, makes coordinate pointsMeet formula V as accurately as possible2=ABT0 1/2/(1+
B·T0 1/2)-V1 SolutionThe extension trend of curve;
(3) the desorption time coordinate that finally obtains optimal fitting result obtains the A of fit equation, B parameter value and under this condition
V1 Solution, V1 SolutionThe then gas loss amount to solve;
The method of the two-point calibration are as follows:
(1) coordinate is establishedWithThrowing is drawn on graph paper,
Wherein V3To crush initial vacuum heat de-airing amount, V4For amount of gas evolved after crushing, coordinateCorresponding V3Depending on
For the component part for adding up desorption quantity, the corresponding parsing time is ti-1, coordinateCorresponding V3+
V4It is considered as the composition portion of accumulative desorption quantity, the corresponding parsing time is ti;
(2) formula V is used2=ABT0 1/2/(1+B·T0 1/2)-V1 SolutionTo coordinate pointsWithFitting optimization repeatedly is carried out, coordinate points are madeWith
Meet formula V as accurately as possible2=ABT0 1/2/(1+B·T0 1/2)-V1 SolutionThe extension trend of curve;
(3) the desorption time coordinate that finally obtains optimal fitting result obtains the A of fit equation, B parameter value and under this condition
V1 Solution, V1 SolutionThe then gas loss amount to solve.
5. a kind of Coal Bed Methane Content Estimation Method according to claim 4, which is characterized in that desorption of mash gas can be calculated
To V2With V3+V4When corresponding theoretical observation interval be tΔ=ti-t.At this point, t is cumulative maximum desorption quantity V2It is corresponding tired
Meter desorption observation time.
6. a kind of error analysis method of any Coal Bed Methane Content Estimation Method, special according to claim 1~5
Sign is, the error analysis of error analysis and gas total content including remaining desorption quantity;
Wherein, the error analysis equation of the remaining desorption quantity are as follows:
RE=100% × (V3 Solution-V3 It is real)/V3 It is real;
The error analysis equation of the gas total content are as follows:
R=100% × (A-VAlways Solution)/VAlways Solution, wherein VAlways Solution=V1 Solution+V2+V3 It is real。
7. a kind of error analysis method of Coal Bed Methane Content Estimation Method according to claim 6, which is characterized in that also
The gas loss amount error analysis for boring the time is proposed including difference;
Wherein, the difference mentions the gas loss amount error analysis equation for boring the time are as follows:
RE=100% × [V0(t1’)-V0(t1)]/V0(t1);
V0(t1) it is that the smallest simulation mentions brill time t1' coal-bed gas loss amount, i.e. actual field desorbs true exposure duration t1
Under the conditions of according to site desorption equation V2=ABT0 1/2/(1+B·T0 1/2)-V1 SolutionThe V of acquisition1 SolutionValue;
V0(t1') it is that different simulations mention brill time point t1' the exposure duration t that acquires1Coal-bed gas loss amount;
Wherein, V0(t1')=V1’(t1’)-V0’(t1');
In formula, V0’(t1') be simulation mention bore the time be t1' when reality coal seams gas discharging quantity, read by desorption data;
V1’(t1') be simulation mention bore the time be t1' when gas loss amount, by site desorption equation V2=ABT0 1/2/(1+B·
T0 1/2)-V1 Solution, mentioned according to simulation and bore time t1' solve;
T0'=t0'+t, T0' it is total desorption time;
t0'=1/2t1’+t2 t0' be coal core sample board design the calculating loss amount time, t2For the tinning time;
t1'=t1+2×T0, t1' be simulation the coal heart mention be drilled into reach the earth's surface time, T0For in t1' under the conditions of the coal heart solution cucurbitula
In time for having desorbed.
8. a kind of error analysis method of Coal Bed Methane Content Estimation Method according to claim 7, which is characterized in that also
Including loss amount caused by different number desorption time point and upper limit desorption quantity relative error analysis;
Wherein, loss amount relative error analysis equation caused by the different number desorption time point are as follows:
REV=100% × (V1 k-V1 n)/V1 n;
Wherein, upper limit desorption quantity relative error analysis equation caused by the different number desorption time point are as follows:
REA=100% × (Ak-An)/An;
In formula, V1 kAnd AkRespectively correspond desorption time point tnGas loss amount and upper limit desorption quantity;All in accordance with equation V2=AB
T0 1/2/(1+B·T0 1/2)-V1 SolutionFitting obtains, n be observation time point, n=1,2,3,4 ... ...;The corresponding cumulative observations time
For tn;K is to calculate time point, k=1,2,3,4 ... ... n, k≤n;
V1 nAnd AnFor the gas loss amount and upper limit desorption quantity of maximum duration observation point, i.e. tnValue when being maximized.
9. a kind of error analysis method of Coal Bed Methane Content Estimation Method according to claim 8, which is characterized in that also
Including different desorption time point tnWhen accumulative desorption of mash gas amount relative error analysis;
Wherein, the different desorption time point tnWhen accumulative desorption of mash gas amount relative error analysis equation are as follows:
RE=100% × (V2 n’-V2 n)/V2 n;
In formula, V2 nFor different number desorption time point tnThe accumulative desorption of mash gas amount of corresponding authentic testing, V2 n' it is time point tn
It is corresponding to be fitted accumulative desorption of mash gas amount, V2 n' by site desorption equation V2=ABT0 1/2/(1+B·T0 1/2)-V1 SolutionFitting is asked
Out, V2 n' i.e. V2, n=1,2,3,4 ....
10. a kind of error analysis method of Coal Bed Methane Content Estimation Method according to claim 9, which is characterized in that
The error analysis equation of the remnants desorption quantity, the error analysis equation of gas total content, the different gas loss for mentioning the brill time
When measuring loss amount relative error analysis equation caused by error analysis equation, different number desorption time point, different number desorption
Between put caused by upper limit desorption quantity relative error analysis equation and different desorption times point tnWhen accumulative desorption of mash gas amount it is opposite
Error analysis equation is used equally for the gas loss amount calculation method such as straight-line method, logarithm method, polynomial method, Amoco method to be missed
Difference analysis.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910956488.8A CN110533349B (en) | 2019-10-08 | 2019-10-08 | Coal seam gas content calculation and error analysis method |
ZA2020/02560A ZA202002560B (en) | 2019-10-08 | 2020-05-08 | Coal seam gas content calculation method and error analysis method thereof |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910956488.8A CN110533349B (en) | 2019-10-08 | 2019-10-08 | Coal seam gas content calculation and error analysis method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110533349A true CN110533349A (en) | 2019-12-03 |
CN110533349B CN110533349B (en) | 2023-03-24 |
Family
ID=68671632
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910956488.8A Active CN110533349B (en) | 2019-10-08 | 2019-10-08 | Coal seam gas content calculation and error analysis method |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN110533349B (en) |
ZA (1) | ZA202002560B (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102128765A (en) * | 2010-12-09 | 2011-07-20 | 中国矿业大学 | Method for directly and rapidly measuring coal seam gas content in underground coal mine |
CN103983534A (en) * | 2014-05-26 | 2014-08-13 | 河南理工大学 | Gas loss calculation method |
RU2601205C1 (en) * | 2015-08-07 | 2016-10-27 | Акционерное общество "Метан Кузбасса" | Method of determining gas content of coal beds |
CN106680451A (en) * | 2015-11-09 | 2017-05-17 | 河南理工大学 | Underground rapid measurement method for coal and gas outburst parameter as well as apparatus thereof |
CN108458951A (en) * | 2018-04-02 | 2018-08-28 | 河南理工大学 | The method that coal seam gas-bearing capacity is quickly determined by desorption efficiency |
CN109406750A (en) * | 2018-12-06 | 2019-03-01 | 西安科技大学 | A kind of coal seam gas-bearing capacity directly measures loss amount calculation method |
-
2019
- 2019-10-08 CN CN201910956488.8A patent/CN110533349B/en active Active
-
2020
- 2020-05-08 ZA ZA2020/02560A patent/ZA202002560B/en unknown
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102128765A (en) * | 2010-12-09 | 2011-07-20 | 中国矿业大学 | Method for directly and rapidly measuring coal seam gas content in underground coal mine |
CN103983534A (en) * | 2014-05-26 | 2014-08-13 | 河南理工大学 | Gas loss calculation method |
RU2601205C1 (en) * | 2015-08-07 | 2016-10-27 | Акционерное общество "Метан Кузбасса" | Method of determining gas content of coal beds |
CN106680451A (en) * | 2015-11-09 | 2017-05-17 | 河南理工大学 | Underground rapid measurement method for coal and gas outburst parameter as well as apparatus thereof |
CN108458951A (en) * | 2018-04-02 | 2018-08-28 | 河南理工大学 | The method that coal seam gas-bearing capacity is quickly determined by desorption efficiency |
CN109406750A (en) * | 2018-12-06 | 2019-03-01 | 西安科技大学 | A kind of coal seam gas-bearing capacity directly measures loss amount calculation method |
Also Published As
Publication number | Publication date |
---|---|
CN110533349B (en) | 2023-03-24 |
ZA202002560B (en) | 2021-04-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102128765B (en) | Method for directly and rapidly measuring coal seam gas content in underground coal mine | |
CN106442938B (en) | A kind of accurate survey calculation method for obtaining shale air content | |
US11435336B2 (en) | Method for determining content of lost gas in shale gas content test | |
CN108982817A (en) | A kind of shale air content appraisal procedure based on Methane Carbon Isotope | |
CN105893679B (en) | Well Test Data Analysis Method is corrected in low yield horizontal well afterflow | |
WO2011133885A1 (en) | Total storage capacity and total porosity of porous media | |
CN104379870A (en) | Method for determining a permeability or mobility of a radial flow response of a reservoir | |
CN107703023A (en) | A kind of surface drilling coal core modeled pressure drop desorption apparatus and test device and method | |
CN108169062A (en) | Simulate the visual test device and method of subterranean coal gas preservation desorption process | |
CN107832574B (en) | Horizontal well water flooded layer interpretation method based on logging while drilling | |
CN105628549A (en) | Method for directly, quickly and full-automatically determining content of desorbing gas in coal seam | |
CN110633451B (en) | Shale gas loss calculation method and system | |
CN108458951A (en) | The method that coal seam gas-bearing capacity is quickly determined by desorption efficiency | |
CN105484735A (en) | Method for evaluating coincidence rate of actual drilling borehole trajectory and design track | |
CN102323181B (en) | Method and apparatus for detecting coalbed gas content based on three-component seismic data | |
CN110173259B (en) | Stoneley wave energy compensation correction method | |
CN106019403B (en) | From being born from reservoir formation porosity measurement method | |
CN110533349A (en) | A kind of coal seam gas-bearing capacity calculates and error analysis method | |
CA2649483A1 (en) | Refined analytical model for formation parameter calculation | |
Nie et al. | Evaluation of the in-place adsorbed gas content of organic-rich shales using wireline logging data: a new method and its application | |
CN109538199A (en) | A kind of coal measure strata air content evaluation method, device and electronic equipment | |
CN115615869A (en) | Loss gas amount testing device and method for simulating rope coring process | |
CN207408265U (en) | A kind of surface drilling coal core modeled pressure drop desorption apparatus and test device | |
CN113775329A (en) | Method and system for acquiring deep carbonate rock mechanical parameters | |
GB2458547A (en) | Earth formation testing by regression analysis of induced flow pressure measurements |
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 |