CN110162851B - Cable formation test pumping numerical simulation and numerical correction method of process thereof - Google Patents
Cable formation test pumping numerical simulation and numerical correction method of process thereof Download PDFInfo
- Publication number
- CN110162851B CN110162851B CN201910382717.XA CN201910382717A CN110162851B CN 110162851 B CN110162851 B CN 110162851B CN 201910382717 A CN201910382717 A CN 201910382717A CN 110162851 B CN110162851 B CN 110162851B
- Authority
- CN
- China
- Prior art keywords
- pumping
- ellipsoid
- formation
- fluid
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a numerical simulation of pumping numerical simulation of a cable formation test and a numerical correction method of the process of the pumping numerical simulation, which comprises the following steps: s1, judging the flow pattern of cable stratum test pressure measurement data and solving the corresponding fluidity; s2, opening a test hole on a well wall by a probe, enabling an invaded zone and formation fluid to start flowing to a well bore, extracting the formation fluid by a sampler, enabling the formation fluid to instantly flow and form a flowing form and continuous conduction to generate a swept range, and simulating the swept range formed by the continuous conduction of the formation fluid flow by an ellipsoid; s3, on the basis of giving reservoir physical property parameters and invasion depth, obtaining the predicted pumping breakthrough time before measurement, pump pure time and corresponding volume according to an algorithm; and S4, testing pumping monitoring parameters by a cable stratum, picking up pumping breakthrough time and volume, correcting a pumping result by optimizing an algorithm fitting function and simulation calculation, obtaining the actual effect of pumping in time, and guiding engineering operation, thereby improving the operation efficiency.
Description
Technical Field
The invention relates to the technical field of petroleum and natural gas exploration and development, in particular to a numerical simulation method for a pumping numerical simulation of a cable stratum test pump and a numerical correction method for a process of the pumping numerical simulation method.
Background
The technology of cable formation testing began in 1950, and with the development of the technology, formation testing logging instruments have been able to meet five major application requirements of formation pressure testing, downhole fluid analysis, downhole fluid sampling, instability testing, micro fracturing and the like under different formation and borehole conditions.
The latest generation of cable formation testers are equipped with some necessary components. For example, an Optical Fluid Analyzer (OFA) can distinguish liquid from gas, can also distinguish water from oil, and the optical response of the OFA has good correlation with crude oil density, saturation pressure, crude oil compressibility, formation volume coefficient and gas-oil ratio.
With the development of formation testing acquisition technology, the schlumberger formation testing well logging interpretation software has also undergone many update upgrades from the original PDPlot software to the InSitu software and the InSitu Pro software, up to now the Techlog software. At present, a comprehensive analysis and application technology of Formation testing logging becomes an application module of Techlog software, namely a Formation testing suite Formation testing logging advanced analysis and application technology. While the development of Downhole fluid analysis technology (Downhole fluid analysis [ DFA ]) is rapid, DFA is mainly used for comprehensively analyzing the following data of the fluid in a pumping pipeline measured by a Downhole fluid identification module of a formation testing logging instrument, and the applied parameters include: resistivity, temperature, pressure, density, viscosity, fluorescence, bubble detection, absorption and diffraction spectra, PH, contamination level by oil-based mud.
However, the application of the current cable stratum test urgently needs to solve the following problems that 1, the mud cake plugging property is poor due to the influence of the performance of the drilling fluid, particularly in a low-porosity reservoir stratum, the invasion process of the drilling fluid is prolonged, the invasion amount is large, the reservoir stratum is seriously polluted, and the pumping effect is difficult to evaluate when the static invasion is continuously generated in the pumping and sampling stage; 2. establishing a design standard of pressure measurement sampling operation, forming a quantitative evaluation means for guiding real-time decision of pumping sampling, and solving the problem of decision standard of pressure measurement and pumping operation in production operation; establishing a pumping simulation model to provide reliable basis for operation key decision points so as to save operation time and reduce operation risk; 3. the method has the advantages that the fluid parameters are monitored by combining the simulation of the pumping process with the actual pumping, so that the representativeness of the long-time pumping sampling sample of the low-permeability reservoir is improved; and the research and application of parameters of each pumping process are deeply excavated, the reliable and effective permeability is obtained, and the DST test operation is guided.
Disclosure of Invention
The invention aims to provide a numerical simulation of pumping numerical simulation of cable formation testing and a numerical correction method of the process thereof.
In order to solve the above problems, the present invention provides a method for simulating pumping values of a cable formation test and correcting the pumping values of the cable formation test, the method comprising the steps of:
s1, obtaining basic parameters of a reservoir stratum, namely, stratum porosity, permeability, reservoir stratum thickness, slurry performance, borehole size, basic pumping parameters of a cable stratum test pump and slurry invasion conditions by using conventional logging information, judging the flow pattern of cable stratum test pressure measurement information and calculating corresponding fluidity;
s2, opening a test hole on a well wall by a probe, enabling invasion zones and formation fluid to start flowing to a well bore, extracting the formation fluid by a sampler, enabling the formation fluid to instantly flow and form a flowing form and continuous conduction to generate a swept range, and simulating the swept range formed by the continuous conduction of the formation fluid flow by an ellipsoid; the formation part of stratum flow outside the well cylinder and the part in the well cylinder are in a complementary relation and integrated into a complete ellipsoid, the size of the major axis and the minor axis of the ellipsoid depends on the permeability of a reservoir and the thickness of the reservoir, the volume relation among a cylinder, the ellipsoid and a cylindrical paraboloid in a volume model is obtained through calculation, the well cylinder is defined as a well cylinder, a pollution zone is the ellipsoid, the original formation fluid is the cylindrical paraboloid, the pressure measurement data of pumping operation is processed, and a pressure gradient curve, a spherical flow gradient curve and a radial flow gradient curve are combined together through analyzing the fluidity information of a pressure measurement point and the possible flow pattern characteristics of the pumping point, so that the flow form of the formation fluid is judged;
s3, inputting reservoir parameters, cable stratum test operation parameters and reservoir invasion depth, and calculating the volume among a cylinder (a shaft), an ellipsoid (a pollution zone) and a cylindrical projectile (undisturbed stratum fluid) of the model; in formation pressure test and fluid sampling operation, oil gas breakthrough time and completion time are main time;
s4, when the pumping operation is tested in the cable formation, the following data of the fluid in the pumping pipeline measured by the underground fluid identification module of the cable formation tester is comprehensively analyzed by utilizing an underground fluid analysis technology (DFA), the pumping result is corrected by optimizing an algorithm fitting function and simulating calculation, and the applied parameters comprise: the resistivity, temperature, pressure, density, water content, viscosity, fluorescence, bubble detection, absorption spectrum and diffraction spectrum, PH value and pollution degree of oil-based mud of the inflow/outflow pipeline; based on the pump breakthrough time, the radius of invasion of the invaded zone can be calculated, which is an important basis for the calibration of the pumping process.
Further, the volume relationship among the cylinder, the ellipsoid and the cylindrical paraboloid of the calculation model in the step S3 is:
x-x 0 ＝a 2 (y 2 +z 2 )，
(x+R) 2 +y 2 ＝R 2 ，
in the formula, x 0 Is the intersection point coordinate of the vertex of the revolution paraboloid on the x axis;
a is a cylinder paraboloid system, and a,
b is the length of the horizontal axis of the ellipsoid, and the unit is meter;
c is the length of the vertical axis of the ellipsoid, and the unit is meter;
r is the radius of the cylinder and is measured in meters;
the region enclosed by the intersecting line of the cylindrical paraboloid and the ellipsoid is V 1 ：
The volume of the ellipsoid cut by the cylinder is V 2 :
Wherein the content of the first and second substances,
a 0 ＝-c,b 0 ＝c，
n is an interval [ -c, c ],
intrusion into the space occupied by the fluid by the volume of the ellipsoid minus V 1 And V 2 Definition of V 3 Is the remaining volume of the ellipsoid cut by the formation virgin fluid and wellbore fluid:
further, the algorithm of the burst time and the completion time in step S3 is:
defining the percentage of the total test fluid volume occupied by the invaded zone fluid is called the contamination rate, during pumping by the cable formation tester, the breakthrough time for pumping is when the contamination rate =95%, and the completion time is when the contamination rate = 5%:
calculating the horizontal and vertical axes (b, c) of the ellipsoid, the units of b, c being m,
b＝c 1 ×Qa×t×DP
c＝c 2 ×Qa×t×DP)
calculating the time (t) of ellipsoid to reach the invaded zone 0 )，t 0 The unit of (a) is s,
b＝Di
calculating the speed of movement of the rotating projectile, sd 0 The unit of (b) is m/s,
sd 0 ＝c 3 ×DP×Qa
calculating the coefficient a of the ellipsoid, wherein the unit of a is m,
a＝c 4 ×K v /K h
in the above formula, the first and second carbon atoms are,
h is the thickness of the stratum in m;
K h is the formation horizontal permeability in md;
K v /K h is the anisotropy coefficient;
di is the depth of invasion in m;
uf is the formation fluid viscosity in cp;
um is the mud filtrate viscosity in cp;
ps is the position of the probe from the top boundary, the default center position is m;
r is the radius of the shaft and is m;
qa is the pump speed in m 3 /s；
t is pump discharge time in units of s;
DP is the differential pressure in MPa;
sd 0 the speed of the rotating projectile moving towards the probe;
t 0 the time for the outer boundary to reach the invaded zone during the expansion process of the ellipsoid;
b, the expression shows that the outer boundary of the ellipsoid reaches the invasion zone, and accordingly t is obtained 0 The expression of (2);
sd 0 means the speed of the rotating projectile moving in the direction of the probe;
Q a the physical meaning of xt is to reflect the expansion velocity of ellipsoids;
c in the above equation 1 、c 2 、c 3 、c 4 The coefficients are based on the relationship between the formation testing parameter Differential Pressure (DP), the testing flow rate (Qa), and the formation anisotropy coefficient (Kv/Kh),
c 1 unit is MPa -1 *m -2 ；
c 2 Unit is MPa -1 *m -2 ；
c 3 Unit is MPa -1 *m -2 ；
c 4 The unit is m;
calculating x 0 Value (distance between the vertex of the cylindrical projectile and the probe point) when the pumping time t is greater than t 0 At the moment, i.e. when the ellipsoid starts to be tangent to the cylindrical surface of the invaded zone, the cylindrical paraboloid (containing the formation virgin fluid) starts to taper towards the inside of the ellipsoid, i.e.:
t≥t 0
x 0 ＝Di-sd 0 (t-t 0 )
calculating the pollution rate C according to the V1, V2 and V3 obtained in the step S2,
t≥t 0
when t is less than or equal to t 0 Then, the contamination rate C (t) =1.0; t from t 0 Then the pollution rate C (t) begins to decrease; when C (t) =0.95 test elapsed time, it is breakthrough time; when C (t) =0.05 test elapsed time, it is completion time.
Further, in the step S4, the pumping result is corrected by using a numerical simulation intelligent optimization algorithm, and the fitting function is established and the simulation calculation is as follows:
and m kinds of logging data are set to participate in optimization calculation, represented by a vector a,
a＝(a 1 ,a 2 ,...,a m ) T
the vector x is m unknown parameters to be solved:
x＝(x 1 ,x 2 ,...,x m ) T
based on the volume model changes employed above, the various response equations are expressed as:
a i ＝f i (x,z),(i＝1,2,...,m)
wherein, the vector z is a regional interpretation parameter (stratum horizontal permeability Kh, anisotropy coefficient Kv/Kh and invasion depth Di), and the fitting function is as follows:
in the formula, ei and fi are the measurement error of the ith actual logging data and the error of the logging response equation respectively.
The beneficial effects of the invention are: according to the volume change of the mud pollution belt of the instrument probe accessory in the pumping process of the cable formation tester, the pumping efficiency of the instrument, the pumping breakthrough time and the pumping pure time are simulated, and meanwhile, the pumping prediction is corrected according to the detection parameters of the instrument in the pumping process, so that the purpose of utilizing the pumping operation prediction and correction is realized, the pumping and sampling operation efficiency of the cable formation tester is improved, and the operation safety of the instrument is guaranteed.
Drawings
FIG. 1 is a flow chart of a numerical correction method in an embodiment of the present invention;
FIG. 2 is a schematic illustration of a coordinate system for ellipsoidal formation fluid flow established in an embodiment of the invention;
FIG. 3 is a sampling graph showing the pump breakthrough time and pump completion time in an embodiment of the present invention;
FIG. 4 is a parameter diagram illustrating the interaction between monitoring parameters of the instrument according to an embodiment of the present invention;
FIG. 5 is a graph of parameters of real-time pump breakthrough time pick-up provided in an embodiment of the present invention;
FIG. 6 is a graph of back calculated mud invasion depth data after pump breakthrough time pickup according to an embodiment of the invention;
fig. 7 is a comparison graph of the calibration of the pumping simulation results calculated in the example of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Referring to fig. 1, a method for simulating pumping values of a cable formation test and correcting the pumping values of the cable formation test comprises the following steps:
s1, obtaining basic parameters (stratum porosity, permeability, reservoir thickness, slurry performance, borehole size, basic pumping parameters of a cable stratum test pump, slurry invasion condition and the like) of a reservoir by using conventional logging information, judging the flow pattern of the cable stratum test pressure measurement information and calculating the corresponding fluidity;
s2, opening a test hole on a well wall by a probe, enabling an invaded zone and formation fluid to start flowing to a well bore, extracting the formation fluid by a sampler, enabling the formation fluid to instantly flow and form a flowing form and continuous conduction to generate a swept range, and simulating the swept range formed by the continuous conduction of the formation fluid flow by an ellipsoid; the part formed by stratum flow outside the well cylinder and the part in the well cylinder are in a complementary relation and integrated into a complete ellipsoid, the size of the major axis and the minor axis of the ellipsoid depends on the permeability of the reservoir and the thickness of the reservoir, the volume relation among a cylinder, the ellipsoid and a cylindrical paraboloid in a volume model is obtained through calculation, the well cylinder is defined as a cylinder, the pollution zone is the ellipsoid, and the undisturbed stratum fluid is the cylindrical paraboloid; the flow form of the formation fluid is judged by combining a pressure gradient curve, a spherical flow gradient curve and a radial flow gradient curve together through analyzing the fluidity information of a pressure measuring point and the possible flow pattern characteristics of a pumping point;
and S3, inputting reservoir parameters, cable formation test operation parameters and reservoir invasion depth, and calculating the volume among a cylinder (a shaft), an ellipsoid (a pollution zone) and a cylindrical paraboloid (undisturbed formation fluid) of the model shown in the figure 2.
The equations for defining the cylindrical projectile (undisturbed formation fluid), ellipsoid (pollution zone) and cylinder (wellbore) are respectively
x-x 0 ＝a 2 (y 2 +z 2 )，
And
(x+R) 2 +y 2 ＝R 2 ，
wherein the parameter x0 is the intersection point coordinate of the vertex of the revolution paraboloid on the x axis, a is a cylindrical paraboloid system, b and c are the horizontal axis length and the vertical axis length of the ellipsoid respectively, and the R is the radius of the cylinder and the meter.
The region enclosed by the intersecting line of the cylindrical paraboloid and the ellipsoid is V 1 ：
The volume of the truncated part of the ellipsoid by the cylinder is V 2 :
Wherein the content of the first and second substances,
a 0 ＝-c,b 0 view interval = c, n [ -c, c]Length-dependent, usually 80, intrusion into the void occupied by the fluidM, i.e. subtracting V from the volume of the ellipsoid 1 And V 2 I.e. the residual volume of the ellipsoid cut by the formation virgin fluid and the wellbore fluid is set as V 3 ，
Calculating the breakthrough time and the pump pure time of pumping of the cable formation test, and defining the volume percentage of the total test fluid occupied by the fluid in the invaded zone as the pollution rate, wherein in the pumping process of the cable formation tester, the moment when the pollution rate is equal to 95% is the breakthrough time of pumping, and the moment when the pollution rate reaches 5% is the pump pure (or completion time). Input parameters for contamination rate calculation:
h is the thickness of the stratum in m;
K h is the formation horizontal permeability in md;
K v /K h is the anisotropy coefficient;
di is the invasion depth in m;
uf is the formation fluid viscosity in cp;
um is the mud filtrate viscosity in cp;
ps is the position of the probe from the top boundary, the default center position is m;
r is the radius of the shaft and is m;
qa is the pump speed in m 3 /s；
t is pump discharge time in units of s;
DP is the differential pressure in MPa;
the horizontal and vertical axes (b, c) of the ellipsoid are calculated, with b, c being in meters.
b＝c 1 ×Qa×t×DP
c＝c 2 ×Qa×t×DP
The time (t 0) for the ellipsoid to reach the invaded zone is calculated, with the unit of t0 being seconds.
The speed of motion of the rotating projectile is calculated, and sd0 is expressed in meters per second.
sd 0 ＝c 3 ×DP×Qa
The coefficients a of the ellipsoid are calculated in meters.
a＝c 4 ×K v /K h
In the formula, t 0 Meaning the time for the outer boundary to reach the invaded zone during the expansion of the ellipsoid. See the expression for b, b = Di, which shows that the outer boundary of the ellipsoid reaches the invasion zone, from which t is obtained 0 Is described in (1).
sd 0 Meaning the speed at which the rotating projectile moves in the direction of the probe.
Q a The physical meaning of x t is: reflecting the propagation velocity of ellipsoids.
C in the above equation 1 、c 2 、c 3 、c 4 The coefficients are the relationships established between formation testing parameters Differential Pressure (DP), testing flow rate (Qa), and formation anisotropy coefficient (Kv/Kh), and are determined according to the existing international system of units, and the units of the coefficients are as follows.
| Parameter(s) | c 1 | c 2 | c 3 | c 4 |
| Dimension line | MPa -1 *m -2 | MPa -1 *m -2 | MPa -1 *m -2 | m |
Calculating x 0 Value (distance between the apex of the cylindrical projectile and the probe point), current pumping time (t) is greater than t 0 At that time, i.e., the point at which the ellipsoid begins to make contact with the invaded cylindrical surface, the cylindrical paraboloid (containing the original fluid of the formation) begins to taper into the interior of the ellipsoid.
t≥t 0
x 0 ＝Di-sd 0 (t-t 0 )
V obtained according to step 3) 1 、V 2 、V 3 And calculating the pollution rate C of the waste water,
t≥t 0
when t is less than or equal to t as shown in FIG. 3 0 The contamination rate C (t) =1.0; t from t 0 The pollution rate C (t) begins to decrease after the increase; when C (t) =95% of the test elapsed time, defined as the breakthrough time; when C (t) =5% of the elapsed time of the test, the completion time is defined.
S4, during the pumping operation of the cable formation testing, comprehensively analyzing the following data of the fluid in the pumping pipeline measured by the Downhole fluid identification module of the cable formation tester by using a Downhole fluid analysis technique (DFA), with reference to fig. 4, the applied parameters include: the resistivity, temperature, pressure, density, water content, viscosity, fluorescence, bubble detection, absorption spectrum and diffraction spectrum, PH value, contamination level of the oil-based mud of the inflow/outflow line, fig. 5 is an example diagram of real-time picking up the breakthrough time of pumping according to the above method, and according to the breakthrough time of pumping, as shown in fig. 6, the intrusion radius of the intrusion zone, which is an important basis for the calibration of the pumping process, can be calculated.
Referring to fig. 7, the intrusion depth of the intrusion zone obtained in the measurement is used, and the pumping result is corrected by using a numerical simulation intelligent optimization algorithm, wherein the main algorithm (establishing a fitting function and simulation calculation) is as follows:
and m kinds of logging data are set to participate in optimization calculation, represented by a vector a,
a＝(a 1 ,a 2 ,...,a m ) T ；
the vector x is m unknown parameters to be solved:
x＝(x 1 ,x 2 ,...,x m ) T
depending on the volumetric model variations employed above, various response equations can be abstractly expressed as:
a i ＝f i (x,z),(i＝1,2,...,m)
wherein the vector z is a regional interpretation parameter (determined in the research as formation horizontal permeability Kh, anisotropy coefficient Kv/Kh and invasion depth Di). We select a fitting function of the form:
in the formula, ei and fi are the measurement error of the ith actual logging data and the error of the logging response equation respectively.
In summary, the invention observes the ratio of the slurry filtrate entering the probe in the pumping stage to the total volume of the sample, the pumping time and the change of the pumping volume according to the volume model numerical simulation and intelligent optimization method on the basis of the classification of the flow pattern of the cable formation testing fluid, and calculates the breakthrough time and the completion time of the pumping of the cable formation testing, thereby realizing the establishment of the numerical simulation model of the pumping, carrying out the simulation and correction of the pumping by comparing the numerical simulation result with the detection parameter, estimating the operation time of the pumping point in advance, providing reference for the real-time decision of the pumping sampling, and improving the operation efficiency and the sampling success rate of the pumping.
Claims (3)
1. A numerical correction method for cable formation test pumping numerical simulation and a process thereof is characterized by comprising the following steps:
s1, judging the flow pattern of the cable stratum test pressure measurement data and calculating the corresponding fluidity;
s2, opening a test hole on a well wall by a probe, enabling invasion zones and formation fluid to start flowing to a well bore, extracting the formation fluid by a sampler, enabling the formation fluid to instantly flow and form a flowing form and continuous conduction to generate a swept range, and simulating the swept range formed by the continuous conduction of the formation fluid flow by an ellipsoid; the formation part of stratum flow outside the well cylinder and the part in the well cylinder are in a complementary relationship and are integrated into a complete ellipsoid, the size of the major axis and the minor axis of the ellipsoid depends on the permeability of a reservoir and the thickness of the reservoir, the volume relationship among a cylinder, the ellipsoid and a cylindrical paraboloid in a volume model is obtained through calculation, the well cylinder is defined as a cylinder, a pollution zone is the ellipsoid, and undisturbed stratum fluid is the cylindrical paraboloid;
s3, on the basis of giving reservoir physical property parameters and invasion depth, obtaining pre-measurement predicted pumping breakthrough time, completion time and corresponding volume according to an algorithm;
s4, testing pumping monitoring parameters by a cable stratum, picking up pumping breakthrough time and volume, and correcting a pumping result by optimizing an algorithm fitting function and simulating calculation;
the algorithm of the burst time and the completion time in the step S3 is as follows:
defining the percentage of the total test fluid volume occupied by the invaded zone fluid is called the contamination rate, during pumping by the cable formation tester, the breakthrough time for pumping is when the contamination rate =95%, and the completion time is when the contamination rate = 5%:
calculating the horizontal and vertical axes (b, c) of the ellipsoid, the units of b, c being m,
b＝c 1 ×Qa×t×DP
c＝c 2 ×Qa×t×DP
calculating the time of arrival of ellipsoid at the invaded zoneTime t 0 ，t 0 The unit of (a) is s,
b＝Di
calculating the speed of movement of the rotating projectile, sd 0 The unit of (a) is m/s,
sd 0 ＝c 3 ×DP×Qa
calculating the coefficient a of the ellipsoid, wherein the unit of a is m,
a＝c 4 ×K v /K h
in the above formula, the first and second carbon atoms are,
h is the thickness of the stratum in m;
K h is the formation horizontal permeability in md;
K v /K h is the anisotropy coefficient;
di is the invasion depth in m;
uf is the formation fluid viscosity in cp;
um is the mud filtrate viscosity in cp;
ps is the position of the probe from the top boundary, and the default center position is m;
r is the radius of the shaft and is m;
qa is the pump speed in m 3 /s；
t is pump discharge time in units of s;
DP is the differential pressure in MPa;
calculating x 0 The distance between the top of the cylindrical projectile and the probe point is obtained, and when the pumping time t is more than t 0 When the ellipsoid starts to be tangent with the cylindrical surface of the invasion zone, the cylindrical paraboloid contains the original fluid of the stratum and starts to be pushed towards the inside of the ellipsoid, namely:
t≥t 0
x 0 ＝Di-sd 0 (t-t 0 )
calculating the pollution rate C according to the V1, V2 and V3 obtained in the step S2,
t≥t 0
when t is less than or equal to t 0 Then, the contamination rate C (t) =1.0; t from t 0 The pollution rate C (t) begins to decrease after the increase; when C (t) =0.95, it is a breakthrough time; when C (t) =0.05, completion time is.
2. A method for simulating pumping values of a cable formation testing system and correcting the pumping values according to claim 1, wherein the volume relationship among the cylinder, ellipsoid and cylindrical paraboloid of the model calculated in step S2 is:
x-x 0 ＝a 2 (y 2 +z 2 )，
(x+R) 2 +y 2 ＝R 2
in the formula, x 0 The coordinate of the intersection point of the vertex of the revolution paraboloid on the x axis is shown;
a is a cylindrical paraboloid system,
b is the length of the horizontal axis of the ellipsoid in meters;
c is the length of the vertical axis of the ellipsoid in meters;
r is the radius of the cylinder and is measured in meters;
the region enclosed by the intersecting line of the cylindrical paraboloid and the ellipsoid is V 1 ：
The volume of the ellipsoid cut by the cylinder is V 2 :
Wherein the content of the first and second substances,
a 0 ＝-c,b 0 ＝c，
n is an interval [ -c, c ],
intrusion into the space occupied by the fluid by the volume of the ellipsoid minus V 1 And V 2 Definition of V 3 Is the remaining volume of ellipsoid cut by formation virgin fluid and wellbore fluid:
3. the method for numerical calibration of cable formation testing pumping numerical simulation and process thereof according to claim 1, wherein the step S4 utilizes a numerical simulation intelligent optimization algorithm to calibrate the pumping result, which establishes a fitting function and the simulation calculation as follows:
and m kinds of logging data are taken into the optimization calculation, represented by a vector a',
a′＝(a 1 ,a 2 ,...,a m ) T
the vector x' is m unknown parameters to be solved:
x′＝(x 1 ,x 2 ,...,x m ) T
depending on the volume model change employed, various response equations are expressed as:
a i ＝f i (x′,z),i＝1,2,…,m
wherein the vector z is a regional interpretation parameter: the fitting function of the horizontal permeability Kh, the anisotropy coefficient Kv/Kh and the invasion depth Di of the stratum is as follows:
in the formula, e i 、f i The measurement error of the ith actual logging data and the error of the logging response equation are respectively.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910382717.XA CN110162851B (en) | 2019-05-09 | 2019-05-09 | Cable formation test pumping numerical simulation and numerical correction method of process thereof |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910382717.XA CN110162851B (en) | 2019-05-09 | 2019-05-09 | Cable formation test pumping numerical simulation and numerical correction method of process thereof |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110162851A CN110162851A (en) | 2019-08-23 |
| CN110162851B true CN110162851B (en) | 2022-11-22 |
Family
ID=67633875
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910382717.XA Active CN110162851B (en) | 2019-05-09 | 2019-05-09 | Cable formation test pumping numerical simulation and numerical correction method of process thereof |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110162851B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110984970B (en) * | 2019-10-09 | 2023-03-24 | 中国海洋石油集团有限公司 | Method for determining starting pressure gradient by utilizing formation test |
| CN110905493B (en) * | 2019-11-21 | 2023-05-12 | 中国海洋石油集团有限公司 | Method for measuring pollution rate of underground stratum fluid |
| CN112147051B (en) * | 2019-12-25 | 2023-04-25 | 中国海洋石油集团有限公司 | Pressure measurement fluidity standardization method based on permeability distribution form |
| CN114382466B (en) * | 2021-12-28 | 2023-06-06 | 中海油田服务股份有限公司 | Method and device for optimizing test parameters |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5337821A (en) * | 1991-01-17 | 1994-08-16 | Aqrit Industries Ltd. | Method and apparatus for the determination of formation fluid flow rates and reservoir deliverability |
| CN1588347A (en) * | 2004-07-13 | 2005-03-02 | 中海油田服务股份有限公司 | Method for detecting analogue cable formation tester |
| CN101906966A (en) * | 2010-07-16 | 2010-12-08 | 中国石油天然气股份有限公司 | Method and device for forecasting reservoir yield |
| CN102073772A (en) * | 2011-01-18 | 2011-05-25 | 大庆油田有限责任公司 | MDT pressure value test simulation method |
| CN107605477A (en) * | 2017-08-17 | 2018-01-19 | 中国海洋石油总公司 | It is a kind of to determine that underground pumping sampling oil gas breaks through the method with moisture content stabilization time |
| CN108982319A (en) * | 2018-07-10 | 2018-12-11 | 中国海洋石油集团有限公司 | A kind of acquisition methods of oil field stratum condition phase percolation curve |
-
2019
- 2019-05-09 CN CN201910382717.XA patent/CN110162851B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5337821A (en) * | 1991-01-17 | 1994-08-16 | Aqrit Industries Ltd. | Method and apparatus for the determination of formation fluid flow rates and reservoir deliverability |
| CN1588347A (en) * | 2004-07-13 | 2005-03-02 | 中海油田服务股份有限公司 | Method for detecting analogue cable formation tester |
| CN101906966A (en) * | 2010-07-16 | 2010-12-08 | 中国石油天然气股份有限公司 | Method and device for forecasting reservoir yield |
| CN102073772A (en) * | 2011-01-18 | 2011-05-25 | 大庆油田有限责任公司 | MDT pressure value test simulation method |
| CN107605477A (en) * | 2017-08-17 | 2018-01-19 | 中国海洋石油总公司 | It is a kind of to determine that underground pumping sampling oil gas breaks through the method with moisture content stabilization time |
| CN108982319A (en) * | 2018-07-10 | 2018-12-11 | 中国海洋石油集团有限公司 | A kind of acquisition methods of oil field stratum condition phase percolation curve |
Non-Patent Citations (2)
| Title |
|---|
| 泵抽式电缆地层测试储层产能预测方法研究;尹邦堂等;《断块油气田》;20110925(第05期);全文 * |
| 泵抽式电缆地层测试复合流形模型研究及应用;张国海等;《中外能源》;20120915(第09期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110162851A (en) | 2019-08-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110162851B (en) | Cable formation test pumping numerical simulation and numerical correction method of process thereof | |
| CN107701180B (en) | Original oil reservoir water saturation calculation method based on closed coring | |
| EP0698722B1 (en) | Method for testing low permeability formations | |
| CN108694264B (en) | Method for determining permeability of shale gas reservoir | |
| CN110043254B (en) | Method for obtaining stratum effective permeability based on cable stratum test data | |
| CN107462936B (en) | Utilize the method for pressure monitoring Data Inversion low permeability reservoir non-Darcy percolation law | |
| US8020437B2 (en) | Method and apparatus to quantify fluid sample quality | |
| WO2006120366A1 (en) | Methods for analysis of pressure response in underground formations | |
| CN109736788A (en) | A kind of experimental method for differentiating chemical flooding leading edge and involving state | |
| CN108825226A (en) | A kind of method and device for assessing gas production after pressure using chemical tracer | |
| CN110552694B (en) | Argillaceous dolomite oil reservoir oil well productivity evaluation method considering multi-factor influence | |
| RU2577865C1 (en) | Method of indicating investigation of wells and interwell space | |
| CN112036097A (en) | Capacity calculation method for water-lock gas well | |
| US11852576B2 (en) | Method for determining the pore size distribution in a reservoir | |
| CN112647930B (en) | Water finding method for oil reservoir engineering of horizontal well | |
| CN114922614A (en) | Formation pressure monitoring method under pressure control drilling working condition | |
| CN112014881B (en) | Water flooding speed prediction method based on time lapse earthquake | |
| Carnegie | Understanding the pressure gradients improves production from oil/water transition carbonate zones | |
| CN110764161B (en) | Comprehensive method for judging and identifying fractured and broken zone of carbonate rock through drilling data | |
| CN108150158B (en) | Deep fractured compact sandstone gas reservoir early water body analysis and prediction method | |
| CN110619490A (en) | Method for evaluating cementing sealing performance of fractured and crushed compact carbonate rock zone | |
| CN113107464B (en) | Horizontal well stepping type water flooded layer identification logging method | |
| RU2768341C1 (en) | Well production rate prediction method taking into account anisotropy of permeability of carbonate rocks | |
| CN113589378B (en) | Fault plugging property evaluation method based on three-dimensional seismic data | |
| CN112147051B (en) | Pressure measurement fluidity standardization method based on permeability distribution form |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |