Method for determining static reserve threshold value of fault-solubility oil reservoir, electronic device and medium

Abstract

The application discloses a method for determining a static reserve threshold value of a faulted-dissolved body reservoir, an electronic device and a medium. The method can comprise: obtaining an actual well test interpretation curve according to well test monitoring data; establishing a well test curve theoretical chart of different threshold values through seismic carving; and comparing and analyzing the well test curve theoretical chart of different threshold values and the actual well test interpretation curve respectively, so as to take the threshold value corresponding to the well test curve with the highest fitting accuracy as an accurate threshold value. The application determines a seismic attribute threshold value in a static reserve calculation process of a faulted-dissolved body reservoir by using numerical well test means based on pressure recovery well test monitoring data, effectively eliminates the influence of working systems, production rates and other factors on the seismic threshold value, reduces the uncertainty of the static reserve calculation of the faulted-dissolved body reservoir, and consolidates the theoretical basis of the reserve calculation of the faulted-dissolved body reservoir.

Description

Technical Field

[0001] This invention relates to the field of oil and gas field development, and more specifically, to a method, electronic equipment, and medium for determining the static reserve threshold value of fractured solution reservoirs. Background Technology

[0002] Fault-dissolved reservoirs are a new type of carbonate reservoir discovered in recent years. Unlike fracture-vuggy carbonate reservoirs, fault-dissolved reservoirs exhibit weaker dissolution, and the main reservoir spaces are cavities or fractures formed by faulting. Furthermore, unlike the planar distribution of conventional sandstone reservoirs, fault-dissolved reservoirs typically present as massive, blocky "plate-like bodies," such as... Figure 1 As shown.

[0003] The calculation of early static reserves in fractured solution reservoirs has traditionally followed the static volumetric method for sandstone reservoirs.

[0004]

[0005] Where, N o For oil geological reserves, 10 4 m 3 V is the volume of the storage mass, in meters. 3 V = Ah, where A is the oil-bearing area in km². 2 h is the average effective thickness, in meters; For the average effective porosity, f; S oi For the average original oil saturation, f; ρ o The average ground crude oil density is expressed in g / cm³. 3 B oi The average crude oil volume coefficient; the reservoir volume V is obtained by carving the fractured cavity; porosity The curves showing the relationship between wave impedance and porosity were obtained through statistical analysis of a large number of oil wells.

[0006] However, because this type of reservoir is mainly controlled by multiple phases of tectonic activity, karstification, and oil-water injection adjustments, its geological boundaries and oil-water interfaces are difficult to determine effectively. This results in significant deviations in reserve calculations. Currently, based on the clastic rock volumetric method for reserve calculation, a preliminary reserve calculation method for fracture-cavity reservoirs—the fracture-cavity carving method—has been developed, which involves "planar unit division, reservoir type division, and vertical segmentation."

[0007] The fractured-vuggy body carving method can be simply described as calculating and determining reserves directly through three-dimensional spatial grid integration based on the results of carving the effective pore volume. By combining well and seismic data and inverting reservoir parameters of a single well, different types of fractured-vuggy body volumes are directly carved in three-dimensional space. The seismic volume of cave-type and porous reservoirs is directly identified, while the volume of fractured reservoirs is calculated using conventional well logging data and the volumetric method.

[0008] Oil saturation was obtained through well logging interpretation and laboratory tests. While this method improves accuracy, the uncertainty in reserve calculation parameters remains due to the ambiguity of seismic attributes and the fact that most wells encountering caverns cannot be logged. The reserve calculation formula shows that the accuracy of the reservoir volume V significantly impacts the precision of the reserve results. Furthermore, in the calculation of fractured-cavity body sculpting, it was found that the calculated reservoir size and volume differed considerably under different seismic attribute threshold values, with the sculpted body size showing a significant difference compared to other values. Figure 2a , Figure 2b As shown, determining the accurate threshold value and thus clarifying the reasonable volume of the sculpted body is the basis for calculating the reserves of fractured solution reservoirs.

[0009] In existing technologies, the determination of seismic attribute thresholds mainly relies on calibrating seismic attributes (amplitude variation rate) and dynamic data (production capacity). However, production capacity is significantly affected by factors such as working conditions and mining rate, resulting in considerable uncertainty in determining the threshold value using this method.

[0010] Therefore, it is necessary to develop a method, electronic equipment, and medium for determining the static reserve threshold value of fault-dissolved reservoirs, effectively eliminating the influence of factors such as working conditions and production rate on the seismic threshold value, reducing the uncertainty in the calculation of static reserves of fault-dissolved reservoirs, and consolidating the theoretical foundation for the calculation of reserves of fault-dissolved reservoirs.

[0011] The information disclosed in the background section of this invention is intended only to enhance the understanding of the general background of this invention, and should not be construed as an admission or in any way implying that such information constitutes prior art known to those skilled in the art. Summary of the Invention

[0012] This invention proposes a method, electronic equipment, and medium for determining the static reserve threshold value of fractured-dissolved reservoirs. It can determine the seismic attribute threshold value in the static reserve calculation process of fractured-dissolved reservoirs by means of numerical well testing and monitoring data of pressure recovery well test. It effectively eliminates the influence of factors such as working system and production rate on the seismic threshold value, reduces the uncertainty of static reserve calculation of fractured-dissolved reservoirs, and consolidates the theoretical basis of reserve calculation of fractured-dissolved reservoirs.

[0013] In a first aspect, embodiments of this disclosure provide a method for determining the static reserve threshold value of a fractured-solution reservoir, including:

[0014] Based on well test monitoring data, obtain the actual well test interpretation curves;

[0015] By using seismic engraving, theoretical charts of well test curves with different threshold values ​​are established;

[0016] The theoretical charts of well test curves with different threshold values ​​are compared and analyzed with the actual well test interpretation curves. The threshold value corresponding to the well test curve with the highest fitting accuracy is the precise threshold value.

[0017] Preferably, obtaining the actual well test interpretation curve based on well test monitoring data includes:

[0018] Based on the well test monitoring data, the derivative of the pressure change was calculated by recovering the pressure change through shut-in pressure testing.

[0019] The relationship curves between pressure change, the derivative of pressure change, and time change are constructed respectively, which are the interpretation curves of the actual well test.

[0020] Preferably, the derivative of the pressure change is calculated using formula (1):

[0021]

[0022] Where △P′ is the derivative of the pressure change, △P is the pressure change, △t is the time change, n is the total number of data points calculated in segments, and i is the number of data points calculated in segments.

[0023] Preferably, the theoretical charts for well test curves with different threshold values ​​are established through seismic engraving, including:

[0024] By using seismic sculpting, we can obtain distribution maps of storage groups with different threshold values;

[0025] Numerical processing was performed on reservoirs with different threshold values ​​to establish theoretical charts of well test curves for different threshold values.

[0026] Preferably, reservoirs with different threshold values ​​are quantified to establish theoretical charts of well test curves for different threshold values:

[0027] A numerical well test model is established, including the seepage equations for porous media and fractured media, and the initial conditions, inner boundary conditions, and outer boundary conditions of the numerical well test model are determined.

[0028] The numerical well test model is used to quantify reservoirs with different threshold values, and the numerical well test model is solved to obtain the bottom hole pressure at different threshold values.

[0029] Establish theoretical charts for well test curves with different threshold values.

[0030] Preferably, the seepage equation of the porous medium is:

[0031]

[0032] Where, k z For vertical penetration rate, C t p is the overall compression coefficient. v Let r be the pressure of the porous medium, r be the pressure propagation radius, and k be the pressure of the porous medium. vr Let z be the planar permeability of the porous medium, and z be the vertical distance of pressure propagation. Porosity of porous media, μ v This refers to the viscosity of crude oil within the porous medium.

[0033] Preferably, the seepage equation of the fractured medium is:

[0034]

[0035] Among them, C t p is the overall compression coefficient. f h is the bottom hole flowing pressure of the oil well. F For reservoir thickness, μ f The viscosity of crude oil within the fractured medium. Porosity of the fractured medium.

[0036] Preferably, the initial conditions are:

[0037]

[0038] The inner boundary conditions are as follows:

[0039]

[0040] The outer boundary conditions are as follows:

[0041]

[0042] Where h is the thickness, P is the pressure, and P i k represents the original formation pressure. r k represents the lateral penetration rate. z φ is the longitudinal permeability, μ is the porosity, and C is the viscosity. t r is the overall compression coefficient. e Let S be the outer boundary distance, C be the skin coefficient, C be the wellbore reservoir coefficient, Q be the production rate, and p be the outer boundary distance. wf p is the bottom hole flowing pressure of the oil well. F For the pressure of the fracture medium, μ F k represents the viscosity of crude oil within the fractured medium. F r represents the crack permeability. w h is the radius of the wellbore. V The thickness of the porous medium.

[0043] As one specific implementation of this disclosure,

[0044] Secondly, embodiments of this disclosure also provide an electronic device, the electronic device comprising:

[0045] Memory, which stores executable instructions;

[0046] A processor that executes the executable instructions in the memory to implement the method for determining the static reserve threshold value of a fractured solution reservoir.

[0047] Thirdly, this disclosure also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the method for determining the static reserve threshold value of a fractured solution reservoir.

[0048] The methods and apparatus of the present invention have other features and advantages that will be apparent from or will be set forth in detail in the accompanying drawings and following detailed description, which together serve to explain the particular principles of the invention. Attached Figure Description

[0049] The above and other objects, features and advantages of the present invention will become more apparent from the more detailed description of exemplary embodiments of the invention in conjunction with the accompanying drawings, wherein the same reference numerals generally represent the same parts.

[0050] Figure 1 A schematic diagram of a typical fractured solution reservoir engraving according to an embodiment of the present invention is shown.

[0051] Figure 2a and Figure 2b The diagrams show a comparison of the volumes of sculpted objects with seismic attribute threshold values ​​of 2.1 and 1 for the same oilfield, according to an embodiment of the present invention.

[0052] Figure 3 A flowchart illustrating the steps of a method for determining the static reserve threshold value of a fractured solution reservoir according to an embodiment of the present invention is shown.

[0053] Figure 4 A schematic diagram of the actual well test interpretation curve of the SHB3 well according to an embodiment of the present invention is shown.

[0054] Figure 5a , Figure 5b , Figure 5c , Figure 5d , Figure 5e The diagrams show the results of carving the grout with discontinuity threshold values ​​of 0.075, 0.103, 0.136, 0.165 and 0.195, respectively, according to an embodiment of the present invention.

[0055] Figure 6 A schematic diagram of a gravity-based multi-medium reservoir well test physical model is shown according to an embodiment of the present invention.

[0056] Figure 7a , Figure 7b , Figure 7c , Figure 7d , Figure 7e Schematic diagrams of numerical well test grid models of fractured and cavitated bodies with discontinuity threshold values ​​of 0.075, 0.103, 0.136, 0.165, and 0.195, respectively, according to an embodiment of the present invention, are shown.

[0057] Figure 8 A schematic diagram showing the fitting of well test curves with different threshold values ​​and actual test curves according to an embodiment of the present invention is illustrated. Detailed Implementation

[0058] Preferred embodiments of the invention will now be described in more detail. While preferred embodiments of the invention are described below, it should be understood that the invention can be implemented in various forms and should not be limited to the embodiments set forth herein.

[0059] Numerical well testing and forward modeling analysis revealed the following patterns in the interpretation curves of fractured-dissolved reservoirs: the greater the reservoir thickness, the deeper and wider the concave area of ​​the well test curve; the larger the reservoir radius, the later the reservoir boundary response; generally, the larger the reservoir size, the deeper the concave area of ​​the well test curve. The reservoir size is significantly affected by the fracture-cavity sculpting threshold. Therefore, based on the forward modeling results, numerical well testing methods, combining static and dynamic approaches, are considered to determine the threshold value.

[0060] This invention provides a method for determining the static reserve threshold value of fractured-solution reservoirs, comprising:

[0061] Based on well test monitoring data, obtain the actual well test interpretation curve; in one example, obtaining the actual well test interpretation curve based on well test monitoring data includes:

[0062] Based on well test monitoring data, the derivative of pressure change is calculated by recovering pressure through shut-in pressure testing.

[0063] The relationship curves between pressure change, the derivative of pressure change, and time change are constructed respectively, which are the interpretation curves of the actual well test.

[0064] In one example, the derivative of the pressure change is calculated using formula (1):

[0065]

[0066] Where △P′ is the derivative of the pressure change, △P is the pressure change, △t is the time change, n is the total number of data points calculated in segments, and i is the number of data points calculated in segments.

[0067] Using seismic engraving, theoretical charts of well test curves with different threshold values ​​are established. In one example, establishing theoretical charts of well test curves with different threshold values ​​using seismic engraving includes:

[0068] By using seismic sculpting, we can obtain distribution maps of storage groups with different threshold values;

[0069] Numerical processing was performed on reservoirs with different threshold values ​​to establish theoretical charts of well test curves for different threshold values.

[0070] In one example, reservoirs with different threshold values ​​are quantified to establish theoretical charts of well test curves for different threshold values:

[0071] Establish a numerical well test model, including the seepage equations for porous media and fractured media, and determine the initial conditions, inner boundary conditions, and outer boundary conditions of the numerical well test model;

[0072] The numerical well test model is used to quantify reservoirs with different threshold values, and the numerical well test model is solved to obtain the bottom hole pressure at different threshold values.

[0073] Establish theoretical charts for well test curves with different threshold values.

[0074] In one example, the seepage equation for porous media is:

[0075]

[0076] Where, k z For vertical penetration rate, C t p is the overall compression coefficient. v Let r be the pressure of the porous medium, r be the pressure propagation radius, and k be the pressure of the porous medium. vr Let z be the planar permeability of the porous medium, and z be the vertical distance of pressure propagation. Porosity of porous media, μ v This refers to the viscosity of crude oil within the porous medium.

[0077] In one example, the seepage equation for fractured media is:

[0078]

[0079] Among them, C t p is the overall compression coefficient. f h is the bottom hole flowing pressure of the oil well. F For reservoir thickness, μ f The viscosity of crude oil within the fractured medium. Porosity of the fractured medium.

[0080] In one example, the initial condition is:

[0081]

[0082] The inner boundary conditions are:

[0083]

[0084] The outer boundary conditions are:

[0085]

[0086] Where h is the thickness, P is the pressure, and P i k represents the original formation pressure. r k represents the lateral penetration rate. z φ is the longitudinal permeability, μ is the porosity, and C is the viscosity. t r is the overall compression coefficient. e Let S be the outer boundary distance, C be the skin coefficient, C be the wellbore reservoir coefficient, Q be the production rate, and p be the outer boundary distance. wf p is the bottom hole flowing pressure of the oil well. F For the pressure of the fracture medium, μ F k represents the viscosity of crude oil within the fractured medium. F r represents the crack permeability. w h is the radius of the wellbore. V The thickness of the porous medium.

[0087] Specifically, a large number of geophysical properties and dynamic comparisons with production wells show that discontinuity properties can effectively characterize the reservoir mass of fractured-dissolved oil reservoirs. For discontinuity properties, reservoir mass distribution maps under different threshold values ​​are extracted at equal intervals to clarify the distribution of reservoir mass and fractures under different threshold values.

[0088] Fault-depression reservoirs are typical "thick plate-shaped" reservoirs, with oil and gas flows mainly vertically. The numerical well testing methods that were previously used for sandstone and were mainly based on planar flow are not applicable. It is necessary to establish numerical well testing methods that are suitable for the characteristics of fault-depression reservoirs.

[0089] In multi-medium reservoirs, the permeability of the horizontal and vertical reservoir layers differs. Considering the influence of fluid gravity, a gravity effect term is added to the vertical seepage velocity term. Then, the reservoir physical parameters are constant, the crude oil is slightly compressible single-phase, the compressibility coefficient is constant, and the reservoir outer boundary can be an infinite boundary, a constant-pressure boundary, or a closed boundary. Skin effect and wellbore storage effect are considered, and the fluid is produced at a constant rate Q. Simultaneously, radial and vertical seepage (three-dimensional unsteady-state seepage) is considered in fractured-solution reservoirs, the fluid is Darcy flow, and the capillary force effect is ignored.

[0090] The numerical well test model includes the seepage equation (2) for porous media and the seepage equation (3) for fractured media. The initial conditions, inner boundary conditions and outer boundary conditions of the numerical well test model are determined, namely formulas (4)-(6).

[0091] A central difference method with an unstructured mesh was employed for numerical solution. To prevent non-convergence, all discretizations used implicit differences, and a virtual mirror mesh was introduced for boundary condition handling. To improve computational efficiency, a radially unstructured mesh was used, with denser early time steps. Due to the use of a non-uniform mesh, to increase computational accuracy and convergence, the upwind model was introduced to discretize the seepage partial differential equations, particularly the first-order partial differential equations. Finally, based on the constitutive matrix, the bottomhole pressure solution was obtained by programming and solving the matrix.

[0092] For reservoir models based on different threshold values, unstructured meshes are selected to quantify the reservoir models and establish quasi-three-dimensional numerical models that conform to the distribution characteristics of reservoirs and fractures.

[0093] For geological models of striped rock reservoirs with different threshold values, numerical well test forward modeling analysis was carried out using a fault-dissolved reservoir numerical well test model to establish well test curves under different threshold values.

[0094] The theoretical charts of well test curves with different threshold values ​​were compared and analyzed with the actual well test interpretation curves. The threshold value corresponding to the well test curve with the highest fitting accuracy was selected as the precise threshold value.

[0095] Specifically, once the precise threshold value is determined, parameters such as the volume of the reservoir can be accurately calibrated in the reserve calculation.

[0096] The present invention also provides an electronic device, comprising: a memory storing executable instructions; and a processor executing the executable instructions in the memory to implement the above-described method for determining the static reserve threshold value of fractured solution reservoirs.

[0097] The present invention also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method for determining the static reserve threshold value of a fractured solution reservoir.

[0098] To facilitate understanding of the solutions and effects of the embodiments of the present invention, three specific application examples are given below. Those skilled in the art should understand that these examples are merely for the purpose of understanding the present invention, and any specific details therein are not intended to limit the present invention in any way.

[0099] Example 1

[0100] Figure 3 A flowchart illustrating the steps of a method for determining the static reserve threshold value of a fractured solution reservoir according to an embodiment of the present invention is shown.

[0101] like Figure 3As shown, the method for determining the static reserve threshold value of the fractured solution reservoir includes: Step 101, obtaining the actual well test interpretation curve based on well test monitoring data; Step 102, establishing theoretical charts of well test curves with different threshold values ​​through seismic engraving; Step 103, comparing and analyzing the theoretical charts of well test curves with different threshold values ​​and the actual well test interpretation curves respectively, and taking the threshold value corresponding to the well test curve with the highest fitting accuracy as the precise threshold value.

[0102] The completed SHB3 well reached a depth of 8342.00m (incline) / 7913.78m (vertical), using acid fracturing completion. 751m³ of drilling mud was lost during drilling. 3 After being put into production, the well produced 53 tons of oil per day, confirming the presence of an oil and gas-bearing type with chaotic and strong reflection characteristics on the fault zone. However, the well has encountered many problems after production, such as low productivity and rapid energy decline, which are inconsistent with the early reserve estimates. Furthermore, due to the short production period and frequent adjustments to the nozzles during the trial production, the operating system was unstable, making it difficult to use seismic attributes and dynamic data to determine the calibration threshold value, resulting in significant uncertainty in the reserve calculation results.

[0103] Figure 4 A schematic diagram of the actual well test interpretation curve of the SHB3 well according to an embodiment of the present invention is shown.

[0104] Figure 5a , Figure 5b , Figure 5c , Figure 5d , Figure 5e The diagrams show the results of carving the grout with discontinuity threshold values ​​of 0.075, 0.103, 0.136, 0.165 and 0.195, respectively, according to an embodiment of the present invention.

[0105] Figure 6 A schematic diagram of a gravity-based multi-medium reservoir well test physical model is shown according to an embodiment of the present invention.

[0106] Well test monitoring data from well SHB3 was selected for analysis to clarify the characteristics of the SHB3 well test curves. The actual well test interpretation curves are shown below. Figure 4 As shown. Five threshold values ​​of 0.075, 0.103, 0.136, 0.165, and 0.195 were selected for discontinuity attribute threshold carving, and the carving results are shown below. Figures 5a-5e As shown. The results of determining the storage space under different threshold values ​​are based on, as follows: Figure 6 The well test physical model shown is used for numerical well test forward modeling studies of five energy attribute threshold values ​​in the SHB3 well area.

[0107] Figure 7a , Figure 7b , Figure 7c , Figure 7d , Figure 7e Schematic diagrams of numerical well test grid models of fractured and cavitated bodies with discontinuity threshold values ​​of 0.075, 0.103, 0.136, 0.165, and 0.195, respectively, according to an embodiment of the present invention, are shown.

[0108] Figure 8 A schematic diagram showing the fitting of well test curves with different threshold values ​​and actual test curves according to an embodiment of the present invention is illustrated.

[0109] According to such Figure 8 The fitting results of the well test curves with different threshold values ​​shown are compared with the actual test curves, and the attribute threshold value is clearly 0.103.

[0110] Using a defined threshold value, reserve calculations were conducted, confirming a geological reserve of 900,000 tons for this well. It was determined that the well controls only one reservoir and has poor communication with deep fractured-vuggy formations, leading to an underestimation of initial production. The reserve calculation results are consistent with dynamic understanding, indicating their reliability.

[0111] This method establishes a reasonable threshold for fracture-cavity sculpting, guiding the sculpting of fracture-cavity bodies and the determination of reservoir volume, thereby guiding oilfield reserve calculation. The final calculation results show significantly improved accuracy in aligning with reservoir geological understanding, effectively guiding the preparation of rolling assessment schemes and reserve declarations. This application has achieved excellent development results and confirms that the method's understanding of the reservoir is consistent with actual reservoir conditions. The research findings have broad application prospects in similar oilfields both domestically and internationally.

[0112] Example 2

[0113] This disclosure provides an electronic device comprising: a memory storing executable instructions; and a processor executing the executable instructions in the memory to implement the aforementioned method for determining the static reserve threshold value of a fractured solution reservoir.

[0114] An electronic device according to an embodiment of the present disclosure includes a memory and a processor.

[0115] This memory is used to store non-transitory computer-readable instructions. Specifically, the memory may include one or more computer program products, which may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include random access memory (RAM) and / or cache memory. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, etc.

[0116] The processor may be a central processing unit (CPU) or other form of processing unit with 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 this disclosure, the processor is used to execute computer-readable instructions stored in the memory.

[0117] Those skilled in the art will understand that, in order to solve the technical problem of how to achieve a good user experience, this embodiment may also include well-known structures such as communication buses and interfaces, and these well-known structures should also be included within the protection scope of this disclosure.

[0118] For a detailed description of this embodiment, please refer to the corresponding descriptions in the foregoing embodiments, which will not be repeated here.

[0119] Example 3

[0120] This disclosure provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the method for determining the static reserve threshold value of a fractured solution reservoir.

[0121] A computer-readable storage medium according to embodiments of the present disclosure stores non-transitory computer-readable instructions. When these non-transitory computer-readable instructions are executed by a processor, all or part of the steps of the methods described in the foregoing embodiments of the present disclosure are performed.

[0122] The aforementioned computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROM and DVD), magneto-optical storage media (e.g., MO), magnetic storage media (e.g., magnetic tape or portable hard drive), media with built-in rewritable non-volatile memory (e.g., memory card), and media with built-in ROM (e.g., ROM cartridge).

[0123] Those skilled in the art should understand that the above description of the embodiments of the present invention is only intended to illustrate the beneficial effects of the embodiments of the present invention, and is not intended to limit the embodiments of the present invention to any of the examples given.

[0124] The various embodiments of the present invention have been described above. These descriptions are exemplary and not exhaustive, nor are they limited to the disclosed embodiments. Many modifications and variations will be apparent to those skilled in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A method for determining the static reserve threshold value of a fractured-solution reservoir, characterized in that, include: Based on well test monitoring data, obtain the actual well test interpretation curves; By using seismic engraving, theoretical charts of well test curves with different threshold values ​​are established; The theoretical charts of well test curves with different threshold values ​​are compared and analyzed with the actual well test interpretation curves. The threshold value corresponding to the well test curve with the highest fitting accuracy is the precise threshold value. Among them, the theoretical charts for well test curves with different threshold values ​​established through seismic engraving include: By using seismic sculpting, we can obtain distribution maps of storage groups with different threshold values; Numerical processing was performed on reservoirs with different threshold values ​​to establish theoretical charts of well test curves for different threshold values. Specifically, reservoirs with different threshold values ​​are quantified to establish theoretical charts of well test curves for different threshold values: A numerical well test model is established, including the seepage equations for porous media and fractured media, and the initial conditions, inner boundary conditions, and outer boundary conditions of the numerical well test model are determined. The numerical well test model is used to quantify reservoirs with different threshold values, and the numerical well test model is solved to obtain the bottom hole pressure at different threshold values. Establish theoretical charts of well test curves with different threshold values; The seepage equation for the fractured medium is as follows: （3） Among them, C t The overall compression coefficient is... This refers to the bottom flow pressure of the oil well. For reservoir thickness, The viscosity of crude oil within the fractured medium. k represents the porosity of the fractured medium. z Where r is the longitudinal permeability and r is the pressure propagation radius. Let be the planar permeability of the fractured medium, and z be the vertical distance of pressure propagation.

2. The method for determining the static reserve threshold value of a fractured-solution reservoir according to claim 1, wherein, Based on well test monitoring data, the actual well test interpretation curves obtained include: Based on the well test monitoring data, the derivative of the pressure change was calculated by recovering the pressure change through shut-in pressure testing. The relationship curves between pressure change, the derivative of pressure change, and time change are constructed respectively, which are the interpretation curves of the actual well test.

3. The method for determining the static reserve threshold value of fractured-solution reservoirs according to claim 2, wherein, The derivative of the pressure change is calculated using formula (1): （1） Among them, △ P ′ is the derivative of the pressure change, Δ P For the change in pressure, Δ t For the time change, n This represents the total number of data points calculated in segments. i This represents the number of data points calculated in segments.

4. The method for determining the static reserve threshold value of a fractured-solution reservoir according to claim 1, wherein, The seepage equation for the porous medium is: （2） Where, k z For vertical penetration rate, C t The overall compression coefficient is... Let r be the pressure of the porous medium, and r be the pressure propagation radius. Let z be the planar permeability of the porous medium, and z be the vertical distance of pressure propagation. Porosity of porous media This refers to the viscosity of crude oil within the porous medium.

5. The method for determining the static reserve threshold value of a fractured solution reservoir according to claim 1, wherein, The initial conditions are: （4） The inner boundary conditions are as follows: （5） The outer boundary conditions are as follows: （6） Where h is the thickness, P is the pressure, and P i k represents the original formation pressure. r k represents the lateral penetration rate. z Φ is the longitudinal permeability, μ is the porosity, and C is the viscosity. t r is the overall compression coefficient. e Let s be the distance to the outer boundary, s be the skin coefficient, C be the wellbore reservoir coefficient, and Q be the production rate. This refers to the bottom flow pressure of the oil well. For the pressure of the fracture medium, The viscosity of crude oil within the fractured medium. For crack permeability, r w Where is the wellbore radius. The thickness of the porous medium.

6. An electronic device, characterized in that, The electronic device includes: Memory, which stores executable instructions; A processor that executes the executable instructions in the memory to implement the method for determining the static reserve threshold value of a fractured solution reservoir as described in any one of claims 1-5.

7. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, implements the method for determining the static reserve threshold value of a fractured solution reservoir as described in any one of claims 1-5.

Images