CN106290058B - A kind of measuring method and device of rock gassiness data - Google Patents
A kind of measuring method and device of rock gassiness data Download PDFInfo
- Publication number
- CN106290058B CN106290058B CN201610725811.7A CN201610725811A CN106290058B CN 106290058 B CN106290058 B CN 106290058B CN 201610725811 A CN201610725811 A CN 201610725811A CN 106290058 B CN106290058 B CN 106290058B
- Authority
- CN
- China
- Prior art keywords
- rock
- measured
- desorption
- volume
- gas
- 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
Links
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/14—Analysing materials by measuring the pressure or volume of a gas or vapour by allowing the material to emit a gas or vapour, e.g. water vapour, and measuring a pressure or volume difference
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
The embodiment of the present application discloses the measuring method and device of a kind of rock gassiness data.The described method includes: desorbing naturally the described method includes: rock to be measured is fitted into desorption apparatus, desorption rate is monitored；Fitting obtains the functional relation of desorption rate and desorption time, according to functional relation, calculates the volume of loss gas；Naturally desorption measures the volume of stripping gas obtained by nature desorption and fast desorption up to monitoring that gas contained by rock to be measured is desorbed completely；The volume for losing gas is added to obtain total volume with the volume of the stripping gas, and then the air content of rock to be measured is calculated.Using each embodiment of the present invention, it can effectively simplify continuous mode, improve the accuracy of measurement result.
Description
Technical field
This application involves oil-gas exploration and development technical field, in particular to a kind of the measuring method and dress of rock gassiness data
It sets.
Background technique
With China's rapid economic development, petroleum resources supply anxiety is increasingly enhanced.Shale gas, coal bed gas and shale oil etc.
Unconventional oil and gas resource exploration, which is developed, to be paid attention to.Rock air content is to evaluate the pass of the tastes of resources such as shale gas, coal bed gas
Bond parameter, how Accurate Determining rock air content is particularly important.
In the prior art, most widely used in the measuring method of rock air content is United States Bureau of Mines's direct method, the party
Gas contained by rock to be measured is divided into three parts by method: the quantity of stripping gas, loss gas and residual gas, three parts gas content is closed
It is that can be indicated with following formula,
In formula, G indicates that rock air content, unit are every gram of cubic centimetre；
GlIndicate loss Gas content, unit is every gram of cubic centimetre；
GdIndicate desorption Gas content, unit is every gram of cubic centimetre；
GrIndicate residual gas content, unit is every gram of cubic centimetre；
VlIndicate loss air volume, unit is cubic centimetre；
VdIndicate stripping gas volume, unit is cubic centimetre；
VrIndicate residual gas volume, unit is cubic centimetre；
M indicates that sample gross mass, unit are gram；
mrIndicate residual gas sample quality, unit is gram.
The volume of three parts gas is determined respectively, so that it may calculate the air content of rock to be measured.It is described in this method
The definition of stripping gas is that rock to be measured is fitted into desorption apparatus after sealing, to be measured under certain temperature, pressure, time conditions
The gas desorbed in rock, volume, which can be determined directly, to be come, but this process needs to carry out one week or more, efficiency
Lowly；The definition of the loss gas is that rock to be measured encounters rock to be measured to rock to be measured from brill and is loaded into this section of desorption apparatus
The volume of the gas that time is desorbed, this portion gas can not be measured directly, be needed through the stripping gas volume at any time
The rule of variation calculates, but since the pressure in atmospheric environment is different from the pressure in stratum, adsorption law is not
Together, the total volume for losing gas in atmosphere and stratum is directly extrapolated by the rule that stripping gas volume changes over time, and can be produced
Raw very big error；The definition of the residual gas is, rock to be measured passes through in desorption apparatus desorb naturally after, by a part to
It surveys rock to take out and crush, remaining gas in the rock to be measured desorbed, the volume of this portion gas can be measured directly
Out, but rock to be measured has gas escape to come out, this part from desorption apparatus is opened to crushing in rock this period to be measured
The volume of gas is difficult to measure, and final measurement result is caused to generate error.
At least there are the following problems in the prior art: the continuous mode time is long, and process is cumbersome, the gas of loss in continuous mode
Body is difficult to measure, and measurement result is caused to generate error.In addition, the calculation method accuracy for losing gas in the prior art is low, loss
The calculated result and true value deviation of gas are big.
Summary of the invention
The purpose of the embodiment of the present application is to provide the measuring method and device of a kind of rock gassiness data, is measured with simplification
Journey shortens minute, improves the accuracy of measurement result.
In order to solve the above technical problems, a kind of measuring method of rock gassiness data provided by the embodiments of the present application is in this way
It realizes:
A kind of measuring method of rock gassiness data, which comprises
Rock to be measured is fitted into progress nature desorption in desorption apparatus, monitor desorption rate and is obtained obtained by nature desorption
Stripping gas volume；
Fitting obtains the functional relation of the desorption rate desorbed naturally and desorption time, according to the functional relation,
The volume in first time period with second time period internal loss gas is calculated separately, total losses air volume is obtained；The first time
Section encounters the rock to be measured to rock to be measured for brill and reaches ground this period, and the second time period is the rock to be measured
Ground is reached to loading desorption apparatus this period；
Naturally after the time desorbed reaches predetermined period, the rock to be measured is carried out in the desorption apparatus fast
Speed desorption, and monitor the fast desorption；Until monitoring that gas contained by rock to be measured is desorbed completely, the quick solution is measured
Inhale the volume of obtained stripping gas；
By desorption obtained by the volume of the loss gas, the volume of the gained stripping gas of desorption naturally, the fast desorption
The volume of gas is added, and obtains the total volume of contained gas in the rock to be measured.
In preferred embodiment, the fitting includes fitting of a polynomial, and the fitting obtains the desorption speed desorbed naturally
The functional relation of rate and desorption time, comprising:
Using the data of the gained desorption rate of desorption naturally, fitting coefficient is calculated, the letter of fitting of a polynomial is obtained
Number relational expression.
In preferred embodiment, the functional relation expression formula of the fitting of a polynomial includes:
Q (t)=antn+an-1tn-1+…+a2t2+a1t1+a0
In formula, n indicates integer；
an, an-1... ... a1, a0Indicate fitting coefficient；
T indicates desorption time, and unit is minute；
Q (t) indicates the desorption rate when desorption t time, and unit is that milliliter is per minute.
In preferred embodiment, the volume for calculating the loss gas in the second time period includes:
Using the functional relation for being fitted and obtaining, it is calculated by the following formula out the second time period internal loss gas
Volume Vl2:
Vl2=qm(tm-tm-1)+qm-1(tm-1-tm-2)+…+q2(t2-t1)+q1t1
In formula, m indicates integer；
tmIndicate that rock to be measured reaches the time experienced behind ground, unit is minute；
qmIt indicates after rock to be measured reaches ground through going through tmDesorption rate after time, unit are that milliliter is per second.
In preferred embodiment, the volume V of the second time period internal loss gasl1Evaluation method include: using described quasi-
The functional relation obtained is closed, the desorption rate q that rock to be measured reaches the ground moment is calculated0, using desorption rate with it is to be measured
The proportional relation of rock inside and outside differential pressure obtains during drilling rock to be measured in the desorption rate q of stratum different depth0i:
pfi=GDp(h-hf)=GDp(h-vt0i)
In formula, GDpIndicate formation pressure gradient, unit is every meter of megapascal；
H indicates prime stratum depth locating for rock to be measured, and unit is rice；
V indicates that the Mean Speed promoted in rock lifting process to be measured, unit are metre per second (m/s)s；
t0iIndicate that rock to be measured was promoted to a certain depth elapsed time, unit is minute；
pfiIndicate strata pressure when rock to be measured is promoted to a certain depth, unit is megapascal；
pfIt indicates that original formation pressure locating for rock to be measured, unit are megapascal, is known quantity；
And then when calculating the brill by following algorithm and encountering the rock to be measured to rock to be measured this section of ground of arrival
Between internal loss gas volume Vl1:
Vl1=q0i(t0i-t0i-1)+q0i-1(t0i-1-t0i-2)+…+q02(t02-t01)+q01t01
In formula, i indicates integer；
t0iIndicate that rock to be measured was promoted to a certain depth elapsed time, unit is minute；
q0iIndicate that rock to be measured is in the desorption rate of stratum different depth during drilling.
In preferred embodiment, it includes being solved using mode of oscillation to rock to be measured that the rock to be measured, which carries out fast desorption,
It inhales.
The method also includes:
The ratio for calculating the total volume of contained gas and the rock quality to be measured in the rock to be measured, obtain it is described to
Survey the air content of rock.
A kind of measurement device of rock gassiness data, described device include:
Cucurbitula is solved, for holding the rock to be measured, the desorption environment of rock to be measured is adjusted, nature is carried out to rock to be measured
Desorption and fast desorption；
With desorb tank connected desorption rate monitoring device, the accumulative body for real-time measurement desorption rate and stripping gas
Product；
The gas collector being connected with monitoring device, for collecting nature desorption and the resulting desorption of fast desorption
Gas.
In preferred embodiment, described device further include:
Calculation processing unit, the desorption rate data reconciliation desorbed naturally obtained for obtaining the monitoring device measurement
The cumulative volume data of air-breathing；
According to the desorption rate data desorbed naturally that measurement obtains, nature desorption speed is obtained using the fitting algorithm of setting
The functional relation of rate and desorption time；
Using the proportional relation of desorption rate and rock inside and outside differential pressure to be measured, rock to be measured during drilling is calculated and exists
The desorption rate q of stratum different depth0i；
Using the calculation formula of setting, the volume of the loss gas in the second time period is calculated；
Using the calculation formula of setting, the volume of the loss gas in the first time period is calculated；
It is also used to calculate the volume for losing gas in the resulting first time period, the damage in the first time period
It loses the volume of gas, the volume of the gained stripping gas of desorption naturally, the volume of fast desorption gained stripping gas to be added, obtain
The total volume of gas contained by the rock to be measured.
In preferred embodiment, the desorption is canned pressure gauge and temp-controlling element, for adjusting the environment in solution cucurbitula.
In preferred embodiment, the desorption is canned oscillation device, for carrying out fast desorption to the rock to be measured.
In preferred embodiment, the fitting algorithm of the setting includes:
Polynomial fitting；
Least square fitting algorithm.
In preferred embodiment, the volume V of the loss gas in the calculating second time period of the settingl2Calculation formula
Include:
Vl2=qm(tm-tm-1)+qm-1(tm-1-tm-2)+…+q2(t2-t1)+q1t1
In formula, m indicates integer；
tmIndicate that rock to be measured reaches the time experienced behind ground, unit is minute；
qmIt indicates after rock to be measured reaches ground through going through tmDesorption rate after time, unit are that milliliter is per second.
In preferred embodiment, the volume V of the loss gas in the calculating second time period of the settingl1Calculation formula
Include:
Vl1=q0i(t0i-t0i-1)+q0i-1(t0i-1-t0i-2)+…+q02(t02-t01)+q01t01
In formula, i indicates integer；
t0iIndicate that rock to be measured was promoted to a certain depth elapsed time, unit is minute；
q0iIndicate that rock to be measured is in the desorption rate of stratum different depth during drilling.
In preferred embodiment, described device further include:
Air content computing module, for the matter by the total volume of gas contained by the rock to be measured divided by the rock to be measured
Amount, calculates the air content of rock to be measured.
The embodiment of the present application eliminates the continuous mode of residual gas in the prior art, and nature is carried out in the desorption apparatus
Desorption and two processes of fast desorption, rock to be measured is desorbed completely in the desorption apparatus, does not need to open desorption apparatus,
It does not need that rock measurement residual gas content to be measured is taken to accelerate to simplify the continuous mode of the prior art in the desorption
The rate desorbed in device, shortens minute, avoids in the prior art because of error brought by residual gas measurement.Separately
Outside, due to considering the difference in atmospheric environment with adsorption law in stratum in the present invention, the plan point of air volume will be lost
Rock to be measured to the volume of rock to be measured arrival internal loss ground this period gas and rock to be measured, which is encountered, for brill reaches ground extremely
It is packed into the volume of solution internal loss cucurbitula this period gas, is calculated two-part according to different theories using different calculation methods
Air volume is lost, the accuracy that loss air volume calculates is improved.Entire technical solution has reached simplified continuous mode, shortens and surveys
It fixes time, improves the effect of the accuracy of measurement result.The measurement device of the rock gassiness data can carry out entirely automatically
The implementation process of measuring method improves user experience without implementing the participation of personnel.
Detailed description of the invention
In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only used for solving
Purpose is released, and is not intended to limit the scope of the present disclosure in any way, in addition, the shape and ratio of each component in figure
Size etc. be only it is schematical, be used to help the understanding of the present invention, be not the specific shape for limiting each component of the present invention and
Proportional sizes.It for those of ordinary skill in the art, without any creative labor, can also be according to this
A little attached drawings obtain other attached drawings.
Fig. 1 is the flow chart of the method for measurement rock gassiness data of the embodiment of the present invention；
Fig. 2 is the structure chart of the device of measurement rock gassiness data of the embodiment of the present invention；
Fig. 3 is the desorption rate data composition scatter plot of the different moments in one embodiment of the present of invention, monitored；
Fig. 4 is the letter of the desorption rate desorbed naturally and desorption time that are fitted in one embodiment of the present of invention
Number relation curve.
Specific embodiment
The embodiment of the present application provides the measuring method and device of a kind of rock gassiness data.
In order to make those skilled in the art better understand the technical solutions in the application, below in conjunction with the application reality
The attached drawing in example is applied, the technical scheme in the embodiment of the application is clearly and completely described, it is clear that described implementation
Example is merely a part but not all of the embodiments of the present application.Based on the embodiment in the application, this field is common
The application protection all should belong in technical staff's every other embodiment obtained without creative efforts
Range.
The embodiment of the invention discloses the measuring method and device of a kind of rock gassiness data, Fig. 1 is a kind of reality of the present invention
The method flow schematic diagram of example is applied, specifically, as shown in Figure 1, the measuring method of the rock gassiness data may include:
S01: rock to be measured is fitted into desorption apparatus, desorbs the predetermined time naturally, desorption rate is monitored, according to described pre-
The desorption rate data that interior monitoring of fixing time obtains, are obtained by polynomial fitting method, desorb desorption rate and desorption naturally
The functional relation of time, the expression formula of the functional relation are as follows:
Q (t)=antn+an-1tn-1+…+a2t2+a1t1+a0 (1)
In formula, n indicates integer；
an, an-1... ... a1, a0Indicate fitting coefficient；
T indicates desorption time, and unit is minute
Q (t) indicates the desorption rate when desorption t time, and unit is that milliliter is per minute.
S02: it using the expression formula (1) for the functional relation that fitting obtains, calculates the rock to be measured and reaches ground
To being packed into solution cucurbitula this period, air volume V is lostl2, the volume of the period internal loss gas is calculated by following algorithm
Vl2:
Vl2=qm(tm-tm-1)+qm-1(tm-1-tm-2)+…+q2(t2-t1)+q1t1 (2)
In formula, m indicates integer；
tmIndicate that rock to be measured reaches the time experienced behind ground, unit is minute；
qmIt indicates after rock to be measured reaches ground through going through tmDesorption rate after time, unit are that milliliter is per second.
The desorption rate q that rock to be measured reaches the ground moment is calculated using the expression formula (1) of the functional relation0, benefit
With the proportional relation of desorption rate and rock inside and outside differential pressure to be measured, rock to be measured is obtained during drilling in stratum different depth
Desorption rate q0i:
pfi=GDp(h-hf)=GDp(h-vt0i) (4)
In formula, GDpIndicate formation pressure gradient, unit is every meter of megapascal, is known parameters；
H indicates prime stratum depth locating for rock to be measured, and unit is rice, is known quantity；
V indicates that the Mean Speed promoted in rock lifting process to be measured, unit are metre per second (m/s)s, is known quantity；
t0iIndicate that rock to be measured was promoted to a certain depth elapsed time, unit is minute；
pfiIndicate strata pressure when rock to be measured is promoted to a certain depth, unit is megapascal；
pfIt indicates that original formation pressure locating for rock to be measured, unit are megapascal, is known quantity.
And then when calculating the brill by following algorithm and encountering the rock to be measured to rock to be measured this section of ground of arrival
Between internal loss gas volume Vl1:
Vl1=q0i(t0i-t0i-1)+q0i-1(t0i-1-t0i-2)+…+q02(t02-t01)+q01t01。(5)
By the Vl1With the Vl2It is added, obtains the volume V of loss gasl。
S03: after the time desorbed naturally reaches predetermined period, by mode of oscillation, fill the rock to be measured in desorption
Interior fast desorption is set, until monitoring that gas contained by rock to be measured is desorbed completely, is measured obtained by nature desorption and fast desorption
The volume V of stripping gasd。
S04: by the volume V of the loss gaslWith the volume V of the stripping gasdIt is added, calculates addition result and rock to be measured
The ratio of stone quality obtains the total volume of gas contained by rock.
As an embodiment of the present invention, the measurement device in conjunction with used in implementation process of the present invention as shown in Fig. 2,
The specific embodiment of the measuring method of the rock gassiness data are as follows:
Rock 1 to be measured is fitted into the desorption apparatus for filling with filler 2, the desorption apparatus includes desorption container 3, temperature control
Element 4, ultrasonic oscillator 5, pressure gauge 6 and valve 7, by temp-controlling element 4 and pressure gauge 6 adjust in desorption apparatus temperature and
Temperature and pressure is all adjusted to identical with atmosphere by pressure, and rock 1 to be measured is made to carry out nature desorption, open simultaneously valve 7,
Valve 8 and valve 10, using gas flowmeter 9, real-time monitoring desorption rate flows out the gas of valve 10, utilizes 12 He of sink
Gas collector composed by gas bottle 11 is collected.
Naturally desorption carries out 900 minutes, obtains in one embodiment of the present of invention, a certain rock to be measured is in desorption apparatus
In desorption process, the desorption rate data of the different moments monitored form scatter plot, and the scatter plot is as shown in figure 3, root
It according to the desorption rate data monitored in the predetermined time, is obtained by polynomial fitting method, the desorption desorbed naturally
The functional relation of rate and desorption time, the expression formula of the functional relation are as follows:
Q (t)=antn+an-1tn-1+…+a2t2+a1t1+a0 (6)
In formula, n indicates integer；
an, an-1... ... a1, a0Indicate fitting coefficient；
T indicates desorption time, and unit is minute
Q (t) indicates the desorption rate when desorption t time, and unit is that milliliter is per minute.
Fig. 4 is to be monitored in desorption process in desorption apparatus in one embodiment of the present of invention according to a certain rock to be measured
The desorption rate data of obtained different moments, the function relation curve of the desorption rate and desorption time that are fitted.
The rock to be measured, which is calculated, using the expression formula (6) of the functional relation for being fitted and obtaining reaches ground to loading
It solves in cucurbitula this period, loses air volume Vl2, the volume V of the period internal loss gas is calculated by following algorithml2:
Vl2=qm(tm-tm-1)+qm-1(tm-1-tm-2)+…+q2(t2-t1)+q1t1 (7)
In formula, m indicates integer；
tmIndicate that rock to be measured reaches the time experienced behind ground, unit is minute；
qmIt indicates after rock to be measured reaches ground through going through tmDesorption rate after time, unit are that milliliter is per second.
The desorption rate q that rock to be measured reaches the ground moment is calculated using the expression formula (6) of the functional relation0, benefit
With the proportional relation of desorption rate and rock inside and outside differential pressure to be measured, rock to be measured is obtained during drilling in stratum different depth
Desorption rate q0i:
pfi=GDp(h-hf)=GDp(h-vt0i) (9)
In formula, GDpIndicate formation pressure gradient, unit is every meter of megapascal, is known parameters；
H indicates prime stratum depth locating for rock to be measured, and unit is rice, is known quantity；
V indicates that the Mean Speed promoted in rock lifting process to be measured, unit are metre per second (m/s)s, is known quantity；
t0iIndicate that rock to be measured was promoted to a certain depth elapsed time, unit is minute；
pfiIndicate strata pressure when rock to be measured is promoted to a certain depth, unit is megapascal；
pfIt indicates that original formation pressure locating for rock to be measured, unit are megapascal, is known quantity.
And then when calculating the brill by following algorithm and encountering the rock to be measured to rock to be measured this section of ground of arrival
Between internal loss gas volume Vl1:
Vl1=q0i(t0i-t0i-1)+q0i-1(t0i-1-t0i-2)+…+q02(t02-t01)+q01t01 (10)
By the Vl1With the Vl2It is added, obtains the volume V of loss gasl。
Naturally after desorbing 900 minutes, opening ultrasonic oscillator 5 makes rock 1 to be measured under the dissemination of filler 2
The fast desorption in desorption container 3 reads gas flowmeter 9 and shows nature until gas flowmeter 9 shows that desorption rate is 0
The cumulative volume V of stripping gas obtained by desorption and fast desorptiond。
By the volume V of the loss gaslWith the cumulative volume V of the stripping gasdIt is added, obtains the total of gas contained by rock
Volume calculates the total volume of gas contained by the rock and the ratio of rock quality to be measured, obtains the gassiness of the rock to be measured
Amount.Then the gas in gas bottle 11 can be subjected to component analysis or exhaust-gas treatment etc..
Supersonic oscillations mode is used in the present embodiment, and the rock to be measured is made to carry out fast desorption in desorption apparatus,
The process for opening desorption apparatus test residual gas content is eliminated, simplifies continuous mode, while fast desorption shortens desorption
Time；In addition, by the environment in the pressure gauge 6 and the adjustment desorption apparatus of temp-controlling element 4 in desorption apparatus, simulated atmosphere ring
Border, and by the improvement of algorithm, improve the accuracy of measurement result.
As another preferred embodiment of the invention, the specific embodiment of the measuring method of the rock gassiness data
Are as follows:
Rock to be measured is fitted into the desorption apparatus of no temp-controlling element and pressure gauge, it is not any in the desorption apparatus
Filler desorbs the rock to be measured naturally 500 minutes, utilizes the desorption speed of gas contained by gas detector monitoring rock
Rate is obtained according to the desorption rate data monitored in the predetermined time by least-square fitting approach, is desorbed naturally
The functional relation of desorption rate and desorption time, functional relation expression formula are as follows:
Q (t)=atn+btn-1+…+ct2+dt1+e (11)
In formula, n indicates integer；
A, b ... ... c, d, e indicate fitting coefficient；
T indicates desorption time, and unit is minute
Q (t) indicates the desorption rate when desorption t time, and unit is that milliliter is per minute.
The rock to be measured, which is calculated, using the expression formula (11) of the functional relation for being fitted and obtaining reaches ground to dress
Enter to solve in cucurbitula this period, loses air volume Vl2, the volume V of the period internal loss gas is calculated by following algorithml2:
Vl2=qm(tm-tm-1)+qm-1(tm-1-tm-2)+…+q2(t2-t1)+q1t1 (12)
In formula, m indicates integer；
tmIndicate that rock to be measured reaches the time experienced behind ground, unit is minute；
qmIt indicates after rock to be measured reaches ground through going through tmDesorption rate after time, unit are that milliliter is per second.
The desorption rate q that rock to be measured reaches the ground moment is calculated using the expression formula (11) of the functional relation0, benefit
With the proportional relation of desorption rate and rock inside and outside differential pressure to be measured, rock to be measured is obtained during drilling in stratum different depth
Desorption rate q0i:
pfi=GDp(h-hf)=GDp(h-vt0i) (14)
In formula, GDpIndicate formation pressure gradient, unit is every meter of megapascal, is known parameters；
H indicates prime stratum depth locating for rock to be measured, and unit is rice, is known quantity；
V indicates the Mean Speed that rock lifting process to be measured is promoted, and unit is metre per second (m/s), is known quantity；
t0iIndicate that rock to be measured was promoted to a certain depth elapsed time, unit is minute；
pfiIndicate strata pressure when rock to be measured is promoted to a certain depth, unit is megapascal；
pfIt indicates that original formation pressure locating for rock to be measured, unit are megapascal, is known quantity.
And then when calculating the brill by following algorithm and encountering the rock to be measured to rock to be measured this section of ground of arrival
Between internal loss gas volume Vl1:
Vl1=q0i(t0i-t0i-1)+q0i-1(t0i-1-t0i-2)+…+q02(t02-t01)+q01t01。(15)
By the Vl1With the Vl2It is added, obtains the volume V of loss gasl。
Naturally after desorbing 500 minutes, by vibrating desorption apparatus Quick mechanical, by rock to be measured in the desorption apparatus
It is interior shatter, so that the rock fast desorption to be measured passes through until gas detector shows that the concentration of gas contained by rock is 0
Gas detector measures the cumulative volume V of stripping gas obtained by nature desorption and fast desorptiond.By the volume V of the loss gaslWith
The cumulative volume V of the stripping gasdIt is added, the total volume of gas contained by rock is calculated, calculate addition result and rock to be measured
The ratio of stone quality obtains the air content of the rock to be measured.The component of the gas as contained by known rock to be measured and determine to
Gas environmental sound contained by rock is surveyed, therefore the gas of eluting gas detector is directly discharged in atmosphere.
In this preferred embodiment, the measurement device is also equipped with calculation processing unit, can be used for obtaining the monitoring
The cumulative volume data of desorption rate data and stripping gas that device measurement obtained desorb naturally；
According to the desorption rate data desorbed naturally that measurement obtains, nature desorption speed is obtained using the fitting algorithm of setting
The functional relation of rate and desorption time；
Using the proportional relation of desorption rate and rock inside and outside differential pressure to be measured, rock to be measured during drilling is calculated and exists
The desorption rate q of stratum different depth0i；
Using the calculation formula of setting, the volume of the loss gas in the second time period is calculated；
Using the calculation formula of setting, the volume of the loss gas in the first time period is calculated；
It can be also used for calculate the loss volume of gas in the resulting first time period, in the first time period
The volume of loss gas, the gained of the desorption naturally volume of stripping gas, stripping gas obtained by the fast desorption volume addition,
Obtain the total volume of gas contained by the rock to be measured.
The fitting algorithm of the setting includes: polynomial fitting；Least square fitting algorithm.
The volume V of loss gas in the calculating second time period of the settingl2Calculation formula include:
Vl2=qm(tm-tm-1)+qm-1(tm-1-tm-2)+…+q2(t2-t1)+q1t1
In formula, m indicates integer；
tmIndicate that rock to be measured reaches the time experienced behind ground, unit is minute；
qmIt indicates after rock to be measured reaches ground through going through tmDesorption rate after time, unit are that milliliter is per second.
The volume V of loss gas in the calculating second time period of the settingl1Calculation formula may include:
Vl1=q0i(t0i-t0i-1)+q0i-1(t0i-1-t0i-2)+…+q02(t02-t01)+q01t01
In formula, i indicates integer；
t0iIndicate that rock to be measured was promoted to a certain depth elapsed time, unit is minute；
q0iIndicate that rock to be measured is in the desorption rate of stratum different depth during drilling.
In this preferred embodiment, described device is further comprised: air content computing module, and being used for will be contained by the rock to be measured
The total volume of gas calculates the air content of rock to be measured divided by the quality of the rock to be measured.
In this preferred embodiment, nature desorption time is shortened, by the mode that rock to be measured is shatter, makes rock to be measured
The rate of fast desorption is higher, further shortens minute；Select another fit approach, using with a upper embodiment
Identical calculation method has equally also achieved the effect that improve measurement result accuracy, simultaneously because including in measurement device
Calculation processing module, the data that can carry out corresponding each process automatically calculate, and do not need implementation personnel and participate in calculating, in this way
Implementation process is simplified, user experience is improved.
Although this application provides the method operating procedure as described in embodiment or flow chart, based on conventional or noninvasive
The means for the property made may include more or less operating procedure.The step of enumerating in embodiment sequence is only numerous steps
One of execution sequence mode, does not represent and unique executes sequence.It, can when device or client production in practice executes
To execute or parallel execute (such as at parallel processor or multithreading according to embodiment or method shown in the drawings sequence
The environment of reason, even distributed data processing environment).The terms "include", "comprise" or its any other variant are intended to contain
Lid non-exclusive inclusion, so that process, method, product or equipment including a series of elements are not only wanted including those
Element, but also including other elements that are not explicitly listed, or further include for this process, method, product or equipment
Intrinsic element.In the absence of more restrictions, be not precluded include the process, method of the element, product or
There is also other identical or equivalent elements in person's equipment.
Unit, device or module that above-described embodiment illustrates etc. can specifically realize by computer chip or entity, or
It is realized by the product with certain function.For convenience of description, various modules point are divided into function when describing apparatus above
It does not describe.It certainly, when implementing the application can the function of each module is real in the same or multiple software and or hardware
It is existing, the module for realizing same function can also be realized by the combination of multiple submodule or subelement etc..Dress described above
Set that embodiment is only schematical, for example, the division of the unit, only a kind of logical function partition, in actual implementation
There may be another division manner, such as multiple units or components can be combined or can be integrated into another system or one
A little features can be ignored, or not execute.Another point, shown or discussed mutual coupling or direct-coupling or communication link
Connecing can be through some interfaces, the indirect coupling or communication connection of device or unit, can be electrical property, mechanical or other shapes
Formula.
It is also known in the art that other than realizing controller in a manner of pure computer readable program code, it is complete
Entirely can by by method and step carry out programming in logic come so that controller with logic gate, switch, specific integrated circuit, programmable
Logic controller realizes identical function with the form for being embedded in microcontroller etc..Therefore this controller is considered one kind
Hardware component, and the structure that the device for realizing various functions that its inside includes can also be considered as in hardware component.Or
Person even, can will be considered as realizing the device of various functions either the software module of implementation method can be hardware again
Structure in component.
The application can describe in the general context of computer-executable instructions executed by a computer, such as program
Module.Generally, program module includes routines performing specific tasks or implementing specific abstract data types, programs, objects, group
Part, data structure, class etc..The application can also be practiced in a distributed computing environment, in these distributed computing environments,
By executing task by the connected remote processing devices of communication network.In a distributed computing environment, program module can
To be located in the local and remote computer storage media including storage equipment.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
It realizes by means of software and necessary general hardware platform.Based on this understanding, the technical solution essence of the application
On in other words the part that contributes to existing technology can be embodied in the form of software products, the computer software product
It can store in storage medium, such as ROM/RAM, magnetic disk, CD, including some instructions are used so that a computer equipment
(can be personal computer, mobile terminal, server or the network equipment etc.) executes each embodiment of the application or implementation
Method described in certain parts of example.
Each embodiment in this specification is described in a progressive manner, the same or similar portion between each embodiment
Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.The application can be used for crowd
In mostly general or special purpose computing system environments or configuration.Such as: personal computer, server computer, handheld device or
Portable device, laptop device, multicomputer system, microprocessor-based system, set top box, programmable electronics set
Standby, network PC, minicomputer, mainframe computer, distributed computing environment including any of the above system or equipment etc..
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application there are many deformation and
Variation is without departing from spirit herein, it is desirable to which the attached claims include these deformations and change without departing from the application's
Spirit.
Claims (14)
1. a kind of measuring method of rock gassiness data, which is characterized in that the described method includes:
Rock to be measured is fitted into progress nature desorption in desorption apparatus, monitor desorption rate and obtains the obtained solution of nature desorption
The volume of air-breathing；
Fitting obtains the functional relation of the desorption rate desorbed naturally and desorption time, according to the functional relation, respectively
The volume in first time period with second time period internal loss gas is calculated, total losses air volume is obtained；The first time period is
Brill encounters the rock to be measured to rock to be measured and reaches ground this period, and the second time period is that the rock to be measured reaches
Ground is extremely packed into desorption apparatus this period；
Naturally after the time desorbed reaches predetermined period, the rock to be measured is quickly solved in the desorption apparatus
It inhales, and monitors the fast desorption；Until monitoring that gas contained by rock to be measured is desorbed completely, the fast desorption institute is measured
The volume of obtained stripping gas；
By stripping gas obtained by the volume of the loss gas, the volume of the gained stripping gas of desorption naturally, the fast desorption
Volume is added, and obtains the total volume of contained gas in the rock to be measured.
2. the measuring method of rock gassiness data as described in claim 1, which is characterized in that the fitting is quasi- including multinomial
It closes, the fitting obtains the functional relation of the desorption rate desorbed naturally and desorption time, comprising:
Using the data of the gained desorption rate of desorption naturally, fitting coefficient is calculated, the function for obtaining fitting of a polynomial closes
It is expression formula.
3. the measuring method of rock gassiness data as claimed in claim 2, which is characterized in that the function of the fitting of a polynomial
Relational expression includes:
Q (t)=antn+an-1tn-1+…+a2t2+a1t1+a0
In formula, n indicates integer；
an,an-1,……a1,a0Indicate fitting coefficient；
T indicates desorption time, and unit is minute；
Q (t) indicates the desorption rate when desorption t time, and unit is that milliliter is per minute.
4. the measuring method of rock gassiness data as described in claim 1, which is characterized in that calculate in the second time period
The volume of loss gas include:
Using the functional relation for being fitted and obtaining, it is calculated by the following formula out the body of the second time period internal loss gas
Product Vl2:
Vl2=qm(tm-tm-1)+qm-1(tm-1-tm-2)+...+q2(t2-t1)+q1t1
In formula, m indicates integer；
tmIndicate that rock to be measured reaches the time experienced behind ground, unit is minute；
qmIt indicates after rock to be measured reaches ground through going through tmDesorption rate after time, unit are that milliliter is per second.
5. the measuring method of rock gassiness data as described in claim 1, which is characterized in that the second time period internal loss
The volume V of gasl1Evaluation method include: to be fitted the functional relation that obtains using described, calculate rock to be measured and reach ground
The desorption rate q at moment0, using the proportional relation of desorption rate and rock inside and outside differential pressure to be measured, obtain to be measured during drilling
Desorption rate q of the rock in stratum different depth0i:
pfi=GDp(h-hf)=GDp(h-vt0i)
In formula, GDpIndicate formation pressure gradient, unit is every meter of megapascal；
H indicates prime stratum depth locating for rock to be measured, and unit is rice；
V indicates that the Mean Speed promoted in rock lifting process to be measured, unit are metre per second (m/s)s；
t0iIndicate that rock to be measured was promoted to a certain depth elapsed time, unit is minute；
pfiIndicate strata pressure when rock to be measured is promoted to a certain depth, unit is megapascal；
pfIt indicates that original formation pressure locating for rock to be measured, unit are megapascal, is known quantity；
And then by following algorithm calculate it is described brill encounter the rock to be measured to rock to be measured arrival ground this period in
Lose the volume V of gasl1:
Vl1=q0i(t0i-t0i-1)+q0i-1(t0i-1-t0i-2)+...+q02(t02-t01)+q01t01
In formula, i indicates integer；
t0iIndicate that rock to be measured was promoted to a certain depth elapsed time, unit is minute；
q0iIndicate that rock to be measured is in the desorption rate of stratum different depth during drilling.
6. the measuring method of rock gassiness data as described in claim 1, which is characterized in that the rock to be measured carries out quick
Desorption includes being desorbed using mode of oscillation to rock to be measured.
7. the measuring method of the rock gassiness data as described in any one of claim 1 to 6, which is characterized in that the side
Method further include:
The ratio for calculating the total volume of contained gas and the rock quality to be measured in the rock to be measured obtains the rock to be measured
The air content of stone.
8. a kind of measurement device of rock gassiness data characterized by comprising
Cucurbitula is solved, for holding rock to be measured, the desorption environment of rock to be measured is adjusted, to rock to be measured progress nature desorption and fastly
Speed desorption；
With desorb tank connected desorption rate monitoring device, the cumulative volume for real-time measurement desorption rate and stripping gas；
The gas collector being connected with monitoring device, for collecting nature desorption and the resulting stripping gas of fast desorption；
Calculation processing unit, for obtaining the desorption rate data desorbed naturally and stripping gas that the monitoring device measurement obtains
Cumulative volume data；
According to the obtained desorption rate data desorbed naturally of measurement, using the fitting algorithm of setting obtain nature desorption rate with
The functional relation of desorption time；
Using the proportional relation of desorption rate and rock inside and outside differential pressure to be measured, rock to be measured is calculated during drilling on stratum
The desorption rate q of different depth0i；
Using the calculation formula of setting, the volume of the loss gas in second time period is calculated；
Using the calculation formula of setting, the volume of the loss gas in first time period is calculated；
It is also used to calculate the volume for losing gas in the resulting first time period, the loss gas in the first time period
Volume, the gained of the desorption naturally volume of stripping gas, the volume of stripping gas obtained by the fast desorption be added, obtain described
The total volume of gas contained by rock to be measured.
9. the measurement device of rock gassiness data as claimed in claim 8, which is characterized in that the desorption is canned pressure gauge
And temp-controlling element, for adjusting the environment in solution cucurbitula.
10. the measurement device of rock gassiness data as claimed in claim 8, which is characterized in that the desorption is canned oscillation
Device, for carrying out fast desorption to the rock to be measured.
11. the measurement device of rock gassiness data as claimed in claim 8, which is characterized in that the fitting algorithm of the setting
It include: polynomial fitting；
Least square fitting algorithm.
12. the measurement device of rock gassiness data as claimed in claim 8, which is characterized in that described in the calculating of the setting
The volume V of loss gas in second time periodl2Calculation formula include:
Vl2=qm(tm-tm-1)+qm-1(tm-1-tm-2)+...+q2(t2-t1)+q1t1
In formula, m indicates integer；
tmIndicate that rock to be measured reaches the time experienced behind ground, unit is minute；
qmIt indicates after rock to be measured reaches ground through going through tmDesorption rate after time, unit are that milliliter is per second.
13. the measurement device of rock gassiness data as claimed in claim 8, which is characterized in that described in the calculating of the setting
The volume V of loss gas in second time periodl1Calculation formula include:
Vl1=q0i(t0i-t0i-1)+q0i-1(t0i-1-t0i-2)+...+q02(t02-t01)+q01t01
In formula, i indicates integer；
t0iIndicate that rock to be measured was promoted to a certain depth elapsed time, unit is minute；
q0iIndicate that rock to be measured is in the desorption rate of stratum different depth during drilling.
14. the measurement device of the rock gassiness data as described in any one of claim 8 to 13, which is characterized in that described
Device further include:
Air content computing module, for the quality by the total volume of gas contained by the rock to be measured divided by the rock to be measured,
Calculate the air content of rock to be measured.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610725811.7A CN106290058B (en) | 2016-08-25 | 2016-08-25 | A kind of measuring method and device of rock gassiness data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610725811.7A CN106290058B (en) | 2016-08-25 | 2016-08-25 | A kind of measuring method and device of rock gassiness data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106290058A CN106290058A (en) | 2017-01-04 |
| CN106290058B true CN106290058B (en) | 2019-05-07 |
Family
ID=57616538
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610725811.7A Active CN106290058B (en) | 2016-08-25 | 2016-08-25 | A kind of measuring method and device of rock gassiness data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106290058B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107014717B (en) * | 2017-03-29 | 2024-04-12 | 上海瑞达峰致能源科技股份有限公司 | Method and device for testing desorption gas amount of lost gas in shale gas well |
| CN108254289B (en) * | 2017-12-18 | 2020-07-10 | 中国石油天然气股份有限公司 | Method and device for determining gas content of shale |
| CN108761033B (en) * | 2018-05-03 | 2023-10-31 | 中国石油天然气股份有限公司 | Method and device for determining total content of shale loss gas and storage medium |
| CN111175471B (en) * | 2018-11-12 | 2022-06-24 | 中石化石油工程技术服务有限公司 | Rock debris logging method |
| CN109765097B (en) * | 2018-12-21 | 2021-07-02 | 中铁三局集团有限公司 | Tunnel surrounding rock rapid classification method based on RPD drilling machine |
| CN112378812B (en) * | 2020-11-06 | 2021-12-14 | 西南石油大学 | Experimental device and method for determining desorption rate of adsorption type shale gas |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2112004C (en) * | 1992-12-21 | 2000-02-22 | Kenneth M. Ledez | Gas amount and solubility investigation apparatus |
| AU2002366235A1 (en) * | 2001-12-14 | 2003-06-30 | Hydro Geo Chem, Inc. | Method and system for estimating gas production by a landfill or other subsurface source |
| CN103063545A (en) * | 2012-12-28 | 2013-04-24 | 恒泰尚合能源技术（北京）有限公司 | Novel shale gas content tester and shale gas content detecting method |
| CN103196702A (en) * | 2013-04-26 | 2013-07-10 | 中联煤层气国家工程研究中心有限责任公司 | Method for measuring gas content of coal bed through airtight coring |
| CN103308634A (en) * | 2013-05-30 | 2013-09-18 | 汪双清 | Device and method for measuring gas content |
| CN103776713A (en) * | 2014-01-09 | 2014-05-07 | 中国石油天然气股份有限公司 | Dense rock gas desorption rate testing device |
| CN104155207A (en) * | 2013-05-14 | 2014-11-19 | 中国石油化工股份有限公司 | Gas content desorption testing device |
| CN105203428A (en) * | 2015-11-04 | 2015-12-30 | 中国地质科学院地质力学研究所 | Method for determining lost gas content in gas content of shale |
| CN105675434A (en) * | 2016-03-16 | 2016-06-15 | 中国地质大学（北京） | System and method for measuring gas content |
| CN206020213U (en) * | 2016-08-23 | 2017-03-15 | 中国石油天然气股份有限公司 | Rock content measurement device |
-
2016
- 2016-08-25 CN CN201610725811.7A patent/CN106290058B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2112004C (en) * | 1992-12-21 | 2000-02-22 | Kenneth M. Ledez | Gas amount and solubility investigation apparatus |
| AU2002366235A1 (en) * | 2001-12-14 | 2003-06-30 | Hydro Geo Chem, Inc. | Method and system for estimating gas production by a landfill or other subsurface source |
| CN103063545A (en) * | 2012-12-28 | 2013-04-24 | 恒泰尚合能源技术（北京）有限公司 | Novel shale gas content tester and shale gas content detecting method |
| CN103196702A (en) * | 2013-04-26 | 2013-07-10 | 中联煤层气国家工程研究中心有限责任公司 | Method for measuring gas content of coal bed through airtight coring |
| CN104155207A (en) * | 2013-05-14 | 2014-11-19 | 中国石油化工股份有限公司 | Gas content desorption testing device |
| CN103308634A (en) * | 2013-05-30 | 2013-09-18 | 汪双清 | Device and method for measuring gas content |
| CN103776713A (en) * | 2014-01-09 | 2014-05-07 | 中国石油天然气股份有限公司 | Dense rock gas desorption rate testing device |
| CN105203428A (en) * | 2015-11-04 | 2015-12-30 | 中国地质科学院地质力学研究所 | Method for determining lost gas content in gas content of shale |
| CN105675434A (en) * | 2016-03-16 | 2016-06-15 | 中国地质大学（北京） | System and method for measuring gas content |
| CN206020213U (en) * | 2016-08-23 | 2017-03-15 | 中国石油天然气股份有限公司 | Rock content measurement device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106290058A (en) | 2017-01-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106290058B (en) | A kind of measuring method and device of rock gassiness data | |
| CN105319603B (en) | The Forecasting Methodology of tight sandstone reservoir complexity chicken-wire cracking | |
| Diallo et al. | Age of stratospheric air in the ERA-Interim | |
| CN106442938B (en) | A kind of accurate survey calculation method for obtaining shale air content | |
| Kummerow et al. | Experimental evaluation of the impact of the interactions of CO2‐SO2, brine, and reservoir rock on petrophysical properties: A case study from the Ketzin test site, Germany | |
| CN108169062B (en) | Simulate the visual test device and method of subterranean coal gas preservation desorption process | |
| CN107121359A (en) | Shaked out mechanics parameter coupling process analogue means and method containing hydrate sediment | |
| CN105089615B (en) | A kind of log data historical regression processing method based on reservoir model | |
| CN103913772B (en) | Micro-seismic event the Forward Modeling based on reservoir geology mechanics parameter | |
| CN105223099B (en) | Shale gas gassiness tester and its test method | |
| Dontsov et al. | A new technique for proppant schedule design | |
| Pitman et al. | A model of granular flows over an erodible surface | |
| CN106842364A (en) | Unconformity biogenic limestone seismic reservoir petrophysical model method for building up | |
| CN105021493B (en) | A kind of absorption analysis method and device of multicomponent gas | |
| CN106153495A (en) | Rock gas absorption phase volume measuring method and device | |
| CN106959464A (en) | The measurement apparatus and measuring method of a kind of precipitation rate of radon | |
| Gavrilov et al. | Verifications of the high-resolution numerical model and polarization relations of atmospheric acoustic-gravity waves | |
| CN106677766B (en) | A kind of experimental provision for measuring Breakup of Liquid Ring air and invading rear low-frequency elastic wave response characteristic | |
| CN105607123B (en) | A kind of method and apparatus for calculating random pore dielectric model seismic wave characteristic information | |
| CN105629319B (en) | The apparatus and method for determining stratum geochemical well logging standard spectrum and relative sensitivity | |
| Pisso et al. | Turbulent vertical diffusivity in the sub-tropical stratosphere | |
| CN209990484U (en) | Giant magnetostrictive forced vibration permeability increasing system | |
| CN104792684B (en) | Sand shale porosity type inversion method based on 3D seismic data | |
| CN111220497A (en) | Shale reservoir fluid utilization range quantitative evaluation device and method after seam network transformation | |
| CN109596497A (en) | Simulate the experimental provision and method of sand body liquefaction intrusion trigger mechanism and forming process |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |