Methods and related equipment for determining oil column height in carbonate rock fractured-collapse reservoirs

Abstract

A method and related apparatus for determining the oil column height in carbonate fault-controlled reservoirs include: collecting temperature and geological data from wells in the carbonate fault-controlled reservoir to obtain static temperature test data at the initial production stage and flow temperature test data for the scheduled production date; establishing a coupled thermodynamic physical model of the wellbore and reservoir in the fault-controlled reservoir; establishing a wellbore temperature distribution model using energy conservation, and calculating the bottom hole temperature from the wellhead temperature; then establishing a one-dimensional vertical thermodynamic model of the reservoir, and iteratively calculating the height from the bottom of the well to the bottom of the reservoir; adding the calculated height to the height of the drilling and development formations to obtain the oil column height in the carbonate fault-controlled reservoir. This invention, based on a coupled thermodynamic model of the wellbore and reservoir, determines a method suitable for calculating the oil column height in carbonate fault-controlled reservoirs. It can directly calculate the oil column height of fault-controlled reservoirs using wellhead temperature data, solving the problem of large calculation errors in conventional static description methods.

Description

Technical Field

[0001] This invention belongs to the field of oil and gas field development and reservoir engineering technology, specifically relating to a method and related apparatus for determining the oil column height in carbonate rock fractured solution reservoirs. Background Technology

[0002] In recent years, the proportion of carbonate oil and gas reservoirs in my country's oil and gas resources has been increasing, and their development potential has been gradually growing. They are a key area for future oil and gas reserve growth. How to develop such reservoirs in a reasonable and effective manner has become an urgent problem to be solved in the petroleum industry.

[0003] Significant oil and gas breakthroughs have been achieved in the ultra-deep carbonate rock field below 7,000 meters. The carbonate rock fault-dissolved oil reservoirs in the Tarim Basin have developed a large number of strike-slip faults. After being transformed by multiple phases of karstification, they have formed fracture-vuggy reservoirs. The reservoir types are mainly large caves, supplemented by fractures and pores. They are characterized by overall oil-bearing along the fault zone and uneven enrichment. The discovery of this type of oil and gas reservoir demonstrates a new model of ultra-deep oil and gas enrichment.

[0004] Oil column height refers to the vertical distance from the oil-water interface to the top boundary of the reservoir. The magnitude of the oil column height is a crucial parameter for determining reservoir size, reservoir research, and production management. For layered oil and gas reservoirs, the oil column height parameter can be calculated by determining structural high points, the oil-water interface, and trap overflow points. However, for carbonate rock fractured-karst reservoirs, the internal heterogeneity is extremely high, and the presence of large caverns that easily lead to air venting and leakage makes it impossible to drill through the reservoir bottom to obtain the true oil column height. Currently, the industry commonly uses static seismic etching to determine the oil column height, but these static methods often have low accuracy and inaccurate results. Summary of the Invention

[0005] The purpose of this invention is to provide a method and related apparatus for determining the oil column height in carbonate rock fractured solution reservoirs, so as to solve the problem that static methods are often not very accurate and the results are inaccurate.

[0006] To achieve the above objectives, the present invention adopts the following technical solution:

[0007] In a first aspect, the present invention provides a method for determining the oil column height in carbonate rock fractured-dissolved oil reservoirs, comprising:

[0008] Collect temperature test data and geological data of oil wells in carbonate rock fault-controlled oil reservoirs to obtain static temperature test data in the early stage of production and flow temperature test data on the scheduled production date;

[0009] Based on the above data, a coupled thermodynamic physical model of wellbore-reservoir in fractured solution reservoirs was established.

[0010] By utilizing the law of conservation of energy, a wellbore temperature distribution model is established, and the bottom hole temperature is calculated from the wellhead temperature.

[0011] Then, a one-dimensional vertical thermodynamic model of the reservoir was established. The height from the bottom of the well to the bottom of the reservoir was obtained by iterative calculation using two temperature test data.

[0012] Finally, the calculated height is added to the height of the drilling and excavation layers to obtain the oil column height of the carbonate rock fractured solution reservoir.

[0013] Furthermore, based on the above data, a coupled thermodynamic physical model of wellbore-reservoir in fractured solution reservoirs is established:

[0014] Based on the fact that the heat transfer from the oil pipe axis to the formation during the crude oil extraction process is a steady-state one-dimensional heat transfer, the amount of heat change per unit length of the micro-element is obtained.

[0015]

[0016] In the formula, dQ is the change in heat along the length of the infinitesimal element, kcal / h; dZ is the unit length, m; To is the temperature of the oil in the pipeline, °C; T r - Formation temperature, °C; R t - Thermal resistance from the tubing axis to the formation, kcal / (h·m·℃).

[0017] Furthermore, by utilizing the law of conservation of energy, a wellbore temperature distribution model is established, and the bottom-hole temperature is calculated from the wellhead temperature.

[0018] Based on the heat loss of crude oil during wellbore production, the total thermal resistance R of crude oil in the wellbore is calculated. t and overall heat transfer coefficient U t ;

[0019] Taking the bottom of the well as the origin of the coordinate axis, and pointing the Z-axis upwards, the infinitesimal element at a distance of (n-1)dZ from the bottom of the well satisfies energy conservation, and the heat Q emanating from the infinitesimal element is transferred out. ndZ It should be equal to the heat Q transferred into the infinitesimal element. (n-1)dZ Subtract the heat lost through radial heat transfer The expression is described as follows:

[0020]

[0021]

[0022]

[0023]

[0024] In the formula, m o(n-1) - Mass of crude oil in micro-elemental volume, kg; c o- Crude oil specific heat capacity, J / (kg·℃); T o(n-1) - Temperature at the inlet of the crude oil micro-element, °C; T on -Temperature at the crude oil micro-element outlet, °C; U t - The total heat transfer coefficient from the tubing axis to the formation, W / (m·℃).

[0025] The expression for the temperature of the stratum where the micro-element is located is obtained from the ground temperature and the geothermal gradient:

[0026]

[0027] In the formula, T d(n-1) - Temperature of the stratum where the micro-element is located, °C; T d - Surface temperature, °C; Z0 - Geothermal gradient, °C / m; h - Wellbore depth, m;

[0028] Substituting the parameters obtained above into the formula for the change in heat of the infinitesimal element, we obtain the change in heat of the infinitesimal element flowing from the wellbore to the bottom of the well with respect to T. o(n) The iterative process is implemented using VB programming to perform iterative calculations, dividing the wellbore depth into n infinitesimal elements, and calculating the bottom hole temperature from the wellhead temperature:

[0029] .

[0030] Furthermore, the total thermal resistance R of crude oil heat loss in the wellbore was calculated. t and overall heat transfer coefficient U t :

[0031] Collect well data to obtain well parameter tables, including the total thermal resistance R from the tubing axis to the formation. t It mainly consists of 5 parts: oil pipe thermal resistance R o Circulatory and radiative thermal resistance R inside the annulus h The thermal resistance R of the bushing g Cement ring thermal resistance R s Formation thermal resistance R d ;

[0032]

[0033] In the formula, λ 管 - Thermal conductivity of the tubing, W / (m·℃); r1 - Inner diameter of the tubing, m; r2 - Outer diameter of the tubing, m; h c - Convection heat transfer coefficient, W / (m·℃); h r- Radiative heat transfer coefficient, W / (m·℃); λ-tube-to-tube thermal conductivity, W / (m·℃); r 3-Casing inner diameter, m; r4 - casing outer diameter, m; λ cement sheath - cement sheath thermal conductivity, W / (m·℃); r5 - cement sheath outer diameter, m; λ formation - formation thermal conductivity, W / (m·℃); r6 - formation outer diameter, m;

[0034] Summing up the thermal resistances in the calculation parameter table yields the complete expression for the thermal resistance from the tubing axis to the formation:

[0035]

[0036]

[0037] The overall heat transfer coefficient U from the tubing axis to the formation is obtained from the above thermal resistance expression. t :

[0038] .

[0039] Furthermore, a one-dimensional vertical thermodynamic model of the reservoir was then established. By substituting temperature test data from two separate tests into the heat transfer model, the height from the bottom of the well to the bottom of the reservoir was calculated iteratively.

[0040] In the vertical direction, considering one-dimensional heat conduction in the Z direction, the energy conservation equation is obtained, leading to the thermodynamic differential equation.

[0041]

[0042] Initial conditions:

[0043] Boundary conditions:

[0044] In the formula, ξ(T) — thermal diffusivity, m 2 / s; T—temperature, °C; Z—distance from any point in the vertical direction to the oil-water interface, m; α—geothermal gradient, °C / m;

[0045] The differential equation is solved by separating variables, introducing eigenvalues, performing inverse Fourier transform, and introducing an error function.

[0046] Solving for the semi-analytical expression of the temperature field, for point A at the bottom of the wellbore, we have:

[0047]

[0048]

[0049] In the formula, h A — Height from the oil-water interface to the bottom of the well (m); T A1 —Initial bottom hole temperature (°C); T A2 —Second test wellbore temperature (°C); TRo —Temperature at the oil-water interface (°C); t—Production time (d);

[0050] The model is programmed in VB to obtain the bottom hole temperature T from the two tests obtained from the wellbore thermodynamic model. A1 T A2 Substitute the values ​​into the formula and perform iterative calculations to obtain the height h from the oil-water interface to the bottom of the well. A .

[0051] Furthermore, the differential equation is solved by separating variables, introducing eigenvalues, performing inverse Fourier transform, and introducing an error function:

[0052] make Separating variables from the differential equation yields:

[0053]

[0054] Introducing eigenvalue -β 2 The above formula is written as

[0055]

[0056] By introducing characteristic functions into the above equation and superimposing them, we obtain the general solution.

[0057]

[0058] After performing an inverse Fourier transform on the above equation, substituting the boundary conditions, we obtain...

[0059]

[0060] Introducing an error function into the above equation, the equation is rewritten as follows:

[0061]

[0062] Series expansion of the error function: .

[0063] Furthermore, the calculated height is added to the height of the drilling and excavation layers to obtain the oil column height of the carbonate rock fractured solution reservoir:

[0064] The obtained height h A Add the known height h of the target formation penetrated from the drilling and completion data. c Adding them together, we get the oil column height h of the reservoir.

[0065] .

[0066] Secondly, the present invention provides a system for determining the oil column height of carbonate rock fractured-dissolved reservoirs, comprising:

[0067] The data acquisition module is used to collect temperature test data and geological data of oil wells in carbonate rock fault-controlled oil reservoirs, and to obtain static temperature test data in the early stage of production and flow temperature test data on the scheduled production date.

[0068] The model building module is used to build a coupled thermodynamic and physical model of wellbore-reservoir in fractured solution reservoirs based on the above data.

[0069] The processing module is used to establish a wellbore temperature distribution model using energy conservation, and calculate the bottom-hole temperature from the wellhead temperature; then, a one-dimensional vertical thermodynamic model of the reservoir is established, and the height from the bottom of the well to the bottom of the reservoir is obtained by iterative calculation using temperature test data from two tests substituted into the heat transfer model; finally, the calculated height is added to the height of the drilling and development formation to obtain the oil column height of the carbonate rock fractured solution reservoir.

[0070] Thirdly, the present invention provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of a method for determining the oil column height of a carbonate rock fractured solution reservoir.

[0071] Fourthly, the present invention provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of a method for determining the oil column height in a carbonate rock fractured solution reservoir.

[0072] Compared with the prior art, the present invention has the following technical effects:

[0073] This invention establishes a method for calculating the oil column height in carbonate rock fractured solution reservoirs based on a thermodynamic model of wellbore-reservoir coupling. This method can directly calculate the oil column height in fractured solution reservoirs using wellhead temperature data, solving the problem of large calculation errors in conventional static description methods. It provides technical support for judging the reserve scale of carbonate rock fractured solution reservoirs and for further development. Attached Figure Description

[0074] Figure 1 This is a flowchart of the present invention;

[0075] Figure 2 This is a diagram of a coupled thermodynamic and physical model of a wellbore-reservoir in a fractured solution reservoir.

[0076] Figure 3 This is a schematic diagram of heat loss during crude oil production in wellbore;

[0077] Figure 4 This is a schematic diagram of energy conservation in a micro-element;

[0078] Figure 5 This is a schematic diagram of a reservoir thermodynamic micro-element heat transfer model. Detailed Implementation

[0079] The present invention will be further described below with reference to the accompanying drawings:

[0080] Please see Figures 1 to 5 Methods for determining the oil column height in carbonate rock fractured solution reservoirs:

[0081] Collect temperature test data and geological data of a certain oil well in a carbonate rock fractured solution reservoir to obtain static temperature test data at the initial stage of production and flow temperature test data at the scheduled production date (generally six months or one year).

[0082] Establish a coupled thermodynamic physical model of wellbore-reservoir in fractured solution reservoirs;

[0083] By utilizing the law of conservation of energy, a wellbore temperature distribution model is established, and the bottom hole temperature can be calculated from the wellhead temperature.

[0084] Then, a one-dimensional vertical heat transfer model of the reservoir was established. The height from the bottom of the well to the bottom of the reservoir was obtained by iterative calculation using two temperature test data as input into the heat transfer model.

[0085] Finally, the calculated height is added to the height of the drilling and excavation layers to obtain the oil column height of the carbonate rock fractured solution reservoir.

[0086] The detailed steps of implementing this invention are as follows:

[0087] Step S10: Collect temperature test data, geological data and drilling and completion data of a certain oil well in a carbonate rock fault-controlled oil reservoir, including static temperature test and flow temperature test, and obtain the wellhead temperature table, as shown in Table 1.

[0088] Step S20: Based on the principle of energy conservation of micro-element, and according to the fact that the heat transfer from the oil pipe axis to the formation during the crude oil extraction process is a steady-state one-dimensional heat transfer, the change in heat per unit length of the micro-element can be obtained.

[0089]

[0090] In the formula, dQ represents the change in heat along the length of the infinitesimal element, in kcal / h; d Z -Unit length, m; T o - Temperature of the oil in the pipeline, °C; T r - Formation temperature, °C; R t - Thermal resistance from the tubing axis to the formation, kcal / (h·m·℃).

[0091] Table 1 Sample table of temperature test data for an oil well in a carbonate rock fractured karst reservoir.

[0092]

[0093] Step S30: Based on the heat loss that occurs during crude oil production in the wellbore, such as... Figure 1 As shown, the total thermal resistance R of crude oil heat loss in the wellbore is calculated. t and overall heat transfer coefficient U t ;

[0094] Step S301: Collect well data to obtain the well parameter table, as shown in Table 2, which shows the total thermal resistance R from the tubing axis to the formation. t It mainly consists of 5 parts: oil pipe thermal resistance (R) o ), convective and radiative thermal resistance within the annulus (R h ), sleeve thermal resistance (R) g ), Cement ring thermal resistance (R) s ), formation thermal resistance (R) d );

[0095]

[0096] In the formula, λ 管 - Thermal conductivity of the tubing, W / (m·℃); r1 - Inner diameter of the tubing, m; r2 - Outer diameter of the tubing, m; h c - Convection heat transfer coefficient, W / (m·℃); h r - Radiative heat transfer coefficient, W / (m·℃); λ 套管 - Thermal conductivity of the casing, W / (m·℃); r3 - Inner diameter of the casing, m; r4 - Outer diameter of the casing, m; λ 水泥环 - Thermal conductivity of the cement ring, W / (m·℃); r5 - Outer diameter of the cement ring, m; λ 地层 - Formation thermal conductivity, W / (m·℃); r6 - Formation outer diameter, m.

[0097] Table 2 Sample Table of Thermal Resistance Calculation Parameters for Oil Wells in a Carbonate Rock Fault-Degradable Reservoir

[0098]

[0099] Step S302: Sum the thermal resistances on the calculation parameter table to obtain the complete expression for the thermal resistance from the tubing axis to the formation.

[0100]

[0101]

[0102] Step S303: The overall heat transfer coefficient from the tubing axis to the formation can be obtained using the above thermal resistance expression. U t ;

[0103]

[0104] Step S40: Taking the bottom of the well as the origin of the coordinate axis, and letting the Z-axis point upwards, the infinitesimal element at a distance of (n-1)dZ from the bottom of the well satisfies the following... Figure 2 According to the energy conservation shown, the heat Q(ndZ) transferred out of the infinitesimal element should be equal to the heat Q[(n-1)dZ] transferred into the infinitesimal element minus the heat lost through radial heat transfer, which can be described by the following expression;

[0105]

[0106]

[0107]

[0108]

[0109] In the formula, m o(n-1) - Mass of crude oil in micro-elemental volume, kg; c o - Crude oil specific heat capacity, J / (kg·℃); T o(n-1) - Temperature at the inlet of the crude oil micro-element, °C; T on -Temperature at the crude oil micro-element outlet, °C; U t - The total heat transfer coefficient from the tubing axis to the formation, W / (m·℃).

[0110] Step S50: Obtain the temperature expression of the stratum where the micro-element is located through the ground temperature and the geothermal gradient;

[0111]

[0112] In the formula, T d(n-1) -Temperature of the stratum where the micro-element is located, °C; T d - Surface temperature, °C; Z0 - Geothermal gradient, °C / m; h - Wellbore depth, m.

[0113] Step S60: Substitute the parameters obtained above into the formula for the change in heat of the infinitesimal element to obtain the flow of the infinitesimal element from the wellbore to the bottom of the well with respect to T. o(n) The iterative process is used, and VB programming is used for iterative calculation. The wellbore depth is divided into n micro-elements, and the bottom hole temperature is calculated from the wellhead temperature.

[0114]

[0115] Step S70: Establish a reservoir thermodynamic model. In the vertical direction, the micro-element considers one-dimensional heat conduction in the Z direction, such as... Figure 3 As shown, the conservation of energy yields the thermodynamic differential equation;

[0116]

[0117] Initial conditions:

[0118] Boundary conditions:

[0119] In the formula, ξ(T) — thermal diffusivity, m 2 / s; T —Temperature, °C; Z—Distance from any point in the vertical direction to the oil-water interface, m; α—Geothermal gradient, °C / m.

[0120] Step S80: Solve the differential equation by separating variables, introducing eigenvalues, performing inverse Fourier transform, and introducing error functions.

[0121] Step S801, let The differential equation is obtained by separating variables.

[0122]

[0123] Step S802, Introduce eigenvalues ​​- β 2 The above formula can be written as:

[0124]

[0125] Step S803: Introduce characteristic functions into the above equation and superimpose them to obtain the general solution;

[0126]

[0127] Step S804: After performing an inverse Fourier transform on the above equation, substitute the boundary conditions to obtain the result.

[0128]

[0129] Step S805: Introduce an error function into the above equation, then the above equation can be written as:

[0130]

[0131] Series expansion of the error function:

[0132] Step S90: Solve for the semi-analytical expression of the temperature field. For the bottom of the well (point A), we have:

[0133]

[0134]

[0135] In the formula, h A— Height from the oil-water interface to the bottom of the well (m); T A1 —Initial bottom hole temperature (°C); T A2 —Second test of bottom hole temperature (°C); T Ro —Temperature at the oil-water interface (°C); t—Production time (d).

[0136] Step S100: Perform VB programming on the model to obtain the bottom hole temperature T from the two tests obtained from the wellbore thermodynamic model. A1 T A2 Substitute the values ​​into the formula and perform iterative calculations to obtain the height h from the oil-water interface to the bottom of the well. A ;

[0137] Step S110: Obtain the height h A In addition to the known height of the target formation penetrated from the drilling and completion data. h c By adding them together, the oil column height of the reservoir can be obtained. h ;

[0138]

[0139] Example:

[0140] Step S10: Collect temperature test data, geological data and drilling and completion data of well XX1 in carbonate rock fractured solution reservoir, obtain relevant production data, and obtain the temperature gradient table of the two tested oil wells, as shown in Table 3.

[0141] Table 3 Temperature test data of well H1 in carbonate rock fractured karst reservoir.

[0142]

[0143] Step S20: Based on the principle of energy conservation of micro-element, and according to the fact that the heat transfer from the oil pipe axis to the formation during the crude oil extraction process is a steady-state one-dimensional heat transfer, the change in heat per unit length of the micro-element can be obtained.

[0144]

[0145] Step S30: Based on the heat loss of crude oil during the wellbore production process, calculate the total thermal resistance R of the crude oil in the wellbore. t and overall heat transfer coefficient U t ;

[0146] Step S301: Collect drilling data and relevant survey data for well XX1 to obtain the oil well thermal resistance calculation parameter table, as shown in Table 4, which shows the total thermal resistance R from the tubing axis to the formation. t It mainly consists of 5 parts: oil pipe thermal resistance (R) o), convective and radiative thermal resistance within the annulus (R h ), sleeve thermal resistance (R) g ), Cement ring thermal resistance (R) s ), formation thermal resistance (R) d );

[0147]

[0148] Table 4. Calculation parameters for thermal resistance of well H1 in carbonate rock fractured karst reservoir.

[0149]

[0150] Step S302: Sum the thermal resistances on the calculation parameter table to obtain the complete expression for the thermal resistance from the tubing axis to the formation.

[0151]

[0152] Step S303: The total heat transfer coefficient from the tubing axis to the formation can be obtained through the complete thermal resistance expression obtained from the above steps. U t ;

[0153]

[0154] Step S40: Divide the wellbore depth of 7400m into n corresponding micro-elements, and consider the relationship between T and... o(n) The iterative process is performed, substituting the wellhead temperatures (35.6℃ and 36.2℃) from the two tests and related parameters into the calculation to obtain the bottom hole temperature T from the two tests. A1 and T A2 .

[0155]

[0156]

[0157] Step S50: Establish a reservoir thermodynamic model. In the vertical direction, the micro-element considers one-dimensional heat conduction in the Z direction to obtain energy conservation and thermodynamic differential equations.

[0158]

[0159] Initial conditions:

[0160] Boundary conditions:

[0161] Step S60: Solve for the semi-analytical expression of the temperature field. For the bottom of the well (point A), we have:

[0162]

[0163] Step S60: Program and calculate the model, with the default oil-water interface temperature. T Ro The bottom hole temperature T, calculated from the wellbore thermodynamic model, is 180℃. A1 T A2 Substitute the values ​​into the formula and perform iterative calculations to obtain the height h from the oil-water interface to the bottom of the well. A .

[0164]

[0165] Step S70: Known drilling height h of the target strata c The value is 12m. Add the calculated height h of the target strata to h. c The oil column height h of well XX1 was obtained;

[0166] .

[0167] In another embodiment of the present invention, an oil column height determination system for carbonate rock fractured-dissolved reservoirs is provided, which can be used to implement the above-mentioned method for determining the oil column height of carbonate rock fractured-dissolved reservoirs. Specifically, the system includes:

[0168] The data acquisition module is used to collect temperature test data and geological data of oil wells in carbonate rock fault-controlled oil reservoirs, and to obtain static temperature test data in the early stage of production and flow temperature test data on the scheduled production date.

[0169] The model building module is used to build a coupled thermodynamic and physical model of wellbore-reservoir in fractured solution reservoirs based on the above data.

[0170] The processing module is used to establish a wellbore temperature distribution model using energy conservation, and calculate the bottom-hole temperature from the wellhead temperature; then, a one-dimensional vertical thermodynamic model of the reservoir is established, and the height from the bottom of the well to the bottom of the reservoir is obtained by iterative calculation using temperature test data from two tests substituted into the heat transfer model; finally, the calculated height is added to the height of the drilling and development formation to obtain the oil column height of the carbonate rock fractured solution reservoir.

[0171] The module division in this embodiment of the invention is illustrative and represents only one logical functional division. In actual implementation, other division methods may be used. Furthermore, the functional modules in the various embodiments of the invention can be integrated into a single processor, exist as separate physical entities, or be integrated into a single module. The integrated modules described above can be implemented in hardware or as software functional modules.

[0172] In another embodiment of the present invention, a computer device is provided, comprising a processor and a memory. The memory stores a computer program, which includes program instructions. The processor executes the program instructions stored in the computer storage medium. The processor may be a Central Processing Unit (CPU), or other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. It is the computing and control core of the terminal, suitable for implementing one or more instructions, specifically suitable for loading and executing one or more instructions in the computer storage medium to achieve a corresponding method flow or corresponding function. The processor described in this embodiment of the present invention can be used in the operation of a method for determining the oil column height in carbonate rock fractured solution reservoirs.

[0173] In another embodiment of the present invention, a storage medium is provided, specifically a computer-readable storage medium (Memory), which is a memory device in a computer device used to store programs and data. It is understood that the computer-readable storage medium here can include both the built-in storage medium in the computer device and extended storage media supported by the computer device. The computer-readable storage medium provides storage space that stores the terminal's operating system. Furthermore, the storage space also stores one or more instructions suitable for loading and execution by a processor. These instructions can be one or more computer programs (including program code). It should be noted that the computer-readable storage medium here can be a high-speed RAM memory or a non-volatile memory, such as at least one disk storage device. The processor can load and execute one or more instructions stored in the computer-readable storage medium to implement the corresponding steps of the method for determining the oil column height of carbonate rock fractured solution reservoirs in the above embodiments.

[0174] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can 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, etc.) containing computer-usable program code.

[0175] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can 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, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0176] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0177] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0178] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit it. Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art should understand that modifications or equivalent substitutions can still be made to the specific implementation of the present invention. Any modifications or equivalent substitutions that do not depart from the spirit and scope of the present invention should be covered within the scope of protection of the claims of the present invention.

Claims

1. A method for determining the oil column height in carbonate rock fault-collapse reservoirs, characterized in that, include: Collect temperature test data and geological data of oil wells in carbonate rock fault-controlled oil reservoirs to obtain static temperature test data in the early stage of production and flow temperature test data on the scheduled production date; Based on the above data, a coupled thermodynamic physical model of wellbore-reservoir in fractured solution reservoirs was established. By utilizing the law of conservation of energy, a wellbore temperature distribution model is established, and the bottom hole temperature is calculated from the wellhead temperature. Then, a one-dimensional vertical thermodynamic model of the reservoir was established. The height from the bottom of the well to the bottom of the reservoir was obtained by iterative calculation using two temperature test data. Finally, the calculated height is added to the height of the drilling and development system to obtain the oil column height of the carbonate rock fractured solution reservoir. Then, a one-dimensional vertical thermodynamic model of the reservoir was established. By substituting the temperature test data from two tests into the heat transfer model, the height from the bottom of the well to the bottom of the reservoir was obtained through iterative calculation. In the vertical direction, the infinitesimal element considers one-dimensional heat conduction in the Z direction, and obtains the thermodynamic differential equation through energy conservation; Initial conditions: Boundary conditions: In the formula, — Thermal diffusivity, m 2 / s; T—temperature, °C; Z—distance from any point in the vertical direction to the oil-water interface, m; α—geothermal gradient, °C / m; The differential equation is solved by separating variables, introducing eigenvalues, performing inverse Fourier transform, and introducing an error function. Solving for the semi-analytical expression of the temperature field, for point A at the bottom of the wellbore, we have: where h A — height of the oil-water interface to the bottom of the well m; T A1 — initial bottom hole temperature °C; T A2 — second test bottom hole temperature °C; T Ro — temperature at the oil-water interface °C; t— production time d; VB programming of the model will be calculated from the wellbore thermodynamic model of the two test bottom hole temperature T A1 , T A2 Substitute formula for iterative calculation, to calculate the height of the oil-water interface to the bottom of the formula h A .

2. The method for determining the oil column height in carbonate rock fractured-dissolved oil reservoirs according to claim 1, characterized in that, Based on the above data, a coupled thermodynamic physical model of wellbore-reservoir in fractured solution reservoirs is established: Based on the fact that the heat transfer from the oil pipe axis to the formation during the crude oil extraction process is a steady-state one-dimensional heat transfer, the amount of heat change per unit length of the micro-element is obtained. In the formula, dQ is the change in heat along the length of the infinitesimal element, kcal / h; dZ is the unit length, m; T o - Temperature of the oil in the pipeline, °C; T r - Formation temperature, °C; R t - Thermal resistance from the tubing axis to the formation, kcal / (h·m·℃).

3. The method for determining the oil column height in carbonate rock fractured-dissolved oil reservoirs according to claim 2, characterized in that, Using the law of conservation of energy, a wellbore temperature distribution model is established, and the bottom hole temperature is calculated from the wellhead temperature. Based on the heat loss of crude oil during wellbore production, the total thermal resistance R of crude oil in the wellbore is calculated. t and overall heat transfer coefficient U t ; Taking the bottom of the well as the origin of the coordinate axis, and pointing the Z-axis upwards, the infinitesimal element at a distance of (n-1)dZ from the bottom of the well satisfies energy conservation, and the heat Q emanating from the infinitesimal element is transferred out. ndZ It should be equal to the heat Q transferred into the infinitesimal element. (n-1)dZ Subtract the heat lost through radial heat transfer The expression is described as follows: In the formula, m o(n-1) - Mass of crude oil in micro-elemental volume, kg; c o - Crude oil specific heat capacity, J / (kg·℃); T o(n-1) - Temperature at the inlet of the crude oil micro-element, °C; T on -Temperature at the crude oil micro-element outlet, °C; U t - The total heat transfer coefficient from the tubing axis to the formation, W / (m·℃); The expression for the temperature of the stratum where the micro-element is located is obtained from the ground temperature and the geothermal gradient: In the formula, T d(n-1) - Temperature of the stratum where the micro-element is located, °C; T d - Surface temperature, °C; Z0 - Geothermal gradient, °C / m; h - Wellbore depth, m; Substituting the parameters obtained above into the formula for the change in heat of the infinitesimal element, we obtain the change in heat of the infinitesimal element flowing from the wellbore to the bottom of the well with respect to T. o(n) The iterative process is implemented using VB programming to perform iterative calculations, dividing the wellbore depth into n infinitesimal elements, and calculating the bottom hole temperature from the wellhead temperature. : 。 4. The method for determining the oil column height in carbonate rock fractured-dissolved oil reservoirs according to claim 3, characterized in that, The total thermal resistance R of crude oil heat loss in the wellbore was calculated. t and overall heat transfer coefficient U t : Collect well data to obtain well parameter tables, including the total thermal resistance R from the tubing axis to the formation. t It mainly consists of 5 parts: oil pipe thermal resistance R o Circulatory and radiative thermal resistance R inside the annulus h The thermal resistance R of the bushing g Cement ring thermal resistance R s Formation thermal resistance R d ; In the formula, λ 管 - Thermal conductivity of the tubing, W / (m·℃); r1 - Inner diameter of the tubing, m; r2 - Outer diameter of the tubing, m; h c - Convection heat transfer coefficient, W / (m·℃); h r- Radiative heat transfer coefficient, W / (m·℃); λ-tube-to-tube thermal conductivity, W / (m·℃); r 3- Casing inner diameter, m; r4 - casing outer diameter, m; λ cement sheath - cement sheath thermal conductivity, W / (m·℃); r5 - cement sheath outer diameter, m; λ formation - formation thermal conductivity, W / (m·℃); r6 - formation outer diameter, m; Summing up the thermal resistances in the calculation parameter table yields the complete expression for the thermal resistance from the tubing axis to the formation: The overall heat transfer coefficient U from the tubing axis to the formation is obtained from the above thermal resistance expression. t : 。 5. The method for determining the oil column height in carbonate rock fractured-dissolved oil reservoirs according to claim 1, characterized in that, The differential equation is solved by separating variables, introducing eigenvalues, performing inverse Fourier transform, and introducing an error function. make Separating variables from the differential equation yields: Introducing eigenvalue -β 2 The above formula is written as By introducing characteristic functions into the above equation and superimposing them, we obtain the general solution. After performing an inverse Fourier transform on the above equation, substituting the boundary conditions, we obtain... Introducing an error function into the above equation, the equation is rewritten as follows: Series expansion of the error function: .

6. The method for determining the oil column height in carbonate rock fractured-dissolved reservoirs according to claim 1, characterized in that, Finally, the calculated height is added to the height of the drilling and excavation layers to obtain the oil column height of the carbonate rock fractured solution reservoir: The obtained height h A Add the known height h of the target formation penetrated from the drilling and completion data. c Adding them together, we get the oil column height h of the reservoir. 。 7. A system for determining the oil column height in carbonate rock fault-collapse reservoirs, characterized in that, include: The data acquisition module is used to collect temperature test data and geological data of oil wells in carbonate rock fault-controlled oil reservoirs, and to obtain static temperature test data in the early stage of production and flow temperature test data on the scheduled production date. The model building module is used to build a coupled thermodynamic and physical model of wellbore-reservoir in fractured solution reservoirs based on the above data. The processing module is used to establish a wellbore temperature distribution model using energy conservation, and calculate the bottom-hole temperature from the wellhead temperature; then, a one-dimensional vertical thermodynamic model of the reservoir is established, and the height from the bottom of the well to the bottom of the reservoir is obtained by iterative calculation using temperature test data from two tests substituted into the heat transfer model; finally, the calculated height is added to the height of the drilling and development formation to obtain the oil column height of the carbonate rock fractured solution reservoir. Then, a one-dimensional vertical thermodynamic model of the reservoir was established. By substituting the temperature test data from two tests into the heat transfer model, the height from the bottom of the well to the bottom of the reservoir was obtained through iterative calculation. In the vertical direction, the infinitesimal element considers one-dimensional heat conduction in the Z direction, and obtains the thermodynamic differential equation through energy conservation; Initial conditions: Boundary conditions: In the formula, — Thermal diffusivity, m 2 / s; T—temperature, °C; Z—distance from any point in the vertical direction to the oil-water interface, m; α—geothermal gradient, °C / m; The differential equation is solved by separating variables, introducing eigenvalues, performing inverse Fourier transform, and introducing an error function. Solving for the semi-analytical expression of the temperature field, for point A at the bottom of the wellbore, we have: In the formula, h A — Height from the oil-water interface to the bottom of the well (m); T A1 —Initial bottom hole temperature (°C); T A2 —Second test wellbore temperature (°C); T Ro —Temperature at the oil-water interface (°C); t—Production time (d); The model is programmed in VB to obtain the bottom hole temperature T from the two tests obtained from the wellbore thermodynamic model. A1 T A2 Substitute the values ​​into the formula and perform iterative calculations to obtain the height h from the oil-water interface to the bottom of the well. A .

8. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the steps of the method for determining the oil column height of a carbonate rock fractured solution reservoir as described in any one of claims 1 to 6.

9. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by the processor, it implements the steps of the method for determining the oil column height of a carbonate rock fractured solution reservoir as described in any one of claims 1 to 6.

Images