Method and device for determining radius range of natural gas filling pore throat in sandstone reservoir

Abstract

The invention discloses a method and device for determining the radius range of a natural gas filling pore throat in a sandstone reservoir, and the method comprises the steps: determining the paleo-power of natural gas filling in the sandstone reservoir; the paleo-resistance of natural gas filling in the sandstone reservoir is determined; determining the maximum filling resistance and the minimum filling resistance during natural gas filling in the sandstone reservoir according to the paleo-power and the paleo-resistance; and correspondingly determining a minimum pore throat radius and a maximum pore throat radius according to the maximum charging resistance and the minimum charging resistance so as to obtain a pore throat radius range. According to the method, the radius range of the natural gas filling pore throat in the sandstone reservoir is accurately calculated by calculating the natural gas filling power and resistance, a basis is provided for research, exploration and development of the natural gas filling process, and the efficiency of the natural gas filling process and exploration and development is improved.

Description

Technical Field

[0001] The present application belongs to the technical field of oil and gas geological exploration, and in particular relates to a method and device for determining the radius range of a natural gas injection pore throat in a sandstone reservoir. Background Art

[0002] The size of pore throats in sandstone reservoirs influences the storage and seepage capacity of natural gas, affecting both its differential enrichment and its recovery. When natural gas is injected into sandstone reservoirs, it experiences capillary resistance of varying magnitude, causing it to fill in different pore throats. Identifying the pore throat radius range within which natural gas fills is crucial for understanding the natural gas accumulation process and developing oil and gas recovery plans based on the pore throat size where natural gas accumulates.

[0003] Limited research exists on the range of pore throat radii for natural gas injection, primarily focusing on characterizing reservoir pore throat size and determining the lower limit of natural gas injection into pore throats. Methods for characterizing reservoir pore throat size include: High-pressure mercury injection (HMI) experiments, which reveal the distribution of different pore throat radii; High-pressure mercury injection (HMI) and nuclear magnetic resonance (NMR) can both reveal reservoir pore structure; the NMR transverse relaxation time (T2) correlates with the throat radius (r), allowing the T2 relaxation time to be converted into pore throat radius for quantitative pore structure assessment. Methods for determining the lower limit of natural gas injection into sandstone reservoirs primarily include empirical statistics, theoretical or numerical calculations, and experimental simulations. Among them, the empirical statistical method has the advantages of being simple and fast, but it is difficult to determine whether it reflects the actual situation of the reservoir due to the lack of sufficient theoretical support; the theoretical calculation is based on mechanical equilibrium, using numerical simulation and establishing a certain geological model to calculate the size of hydrocarbon generation pressure increase, and using the Young-Laplace equation to calculate the lower limit of the filling pore throat, but different scholars have proposed different calculation formulas and obtained different results; there are also various data such as high-pressure mercury injection, nuclear magnetic resonance, logging, and oil testing to calculate the thickness of the bound water film under reservoir conditions, compare the differences between experimental and reservoir conditions, and comprehensively determine the lower limit of filling under tight sandstone reservoir conditions; the experimental simulation method mainly uses core displacement experiments, the key of which is the determination of various parameters and the selection of samples. The experimental process and results are relatively ideal.

[0004] While previous research has focused on reservoir pore structure or pore throat size characterization and the lower limit of natural gas injection pore throats, there has been limited research on the pore throat radius range for natural gas injection. Therefore, accurately determining the pore throat radius range for natural gas injection in sandstone reservoirs to improve the efficiency of the natural gas injection process and exploration and development is an urgent issue. Summary of the Invention

[0005] The embodiments of the present application provide a method and apparatus for determining the pore throat radius range for natural gas injection in sandstone reservoirs, thereby accurately determining the pore throat radius range for natural gas injection, at least to a certain extent, providing a basis for natural gas injection process research and exploration and development, and improving the efficiency of the natural gas injection process and exploration and development.

[0006] Other features and advantages of the present application will become apparent from the following detailed description, or may be learned in part by practice of the present application.

[0007] According to a first aspect of an embodiment of the present application, a method for determining a natural gas charging pore throat radius range in a sandstone reservoir is provided, comprising:

[0008] Determine the paleodynamics of natural gas charge in sandstone reservoirs;

[0009] determining paleo-resistance to natural gas charge in said sandstone reservoir;

[0010] determining a maximum filling resistance and a minimum filling resistance when natural gas is filled into the sandstone reservoir according to the paleodynamic force and the paleoresistance;

[0011] The minimum pore throat radius and the maximum pore throat radius are correspondingly determined according to the maximum filling resistance and the minimum filling resistance to obtain the pore throat radius range.

[0012] In some embodiments of the present application, based on the above solution, determining the paleodynamic force of natural gas charging in the sandstone reservoir includes:

[0013] The residual pressure evolution history of the sandstone reservoir was simulated using Petromod basin simulation software to obtain the paleodynamics of natural gas charging in the sandstone reservoir.

[0014] In some embodiments of the present application, based on the above solution, determining the paleo-resistance of natural gas charging in the sandstone reservoir includes:

[0015] Establishing a diagenetic evolution sequence of the sandstone reservoir;

[0016] determining the porosity evolution history of the sandstone reservoir according to the diagenetic evolution sequence;

[0017] Determining the maximum connected pore throat radius evolution history of the sandstone reservoir at different geological times based on the porosity evolution history;

[0018] The paleo-resistance to natural gas filling in the sandstone reservoir is determined according to the evolution history of the maximum connected pore throat radius.

[0019] In some embodiments of the present application, based on the aforementioned solution, determining the evolution history of the maximum connected pore throat radius of the sandstone reservoir at different geological times according to the porosity evolution history includes:

[0020] Based on the high-pressure pump test results, a positive correlation formula between porosity and the maximum connected pore throat radius was determined;

[0021] The evolution history of the maximum connected pore throat radius is determined based on the porosity evolution history and the positive correlation formula.

[0022] In some embodiments of the present application, based on the aforementioned solution, determining the paleo-resistance to natural gas charging in the sandstone reservoir according to the evolution history of the maximum connected pore throat radius includes:

[0023] The paleo-resistance to natural gas filling in the sandstone reservoir is calculated based on the evolution history of the maximum connected pore throat radius, the gas-water interfacial tension and the wetting contact angle.

[0024] In some embodiments of the present application, based on the aforementioned solution, determining the maximum filling resistance and the minimum filling resistance when natural gas is filled in the sandstone reservoir according to the paleodynamic force and the paleoresistance includes:

[0025] The maximum value of the paleodynamic force is determined as the maximum filling resistance;

[0026] The paleodynamic force equal to the paleoresistance within the same geological time is determined as the minimum filling resistance.

[0027] In some embodiments of the present application, based on the aforementioned solution, determining the minimum pore throat radius and the maximum pore throat radius according to the maximum filling resistance and the minimum filling resistance includes:

[0028] determining the minimum pore throat radius according to the maximum filling resistance, the gas-water interfacial tension, and the wetting contact angle;

[0029] The maximum pore throat radius is determined according to the minimum filling resistance, the air-water interfacial tension and the wetting contact angle.

[0030] According to a second aspect of an embodiment of the present application, a device for determining a natural gas injection pore throat radius range in a sandstone reservoir is provided, comprising:

[0031] Paleodynamic determination module, used to determine the paleodynamic forces of natural gas injection in sandstone reservoirs;

[0032] a paleo-resistance determination module, configured to determine paleo-resistance to natural gas charging in the sandstone reservoir;

[0033] a filling resistance determination module, configured to determine a maximum filling resistance and a minimum filling resistance when natural gas is filled in the sandstone reservoir according to the paleodynamics and the paleoresistance;

[0034] The pore throat range determination module is used to determine the minimum pore throat radius and the maximum pore throat radius according to the maximum filling resistance and the minimum filling resistance to obtain the pore throat radius range.

[0035] According to a third aspect of an embodiment of the present application, a device for determining the radius range of a natural gas injection pore throat in a sandstone reservoir is provided, comprising a processor and a memory, wherein the memory stores computer program instructions that can be executed by the processor, and when the processor executes the computer program instructions, the steps of the method described in any one of the above-mentioned first aspects are implemented.

[0036] According to a fourth aspect of an embodiment of the present application, a computer-readable storage medium is provided, in which computer program instructions are stored. When the computer program instructions are executed by a processor, the processor is prompted to implement the steps of the method described in any one of the first aspects above.

[0037] In this application, the paleodynamic force of natural gas filling in a sandstone reservoir is determined; the paleoresistance of natural gas filling in the sandstone reservoir is determined; the maximum and minimum filling resistances during natural gas filling in the sandstone reservoir are determined based on the paleodynamic force and paleoresistance; and the minimum and maximum pore throat radii are correspondingly determined based on the maximum and minimum filling resistances to obtain a pore throat radius range. By calculating the dynamic force and resistance of natural gas filling, the radius range of the pore throats for natural gas filling in sandstone reservoirs is accurately calculated, providing a basis for research on the natural gas filling process and exploration and development, thereby improving the efficiency of the natural gas filling process and exploration and development.

[0038] It should be understood that the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the present application. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] The accompanying drawings are incorporated into and constitute a part of the specification, illustrating embodiments consistent with the present application and, together with the specification, explaining the principles of the present application. Obviously, the drawings described below are only some embodiments of the present application, and those skilled in the art can derive other drawings based on these drawings without inventive effort. In the drawings:

[0040] Figure 1 A schematic flow chart of a method for determining a pore throat radius range for natural gas injection in a sandstone reservoir according to one embodiment is shown;

[0041] Figure 2 A detailed schematic diagram of the evolution history of natural gas charging dynamics in a sandstone reservoir according to one embodiment is shown;

[0042] Figure 3A schematic diagram showing the relationship between porosity and maximum connected pore throat radius in one embodiment is shown;

[0043] Figure 4 A detailed schematic diagram of the evolution history of natural gas charging resistance in a sandstone reservoir in one embodiment is shown;

[0044] Figure 5 A schematic diagram showing the relationship between ancient power and ancient resistance in one embodiment is shown;

[0045] Figure 6 A block diagram of a device for determining a natural gas injection pore throat radius range in a sandstone reservoir according to one embodiment is shown;

[0046] Figure 7 A schematic structural diagram of an apparatus for determining the range of pore throat radius for natural gas injection in a sandstone reservoir is shown in one embodiment. DETAILED DESCRIPTION

[0047] The following will be combined with the drawings in the embodiments of this application to clearly and completely describe the technical solutions in the embodiments of this application. Obviously, the embodiments described are only part of the embodiments of this application, not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of this application.

[0048] In addition, described feature, structure or characteristic can be combined in one or more embodiments in any suitable manner.In the following description, many specific details are provided so as to provide a full understanding of the embodiments of the present application. However, it will be appreciated by those skilled in the art that the technical scheme of the present application can be put into practice without one or more of the specific details, or other methods, components, devices, steps etc. can be adopted. In other cases, known methods, devices, implementations or operations are not shown or described in detail to avoid blurring the various aspects of the application.

[0049] The block diagrams shown in the accompanying drawings are merely functional entities and do not necessarily correspond to physically separate entities. That is, these functional entities may be implemented in software, in one or more hardware modules or integrated circuits, or in different networks and / or processor devices and / or microcontroller devices.

[0050] The flowcharts shown in the accompanying drawings are for illustrative purposes only and do not necessarily include all contents and operations / steps, nor must they be executed in the order described. For example, some operations / steps may be decomposed, while others may be combined or partially combined. Therefore, the actual execution order may vary depending on the actual situation.

[0051] It should also be noted that the terms "first," "second," and the like in the specification and claims of this application and the accompanying drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or precedence. It should be understood that the terms used in this manner are interchangeable where appropriate, so that the embodiments of the present application described herein can be implemented in an order other than that shown or described.

[0052] Figure 1 FIG1 shows a flow chart of a method for determining the range of pore throat radius for natural gas injection in a sandstone reservoir in one embodiment. Figure 1 As shown, a method for determining the range of the pore throat radius for natural gas injection in a sandstone reservoir is provided. The method may include the following steps 101 to 104.

[0053] In step 101 , the paleodynamics of natural gas charge in a sandstone reservoir is determined.

[0054] It should be noted that whether natural gas can be filled into pore throats in sandstone reservoirs depends on the coupling relationship between dynamics and resistance during filling, and the evolution of dynamics and resistance must be considered. However, relevant technologies rarely consider the lower and upper limits of the natural gas filling pore throat radius from the perspectives of filling dynamics, filling resistance, and filling process. This makes it difficult to accurately determine the radius range of the pore throats filled with natural gas in the reservoir, affecting the efficiency of the natural gas filling process and exploration and development.

[0055] It can be understood that paleodynamics refers to the dynamics of natural gas charging in sandstone reservoirs at different geological times.

[0056] In some embodiments, the Petromod basin simulation software can be used to simulate the residual pressure evolution history of the sandstone reservoir to obtain the paleodynamics of natural gas injection in the sandstone reservoir.

[0057] During the implementation process, for the sandstone reservoir to be studied, based on drilling, oil and gas resource evaluation parameters and other data, the stratigraphic stratification of the well and the geological age, lithology, thickness, erosion thickness of each stratum, as well as parameters such as paleoheat flow, paleosurface temperature, and paleowater depth required for basin simulation can be obtained. Then, the Petromod basin simulation software can be used to restore the evolution history of the residual pressure of a single well, and the evolution history of the natural gas charging power in the sandstone reservoir can be obtained, and then the natural gas charging power (i.e., paleopower) at different geological times can be obtained. Figure 2 FIG. 1 shows a detailed schematic diagram of the evolution history of natural gas injection dynamics in a sandstone reservoir in one embodiment. Figure 2 As shown in the figure, the natural gas charging power (i.e. paleopower) was different at different geological times.

[0058] In step 102 , the paleo-resistance to natural gas charging in a sandstone reservoir is determined.

[0059] It can be understood that paleoresistance refers to the resistance to natural gas charging in sandstone reservoirs at different geological times.

[0060] In some embodiments, a diagenetic evolution sequence of a sandstone reservoir can be established; the porosity evolution history of the sandstone reservoir can be determined based on the diagenetic evolution sequence; the maximum connected pore throat radius evolution history of the sandstone reservoir at different periods can be determined based on the porosity evolution history; and the paleo-resistance to natural gas filling in the sandstone reservoir can be determined based on the maximum connected pore throat radius evolution history.

[0061] It can be understood that the porosity evolution history refers to the porosity of sandstone reservoirs at different geological times, and the maximum connected pore throat radius evolution history refers to the maximum connected pore throat radius of sandstone reservoirs at different geological times.

[0062] During the implementation process, the positive correlation formula between porosity and maximum connected pore throat radius can be determined based on the high-pressure pump test results; and the evolution history of the maximum connected pore throat radius can be determined according to the porosity evolution history and the positive correlation formula.

[0063] After obtaining the evolution history of the maximum connected pore throat radius, the paleo-resistance to natural gas filling in sandstone reservoirs can be calculated based on the evolution history of the maximum connected pore throat radius, gas-water interfacial tension, and wetting contact angle.

[0064] Specifically, based on observations of the reservoir space and diagenetic analysis of the sandstone reservoir under study, a diagenetic evolution sequence can be established. The porosity evolution history of the sandstone reservoir can then be reconstructed based on the porosity increase or decrease caused by different diagenetic processes. Then, based on porosity and high-pressure mercury injection analysis results, a formula for the positive correlation between porosity and the maximum connected pore throat radius can be established. Figure 3 A schematic diagram showing the relationship between porosity and maximum connected pore throat radius in one embodiment is shown. Figure 3 As shown in the figure, there is a positive correlation between porosity and the maximum connected pore throat radius. By fitting the porosity and the maximum connected pore throat radius, a positive correlation formula can be obtained. Furthermore, by substituting the porosity in the porosity evolution history into the positive correlation formula, the maximum connected pore throat radius of sandstone reservoirs at different geological times can be calculated. Finally, combined with the test data of gas-water interfacial tension and wetting contact angle, according to the classic capillary pressure formula Among them, P c The evolution history of the natural gas charging resistance in sandstone reservoirs, i.e., the paleo-resistance of natural gas charging, is calculated using the following equations: σ is the gas-water interfacial tension, θ is the wetting contact angle, and r is the maximum connected pore throat radius. Figure 4 FIG. 1 shows a detailed schematic diagram of the evolution history of the natural gas filling resistance in a sandstone reservoir in one embodiment. Figure 4 As shown in Figure 2, the natural gas filling resistance (i.e., paleo-resistance) is different at different geological times.

[0065] In step 103, the maximum filling resistance and the minimum filling resistance during natural gas filling in the sandstone reservoir are determined based on the paleodynamics and paleoresistance.

[0066] In some embodiments, the maximum value of the paleodynamic force can be determined as the maximum filling resistance; and the paleodynamic force equal to the paleoresistance in the same period can be determined as the minimum filling resistance.

[0067] During the implementation process, the coupling relationship between the evolution of paleodynamics and paleoresistance can be analyzed along the time axis based on the evolution history of natural gas charging power and the evolution history of natural gas charging resistance. Figure 5 FIG. 1 shows a schematic diagram of the relationship between the ancient power (charging power) and the ancient resistance (charging resistance) in one embodiment. Figure 5 As shown in Figure 2, at the same geological time, the paleo-resistance and paleo-dynamic force may be the same or different. The paleo-dynamic force when the paleo-dynamic force and paleo-resistance are equal can be determined based on this relationship diagram, that is, Figure 5 The paleodynamic force corresponding to the intersection of the two is used as the minimum charging resistance for natural gas charging in sandstone reservoirs. Based on this relationship diagram, the maximum charging dynamic force during the evolution process can also be determined and used as the maximum charging resistance for natural gas charging in sandstone reservoirs.

[0068] In step 104, the minimum pore throat radius and the maximum pore throat radius are correspondingly determined according to the maximum filling resistance and the minimum filling resistance to obtain a pore throat radius range.

[0069] In some embodiments, the minimum pore throat radius can be determined based on the maximum filling resistance, air-water interfacial tension, and wetting contact angle; the maximum pore throat radius can be determined based on the minimum filling resistance, air-water interfacial tension, and wetting contact angle.

[0070] During the implementation process, the maximum filling resistance, air-water interfacial tension, and wetting contact angle can be substituted into the above capillary pressure formula to calculate the minimum pore throat radius. Similarly, the minimum filling resistance, air-water interfacial tension, and wetting contact angle can be substituted into the above capillary pressure formula to calculate the maximum pore throat radius.

[0071] This application determines the paleodynamic force of natural gas injection in sandstone reservoirs; determines the paleoresistance of natural gas injection in sandstone reservoirs; determines the maximum and minimum filling resistances during natural gas injection in sandstone reservoirs based on the paleodynamic force and paleoresistance; and determines the minimum and maximum pore throat radii based on the maximum and minimum filling resistances to obtain the pore throat radius range. By calculating the dynamic force and resistance of natural gas injection, the radius range of the natural gas injection pore throat in sandstone reservoirs can be accurately calculated, providing a basis for natural gas injection process research and exploration and development, and improving the efficiency of natural gas injection and exploration and development.

[0072] In other embodiments, the method for determining the range of the pore throat radius for natural gas injection in a sandstone reservoir may further include steps 201 to 204:

[0073] Step 201: simulating the residual pressure evolution history of the sandstone reservoir using Petromod basin simulation software to obtain the paleodynamics of natural gas filling in the sandstone reservoir;

[0074] Step 202: Establish a diagenetic evolution sequence for the sandstone reservoir, determine the porosity evolution history of the sandstone reservoir based on the diagenetic evolution sequence, determine a positive correlation formula between porosity and maximum connected pore throat radius based on high-pressure pump test results, determine the evolution history of the maximum connected pore throat radius based on the porosity evolution history and the positive correlation formula, and calculate the paleo-resistance to natural gas filling in the sandstone reservoir based on the evolution history of the maximum connected pore throat radius, the gas-water interfacial tension, and the wetting contact angle.

[0075] Step 203: The maximum value of the paleodynamic force is determined as the maximum filling resistance of natural gas filling in the sandstone reservoir, and the paleodynamic force equal to the paleoresistance within the same geological time is determined as the minimum filling resistance of natural gas filling in the sandstone reservoir;

[0076] Step 204 : determining the minimum pore throat radius based on the maximum filling resistance, the air-water interfacial tension, and the wetting contact angle, and determining the maximum pore throat radius based on the minimum filling resistance, the air-water interfacial tension, and the wetting contact angle, thereby obtaining a pore throat radius range.

[0077] Taking the determination of the natural gas charging pore throat radius range in the fine sandstone reservoir at 3820.5 m in the upper Sha 3 subsection of the Yuetan 1 well in the eastern sag of the Liaohe Depression as an example, in step 201, for the sandstone reservoir to be studied, based on the drilling data of the Yuehai 1 well and the oil and gas resource evaluation parameter data of the Liaohe Depression, the stratigraphic layers from the upper Sha 3 subsection of the well to the Quaternary and the geological age, lithology, thickness, and denudation thickness parameters of each stratigraphy, as well as the paleoheat flow, paleosurface temperature, paleowater depth and other parameters required for basin simulation are obtained. The Petromod basin simulation software is used to restore the formation residual pressure evolution history at 3820.5 m in the upper Sha 3 subsection of the well, thereby obtaining the evolution history of the natural gas charging power in the sandstone reservoir ( Figure 2 ), which is the ancient power of natural gas injection.

[0078] In step 202, based on the reservoir space observation and diagenetic analysis of the sandstone reservoir in the upper Sha 3 sub-member of the well, the diagenetic evolution sequence of the sandstone reservoir is established, and the porosity evolution history of the sandstone reservoir is restored according to the porosity increase or decrease caused by different diagenesis. Then, based on the existing porosity and high-pressure mercury injection test analysis results, the positive correlation relationship between porosity and maximum connected pore throat radius and the formula ( Figure 3); then, based on the porosity evolution history, the evolution history of the maximum connected pore throat radius of the sandstone reservoir at different geological times was calculated; and combined with the test data of gas-water interfacial tension and wetting contact angle, the evolution history of natural gas filling resistance in the sandstone reservoir was calculated according to the classic capillary pressure formula ( Figure 4 ), which is the ancient resistance to natural gas charging.

[0079] In step 203, based on the paleodynamics and paleoresistance evolution of natural gas charging restored in steps 201 and 202, the coupling relationship between the paleodynamics and paleoresistance evolution is analyzed along the time axis to determine the paleodynamics when the paleodynamics begins to equal the paleoresistance (equivalent to the minimum charging resistance) and the maximum paleodynamics during the evolution process (equivalent to the maximum charging resistance), which occurred at 38 Ma and 0 Ma, respectively. Figure 5 ), thereby obtaining the maximum and minimum filling resistances of natural gas in sandstone reservoirs, which are 0.0094MPa and 8.94MPa, respectively.

[0080] In step 204, based on the minimum and maximum filling resistances when natural gas is injected into the sandstone reservoir of the upper Shahejie Formation 3, the corresponding maximum and minimum pore throat radii are calculated according to the capillary pressure formula, thereby obtaining the pore throat range for natural gas injection, i.e., between 0.0052 and 12.3587 μm.

[0081] The embodiments of the present application restore the paleodynamics and paleoresistance of natural gas charging, analyze the coupled relationship between the evolution of natural gas charging dynamics and resistance, determine the maximum and minimum charging resistances during natural gas charging, and then calculate the minimum and maximum pore throat radii for natural gas charging. This allows the pore throat range of natural gas charging to be determined, providing a basis for research, exploration, and development of the natural gas charging process.

[0082] The following describes an embodiment of the apparatus of the present application, which can be used to implement the method for determining the pore throat radius range for natural gas injection into a sandstone reservoir described in the aforementioned embodiment of the present application. For details not disclosed in the apparatus embodiment of the present application, please refer to the embodiment of the method for determining the pore throat radius range for natural gas injection into a sandstone reservoir described in the aforementioned embodiment of the present application.

[0083] Figure 6 The block diagram of the device for determining the radius range of the pore throat for natural gas filling in the sandstone reservoir in the embodiment of the present application is shown. Figure 6As shown, the device for determining the pore throat radius range of natural gas charging in a sandstone reservoir according to an embodiment of the present application includes: a paleodynamic force determination module 601, a paleoresistance determination module 602, a charging resistance determination module 603 and a pore throat range determination module 604, wherein the paleodynamic force determination module 601 is used to determine the paleodynamic force of natural gas charging in the sandstone reservoir; the paleoresistance determination module 602 is used to determine the paleoresistance of natural gas charging in the sandstone reservoir; the charging resistance determination module 603 is used to determine the maximum charging resistance and the minimum charging resistance when natural gas is charged in the sandstone reservoir according to the paleodynamic force and the paleoresistance; the pore throat range determination module 604 is used to determine the minimum pore throat radius and the maximum pore throat radius according to the maximum charging resistance and the minimum charging resistance, so as to obtain the pore throat radius range.

[0084] In some embodiments of the present application, based on the aforementioned solution, the paleodynamics determination module 601 is further configured to simulate the residual pressure evolution history of the sandstone reservoir using the Petromod basin simulation software to obtain the paleodynamics of natural gas injection in the sandstone reservoir.

[0085] In some embodiments of the present application, based on the aforementioned scheme, the paleo-resistance determination module 602 is also used to establish a diagenetic evolution sequence of the sandstone reservoir; determine the porosity evolution history of the sandstone reservoir according to the diagenetic evolution sequence; determine the maximum connected pore throat radius evolution history of the sandstone reservoir at different geological times according to the porosity evolution history; and determine the paleo-resistance to natural gas filling in the sandstone reservoir according to the maximum connected pore throat radius evolution history.

[0086] In some embodiments of the present application, based on the aforementioned solution, the paleo-resistance determination module 602 is further configured to determine a positive correlation formula between porosity and maximum connected pore throat radius based on high-pressure pump test results;

[0087] According to the porosity evolution history and the positive correlation formula, the evolution history of the maximum connected pore throat radius is determined.

[0088] In some embodiments of the present application, based on the aforementioned scheme, the paleo-resistance determination module 602 is further used to calculate the paleo-resistance of natural gas filling in the sandstone reservoir according to the evolution history of the maximum connected pore throat radius, the gas-water interfacial tension and the wetting contact angle.

[0089] In some embodiments of the present application, based on the aforementioned scheme, the filling resistance determination module 603 is further used to determine the maximum value of the paleodynamics as the maximum filling resistance; and to determine the paleodynamics equal to the paleoresistance within the same geological time as the minimum filling resistance.

[0090] In some embodiments of the present application, based on the aforementioned solution, the pore throat range determination module 604 is further configured to determine the minimum pore throat radius according to the maximum filling resistance, the gas-water interfacial tension, and the wetting contact angle;

[0091] The maximum pore throat radius is determined based on the minimum filling resistance, gas-water interfacial tension and wetting contact angle.

[0092] Based on the same inventive concept, the embodiment of the present application also provides a device for determining the radius range of the pore throat of natural gas injection in a sandstone reservoir, referring to Figure 7 , showing a schematic structural diagram of a device for determining the radius range of pore throats for natural gas injection in sandstone reservoirs in an embodiment of the present application, wherein the device for determining the radius range of pore throats for natural gas injection in sandstone reservoirs comprises one or more memories 704, one or more processors 702, and at least one computer program (computer program instruction) stored on the memories 704 and executable on the processors 702, and when the processors 702 execute the computer programs, the method described above is implemented.

[0093] Among them, Figure 7 In the embodiment of the present invention, a bus architecture (represented by bus 700) is shown. Bus 700 may include any number of interconnected buses and bridges, and bus 700 links together various circuits including one or more processors represented by processor 702 and memory represented by memory 704. Bus 1200 may also link together various other circuits such as peripherals, voltage regulators, and power management circuits, which are well known in the art and therefore will not be described further herein. Bus interface 705 provides an interface between bus 700 and receiver 701 and transmitter 703. Receiver 701 and transmitter 703 may be the same component, namely a transceiver, which provides a unit for communicating with various other devices over a transmission medium. Processor 702 is responsible for managing bus 700 and general processing, while memory 704 may be used to store data used by processor 702 when performing operations.

[0094] Based on the same inventive concept, an embodiment of the present application provides a computer-readable storage medium, in which computer program instructions are stored. When the computer program instructions are executed by a processor, the processor is prompted to implement the steps of the method as described above.

[0095] The functions described herein may be implemented in hardware, software executed by a processor, firmware, or any combination thereof. If implemented in software executed by a processor, the functions may be stored as one or more instructions or codes on or transmitted via a computer-readable medium. Other examples and implementations are within the scope and spirit of this application and the appended claims. For example, due to the nature of software, the functions described above may be implemented using software executed by a processor, hardware, firmware, hardwiring, or a combination of any of these. Furthermore, the functional units may be integrated into a single processing unit, each unit may exist physically separately, or two or more units may be integrated into a single unit.

[0096] In the several embodiments provided in this application, it should be understood that the disclosed technical content can be implemented in other ways. Among them, the device embodiments described above are only exemplary. For example, the division of the units can be a logical function division. In actual implementation, there may be other division methods, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of units or modules, which can be electrical or other forms.

[0097] The units described as separate components may or may not be physically separate, and the components of the control device may or may not be physical units, that is, they may be located in one place or distributed across multiple units. Some or all of the units may be selected according to actual needs to achieve the purpose of the present embodiment.

[0098] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present application is essentially or the part that contributes to the prior art or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions for enabling a computer device (which can be a personal computer, a server or a network device, etc.) to execute all or part of the steps of the method described in each embodiment of the present application. The aforementioned storage medium includes: U disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), mobile hard disk, magnetic disk or optical disk, etc. Various media that can store computer program instructions.

[0099] The foregoing is merely an embodiment of the present application and is not intended to limit the present application. Various modifications and variations are possible for those skilled in the art. Any modifications, equivalent substitutions, or improvements made within the spirit and principles of the present application shall be included within the scope of the claims of the present application.

Claims

1. A method for determining the range of pore throat radius for natural gas injection in sandstone reservoirs, characterized in that: include: Determine the paleodynamics of natural gas charge in sandstone reservoirs; determining paleo-resistance to natural gas charge in said sandstone reservoir; determining a maximum filling resistance and a minimum filling resistance when natural gas is filled into the sandstone reservoir according to the paleodynamic force and the paleoresistance; The minimum pore throat radius and the maximum pore throat radius are correspondingly determined according to the maximum filling resistance and the minimum filling resistance to obtain the pore throat radius range.

2. The method according to claim 1, characterized in that The determination of the paleodynamics of natural gas charging in sandstone reservoirs comprises: The residual pressure evolution history of the sandstone reservoir was simulated using Petromod basin simulation software to obtain the paleodynamics of natural gas charging in the sandstone reservoir.

3. The method according to claim 1, characterized in that Determining the paleo-resistance to natural gas filling in the sandstone reservoir comprises: Establishing a diagenetic evolution sequence of the sandstone reservoir; determining the porosity evolution history of the sandstone reservoir according to the diagenetic evolution sequence; Determining the maximum connected pore throat radius evolution history of the sandstone reservoir at different geological times based on the porosity evolution history; The paleo-resistance to natural gas filling in the sandstone reservoir is determined according to the evolution history of the maximum connected pore throat radius.

4. The method according to claim 3, characterized in that Determining the maximum connected pore throat radius evolution history of the sandstone reservoir at different geological times based on the porosity evolution history includes: Based on the high-pressure pump test results, a positive correlation formula between porosity and the maximum connected pore throat radius was determined; The evolution history of the maximum connected pore throat radius is determined based on the porosity evolution history and the positive correlation formula.

5. The method according to claim 3, characterized in that Determining the paleo-resistance to natural gas filling in the sandstone reservoir based on the evolution history of the maximum connected pore throat radius includes: The paleo-resistance to natural gas filling in the sandstone reservoir is calculated based on the evolution history of the maximum connected pore throat radius, the gas-water interfacial tension and the wetting contact angle.

6. The method according to claim 1, characterized in that The determining of the maximum filling resistance and the minimum filling resistance when natural gas is filled in the sandstone reservoir according to the paleodynamics and the paleoresistance includes: The maximum value of the paleodynamic force is determined as the maximum filling resistance; The paleodynamic force equal to the paleoresistance within the same geological time is determined as the minimum filling resistance.

7. The method according to claim 1, characterized in that The determining of the minimum pore throat radius and the maximum pore throat radius according to the maximum filling resistance and the minimum filling resistance includes: determining the minimum pore throat radius according to the maximum filling resistance, the gas-water interfacial tension, and the wetting contact angle; The maximum pore throat radius is determined according to the minimum filling resistance, the air-water interfacial tension and the wetting contact angle.

8. A device for determining the radius range of pore throats for natural gas injection in sandstone reservoirs, characterized in that: include: Paleodynamic determination module, used to determine the paleodynamic forces of natural gas injection in sandstone reservoirs; a paleo-resistance determination module, configured to determine paleo-resistance to natural gas charging in the sandstone reservoir; a filling resistance determination module, configured to determine a maximum filling resistance and a minimum filling resistance when natural gas is filled in the sandstone reservoir according to the paleodynamics and the paleoresistance; The pore throat range determination module is used to determine the minimum pore throat radius and the maximum pore throat radius according to the maximum filling resistance and the minimum filling resistance to obtain the pore throat radius range.

9. A device for determining the radius range of a natural gas injection pore throat in a sandstone reservoir, comprising a processor and a memory, characterized in that: The memory stores computer program instructions that can be executed by the processor, and when the processor executes the computer program instructions, the steps of the method according to any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium, characterized in that The computer-readable storage medium stores computer program instructions, which, when executed by a processor, prompt the processor to implement the steps of the method according to any one of claims 1 to 7.

Images