Pressure loss degree prediction method suitable for pressure maintaining coring technology of submarine drilling machine
Technical Field
The invention relates to a pressure loss degree prediction method suitable for a pressure maintaining and coring technology of a submarine drilling machine.
Background
The natural gas hydrate is taken as an efficient and clean new energy source, is regarded as an ideal substitute energy source for future petroleum and coal by countries around the world, is mainly endowed with submarine sediments and land permafrost zones, and has reserves which are 2 times of the reserves of known coal, petroleum and natural gas worldwide. Sea natural gas hydrates exist in solid form in high pressure and low temperature environments, and the water depth of the sea natural gas hydrates is generally between 1200 meters and 2000 meters. Submarine natural gas hydrates have instability at normal temperature and pressure, and are usually explored by adopting a pressure maintaining and coring mode. When the submarine drilling machine is used for working, all drilling pipes required for the whole drilling are stored in a drilling tool warehouse on the drilling machine and are lowered to the seabed together with the drilling machine. After the back drilling is finished, the drill pipe filled with the core is temporarily stored in a drilling tool warehouse until all the operation tasks are finished, and then the drill pipe is lifted to the sea surface along with a submarine drilling machine. Compared with the recovery mode that the traditional drilling ship can rapidly lift a single drilling tool out of the water surface after the back drilling is finished, the lifting speed of the submarine drilling machine is changed in different sea water intervals in the whole lifting recovery process, and the duration of the submarine drilling machine in the sea water in the whole lifting process is longer. The submarine drilling machine is provided with a small-caliber thin-wall pressure-maintaining coring drilling tool based on a thin-wall core tube direct sealing principle. The existing pressure maintaining performance calculating method is mainly developed around a thick-wall pressure maintaining cylinder with the ratio of inner diameter to outer diameter being larger than 1.2. In the conventional pressure maintaining performance prediction process, the pressure maintaining performance of the thin-wall pressure maintaining core drilling tool is usually qualitatively predicted by directly adopting a pressure predicting method of a thick-wall cylinder. The influence of temperature change factors is not considered in the calculation process, or the influence of temperature change in the whole recovery process is considered integrally in a constant value coefficient mode. But the small-caliber thin-wall pressure maintaining drilling tool suitable for the submarine drilling machine needs to pass through seawater layers with different environmental temperatures in the recovery process and has longer duration. Therefore, when the traditional pressure maintaining performance prediction method predicts the pressure maintaining performance of the small-caliber thin-wall pressure maintaining drilling tool, the actual problems of larger deviation degree between a theoretical prediction result and an actual test result after water is discharged under the actual recovery working condition, unstable prediction accuracy under different working water depths and the like can be caused.
Disclosure of Invention
In order to solve the technical problems, the invention provides the pressure loss degree prediction method which can reduce the error of taking the temperature into consideration in the whole drilling tool lifting process, so that the calculation result is more accurate and is suitable for the pressure maintaining and coring technology of the submarine drilling machine.
In order to achieve the above object, the present invention provides the following solutions:
the pressure loss degree prediction method suitable for the pressure maintaining and coring technology of the submarine drilling machine comprises the following steps of:
1) Performing interval segmentation on the working water depth of the submarine drilling machine;
2) Calculating the volume change delta V of each interval under the action of internal and external pressure difference Pi I is a section sequence number, and is numbered from bottom to top;
3) Calculating the volume expansion delta V of the core tube due to temperature in each interval Ti ,
4) The total volume change DeltaVi of the core tube in each interval is calculated by utilizing the superposition principle,
5) The total pressure drop value delta Pi in the core tube of each interval is calculated,
6) And adding the total pressure drop values delta Pi of the core tubes in each interval to obtain a final total pressure drop value delta P of the core tubes.
In the pressure loss degree prediction method suitable for the pressure maintaining and coring technology of the submarine drilling machine, the specific operation of the step 2) is as follows:
σ ri =P i R i 2 (1-R 0 2 /r 2 )/(R 0 2 -R i 2 )
σ θi =P i R i 2 (1+R 0 2 /r 2 )/(R 0 2 -R i 2 )
σ Zi =P i R i 2 /(R 0 2 -R i 2 )
wherein sigma ri Is the radial stress and sigma of the inner core tube in the ith section θi Is tangential stress and sigma of the inner core tube in the ith interval Zi Is the axial stress in the inner core tube of the ith section, P i Is the internal pressure of the core tube in the ith section, and the internal pressure P of the core tube in the (i+1) th section i+1 =P i +ΔPi,P 1 Designing pressure for pressure maintaining and coring; r is R i Is the inner diameter of the core tube in the ith section, and the inner diameter of the core tube in the (i+1) th section is R i+1 =R i -ΔR Pi ,R 0 Is the outer diameter of the core tube, R is the radius of the stress analysis place, and r=r is taken as i ;
ε ri =[σ ri -μ(σ zi +σ θi )]/E
ε θi =[σ θi -μ(σ ri +σ Zi )]/E
ε Zi =[σ Zi -μ(σ ri +σ θi )]/E
Wherein ε ri Is the radial strain of the inner core tube in the ith section epsilon θi Is tangential strain epsilon of the inner core tube in the ith section zi Is the axial strain of the inner core tube of the ith interval; mu is Poisson's ratio, E is elastic modulus,
the change amount of the inner core pipe inner diameter in the ith section is obtained by the change amount of the thickness of the pipe body under the action of the internal and external pressure difference:
ΔR pi =δε ri
Δh pi =hε Zi
ΔV pi =π(R i +ΔR pi ) 2 (h+Δh pi )-πR i 2 h
wherein delta is the thickness of the core tube, h is the total length of the core tube, and DeltaR Pi Is the change quantity of the inner diameter of the core tube in the ith section, delta h Pi Is the axial variation of the core tube in the ith section.
In the pressure loss degree prediction method suitable for the pressure maintaining and coring technology of the submarine drilling machine, the specific operation of the step 3) is as follows:
ΔL i =α(2πR i ×ΔT i )
ΔL 0 =α(2πR 0 ×ΔT i )
Δδ Ti =α(R 0 -R i )ΔT i
Δh Ti =αhΔT i
ΔV Ti =(L i +ΔL i +L 0 +ΔL 0 )(h+Δh Ti )(δ+Δδ Ti )/2-π(R 0 2 -R i 2 )h
wherein: l (L) i Is the inner circumference of the core tube, L o Is the outer circumference of the core tube, deltaL i Is the variation of the inner circumference of the core tube in the ith section, deltaL o Delta T is the variation of the outer circumference of the core tube i For the water temperature change in the i-th zone, alpha is the temperature coefficient of linear expansion of the core tube, 1.2X10 -5 K,Δδ Ti For temperature-induced core tube thickness variation Δh Ti Core tube for temperature induced ith intervalThe length varies.
In the pressure loss degree prediction method suitable for the pressure maintaining and coring technology of the submarine drilling machine, the specific operation of the step 4) is as follows:
ΔVi=ΔV Pi +ΔV Ti 。
in the pressure loss degree prediction method suitable for the pressure maintaining and coring technology of the submarine drilling machine, the specific operation of the step 5) is as follows:
ΔP=-KΔVi/V
wherein K is the bulk modulus of the pattern; v is the total volume of the core tube.
Compared with the prior art, the invention has the beneficial effects that:
(1) According to the invention, the working water depth of the submarine drilling machine is segmented in intervals, so that the error of integrally considering the temperature to the drilling tool lifting process can be reduced, and the calculation result is more accurate.
(2) According to the invention, the influence of temperature change on the pressure change in the pipe is considered, and the internal volume change caused by the temperature change is replaced by calculating the pipe body expansion caused by the temperature change, so that the internal pressure error caused by the temperature change is reduced.
Drawings
FIG. 1 is a flow chart of the present invention.
Detailed Description
The technical scheme of the present invention will be clearly and completely described in the following in conjunction with the accompanying drawings and examples.
Taking the drilling depth of 2000m as an example, the pressure drop in the core pipe lifted from the seabed to the sea surface is calculated, and the method comprises the following steps of:
1) The working water depth of 2000 meters is divided into three parts according to the temperature distribution of the seawater: the first interval is 2000m-350m, the second interval is 350m-30m, and the third interval is 30m-0m. The first interval pressure is 20MPa, the second interval pressure is 3.5MPa, and the third interval pressure is 0.3MPa. The first interval water temperature is changed to 1-5 ℃, the first interval water temperature is changed to 5-30 ℃, and the first interval water temperature is changed to 30-30 ℃.
2) Calculating the volume change quantity of the inside of the pipe body under the action of the internal and external pressure difference of the core pipe in each interval, wherein the calculation formula is as follows:
from the elastic mechanics in the material mechanics, it is known that: radial stress sigma of core tube in ith section ri Tangential stress sigma θi And axial stress sigma Zi The method comprises the following steps:
σ ri =P i R i 2 (1-R 0 2 /r 2 )/(R 0 2 -R i 2 )
σ θi =P i R i 2 (1+R 0 2 /r 2 )/(R 0 2 -R i 2 )
σ Zi =P i R i 2 /(R 0 2 -R i 2 )
wherein sigma ri Is the radial stress and sigma of the inner core tube in the ith section θi Is tangential stress and sigma of the inner core tube in the ith interval Zi Is the axial stress in the core tube in the ith section, i is the section serial number, and the number is from bottom to top, P i Is the internal pressure of the core tube in the ith zone, R i Is the inner diameter of the core tube in the ith section, so the inner diameter of the core tube in the (i+1) th section is R i+1 =R i -ΔR Pi ,R 0 Is the outer diameter of the core tube, R is the radius of the stress analysis place, and r=r is taken as i ;P i =20MPa,R i =33.35mm,R 0 =36.6mm,r=R i 。
The strain epsilon of the core tube in the ith zone corresponding to the stress can be obtained from the stress ri 、ε θi And epsilon Zi The method comprises the following steps:
ε ri =[σ ri -μ(σ zi +σ θi )]/E
ε θi =[σ θi -μ(σ ri +σ Zi )]/E
ε Zi =[σ Zi -μ(σ ri +σ θi )]/E
wherein: epsilon ri Is the radial strain of the inner core tube in the ith section epsilon θi Is tangential strain epsilon of the inner core tube in the ith section zi Is the ith regionAxial strain of the inner core tube; μ is poisson's ratio, E is elastic modulus, e=2.01X1011 Pa.
The change of the inner diameter of the core tube is obtained by the change of the thickness of the tube body under the action of the internal and external pressure difference:
ΔR pi =δε ri ;
Δh pi =hε Zi ;
wherein: delta is the thickness of the core tube, delta=3.25 mm, h is the total length of the core tube, h=2670 mm, Δr Pi Is the change quantity of the inner diameter of the core tube in the ith section, delta h Pi Is the axial variation of the core tube in the ith section.
ΔV pi =π(R i +ΔR pi ) 2 (h+Δh pi )-πR i 2 h。
3) Calculating the volume expansion quantity of the pipe body caused by temperature, fully considering the expansion of the inner pipe and the outer pipe, replacing the volume change quantity in the pipe with the volume expansion quantity of the core pipe, and adopting the calculation formula as follows:
ΔL i =α(2πR i ×ΔT i );
ΔL 0 =α(2πR 0 ×ΔT i );
Δδ Ti =α(R 0 -R i )ΔT i ;
Δh Ti =αhΔT i ;
ΔV Ti =(L i +ΔL i +L 0 +ΔL 0 )(h+Δh Ti )(δ+Δδ Ti )/2-π(R 0 2 -R i 2 )h;
wherein L is i Is the inner circumference of the core tube, L o Is the outer circumference of the core tube, deltaL i Is the variation of the inner circumference of the core tube in the ith section, deltaL o Delta T is the variation of circumference of the outer circumference of the core tube i For the water temperature change in the i-th zone, alpha is the temperature coefficient of linear expansion of the core tube, 1.2X10 -5 K,Δδ Ti For temperature-induced core tube thickness variation Δh Ti The core tube length in the i-th zone varies due to temperature.
4) The total volume change of each interval is calculated by utilizing the superposition principle, and the calculation formula is as follows:
ΔVi=ΔV Pi +ΔV Ti 。
5) The integral pressure drop value in the pipe body of the single interval is calculated, and the calculation formula is as follows:
P i+1 =P i +ΔPi;P 1 designing pressure for pressure maintaining and coring;
wherein K is the bulk modulus of elasticity of the sample; v is the total volume of the core tube.
6) And adding the total pressure drop values delta Pi of the core tubes in each interval to obtain a final total pressure drop value delta P of the core tubes.
The water depth of the drilling tool working is 2000m for segment analysis:
a:2000-350m depth of water
In the area, the water temperature is increased from 1 ℃ to 5 ℃, and the pressure outside the sampling core tube is linearly reduced from 20MPa to 3.5MPa. The initial pressure of the inner core tube is 20MPa, the highest temperature in the area is changed to 4 ℃, and the maximum pressure difference between the inside and the outside of the pipeline is 17.5MPa. The above steps show that the pressure drop at 350m and the pressure in the tube are:
ΔP1=-0.2837MPa,P 2 =P 1 +ΔP1=19.7163MPa
b: water depth of 350-30m
In the area, the water temperature is increased from 5 ℃ to 30 ℃, and the pressure outside the sampling core tube is linearly reduced from 3.5MPa to 0.3MPa. The initial pressure of the inner core tube is 19.7163MPa. The maximum temperature change in this region was 25℃and the maximum pressure difference between the inside and outside of the pipe was 19.4163MPa. The calculation mode is the same as that above, and the pressure drop at 30 meters and the pressure in the pipeline are as follows:
ΔP2=-0.8141MPa,P 3 =P 2 +ΔP2=18.9022MPa
c:30-0m (sea level)
In the area, the water temperature is constant at 30 ℃, and the pressure outside the sampling core tube is linearly reduced from 0.3MPa to 0.1MPa. In addition, the initial pressure of the inner core tube was 18.9022MPa. The temperature change of the area is 0 ℃, and the maximum pressure difference between the inside and the outside of the pipeline is 18.8022MPa. The sea surface pressure drop and the in-pipe pressure were calculated as:
ΔP3=-0.4297MPa,
8) From the above calculation, it is possible to: final core tube total pressure drop value Δp=Δp1+Δp2+Δp3= -1.5275MPa; the final pressure in the core tube is P 1 +ΔP=18.4725MPa。