CN115344821A - Method for determining static reserve threshold of dissolved-solution reservoir, electronic device and medium - Google Patents
Method for determining static reserve threshold of dissolved-solution reservoir, electronic device and medium Download PDFInfo
- Publication number
- CN115344821A CN115344821A CN202110518617.2A CN202110518617A CN115344821A CN 115344821 A CN115344821 A CN 115344821A CN 202110518617 A CN202110518617 A CN 202110518617A CN 115344821 A CN115344821 A CN 115344821A
- Authority
- CN
- China
- Prior art keywords
- well testing
- curve
- medium
- well
- pressure
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 230000003068 static effect Effects 0.000 title claims abstract description 34
- 238000012360 testing method Methods 0.000 claims abstract description 100
- 238000012544 monitoring process Methods 0.000 claims abstract description 14
- 238000004519 manufacturing process Methods 0.000 claims abstract description 10
- 238000010835 comparative analysis Methods 0.000 claims abstract description 5
- 238000011084 recovery Methods 0.000 claims abstract description 5
- 238000000926 separation method Methods 0.000 claims abstract 2
- 239000003921 oil Substances 0.000 claims description 36
- 238000003860 storage Methods 0.000 claims description 17
- 230000035699 permeability Effects 0.000 claims description 16
- 239000010779 crude oil Substances 0.000 claims description 12
- 239000011148 porous material Substances 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 12
- 230000006835 compression Effects 0.000 claims description 10
- 238000007906 compression Methods 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 9
- 238000009826 distribution Methods 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 4
- 238000004364 calculation method Methods 0.000 description 22
- 238000010586 diagram Methods 0.000 description 12
- 239000012530 fluid Substances 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 230000005484 gravity Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000005553 drilling Methods 0.000 description 3
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 239000002904 solvent Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 108010001267 Protein Subunits Proteins 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000004090 dissolution Methods 0.000 description 1
- 230000002500 effect on skin Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- 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
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Data Mining & Analysis (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- General Engineering & Computer Science (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Educational Administration (AREA)
- Algebra (AREA)
- Agronomy & Crop Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Game Theory and Decision Science (AREA)
- Computer Hardware Design (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Quality & Reliability (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The application discloses a method for determining a static reserve threshold value of an oil reservoir of an oil-gas separation system, electronic equipment and a medium. The method can comprise the following steps: acquiring an actual well testing interpretation curve according to the well testing monitoring data; building well testing curve theoretical plates with different threshold values through seismic carving; and respectively carrying out comparative analysis on the well testing curve theoretical plate and the actual well testing interpretation curve with different threshold values, and taking the threshold value corresponding to the well testing curve with the highest fitting precision as an accurate threshold value. The method determines the seismic attribute threshold value in the process of calculating the static reserves of the dissolved oil reservoir by means of pressure recovery well testing monitoring data and numerical well testing, effectively eliminates the influence of factors such as a working system, a production speed and the like on the seismic threshold value, reduces the uncertainty of calculating the static reserves of the dissolved oil reservoir, and tamps the theoretical basis of calculating the reserves of the dissolved oil reservoir.
Description
Technical Field
The invention relates to the field of oil and gas field development, in particular to a method for determining a static reserve threshold value of an oil reservoir of an oil-gas dissolving body, electronic equipment and a medium.
Background
The fractured-solvent oil reservoir is a new carbonate reservoir type discovered in recent years, different from a fractured-vuggy carbonate reservoir, the fractured-solvent oil reservoir has weaker corrosion action, and the main reservoir space type is a cave or a section cavity or a crack formed by transformation under the fracture action. In addition, unlike the reservoir plane spread of conventional sandstone reservoirs, the solution-fractured reservoir is typically a massive "platelike body" as shown in fig. 1.
The calculation of early static reserves of an oil reservoir of an interrupted solution is continuously carried out by using a static volume method of a sandstone oil reservoir:
wherein N is o For petroleum geological reserves, 10 4 m 3 (ii) a V is the reservoir volume, m 3 V = Ah, A is the oil-containing area, km 2 H is the average effective thickness, m;
is the average effective porosity, f; s oi Is the average original oil saturation, f; rho o Is the average ground crude oil density, g/cm 3 ；B oi Is the average crude oil volume coefficient; the volume V of the reservoir body is obtained by carving the slot body; porosity of
The wave impedance property and porosity relation curve obtained through analysis and statistics of a large number of oil wells.
However, the oil deposit of this type is mainly controlled by multi-phase tectonic activity, karst action and oil-water filling adjustment, and the geological boundary and the oil-water interface are difficult to be effectively determined. The reserves calculation results have severe deviations. At present, on the basis of a clastic rock volume method reserve calculation method, a fracture-cave body carving method which is a plane sub-unit, sub-reservoir type and longitudinal segmentation solution reservoir reserve calculation method is preliminarily formed.
The carving method of the slot body can be simply described as that the actual storage amount is directly calculated by three-dimensional grid integration according to the carving effective pore volume result. Through well-seismic combination and inversion of reservoir parameters of a single well, different types of fracture and hole body volumes are directly carved in a three-dimensional space, seismic volumes of cave type reservoirs and hole type reservoirs are directly identified, the fracture type reservoir volumes are calculated by using conventional well logging data and adopting a volume method.
The oil saturation is obtained by well logging interpretation and laboratory tests. Although the accuracy of the method is improved, the reservoir calculation parameters have certain uncertainty due to the multi-solution of seismic attributes and the fact that most karst cave drilling wells cannot be logged. As can be seen from the reserve calculation formula, the accuracy of the reservoir volume V has a large impact on the accuracy of the reserve results. In the fracture-cave body carving calculation process, the calculated sizes and volumes of reservoirs are far different under the condition of different seismic attribute threshold values, for example, as shown in fig. 2a and fig. 2b, the accurate threshold value is determined by the scale pairs of different threshold value carving bodies, and therefore, the clear and reasonable carving body volume is the basis of the calculation of the reserves of the fractured-solution oil reservoir.
In the prior art, the seismic attribute threshold value is determined by mainly utilizing the seismic attribute (amplitude change rate) and dynamic data (capacity) for calibration. However, the productivity is obviously influenced by factors such as a working system, a sampling speed and the like, so that the method has larger uncertainty in determining the threshold value.
Therefore, it is necessary to develop a method, an electronic device and a medium for determining the threshold value of the static reserve of the dissolved oil reservoir, effectively eliminate the influence of factors such as a working system and a sampling speed on the seismic threshold value, reduce the uncertainty of the calculation of the static reserve of the dissolved oil reservoir, and tamp the theoretical basis of the calculation of the reserve of the dissolved oil reservoir.
The information disclosed in this background section is only for enhancement of understanding of the general background of the invention and should not be taken as an acknowledgement or any form of suggestion that this information forms the prior art already known to a person skilled in the art.
Disclosure of Invention
The invention provides a method for determining a threshold value of a static reserve of an oil reservoir of an fractured-solution body, electronic equipment and a medium, which can determine a threshold value of an earthquake attribute in the calculation process of the static reserve of the oil reservoir of the fractured-solution body by means of pressure recovery well testing monitoring data and numerical well testing, effectively eliminate the influence of factors such as a working system, a mining speed and the like on the earthquake threshold value, reduce the uncertainty of the calculation of the static reserve of the oil reservoir of the fractured-solution body and tamp the theoretical basis of the calculation of the reserve of the oil reservoir of the fractured-solution body.
In a first aspect, an embodiment of the present disclosure provides a method for determining a threshold of a static reserve of an oil reservoir, including:
acquiring an actual well testing interpretation curve according to the well testing monitoring data;
building well testing curve theoretical charts with different threshold values through seismic carving;
and respectively carrying out comparative analysis on the well testing curve theoretical plates with different threshold values and the actual well testing interpretation curve, and taking the threshold value corresponding to the well testing curve with the highest fitting precision as an accurate threshold value.
Preferably, obtaining an actual well test interpretation curve according to the well test monitoring data comprises:
according to the well testing monitoring data, calculating a pressure variation derivative through well shut-in pressure recovery;
and respectively constructing a relation curve of the pressure variation, the pressure variation derivative and the time variation, namely the actual well testing interpretation curve.
Preferably, the pressure variation derivative is calculated by equation (1):
wherein, Δ P' is a derivative of pressure variation, Δ P is pressure variation, Δ t is time variation, n is a total number of data points calculated in segments, and i is a number of data points calculated in segments.
Preferably, the building of the theoretical plate of the well test curve with different threshold values through seismic carving comprises the following steps:
acquiring distribution maps of reservoirs with different threshold values through seismic carving;
and carrying out numerical processing on the reservoirs with different threshold values, and establishing well testing curve theoretical plates with different threshold values.
Preferably, the reservoirs with different threshold values are subjected to numerical processing, and well testing curve theoretical plates with different threshold values are established:
establishing a numerical well testing model, including a hole medium seepage equation and a fracture medium seepage equation, and determining an initial condition, an inner boundary condition and an outer boundary condition of the numerical well testing model;
performing numerical processing on the reservoirs with different threshold values through the numerical well testing model, and solving the numerical well testing model to obtain the bottom hole pressures with different threshold values;
and establishing a well testing curve theoretical plate with different threshold values.
Preferably, the pore medium seepage equation is as follows:
wherein k is z For longitudinal permeability, C t To synthesize the compression factor, p v Is the pressure of the pore medium, r is the pressure propagation radius, k vr Is the permeability of the medium in the hole plane, z is the vertical distance of pressure propagation,
is the porosity of the porous medium, mu v Is the viscosity of the crude oil in the pore medium.
Preferably, the fracture medium seepage equation is as follows:
wherein, C t To synthesize the compression factor, p f For the well bottom flowing pressure, h F Is the thickness of the reservoir, mu f Is the viscosity of the crude oil in the fracture medium,
is the fracture medium porosity.
Preferably, the initial conditions are:
the inner boundary conditions are as follows:
the outer boundary conditions are as follows:
wherein h is thickness, P is pressure, P is i Is the original formation pressure, k r As transverse permeability, k z Is the longitudinal permeability,. Phi.is the porosity,. Mu.is the viscosity,. C t To synthesize the compression factor, r e Is an outer boundary distance record, S is a skin coefficient, C is a wellbore storage coefficient, Q is a production, p wf For well bottom flowing pressure, p F Pressure of the fracture medium, mu F Is the crude oil viscosity, k, in the fracture medium F Is the permeability of the crack, r w Is the radius of the wellbore, h V Is the pore media thickness.
As a specific implementation of the embodiments of the present disclosure,
in a second aspect, an embodiment of the present disclosure further provides an electronic device, including:
a memory storing executable instructions;
and the processor runs the executable instructions in the memory to realize the method for determining the static reserve threshold of the solution-breaking reservoir.
In a third aspect, an embodiment of the present disclosure further provides a computer-readable storage medium, where a computer program is stored, and when the computer program is executed by a processor, the method for determining a static reserve threshold of an oil-gas reservoir is implemented.
The methods and apparatus of the present invention have other features and advantages which will be apparent from or are set forth in detail in the accompanying drawings and the following detailed description, which are incorporated herein, and which together serve to explain certain principles of the invention.
Drawings
The above and other objects, features and advantages of the present invention will become more apparent by describing in more detail exemplary embodiments thereof with reference to the attached drawings, in which like reference numerals generally represent like parts.
FIG. 1 illustrates a schematic diagram of an exemplary dissolved-fluid reservoir sculpting, according to one embodiment of the present invention.
Fig. 2a and 2b show schematic diagrams comparing volumes of sculpted volumes with thresholds of 2.1 and 1, respectively, for the same oilfield seismic attribute, according to an embodiment of the invention.
FIG. 3 is a flow chart illustrating the steps of a method for determining a static reserve threshold for an interrupted reservoir in accordance with one embodiment of the present invention.
Fig. 4 shows a schematic diagram of an actual well test interpretation curve of an SHB3 well according to an embodiment of the invention.
Fig. 5a, 5b, 5c, 5d, and 5e are diagrams illustrating the result of carving a slot with discontinuity thresholds of 0.075, 0.103, 0.136, 0.165, and 0.195, respectively, according to an embodiment of the present invention.
FIG. 6 shows a schematic diagram of a physical model for a well test considering a gravity multi-media reservoir according to one embodiment of the invention.
Fig. 7a, 7b, 7c, 7d, and 7e respectively show schematic diagrams of a fracture cave carving volume numerical well testing mesh model with discontinuity threshold values of 0.075, 0.103, 0.136, 0.165, and 0.195, respectively, according to an embodiment of the invention.
FIG. 8 is a diagram illustrating fitting of a well testing profile to a true test curve for different threshold values according to one embodiment of the present invention.
Detailed Description
Preferred embodiments of the present invention will be described in more detail below. While the following describes preferred embodiments of the present invention, it should be understood that the present invention may be embodied in various forms and should not be limited by the embodiments set forth herein.
The numerical well testing forward analysis shows that the well testing interpretation curve of the dissolved oil reservoir has the following regularity: the larger the thickness of the reservoir body is, the deeper and wider the well testing curve concave, the larger the radius of the reservoir body is, the later the response of the boundary of the reservoir body is, and the overall expression is that the larger the scale of the reservoir body is, the deeper the well testing curve is concave. The size of the reservoir is obviously influenced by the carving threshold value of the slot hole body. Therefore, based on the forward model result, the threshold value is determined by dynamic-static combination by using a numerical well testing method.
The invention provides a method for determining a static reserve threshold value of a solution reservoir, which comprises the following steps:
acquiring an actual well testing interpretation curve according to the well testing monitoring data; in one example, obtaining an actual well test interpretation curve from the well test monitoring data comprises:
according to the well testing monitoring data, calculating a pressure variation derivative through well shut-in pressure recovery;
and respectively constructing a relation curve of the pressure variation, the pressure variation derivative and the time variation, namely an actual well testing interpretation curve.
In one example, the pressure change amount derivative is calculated by equation (1):
wherein, Δ P' is a derivative of pressure variation, Δ P is pressure variation, Δ t is time variation, n is the total number of data points calculated in segments, and i is the number of data points calculated in segments.
Building well testing curve theoretical plates with different threshold values through seismic carving; in one example, creating a theoretical plate of well test curves for different thresholds by seismic sculpting includes:
acquiring the distribution maps of the reservoirs with different threshold values through seismic carving;
and carrying out numerical processing on the reservoirs with different threshold values, and establishing well testing curve theoretical plates with different threshold values.
In one example, reservoirs with different thresholds are numerically processed, and well testing curve theoretical charts with different thresholds are established:
establishing a numerical well testing model, including a hole medium seepage equation and a fracture medium seepage equation, and determining an initial condition, an inner boundary condition and an outer boundary condition of the numerical well testing model;
carrying out numerical processing on the reservoirs with different threshold values through a numerical well testing model, solving the numerical well testing model, and obtaining bottom hole pressures with different threshold values;
and establishing well testing curve theoretical plates with different threshold values.
In one example, the void medium seepage equation is:
wherein k is z For longitudinal permeability, C t To synthesize the compression factor, p v Is the pressure of the pore medium, r is the pressure propagation radius, k vr Is the permeability of the medium in the hole plane, z is the vertical distance of pressure propagation,
is the porosity of the porous medium, mu v Is the viscosity of the crude oil in the pore medium.
In one example, the fracture medium seepage equation is:
wherein, C t To synthesize the compression factor, p f For the well bottom flowing pressure, h F Is the thickness of the reservoir, mu f Is the viscosity of the crude oil in the fracture medium,
is the porosity of the fracture medium.
In one example, the initial conditions are:
the inner boundary conditions were:
the outer boundary conditions are as follows:
wherein h is thickness, P is pressure, P is i Is the original formation pressure, k r As transverse permeability, k z Is the longitudinal permeability, phi is the porosity, mu is the viscosity, C t To synthesize the compression factor, r e Recording the distance of the outer boundary, S is the skin coefficient, C is the well bore storage coefficient, Q is the production, p wf For well bottom flow pressure, p F Pressure of the fracture medium, mu F Is the crude oil viscosity, k, in the fracture medium F Is the permeability of the crack, r w Is the radius of the wellbore, h V Is the pore media thickness.
Specifically, a large number of geophysical attributes and production well dynamic comparison results show that the discontinuous attributes can effectively represent reservoirs of the dissolved-fluid reservoir, reservoir distribution graphs under different threshold values are extracted at equal intervals according to the discontinuous attributes, and the spreading of the reservoirs and fractures under the different threshold values is determined.
The fractured solution reservoir is a typical 'huge thick plate-shaped' reservoir, the oil-gas flow mainly adopts vertical flow, the traditional numerical well testing method which mainly adopts planar flow and aims at sandstone is not suitable, and a numerical well testing method which is suitable for the characteristics of the fractured solution reservoir needs to be established.
The permeability of the horizontal reservoir and the longitudinal reservoir of the multi-medium oil reservoir is different, the influence of fluid gravity is considered, and the influence term of gravity is added into the longitudinal seepage velocity term. The physical property parameters of the oil deposit are constant, the single-phase crude oil is slightly compressible, the compression coefficient is constant, the outer boundary of the oil deposit can be an infinite boundary, a constant pressure boundary or a closed boundary, the skin effect and the shaft storage effect are considered, the fluid is produced with the constant yield Q, the radial and longitudinal seepage (three-dimensional unsteady seepage) of the solution oil deposit is considered, the fluid is Darcy seepage, and the influence of capillary force is ignored.
The numerical well testing model comprises a hole medium seepage equation (2) and a fracture medium seepage equation (3), and initial conditions, inner boundary conditions and outer boundary conditions of the numerical well testing model are determined, namely formulas (4) - (6).
And (4) carrying out numerical solution on the unstructured grid by adopting a central difference method. In order to prevent the unconvergence condition, all the discrepancies adopt implicit difference, and the boundary condition processing introduces a virtual mirror image grid. In order to improve the calculation efficiency, the space grid adopts a radial unstructured grid, and the early time step is dense. Due to the adoption of the non-uniform grids, in order to increase the calculation accuracy and the convergence, an upwind model discrete seepage partial differential equation, in particular a first-order partial differential equation, is introduced. And finally, programming and solving the matrix according to the composition matrix to obtain a bottom hole pressure solution.
And selecting an unstructured grid for a reservoir body model based on different threshold values, digitizing the reservoir body model, and establishing a quasi-three-dimensional numerical model which accords with reservoir body and crack distribution characteristics.
Aiming at geological models under different threshold value oil reservoir condition, a numerical value well testing forward analysis is carried out by applying a numerical value well testing model of the dissolution reservoir, and well testing curve charts under different threshold value conditions are established.
And respectively carrying out comparative analysis on the well testing curve theoretical plate and the actual well testing interpretation curve with different threshold values, and taking the threshold value corresponding to the well testing curve with the highest fitting precision as an accurate threshold value.
Specifically, after the accurate threshold value is determined, parameters such as the volume of the reservoir body in the reserve calculation can be accurately calibrated.
The present invention also provides an electronic device, comprising: a memory storing executable instructions; and the processor runs the executable instructions in the memory to realize the method for determining the static reserve threshold of the solution-breaking oil reservoir.
The invention further provides a computer-readable storage medium, which stores a computer program, and the computer program is executed by a processor to implement the method for determining the threshold value of the static reserve of the solution-breaking reservoir.
To facilitate understanding of the solution of the embodiments of the present invention and the effects thereof, three specific application examples are given below. It will be understood by those skilled in the art that this example is merely for the purpose of facilitating an understanding of the present invention and that any specific details thereof are not intended to limit the invention in any way.
Example 1
FIG. 3 is a flow chart illustrating the steps of a method for determining a static reserve threshold for an interrupted reservoir in accordance with one embodiment of the present invention.
As shown in fig. 3, the method for determining the threshold value of the static reserve of the solution-cutoff reservoir includes: step 101, acquiring an actual well testing interpretation curve according to well testing monitoring data; 102, building a well testing curve theoretical chart with different threshold values through seismic carving; and 103, respectively carrying out comparative analysis on the well testing curve theoretical plate and the actual well testing interpretation curve with different threshold values, and taking the threshold value corresponding to the well testing curve with the highest fitting precision as an accurate threshold value.
The drilling completion depth of the finished SHB3 well is 8342.00m (inclined)/7913.78 m (vertical), and the completion is completed by acid fracturing. 751m of lost mud in the process of drilling 3 . 53t/d oil is produced at the day after the well is put into production, and the oil-gas containing type with the disordered strong reflection characteristics on the fracture zone is verified. However, the well has problems of low productivity, fast energy reduction, etc. after production, which are not consistent with the results of early reserve estimation. And because the production time is short, the oil nozzle is frequently adjusted in the pilot production process, the working system is unstable, the threshold value is difficult to calibrate by using the seismic attribute and the dynamic data, and the storage calculation result has great uncertainty.
Fig. 4 shows a schematic diagram of an actual well test interpretation curve of an SHB3 well according to an embodiment of the invention.
Fig. 5a, 5b, 5c, 5d, and 5e are diagrams illustrating the result of carving a slot with discontinuity thresholds of 0.075, 0.103, 0.136, 0.165, and 0.195, respectively, according to an embodiment of the present invention.
FIG. 6 shows a schematic diagram of a physical model for a well test considering a gravity multi-media reservoir according to one embodiment of the invention.
Well testing monitoring data of the SHB3 well is selected, well testing data analysis work is carried out, well testing curve characteristics of the SHB3 are determined, and an actual well testing interpretation curve is shown in figure 4. Preferably 5 thresholds, such as discontinuity attribute thresholds 0.075, 0.103, 0.136, 0.165, and 0.195, are used for the engraving of the slot, as shown in fig. 5 a-5 e, respectively. Determining the results of the reservoir space under different threshold values, and performing forward research on the well testing of 5 energy attribute threshold values of the SHB3 well region according to the well testing physical model shown in FIG. 6.
Fig. 7a, 7b, 7c, 7d, and 7e respectively show schematic diagrams of a fracture cave carving volume numerical well testing mesh model with discontinuity threshold values of 0.075, 0.103, 0.136, 0.165, and 0.195, respectively, according to an embodiment of the invention.
FIG. 8 is a diagram illustrating fitting of a well testing profile to a true test curve for different threshold values according to one embodiment of the present invention.
According to the fitting results of the test well curve plates with different threshold values and the real test curve shown in fig. 8, the threshold value of the definite attribute is 0.103.
And carrying out reserve calculation by utilizing the determined threshold value, and determining the geological reserve of the well by 90 ten thousand tons. It is clear that the well only controls one set of reservoirs and communicates poorly with deep slots, resulting in low initial production. The calculation result of the reserves is consistent with the dynamic knowledge, and the calculation result is reliable.
The method determines a reasonable threshold value of the carving of the fracture-cavity body, guides the carving of the fracture-cavity body and the determination of the volume of a reservoir body, and further guides the calculation of the reserves of the oil field. The final calculation result is in accordance with the understanding of the oil reservoir geology, the accuracy is greatly improved, the compiling and reserve declaration of the rolling evaluation scheme is effectively guided, the application has a good development effect, the understanding of the method on the oil reservoir is proved to be in accordance with the oil reservoir practice, and the research result has a wide application prospect in similar oil reservoirs at home and abroad.
Example 2
The present disclosure provides an electronic device including: a memory storing executable instructions; and the processor runs the executable instructions in the memory to realize the method for determining the static reserve threshold of the solution-breaking oil reservoir.
An electronic device according to an embodiment of the present disclosure includes a memory and a processor.
The memory is to store non-transitory computer readable instructions. In particular, the memory may include one or more computer program products that may include various forms of computer-readable storage media, such as volatile memory and/or non-volatile memory. The volatile memory may include, for example, random Access Memory (RAM), cache memory (cache), and/or the like. The non-volatile memory may include, for example, read Only Memory (ROM), hard disk, flash memory, etc.
The processor may be a Central Processing Unit (CPU) or other form of processing unit having data processing capabilities and/or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In one embodiment of the disclosure, the processor is configured to execute the computer readable instructions stored in the memory.
Those skilled in the art should understand that, in order to solve the technical problem of how to obtain a good user experience, the present embodiment may also include well-known structures such as a communication bus, an interface, and the like, and these well-known structures should also be included in the protection scope of the present disclosure.
For the detailed description of the present embodiment, reference may be made to the corresponding descriptions in the foregoing embodiments, which are not repeated herein.
Example 3
The embodiment of the disclosure provides a computer-readable storage medium, which stores a computer program, and when the computer program is executed by a processor, the computer program realizes the method for determining the static reserve threshold of the solution-breaking reservoir.
A computer-readable storage medium according to an embodiment of the present disclosure has non-transitory computer-readable instructions stored thereon. The non-transitory computer readable instructions, when executed by a processor, perform all or a portion of the steps of the methods of the embodiments of the disclosure previously described.
The computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROMs and DVDs), magneto-optical storage media (e.g., MOs), magnetic storage media (e.g., magnetic tapes or removable disks), media with built-in rewritable non-volatile memory (e.g., memory cards), and media with built-in ROMs (e.g., ROM cartridges).
It will be appreciated by persons skilled in the art that the above description of embodiments of the invention is for the purpose of illustrating the benefits of embodiments of the invention only, and is not intended to limit embodiments of the invention to any examples given.
Having described embodiments of the present invention, the foregoing description is intended to be exemplary, not exhaustive, and not limited to the embodiments disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments.
Claims (10)
1. A method for determining a static reserve threshold of an oil reservoir of an oil-water-gas separation system is characterized by comprising the following steps:
acquiring an actual well testing interpretation curve according to the well testing monitoring data;
building well testing curve theoretical charts with different threshold values through seismic carving;
and respectively carrying out comparative analysis on the well testing curve theoretical plates with different threshold values and the actual well testing interpretation curve, and taking the threshold value corresponding to the well testing curve with the highest fitting precision as an accurate threshold value.
2. The method for determining the static reserve threshold of the solution-breaking reservoir according to claim 1, wherein the step of obtaining an actual well test interpretation curve according to the well test monitoring data comprises the steps of:
according to the well testing monitoring data, calculating a pressure variation derivative through well shut-in pressure recovery;
and respectively constructing a relation curve of the pressure variation, the pressure variation derivative and the time variation, namely the actual well testing interpretation curve.
3. The method of determining a static reserve threshold for an fractured-solution reservoir according to claim 2, wherein the derivative of the pressure change is calculated by the formula (1):
wherein, Δ P' is a derivative of pressure variation, Δ P is pressure variation, Δ t is time variation, n is the total number of data points calculated in segments, and i is the number of data points calculated in segments.
4. The method for determining the threshold value of the static reserve of the fractured-solution reservoir according to claim 1, wherein the step of establishing the theoretical plate of the well testing curve with different threshold values through seismic carving comprises the following steps:
acquiring the distribution maps of the reservoirs with different threshold values through seismic carving;
and carrying out numerical processing on the reservoirs with different threshold values, and establishing well testing curve theoretical plates with different threshold values.
5. The method for determining the threshold value of the static reserve of the solution reservoir according to claim 4, wherein the reservoirs with different threshold values are subjected to numerical processing, and well testing curve theoretical plates with different threshold values are established:
establishing a numerical well testing model, including a hole medium seepage equation and a fracture medium seepage equation, and determining an initial condition, an inner boundary condition and an outer boundary condition of the numerical well testing model;
carrying out numerical processing on the reservoirs with different threshold values through the numerical well testing model, solving the numerical well testing model, and obtaining bottom hole pressures with different threshold values;
and establishing a well testing curve theoretical plate with different threshold values.
6. The method for determining the static reserve threshold of the solution-breaking reservoir according to claim 5, wherein the pore medium seepage equation is as follows:
wherein k is z For longitudinal permeability, C t To synthesize the compression factor, p v Is the pressure of the pore medium, r is the pressure propagation radius, k vr Is the permeability of the medium in the hole plane, z is the vertical distance of pressure propagation,
is the porosity of the porous medium, mu v Is the viscosity of the crude oil in the pore medium.
7. The method for determining the static reserve threshold of the fractured-solution reservoir according to claim 5, wherein the fracture medium seepage equation is as follows:
wherein, C t To synthesize the compression factor, p f For the well bottom flowing pressure, h F Is the thickness of the reservoir, mu f Is the viscosity of the crude oil in the fracture medium,
is the fracture medium porosity.
8. The method for determining the static reserve threshold of an oil-gas reservoir according to claim 5, wherein the initial conditions are as follows:
the inner boundary conditions are as follows:
the outer boundary conditions are as follows:
wherein h is thickness, P is pressure, P is i Is the original formation pressure, k r As transverse permeability, k z Is the longitudinal permeability, phi is the porosity, mu is the viscosity, C t To synthesize the compression factor, r e Recording the distance of the outer boundary, S is the skin coefficient, C is the well bore storage coefficient, Q is the production, p wf For well bottom flow pressure, p F Pressure of the fracture medium, μ F Is the viscosity of crude oil in the fracture medium, k F Is the permeability of the crack, r w Is the radius of the wellbore, h V Is the pore media thickness.
9. An electronic device, characterized in that the electronic device comprises:
a memory storing executable instructions;
a processor executing the executable instructions in the memory to implement the method of determining a static reserve threshold for an interrupted solution reservoir of any of claims 1-8.
10. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program which, when executed by a processor, implements the method for determining a static reserve threshold for an oil-solution reservoir according to any of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110518617.2A CN115344821A (en) | 2021-05-12 | 2021-05-12 | Method for determining static reserve threshold of dissolved-solution reservoir, electronic device and medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110518617.2A CN115344821A (en) | 2021-05-12 | 2021-05-12 | Method for determining static reserve threshold of dissolved-solution reservoir, electronic device and medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115344821A true CN115344821A (en) | 2022-11-15 |
Family
ID=83947273
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110518617.2A Pending CN115344821A (en) | 2021-05-12 | 2021-05-12 | Method for determining static reserve threshold of dissolved-solution reservoir, electronic device and medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115344821A (en) |
-
2021
- 2021-05-12 CN CN202110518617.2A patent/CN115344821A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11591900B2 (en) | Method to predict overpressure uncertainty from normal compaction trendline uncertainty | |
| EP2497900B1 (en) | Modeling hydraulic fractures | |
| US20150370934A1 (en) | Completion design based on logging while drilling (lwd) data | |
| EP2904530B1 (en) | System, method and computer program product for determining placement of perforation intervals using facies, fluid boundaries, geobodies and dynamic fluid properties | |
| CN110096718B (en) | Method for obtaining volume of karst cave in carbonate reservoir | |
| CN110107277B (en) | Method for obtaining volume of karst cave encountered by well drilling in carbonate reservoir | |
| US10337294B2 (en) | Reservoir permeability upscaling | |
| CN107861162B (en) | Microelectrode logging data-based natural crack identification method and system | |
| CN105911592A (en) | Practical drilling constraint three-dimensional seismic data long wavelength static correction method | |
| McConnell | Double porosity well testing in the fractured carbonate rocks of the Ozarks | |
| CN115344821A (en) | Method for determining static reserve threshold of dissolved-solution reservoir, electronic device and medium | |
| CN113742962B (en) | Shale reservoir attribute three-dimensional modeling method based on horizontal well | |
| CN111963148B (en) | Method for determining pressure coefficient and drilling mud density of carbonate fracture-cavity type oil and gas reservoir | |
| CN112394416A (en) | Heterogeneous fracture control reservoir prediction method and device | |
| CN113738344B (en) | Oil-water interface depth determining method and early warning method for preventing production well water channeling | |
| US10460051B2 (en) | Computationally-efficient modeling of viscous fingering effect for enhanced oil recovery (EOR) agent injected at multiple injection concentrations | |
| CN110596760A (en) | Method for determining pore roar radius in oil reservoir development process | |
| CN117310799B (en) | Mine floor limestone aquifer identification method based on hole-slot-hole-matrix multiple medium | |
| CN110633904A (en) | Metamorphic rock reservoir classification method and device | |
| CN113803055B (en) | Oil-water interface depth determining method and early warning method for preventing production well water channeling | |
| CN113537668B (en) | Method and system for evaluating crude oil geological reserves in unfilled karst cave | |
| Lamberghini et al. | Performance Based Reservoir Characterization in a Tight Gas Reservoir System-Case Study from Lajas and Punta Rosada Formations in the Neuquen Basin, Argentina | |
| CN116629140A (en) | Stratum four-pressure prediction method and device based on artificial intelligence and real drilling data | |
| CN117386356A (en) | Geosteering method and device based on logging data | |
| AU2012396845B2 (en) | System, method and computer program product for determining placement of perforation intervals using facies, fluid boundaries, geobodies and dynamic fluid properties |
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 |