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 PDF

Info

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
gas
volume
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610725811.7A
Other languages
Chinese (zh)
Other versions
CN106290058A (en
Inventor
薛华庆
周尚文
郭伟
申卫兵
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Natural Gas Co Ltd
Original Assignee
China Petroleum and Natural Gas Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN201610725811.7A priority Critical patent/CN106290058B/en
Publication of CN106290058A publication Critical patent/CN106290058A/en
Application granted granted Critical
Publication of CN106290058B publication Critical patent/CN106290058B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N7/00Analysing materials by measuring the pressure or volume of a gas or vapour
    • G01N7/14Analysing 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

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

A kind of measuring method and device of rock gassiness data
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.
CN201610725811.7A 2016-08-25 2016-08-25 A kind of measuring method and device of rock gassiness data Active CN106290058B (en)

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 (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108254289B (en) * 2017-12-18 2020-07-10 中国石油天然气股份有限公司 Method and device for determining gas content of shale
CN109765097B (en) * 2018-12-21 2021-07-02 中铁三局集团有限公司 Tunnel surrounding rock rapid classification method based on RPD drilling machine

Citations (10)

* Cited by examiner, † Cited by third party
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

Patent Citations (10)

* Cited by examiner, † Cited by third party
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
Sobieski et al. Darcy’s and Forchheimer’s laws in practice. Part 1. The experiment
CN103913772B (en) Micro-seismic event the Forward Modeling based on reservoir geology mechanics parameter
Hall et al. Timescales for the stratospheric circulation derived from tracers
Dontsov et al. A new technique for proppant schedule design
CN106124349A (en) Determine the critical erosion rate of gas well, the method and device of flow
CN106355571A (en) Determination method and device for dolostone reservoir quality
CN104963657B (en) Numerical reservoir simulation method and device based on the constraint of high water-cut stage Monitoring Data
CN105549088B (en) The recognition methods of Fractured tight sand Mesosphere and device
CN107121354A (en) A kind of device and its application method for testing Karst grouting material dispersion resistance matter under water
CN106842364A (en) Unconformity biogenic limestone seismic reservoir petrophysical model method for building up
CN102654045B (en) Top air injection gravity-assisted flooding oil-extraction one-dimensional physical simulation experimental system
CN108169062B (en) Simulate the visual test device and method of subterranean coal gas preservation desorption process
CN106597542A (en) Method and device for predicting reservoir characteristic parameters
CN106153495A (en) Rock gas absorption phase volume measuring method and device
CN105629319B (en) The apparatus and method for determining stratum geochemical well logging standard spectrum and relative sensitivity
CN105223099B (en) Shale gas gassiness tester and its test method
CN106842316A (en) Crack determines method and apparatus
CN106338768A (en) Processing method, device and system for generating reservoir prediction attribute data
Michaels et al. Natural signals in the MSU lower tropospheric temperature record
CN208505928U (en) Microcosmic detection simulation device is evaluated in gravel packing zone blocking
CN109541667A (en) A kind of Radon eduction simulator and radon release rate method
CN106677766B (en) A kind of experimental provision for measuring Breakup of Liquid Ring air and invading rear low-frequency elastic wave response characteristic
CN107167844B (en) A kind of thin sand-mud interbed fluid substitution method
CN207396259U (en) A kind of coal bed gas electronics desorption apparatus

Legal Events

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