CN113673101A - Method and processor for calculating volume of beaded fracture-cavity karst cave - Google Patents
Method and processor for calculating volume of beaded fracture-cavity karst cave Download PDFInfo
- Publication number
- CN113673101A CN113673101A CN202110949690.5A CN202110949690A CN113673101A CN 113673101 A CN113673101 A CN 113673101A CN 202110949690 A CN202110949690 A CN 202110949690A CN 113673101 A CN113673101 A CN 113673101A
- Authority
- CN
- China
- Prior art keywords
- fracture
- beaded
- cavity
- pressure
- coefficient
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000003860 storage Methods 0.000 claims abstract description 47
- 230000006835 compression Effects 0.000 claims description 23
- 238000007906 compression Methods 0.000 claims description 23
- 239000010779 crude oil Substances 0.000 claims description 14
- 239000012530 fluid Substances 0.000 claims description 14
- 239000011148 porous material Substances 0.000 claims description 14
- 210000003491 skin Anatomy 0.000 claims description 14
- 239000003129 oil well Substances 0.000 claims description 13
- 238000004519 manufacturing process Methods 0.000 claims description 10
- 230000035699 permeability Effects 0.000 claims description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 230000005465 channeling Effects 0.000 claims description 7
- 239000000126 substance Substances 0.000 claims description 6
- 238000011156 evaluation Methods 0.000 claims description 4
- 210000002615 epidermis Anatomy 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 10
- 238000011161 development Methods 0.000 abstract description 7
- 230000000694 effects Effects 0.000 abstract description 5
- 230000008901 benefit Effects 0.000 abstract description 2
- 238000009472 formulation Methods 0.000 abstract 1
- 239000000203 mixture Substances 0.000 abstract 1
- 239000003921 oil Substances 0.000 description 16
- 238000012360 testing method Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 9
- 238000004590 computer program Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 5
- 238000004891 communication Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 235000019994 cava Nutrition 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
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)
- Measuring Fluid Pressure (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to a method and a processor for calculating the volume of a beaded fracture-cavity karst cave, wherein the method comprises the following steps: acquiring actual pressure data inside the beaded fracture-cavity body, a first initial value of a karst cave storage coefficient of the beaded fracture-cavity body and a first fracture-cavity body bottom pressure solution model; inputting a first initial value into a first fracture-cavity body bottom pressure solution model to obtain a first bottom pressure solution, and inputting the first bottom pressure solution into a second fracture-cavity body bottom pressure solution model to obtain theoretical pressure data of the beaded fracture-cavity body; determining a solved value of a karst cave storage coefficient of the beaded fracture-cavity body according to actual pressure data and theoretical pressure data; the karst cave volume in the beaded seam body is determined according to the solved value, the method has the advantages of simplicity, understandability and simplicity in operation, the calculation result of the karst cave volume of the beaded seam body is more in line with the actual situation, and the method has an important technical support effect on formulation of a reasonable oil field development scheme and adjustment measures.
Description
Technical Field
The invention relates to the field of oil reservoir engineering, in particular to a method for calculating the volume of a beaded fracture-cavity karst cave, a processor and a machine-readable storage medium.
Background
The oil gas external dependence degree of China climbs year by year, and once exceeds the national energy external dependence degree safety line. In recent years, in response to national calls, oil field companies have increased investment, increased production and stocked, and several hundred million-ton carbonate oil fields have been successively discovered in Tarim basins. Such newly discovered fields exhibit some characteristics that are different from conventional carbonate reservoirs: the oil reservoir is formed by large-scale structure fracture activity and multi-phase karst action, the buried depth of the reservoir generally exceeds 6km, 4 media mainly develop and comprise a large karst cave, small holes, structural fractures and rock matrixes, and oil gas is not stored in the matrixes generally but stored in the karst cave and natural fractures due to the large buried depth of the reservoir; the communication relation between the karst caves and the natural fractures is very complex, wherein the bead-shaped fracture-cavity body is the most common fracture-cavity communication mode in the oil reservoir; production wells are typically deployed around the beaded fracture body at development time to obtain industrial production. Therefore, reasonably evaluating parameters of the beaded fracture-cavity body so as to improve the reservoir development effect becomes a hot problem in carbonate reservoir development.
Well test analysis is the most commonly used technology for inverting reservoir parameters at present, however, the existing well test analysis method is based on the classical Darcy seepage theory and cannot accurately describe the complex fluid flow behavior in a large karst cave; in addition, the existing well testing analysis method only considers 3 media of small holes, structural cracks and rock matrixes, and has deviation from the actual condition of the oil field; more importantly, the existing well testing analysis method does not consider the specific communication mode of the beaded fracture-hole body. If the newly discovered carbonate reservoir is explained by adopting a traditional well testing analysis method, the result of explained parameters (such as permeability, karst cave radius, karst cave volume and the like) is distorted, and the development effect of the reservoir is influenced.
Disclosure of Invention
The invention aims to provide a method, a processor and a machine-readable storage medium for calculating the volume of the beaded cave body cavern, and the method, the processor and the machine-readable storage medium have the advantages of simplicity and understandability, simple operation and capability of enabling the calculation result of the volume of the beaded cave body cavern to better meet the actual situation.
In order to achieve the above object, a first aspect of the present invention provides a method for calculating a karst cave volume of a beaded fracture-cave body, the method comprising:
acquiring actual pressure data inside the beaded fracture-cavity body, a first initial value of a karst cave storage coefficient of the beaded fracture-cavity body and a first fracture-cavity body bottom pressure solution model;
inputting the first initial value into a first fracture-cavity bottom pressure solution model to obtain theoretical pressure data of the beaded fracture-cavity body;
determining a solved value of a karst cave storage coefficient of the beaded fracture-cavity body according to actual pressure data and theoretical pressure data;
and determining the volume of the karst cave in the beaded seam-hole body according to the solved value.
In an embodiment of the present invention, inputting a first initial value into a first fracture-cavity bottom pressure solution model to obtain theoretical pressure data of a beaded fracture-cavity comprises:
inputting a first initial value into a first fracture-cave body bottom pressure solution model to obtain a first bottom pressure solution, wherein the first bottom pressure solution means that dimensionless bottom pressure of a beaded fracture-cave body in a Laplace space does not consider a well reservoir skin;
acquiring a second fracture-cavity bottom pressure solution model;
inputting the first bottom hole pressure solution into a second fracture-cavity bottom pressure solution model to obtain a second bottom hole pressure solution, wherein the second bottom hole pressure solution refers to the pressure drop of a beaded fracture-cavity body in a Laplace space under consideration of a well reservoir skin;
and obtaining theoretical pressure data of the beaded fracture-hole body according to the second bottom hole pressure solution.
In an embodiment of the invention, inputting the first initial value into the first fracture-cavity bottom-hole pressure solution model to obtain the first bottom-hole pressure solution comprises:
acquiring a karst cave equivalent radius, a fracture storage capacity ratio, a karst cave channeling coefficient, a fracture medium porosity, a fracture medium comprehensive compression coefficient, a pore medium porosity and a pore medium comprehensive compression coefficient;
calculating a first bottom hole pressure solution according to the following equation:
wherein the content of the first and second substances,
a first bottom hole pressure solution; s is a laplace variable; omega is a fracture storage-capacity ratio; lambda is a karst cave channeling coefficient;
porosity of the fracture medium; c. CtfThe comprehensive compression coefficient of the crack medium is obtained;
porosity of the porous medium; c. CtvThe comprehensive compression coefficient of the pore medium is; cvIs a first initial value; r isfIs the equivalent radius of the karst cave.
In an embodiment of the invention, inputting the first bottom hole pressure solution into a second fracture cavity bottom hole pressure solution model to obtain a second bottom hole pressure solution comprises:
acquiring the viscosity of crude oil, the formation permeability, the volume coefficient of the crude oil, the yield of an oil well, the storage coefficient of a shaft and the skin coefficient;
calculating a second bottom hole pressure solution according to the following equation:
wherein the content of the first and second substances,
a second bottom hole pressure solution; is the crude oil viscosity; b is the volume coefficient of crude oil; q is the oil well production; k is a radical offIs the formation permeability; s is the epidermis coefficient; c is the wellbore reservoir coefficient.
In the embodiment of the invention, determining the solved value of the beaded fracture-cavity cavern coefficient according to the actual pressure data and the theoretical pressure data comprises the following steps:
drawing an actual pressure double-logarithmic curve according to the actual pressure data, wherein the actual pressure double-logarithmic curve comprises an actual pressure value curve and an actual pressure value derivative curve;
drawing a theoretical pressure double-logarithmic curve according to theoretical pressure data, wherein the theoretical pressure double-logarithmic curve comprises a theoretical pressure value curve and a theoretical pressure value derivative curve;
and determining the calculated value of the storage coefficient of the beaded fracture-cavity karst cave according to the actual pressure double-logarithmic curve and the theoretical pressure double-logarithmic curve.
In the embodiment of the invention, determining the solved value of the karst cave storage coefficient of the beaded fracture-cave body according to the actual pressure double-logarithmic curve and the theoretical pressure double-logarithmic curve comprises the following steps:
calculating a first coincidence rate between the actual pressure value curve and the theoretical pressure value curve and a second coincidence rate between the actual pressure value derivative curve and the theoretical pressure value derivative curve;
and determining the calculated value of the karst cave storage coefficient of the beaded fracture-cave body according to the first coincidence rate and the second coincidence rate.
In an embodiment of the present invention, determining a derived value of a beaded-shaped fracture-cavity cavern reservoir coefficient as a function of the first and second rates of coincidence comprises:
comparing the first coincidence rate with a first preset coincidence rate range to judge whether the first coincidence rate is within the first preset coincidence rate range;
comparing the second coincidence rate with a second preset coincidence rate range to judge whether the second coincidence rate is within the second preset coincidence rate range;
and determining the solved value of the karst cave storage coefficient as a first initial value under the condition that the first coincidence rate is within a first preset coincidence rate range and the second coincidence rate is within a second preset coincidence rate range.
In an embodiment of the present invention, determining the cavern volume in the beaded fracture body from the solved value comprises:
acquiring a fluid compression coefficient;
the cavern volume is calculated according to the following formula:
wherein V is the volume of the karst cave; c. CfIs the fluid compressibility.
A second aspect of the invention provides a processor configured to perform the above-described method of calculating a cavern volume of a beaded fracture body.
A third aspect of the invention provides a machine-readable storage medium having instructions stored thereon for causing a machine to perform the above-described method of calculating a cavern volume of a beaded fracture cavity.
By the technical scheme, actual pressure data inside the beaded fracture-cavity body, a first initial value of the karst cave storage coefficient of the beaded fracture-cavity body and a first fracture-cavity body well bottom pressure solution model are obtained; inputting the first initial value into a first fracture-cavity bottom pressure solution model to obtain theoretical pressure data of the beaded fracture-cavity body; determining a solved value of a karst cave storage coefficient of the beaded fracture-cavity body according to actual pressure data and theoretical pressure data; the method is simple and easy to understand and is simple to operate, the well bottom pressure solution model of the first fracture-cavity body is more suitable for the actual conditions of the oil field, interpretation efficiency can be greatly improved when the method is applied to the oil field, and oil field production is facilitated.
Drawings
The accompanying drawings, which are included to provide a further understanding of the embodiments of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the embodiments of the invention without limiting the embodiments of the invention. In the drawings:
FIG. 1 is a schematic flow chart of a method for calculating the volume of a beaded-shaped cavity solution cavern according to an embodiment of the invention;
FIG. 2 is a schematic view of a physical model of a beaded slot body in an embodiment of the present invention;
FIG. 3 is a logarithmic graph of the derivative of the theoretical pressure value and the theoretical pressure value according to an embodiment of the present invention;
FIG. 4 is a diagram illustrating a fitting result between a theoretical log-log curve and an actual log-log curve according to an embodiment of the present invention.
Description of the reference numerals
1 theoretical pressure value curve 2 theoretical pressure value derivative curve
3 actual pressure value curve 4 actual pressure value derivative curve
Detailed Description
The following detailed description of embodiments of the invention refers to the accompanying drawings. It should be understood that the detailed description and specific examples, while indicating embodiments of the invention, are given by way of illustration and explanation only, not limitation.
Note that, if directional indications (such as up, down, left, right, front, and rear … …) are referred to in the embodiments of the present application, the directional indications are merely used to explain the relative positional relationship between the components, the movement situation, and the like in a specific posture, and if the specific posture is changed, the directional indications are changed accordingly.
In addition, if there is a description of "first", "second", etc. in the embodiments of the present application, the description of "first", "second", etc. is for descriptive purposes only and is not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one such feature. In addition, technical solutions between the various embodiments can be combined with each other, but must be realized by a person skilled in the art, and when the technical solutions are contradictory or cannot be realized, such a combination should not be considered to exist, and is not within the protection scope of the present application.
In an embodiment of the present invention, a novel method for calculating a volume of a beaded cave body cavern is provided, which is suitable for a control device for calculating a volume of a beaded cave body cavern, and specifically, as shown in fig. 1, the method includes the following steps:
step S101: actual pressure data inside the beaded fracture-cavity body, a first initial value of a karst cave storage coefficient of the beaded fracture-cavity body and a first fracture-cavity body bottom pressure solution model are obtained.
It is understood that the control device for calculating the volume of the karst cave of the beaded shape caves in the embodiment comprises a pressure gauge and a processor, wherein the pressure gauge is electrically connected with the processor. Before calculating the volume of the beaded fracture-cavity karst cave, an operator needs to select a beaded fracture-cavity oil well (in the embodiment, the selected oil well is an oil well of a northeast oil field of a Tarim basin, and the well number of the oil well is SHB A), and after closing the well, the pressure recovery test needs to be performed on the oil well, wherein the pressure recovery test at least comprises one of a bottom hole pressure test or a wellhead pressure test, wherein the bottom hole pressure test refers to that a pressure gauge is put right above a production zone of the oil well, and the pressure gauge and upper fluid are sealed by a packer so as to monitor the change relationship of the bottom hole pressure along with time after closing the well in real time; the wellhead pressure test means that the pressure gauge is placed at a wellhead so as to monitor the time-varying relation of wellhead pressure after well shut-in real time, in addition, data obtained by the wellhead pressure test in the embodiment need to be converted into bottom pressure and then subsequent fitting calculation can be carried out, and influence factors such as shaft friction, liquid level height and the like need to be considered in the conversion process so as to reduce the converted pressure error.
And the pressure gauge transmits the actual pressure value to the processor after detection is finished, and the processor acquires actual pressure data in the fracture body based on the actual pressure value. The actual pressure data in this embodiment includes an actual pressure value and an actual pressure value derivative, and the processor performs calculation processing on the actual pressure value to obtain the actual pressure value derivative.
In an embodiment of the present invention, the control device further comprises a memory electrically connected to the processor, wherein the first initial value of the beaded-shaped fracture-cavity karst-cave reservoir coefficient and the first fracture-cavity bottom-of-well pressure solution model are pre-stored in the memory, and the processor can retrieve the initial value and the first fracture-cavity bottom-of-well pressure solution model when in use.
In the embodiment, the processor appropriately simplifies the actual beaded fracture body according to the specific communication mode of the beaded fracture body to establish a physical model in which the karst cave region and the fracture region are connected in series, and after the physical model is established, the processor establishes a first fracture-cave body well bottom pressure solution model based on the physical model and stores the first fracture-cave body well bottom pressure solution model in the memory. As shown in fig. 2, the cavern region in the physical model only contains the cavern medium, and the fluid flow of the cavern medium complies with the free flow (i.e., the flow in the cavern region is free flow); the fracture zone contains 3 media of pores, formation fractures and rock matrix, the fluid flow of which obeys darcy's law (i.e. the flow in the fracture zone is darcy flow). Because the physical model considers the specific communication mode of 4 media and the beaded fracture-cavity body and couples the free flow and the Darcy flow, compared with the traditional Darcy flow theory, the physical model adopts a more advanced theory to establish the first fracture-cavity body bottom pressure solution model, so that the established first fracture-cavity body bottom pressure solution model is more in line with the actual situation.
Step S102: inputting the first initial value into the first fracture-cavity bottom pressure solution model to obtain theoretical pressure data of the beaded fracture-cavity body.
In an embodiment of the present invention, the step S102 of inputting the first initial value into the first fracture-cavity bottom pressure solution model to obtain the theoretical pressure data of the beaded fracture-cavity includes the following steps:
step S201: inputting the first initial value into the first fracture-cavern body bottom-hole pressure solution model to obtain a first bottom-hole pressure solution, wherein the first bottom-hole pressure solution refers to a dimensionless bottom-hole pressure of the beaded fracture-cavern body without considering a well reservoir skin under Laplace space.
It is to be understood that the step S201 of inputting the first initial value into the first fracture-cavity bottom pressure solution model to obtain a first bottom pressure solution includes the following steps:
step S301: acquiring a karst cave equivalent radius, a fracture storage capacity ratio, a karst cave channeling coefficient, a fracture medium porosity, a fracture medium comprehensive compression coefficient, a pore medium porosity and a pore medium comprehensive compression coefficient;
step S302: calculating a first bottom hole pressure solution according to the following equation:
wherein the content of the first and second substances,
a first bottom hole pressure solution; s is a laplace variable; omega is a fracture storage-capacity ratio; lambda is a karst cave channeling coefficient;
porosity of the fracture medium; c. CtfThe comprehensive compression coefficient of the crack medium is obtained;
porosity of the porous medium; c. CtvThe comprehensive compression coefficient of the pore medium is; cvIs a first initial value in m3/Pa；rfIs the equivalent radius of the karst cave, and the unit is m.
In the embodiment of the invention, the equivalent radius of the karst cave, the fracture volume ratio, the blow-by coefficient of the karst cave, the porosity of the fracture medium, the comprehensive compression coefficient of the fracture medium, the porosity of the pore medium and the comprehensive compression coefficient of the pore medium are all stored in a memory in advance, and the processor can call the equivalent radius, the fracture volume ratio, the blow-by coefficient of the karst cave, the porosity of the pore medium and the comprehensive compression coefficient of the pore medium when in application. Further, the fracture medium porosity and poresThe sum of the porosity of the hole medium is the average porosity of the stratum (in this embodiment, the average porosity of the stratum is 0.1, i.e., + -. phi)f+φm0.1); furthermore, ctf＝ctm＝8.36×10-4MPa-1。
Before calculating the first bottom hole pressure solution according to the formula (1), the processor needs to assign a first initial value to the cavern body karst cave storage coefficient, and needs to assign the first initial value according to geological data such as well logging, well completion, earthquake, core and the like, so that the calculation result of the formula (1) is more consistent with the actual situation of the stratum, wherein the first initial value is 60m in the embodiment3(ii) MPa; the first initial value assignment is substituted into equation (1) after it is completed to calculate a first bottom hole pressure solution (i.e., the dimensionless bottom hole pressure of the beaded fracture body in laplace space without considering the reservoir skin of the well).
Step S202: and obtaining a second fracture-cavity bottom pressure solution model.
In an embodiment of the invention, the second fracture-cavity bottom-of-well pressure solution model is pre-stored in the memory and can be retrieved by the processor when applied.
Step S203: inputting the first bottom hole pressure solution into a second fracture-cavity bottom pressure solution model to obtain a second bottom hole pressure solution, wherein the second bottom hole pressure solution refers to a pressure drop of the beaded fracture body under the Laplace space with the well reservoir skin taken into consideration.
In one embodiment of the present invention, the step S203 of inputting the first bottom hole pressure solution into the second fracture cavity bottom hole pressure solution model to obtain the second bottom hole pressure solution comprises the following steps:
step S401: acquiring the viscosity of crude oil, the formation permeability, the volume coefficient of the crude oil, the yield of an oil well, the storage coefficient of a shaft and the skin coefficient;
step S402: calculating a second bottom hole pressure solution according to the following equation:
wherein the content of the first and second substances,
is the second bottom hole pressure solution in Pa; mu is the viscosity of the crude oil in mPa s; b is the volume coefficient of crude oil; q is the oil well production in m3/d；kfIs the formation permeability in m2(ii) a S is the epidermis coefficient; c is the wellbore storage coefficient in m3/Pa。
Further, in the present embodiment, the crude oil viscosity, the formation permeability, the crude oil volume coefficient, the oil well production, the wellbore storage coefficient, and the skin coefficient are all stored in the memory in advance, and μ ═ 2.242mPa · s; b-1.144; q is 70m3And d. The processor retrieves the parameters from the memory and then substitutes them into equation (2) to calculate a second bottom hole pressure solution.
Step S204: and obtaining theoretical pressure data of the beaded fracture-hole body according to the second bottom hole pressure solution.
It can be understood that the theoretical pressure data in this embodiment includes a theoretical pressure value and a theoretical pressure value derivative, the second bottom hole pressure solution obtained by the processor is the theoretical pressure value (that is, the pressure drop of the reservoir skin of the well is considered in the laplace space by the beaded fracture body), after the processor obtains the theoretical pressure value, the processor performs a stepfest numerical inversion on the theoretical pressure value to further obtain data of the change of the bottom hole pressure in the real space with time, and then obtains the theoretical pressure value derivative by calculating the logarithmic time derivative of the data.
Step S103: and determining the solved value of the karst cave storage coefficient of the beaded fracture-cavity body according to the actual pressure data and the theoretical pressure data.
In one embodiment of the present invention, the step S103 of determining the solved value of the karst cave reservoir coefficient of the beaded fracture-cavity based on the actual pressure data and the theoretical pressure data comprises the following steps:
step S501: drawing an actual pressure double-logarithmic curve according to the actual pressure data, wherein the actual pressure double-logarithmic curve comprises an actual pressure value curve 3 and an actual pressure value derivative curve 4;
step S502: drawing a theoretical pressure double-logarithmic curve according to theoretical pressure data, wherein the theoretical pressure double-logarithmic curve comprises a theoretical pressure value curve 1 and a theoretical pressure value derivative curve 2;
step S503: and determining the calculated value of the storage coefficient of the beaded fracture-cavity karst cave according to the actual pressure double-logarithmic curve and the theoretical pressure double-logarithmic curve.
It can be understood that after the processor obtains the actual pressure value, the derivative of the actual pressure value to logarithmic time is calculated to obtain the derivative of the actual pressure value, and then the processor simultaneously draws the actual pressure value and the derivative of the actual pressure value on a double-logarithm chart to obtain a double-logarithm actual curve chart of the beaded fracture-shaped body; after obtaining the theoretical pressure value and the theoretical pressure value derivative, the processor simultaneously draws the theoretical pressure value and the theoretical pressure value derivative on a double-logarithm chart to obtain a double-logarithm theoretical curve chart of the beaded fracture-cavity body, wherein the double-logarithm theoretical curve chart of the beaded fracture-cavity body is shown in FIG. 3 (in the chart, a vertical axis represents the pressure value and the pressure value derivative, and a horizontal axis represents time); it can be understood that, after the processor acquires the actual pressure log-log curve and the theoretical pressure log-log curve, the two pairs of curves are fitted, as shown in fig. 4 (in the figure, the vertical axis represents the pressure value and the pressure value derivative, and the horizontal axis represents time), that is, the actual pressure value curve 3 and the theoretical pressure value curve 1 are fitted, and the actual pressure value derivative curve 4 and the theoretical pressure value derivative curve 2 are fitted, so as to determine the evaluation value of the karst cave storage coefficient of the karst cave body according to the fitting result, wherein the fitting algorithm includes at least one of a least square method or a genetic algorithm.
In an embodiment of the present invention, the step S503 of determining the solved value of the karst cave storage coefficient of the beaded fracture-cavity body according to the actual pressure log-log curve and the theoretical pressure log-log curve includes the following steps:
step S601: calculating a first coincidence rate between the actual pressure value curve 3 and the theoretical pressure value curve 1 and a second coincidence rate between the actual pressure value derivative curve 4 and the theoretical pressure value derivative curve 2;
step S602: and determining the calculated value of the karst cave storage coefficient of the beaded fracture-cave body according to the first coincidence rate and the second coincidence rate.
It can be understood that after the processor fits the actual pressure log-log curve and the theoretical pressure log-log curve, the coincidence rate between the actual pressure log-log curve and the theoretical pressure log-log curve is calculated, that is, the processor calculates a first coincidence rate after the actual pressure value curve 3 and the theoretical pressure value curve 1 are fitted and a second coincidence rate after the actual pressure value derivative curve 4 and the theoretical pressure value derivative curve 2 are fitted, so as to determine the evaluation value of the reservoir parameters of the fracture-cave reservoir according to the first coincidence rate and the second coincidence rate.
In one embodiment of the present invention, the step S602 of determining a calculated value of the string-bead-like cavity cavern reservoir coefficient according to the first coincidence rate and the second coincidence rate includes the following steps:
step S701: comparing the first coincidence rate with a first preset coincidence rate range to judge whether the first coincidence rate is within the first preset coincidence rate range;
step S702: comparing the second coincidence rate with a second preset coincidence rate range to judge whether the second coincidence rate is within the second preset coincidence rate range;
step S703: and determining the solved value of the string-bead-shaped slot cavity karst cave storage coefficient as a first initial value under the condition that the first coincidence rate is within a first preset coincidence rate range and the second coincidence rate is within a second preset coincidence rate range.
It can be understood that after obtaining the first coincidence rate and the second coincidence rate, the processor needs to determine whether the first coincidence rate is within a first preset coincidence rate range and whether the second coincidence rate is within a second preset coincidence rate range, where the first preset coincidence rate range is 0.7-1; the second preset coincidence rate range is 0.7-1. If the first coincidence rate is within a first preset coincidence rate range and the second coincidence rate is within a second preset coincidence rate range, the actual pressure log-log curve and the theoretical pressure log-log curve reach the standard of successful fitting, and then the solved value of the fracture-cave reservoir layer parameter is determined to be a first initial value; if the first coincidence rate is not within the first preset coincidence rate range, and/or the second coincidence rate is notAnd in a second preset coincidence rate range, the actual pressure log-log curve and the theoretical pressure log-log curve do not reach the standard of successful fitting, and further the obtained value of the fracture-cavity reservoir parameter is determined to be not the first initial value, at the moment, the processor needs to re-assign the first initial value of the fracture-cavity reservoir parameter, and continues to execute subsequent steps after assigning until the actual pressure log-log curve and the theoretical pressure log-log curve reach the standard of successful fitting (namely the first coincidence rate is in the first preset coincidence rate range, and the second coincidence rate is in the second preset coincidence rate range), and further the obtained value of the fracture-cavity body karst-cavity reservoir coefficient is determined to be the first initial value. In an embodiment of the invention, the first initial value is reassigned to 62.8m3And at the time of/MPa, the actual pressure double-logarithmic curve and the theoretical pressure double-logarithmic curve reach the standard of successful fitting.
Step S104: and determining the volume of the karst cave in the fracture and cave body according to the obtained value.
In one embodiment of the present invention, the step S104 of determining the cavern volume in the beaded fracture body according to the calculation value includes the following steps:
step S801: acquiring a fluid compression coefficient;
step S802: the cavern volume is calculated according to the following formula:
wherein V is the volume of the karst cave and the unit is m3；cfIs the coefficient of compression of the fluid, in Pa-1。
It can be understood that the fluid compression coefficient is stored in the memory in advance, the processor retrieves the fluid compression coefficient from the memory after determining the solved value of the cavern body cavern coefficient, and then the fluid compression coefficient and the fluid compression coefficient are substituted into the formula (3) to calculate the cavern volume in the beaded cavern body.
Furthermore, in the embodiment of the invention, the calculated value of the beaded cave body karst cave reservoir parameters can be determined according to actual pressure data and theoretical pressure data, wherein the beaded cave body karst cave reservoir parameters at least comprise a shaft reservoir coefficient, a skin coefficient, a formation permeability, a karst cave equivalent radius, a fracture volume ratio and a karst cave channeling coefficient.
Specifically, the parameters are given initial values according to geological data such as well logging, well completion, earthquake, core and the like, wherein the initial value given to the well bore storage coefficient is a second initial value, the initial value given to the skin coefficient is a third initial value, the initial value given to the stratum permeability is a fourth initial value, the initial value given to the karst cave equivalent radius is a fifth initial value, the initial value given to the fracture volume ratio is a sixth initial value, the initial values are stored in a memory in advance, and the second initial value is 6m3The second initial value is 0.1, the third initial value is 100mD, the fifth initial value is 25m, and the sixth initial value is 1; when the processor calculates the first bottom hole pressure solution according to the formula (1), besides the first initial value, a fifth initial value needs to be called and substituted into the formula (1) for calculation; when the processor calculates a second bottom hole pressure solution according to the formula (2), a second initial value, a third initial value and a fourth initial value need to be called out and substituted into the formula (2) for calculation, and theoretical pressure data of the fracture-cavity body are obtained according to the calculation result;
after the processor obtains theoretical pressure data, further determining a calculation value of the parameters of the beaded cave body karst cave reservoir according to the actual pressure data and the theoretical pressure data, namely drawing an actual pressure double-logarithm curve according to the actual pressure data; drawing a theoretical pressure double logarithmic curve according to theoretical pressure data; after fitting the actual pressure double-logarithmic curve and the theoretical pressure double-logarithmic curve, calculating a first coincidence rate and a second coincidence rate, if the first coincidence rate is within a first preset coincidence rate range and the second coincidence rate is within a second preset coincidence rate range, indicating that the actual pressure double-logarithmic curve and the theoretical pressure double-logarithmic curve reach the standard of successful fitting, and further determining that the solved value of each parameter is the corresponding initial value; otherwise, the initial values of the parameters are endowed again and the subsequent steps are executed until the actual pressure double-logarithm curve and principleThe theoretical pressure double logarithmic curve reaches the standard of successful fitting, the solved value of each parameter is determined to be the corresponding newly-assigned initial value, and in the embodiment of the invention, the second initial value is re-assigned to be 5.1m3and/MPa, the third initial value is reassigned to 0.046, the fourth initial value is reassigned to 100.86mD, the fifth initial value is reassigned to 26m, and the sixth initial value is assigned to 1, so that the actual pressure double-logarithmic curve and the theoretical pressure double-logarithmic curve reach the standard of successful fitting. In this embodiment, the above parameters of the beaded fracture-cavity karst cave reservoir may be used to determine whether a stimulation modification measure is required to improve the development effect, so as to guide the development of the oil and gas field more accurately.
The invention provides a method for calculating the volume of a beaded cave body karst cave, which comprises the steps of obtaining actual pressure data inside a beaded cave body, a first initial value of a beaded cave body karst cave storage coefficient and a first cave body well bottom pressure solution model; inputting the first initial value into a first fracture-cavity bottom pressure solution model to obtain theoretical pressure data of the beaded fracture-cavity body; determining a solved value of a karst cave storage coefficient of the beaded fracture-cavity body according to actual pressure data and theoretical pressure data; the method is simple and easy to understand and is simple to operate, the well bottom pressure solution model of the first fracture-cavity body is more suitable for the actual conditions of the oil field, interpretation efficiency can be greatly improved when the method is applied to the oil field, and oil field production is facilitated.
Another embodiment of the present invention provides a novel processor configured to perform the above-described method for calculating a cavern volume of a beaded fracture cavity.
Another embodiment of the present invention provides a novel machine-readable storage medium having stored thereon instructions for causing a machine to perform the above-described method for calculating a cavern volume of a beaded fracture cavity.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In a typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include forms of volatile memory in a computer readable medium, Random Access Memory (RAM) and/or non-volatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). The memory is an example of a computer-readable medium.
Computer-readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), Static Random Access Memory (SRAM), Dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), Read Only Memory (ROM), Electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), Digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, a computer readable medium does not include a transitory computer readable medium such as a modulated data signal and a carrier wave.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other identical elements in the process, method, article, or apparatus that comprises the element.
The above are merely examples of the present application and are not intended to limit the present application. Various modifications and changes may occur to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the scope of the claims of the present application.
Claims (10)
1. A method for calculating the volume of a beaded fracture-cavity karst cave, which is characterized by comprising the following steps:
acquiring actual pressure data inside the beaded fracture-cavity body, a first initial value of a karst cave storage coefficient of the beaded fracture-cavity body and a first fracture-cavity body bottom pressure solution model;
inputting the first initial value into the first fracture-cave body well bottom pressure solution model to obtain theoretical pressure data of the beaded fracture-cave body;
determining a solved value of the string-bead-shaped slot cavity karst cave storage coefficient according to the actual pressure data and the theoretical pressure data;
and determining the volume of the karst cave in the beaded seam hole body according to the solved value.
2. The method for calculating a beaded fracture cavity cavern volume of claim 1, wherein the inputting the first initial value into the first fracture-cavity bottom pressure solution model to obtain theoretical pressure data of the beaded fracture cavity comprises:
inputting the first initial value into the first fracture-cavern body bottom-hole pressure solution model to obtain a first bottom-hole pressure solution, wherein the first bottom-hole pressure solution refers to a dimensionless bottom-hole pressure of the beaded fracture-cavern body without considering a well reservoir skin under Laplace space;
acquiring a second fracture-cavity bottom pressure solution model;
inputting the first bottom hole pressure solution into the second fracture-cavity bottom pressure solution model to obtain a second bottom hole pressure solution, wherein the second bottom hole pressure solution refers to a pressure drop of the beaded fracture in the Laplace space under consideration of the well reservoir skin;
and obtaining theoretical pressure data of the beaded fracture-hole body according to the second bottom hole pressure solution.
3. The method of calculating a beaded fracture-cavity cavern volume of claim 2, wherein the inputting the first initial value into the first fracture-cavity bottom pressure solution model to obtain a first bottom pressure solution comprises:
acquiring a karst cave equivalent radius, a fracture storage capacity ratio, a karst cave channeling coefficient, a fracture medium porosity, a fracture medium comprehensive compression coefficient, a pore medium porosity and a pore medium comprehensive compression coefficient;
calculating the first bottom hole pressure solution according to the following equation:
wherein the content of the first and second substances,
is the first bottom hole pressure solution; s is a laplace variable; omega is the fracture storage-capacity ratio; lambda is the karst cave channeling coefficient;
is the fracture medium porosity; c. CtfThe comprehensive compression coefficient of the fracture medium is obtained;
porosity of the porous medium; c. CtvThe comprehensive compression coefficient of the pore medium is obtained; cvIs the first initial value; r isfIs the equivalent radius of the karst cave.
4. The method of calculating a beaded fracture-cavity cavern volume of claim 2, wherein the inputting the first bottom hole pressure solution into the second fracture-cavity bottom hole pressure solution model to obtain a second bottom hole pressure solution comprises:
acquiring the viscosity of crude oil, the formation permeability, the volume coefficient of the crude oil, the yield of an oil well, the storage coefficient of a shaft and the skin coefficient;
calculating the second bottom hole pressure solution according to the following equation:
wherein the content of the first and second substances,
is the second bottom hole pressure solution; μ is the crude oil viscosity; b is the crude oil volume coefficient; q is the oil well production; k is a radical offIs the formation permeability; s is the epidermis coefficient; c is the wellbore reservoir coefficient.
5. The method for calculating the volume of a beaded cave body cavern according to claim 1, wherein the determining an evaluation value of the beaded cave body cavern coefficient according to the actual pressure data and the theoretical pressure data comprises:
drawing an actual pressure double-logarithmic curve according to the actual pressure data, wherein the actual pressure double-logarithmic curve comprises an actual pressure value curve and an actual pressure value derivative curve;
drawing a theoretical pressure double-logarithmic curve according to the theoretical pressure data, wherein the theoretical pressure double-logarithmic curve comprises a theoretical pressure value curve and a theoretical pressure value derivative curve;
and determining the solved value of the storage coefficient of the beaded cave body cavern according to the actual pressure log curve and the theoretical pressure log curve.
6. The method of calculating a beaded fracture-cavity cavern volume of claim 5, wherein the determining the derived value of the beaded fracture-cavity reservoir coefficient from the actual pressure log-log curve and the theoretical pressure log-log curve comprises:
calculating a first coincidence rate between the actual pressure value curve and the theoretical pressure value curve and a second coincidence rate between the actual pressure value derivative curve and the theoretical pressure value derivative curve;
and determining the calculated value of the string-bead-shaped fracture-cavity karst cave storage coefficient according to the first coincidence rate and the second coincidence rate.
7. The method of calculating a beaded fracture-cavity cavern volume of claim 6, wherein the determining the derived value of the beaded fracture-cavity reservoir coefficient as a function of the first and second rates of coincidence comprises:
comparing the first coincidence rate with the first preset coincidence rate range to judge whether the first coincidence rate is within the first preset coincidence rate range;
comparing the second coincidence rate with the second preset coincidence rate range to judge whether the second coincidence rate is within the second preset coincidence rate range;
and determining the solved value of the karst cave reservoir coefficient as the first initial value under the condition that the first coincidence rate is within the first preset coincidence rate range and the second coincidence rate is within the second preset coincidence rate range.
8. The method of calculating a karst cave volume of a beaded fracture cavity as claimed in claim 7, wherein said determining the karst cave volume in the beaded fracture cavity according to the evaluation value comprises:
acquiring a fluid compression coefficient;
calculating the cavern volume according to the following formula:
wherein V is the volume of the karst cave; c. CfIs the fluid compressibility factor.
9. A processor, characterized in that the processor is configured to perform the method of calculating a beaded fracture cavity cavern volume of any one of claims 1 to 8.
10. A machine-readable storage medium having instructions stored thereon for causing a machine to perform the method of calculating a beaded-cavity body cavern volume of any of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110949690.5A CN113673101B (en) | 2021-08-18 | 2021-08-18 | Calculation method and processor for karst cave volume of beaded type fracture cave body |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110949690.5A CN113673101B (en) | 2021-08-18 | 2021-08-18 | Calculation method and processor for karst cave volume of beaded type fracture cave body |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113673101A true CN113673101A (en) | 2021-11-19 |
| CN113673101B CN113673101B (en) | 2024-03-29 |
Family
ID=78543668
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110949690.5A Active CN113673101B (en) | 2021-08-18 | 2021-08-18 | Calculation method and processor for karst cave volume of beaded type fracture cave body |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113673101B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108798634A (en) * | 2018-07-09 | 2018-11-13 | 合肥辰工科技有限公司 | A kind of fractured-vuggy reservoir list hole-mean value of single layer Well Test Data Analysis Method |
| CN109138974A (en) * | 2017-06-19 | 2019-01-04 | 中国石油化工股份有限公司 | A kind of fractured-cavernous carbonate reservoir discrete values well test analysis method and system |
| CN110593865A (en) * | 2019-09-29 | 2019-12-20 | 中国石油集团川庆钻探工程有限公司 | Well testing interpretation method for characteristic parameters of oil reservoir fracture hole |
| CN111444462A (en) * | 2020-04-03 | 2020-07-24 | 中国石油大学(北京) | Method and equipment for measuring and calculating bead body data according to unstable well testing |
| CN111581786A (en) * | 2020-04-19 | 2020-08-25 | 东北石油大学 | Well testing interpretation model for analyzing fracture-cave series-connection mode double-hole composite reservoir |
| CN112377184A (en) * | 2020-11-30 | 2021-02-19 | 中国石油天然气股份有限公司 | Physical property parameter analysis method and device for beaded fracture-cave carbonate reservoir |
| US20210164345A1 (en) * | 2019-05-06 | 2021-06-03 | Xi'an Huaxian Petroleum Technology Co., Ltd | A Flow Simulation and Transient Well Analysis Method Based on Generalized Tube Flow and Percolation Coupling |
-
2021
- 2021-08-18 CN CN202110949690.5A patent/CN113673101B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109138974A (en) * | 2017-06-19 | 2019-01-04 | 中国石油化工股份有限公司 | A kind of fractured-cavernous carbonate reservoir discrete values well test analysis method and system |
| CN108798634A (en) * | 2018-07-09 | 2018-11-13 | 合肥辰工科技有限公司 | A kind of fractured-vuggy reservoir list hole-mean value of single layer Well Test Data Analysis Method |
| US20210164345A1 (en) * | 2019-05-06 | 2021-06-03 | Xi'an Huaxian Petroleum Technology Co., Ltd | A Flow Simulation and Transient Well Analysis Method Based on Generalized Tube Flow and Percolation Coupling |
| CN110593865A (en) * | 2019-09-29 | 2019-12-20 | 中国石油集团川庆钻探工程有限公司 | Well testing interpretation method for characteristic parameters of oil reservoir fracture hole |
| CN111444462A (en) * | 2020-04-03 | 2020-07-24 | 中国石油大学(北京) | Method and equipment for measuring and calculating bead body data according to unstable well testing |
| CN111581786A (en) * | 2020-04-19 | 2020-08-25 | 东北石油大学 | Well testing interpretation model for analyzing fracture-cave series-connection mode double-hole composite reservoir |
| CN112377184A (en) * | 2020-11-30 | 2021-02-19 | 中国石油天然气股份有限公司 | Physical property parameter analysis method and device for beaded fracture-cave carbonate reservoir |
Non-Patent Citations (5)
| Title |
|---|
| RUIZHONG JIANG等: "Characteristics of transient pressure performance of horizontal wells in fractured vuggy tight fractal reservoirs considering nonlinear seepage", 《IFP ENERGIES NOUVELLES》, pages 1 - 14 * |
| WENYANG SHI等: "Transient Preussure Behavior of Acid Fracturing Well in Tidal Flat Carbonate Reservoirs", 《SPE KINGDOM OF SAUDI ARABIA ANNUAL TECHNICAL SYMPOSIUM AND EXHIBITION》, pages 1 - 35 * |
| 史文洋等: "川西潮坪相裂缝型碳酸盐岩分层酸压井压力动态分析", 《岩性油气藏》, vol. 32, no. 1, pages 152 - 160 * |
| 林永学等: "考虑应力敏感的缝洞型油藏三重介质试井分析理论模型", 《科学技术与工程》, vol. 21, no. 11, pages 4345 - 4351 * |
| 蔡明金等: "碳酸盐岩洞穴型储层试井解释新模型", 《特种油气藏》, vol. 21, no. 2, pages 98 - 101 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113673101B (en) | 2024-03-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107701180B (en) | Original oil reservoir water saturation calculation method based on closed coring | |
| CN106777651B (en) | Oil-water well yield splitting method based on injection-production balance principle | |
| CN105893679B (en) | Well Test Data Analysis Method is corrected in low yield horizontal well afterflow | |
| CN110472372B (en) | Dual-medium-based permeability prediction method and system | |
| CN102748007B (en) | Well testing analytical method and device | |
| CN107038290A (en) | A kind of computational methods for the cave-in pressure for considering pressurized effect power | |
| Permadi et al. | An improved water coning calculation for horizontal wells | |
| CN111222252A (en) | Method and system for predicting oil-water two-phase post-pressure productivity of low-saturation oil reservoir | |
| CN110735636A (en) | Method and system for measuring and calculating equivalent boundary distance of reservoir of multi-layer commingled production well | |
| CN111582532B (en) | Method and device for predicting liquid production capacity of horizontal well of stress sensitive oil reservoir | |
| CN108920824B (en) | Yield decrement analysis method for narrow river channel reservoir | |
| CN113673101B (en) | Calculation method and processor for karst cave volume of beaded type fracture cave body | |
| CN110309611A (en) | Air water two phase fluid flow law forecasting method and system based on air water thickness distribution | |
| CN111577264A (en) | Method and device for predicting capacity of fractured-pore oil reservoir horizontal well | |
| CN113553746B (en) | Method and processor for rapidly diagnosing fracture-cavity oil reservoir parameters | |
| Omagbon et al. | Experiences in developing a dual porosity model of the Leyte geothermal production field | |
| CN114169204B (en) | Sand control opportunity determination method for offshore oil and gas field development and production | |
| CN111444462B (en) | Method and equipment for measuring and calculating bead body data according to unstable well testing | |
| CN113107474B (en) | Method and device for predicting reservoir parameters of hypotonic gas reservoirs by using short-term well closing data | |
| CN106285657B (en) | For determining that the method and system of collapsing event occur for fracture-pore reservoir solution cavity | |
| CN106677764B (en) | Method for calculating pressure difference in production of stress sensitivity gas reservoir test | |
| CN108868724B (en) | Method and device for determining gas lift production increasing oil and gas quantity of condensate gas well | |
| WO2021231275A1 (en) | Construction of a high-resolution advanced 3d transient model with multiple wells by integrating pressure transient data into static geological model | |
| CN115522909A (en) | Single well model size determination method and device based on numerical simulation | |
| CN109658016B (en) | Identification method and system for tight gas reservoir supply boundary |
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 |