CN112698401A - Rock physical template construction method, device, equipment and storage medium - Google Patents

Rock physical template construction method, device, equipment and storage medium Download PDF

Info

Publication number
CN112698401A
CN112698401A CN202011418723.5A CN202011418723A CN112698401A CN 112698401 A CN112698401 A CN 112698401A CN 202011418723 A CN202011418723 A CN 202011418723A CN 112698401 A CN112698401 A CN 112698401A
Authority
CN
China
Prior art keywords
reservoir
rock
density
wave velocity
petrophysical
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011418723.5A
Other languages
Chinese (zh)
Other versions
CN112698401B (en
Inventor
王磊
陈彬滔
白洁
何世琦
薛罗
杜炳毅
潘树新
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN202011418723.5A priority Critical patent/CN112698401B/en
Publication of CN112698401A publication Critical patent/CN112698401A/en
Application granted granted Critical
Publication of CN112698401B publication Critical patent/CN112698401B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/282Application of seismic models, synthetic seismograms
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/622Velocity, density or impedance
    • G01V2210/6222Velocity; travel time
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/622Velocity, density or impedance
    • G01V2210/6224Density
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/624Reservoir parameters
    • G01V2210/6242Elastic parameters, e.g. Young, Lamé or Poisson
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/624Reservoir parameters
    • G01V2210/6246Permeability
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention provides a petrophysical template construction method, a petrophysical template construction device, a petrophysical template construction equipment and a storage medium, wherein the method comprises the following steps: calculating to obtain rock reservoir elastic parameters according to reservoir rock physical parameters of a research work area and a two-phase medium rock physical theoretical model, wherein the rock reservoir elastic parameters comprise reservoir longitudinal wave velocity, transverse wave velocity and density under the conditions of different porosities and different water saturation; constructing a first reservoir sensitive parameter and a second reservoir sensitive parameter which are influenced by eliminating a rock framework and retaining pores and pore fluid according to the rock reservoir elastic parameters; and constructing a rock physical template according to the first reservoir sensitive parameter and the second reservoir sensitive parameter.

Description

Rock physical template construction method, device, equipment and storage medium
Technical Field
The invention relates to the field of petroleum geophysical exploration, in particular to a method, a device, equipment and a storage medium for constructing a petrophysical template.
Background
With the deep development of petroleum exploration and development, reservoir prediction technology is continuously updated and iterated, and exploration means are developed into methods such as prestack inversion, elastic parameter inversion and high-precision geostatistics inversion from poststack attribute analysis and poststack wave impedance inversion. Rock physics plays an important role in a bridge in the process of oil exploration and development, the relationship between physical property parameters of a reservoir and seismic response characteristics is established, the rock physics template established based on a rock physics theoretical model describes the distribution rule of reservoirs with different lithologies, different porosities and different water saturation in the rock physics template, and theoretical basis is provided for quantitative reservoir prediction. The rock physical analysis at the present stage generally carries out quantitative reservoir prediction based on a physical template constructed based on a longitudinal wave impedance-longitudinal wave velocity ratio, and the rock physical theory considers that when a reservoir contains oil and gas, the corresponding longitudinal wave velocity and density are reduced, and the transverse wave velocity only reflects the rock skeleton property and has no obvious transformation, so that the oil and gas reservoir has independent distribution areas in the rock physical template constructed based on the longitudinal wave impedance-longitudinal wave velocity ratio, and the reservoir porosity and the water saturation can be effectively distinguished. The reservoir prediction factor constructed by combining the rock physical template based on longitudinal wave impedance-longitudinal and transverse wave velocity ratio with the prestack inversion parameters is used for quantitative detection of the oil-gas reservoir, and a good application effect is achieved in the practical application process. With the continuous deepening of exploration, a high-porosity high-hydrocarbon-content high-quality reservoir is gradually explored and discovered and is transferred to a development stage, an exploration target is gradually converted into a water-containing oil layer with medium and low porosity, the application of a conventional longitudinal wave impedance-longitudinal wave velocity ratio rock physical template is greatly limited due to the fact that oil and water are distributed in the same layer and the porosity is medium, the distribution of the water-containing oil layer and the oil-containing oil layer cannot be effectively distinguished, the phenomenon is caused because the distribution area of the water-containing sandstone and the oil-containing sandstone in the conventional longitudinal wave impedance-longitudinal wave velocity ratio rock physical template is relatively concentrated, when the porosity is gradually reduced, the difference between the water-containing sandstone and the oil-containing sandstone is not obvious, the conventional rock physical template is difficult to effectively distinguish the distribution of the water-containing sandstone and the oil-containing sandstone under the condition of medium and.
Disclosure of Invention
The method is used for solving the problem that when the petrophysical template used in the prior art is applied to a water-bearing oil layer with medium and low porosity, the distribution of the water-bearing oil layer and the oil-bearing oil layer cannot be effectively distinguished, so that the risk of reservoir prediction is increased.
In order to solve the above technical problem, a first aspect herein provides a petrophysical template construction method, including:
calculating to obtain rock reservoir elastic parameters according to reservoir rock physical parameters of a research work area and a two-phase medium rock physical theoretical model, wherein the rock reservoir elastic parameters comprise reservoir longitudinal wave velocity, transverse wave velocity and density under the conditions of different porosities and different water saturation;
constructing a first reservoir sensitive parameter and a second reservoir sensitive parameter which are influenced by eliminating a rock framework and retaining pores and pore fluid according to the rock reservoir elastic parameters;
and constructing a rock physical template according to the first reservoir sensitive parameter and the second reservoir sensitive parameter.
In further embodiments herein, the investigating the reservoir petrophysical parameters of the work area comprises: matrix bulk modulus, matrix shear modulus, matrix density, water bulk modulus, oil bulk modulus, water density, oil density, porosity, and water saturation.
In a further embodiment of the present disclosure, the calculating the elastic parameter of the rock reservoir according to the petrophysical parameter of the reservoir in the research work area and the two-phase medium petrophysical theoretical model includes: based on the two-phase medium rock physical theory model, the following calculation is performed:
calculating the volume modulus and the shear modulus of the framework according to the volume modulus, the shear modulus and the porosity of the matrix;
calculating the bulk modulus and the density of the mixed fluid according to the water saturation, the bulk modulus of water, the bulk modulus of oil, the water density and the oil density;
calculating the rock equivalent volume modulus of the saturated two-phase fluid according to the skeleton volume modulus, the porosity, the matrix volume modulus and the mixed fluid volume modulus;
equating the skeleton shear modulus to the rock equivalent shear modulus of saturated two-phase fluid;
calculating the rock equivalent density of the saturated two-phase fluid according to the matrix density, the porosity and the mixed fluid density;
calculating longitudinal wave velocity and transverse wave velocity of the reservoir according to the rock equivalent volume modulus of saturated two-phase fluid, the rock equivalent shear modulus of saturated two-phase fluid and the rock equivalent density of saturated two-phase fluid;
the rock equivalent density of the saturated two-phase fluid is taken as the reservoir density.
In further embodiments herein, calculating reservoir compressional and shear velocities from the petrophysical equivalent bulk modulus of the saturated two-phase fluid, the petrophysical equivalent shear modulus of the saturated two-phase fluid, and the petrophysical equivalent density of the saturated two-phase fluid comprises calculating reservoir compressional and shear velocities using the following formulas:
Figure BDA0002821276910000021
wherein v ispIs the reservoir longitudinal wave velocity, vsIs the reservoir shear wave velocity, KsatRock equivalent bulk modulus, μ, for saturated two-phase fluidssatRock equivalent shear modulus, ρ, for saturated two-phase fluidssatIs the rock equivalent density of a two-phase fluid.
In a further embodiment of this document, constructing a first reservoir sensitive parameter and a second reservoir sensitive parameter that eliminate the influence of rock skeleton and retain the influence of pores and pore fluids according to the rock reservoir elastic parameters includes:
determining a first reservoir sensitive parameter according to the ratio of the reservoir density to the reservoir shear wave velocity;
and determining a second reservoir sensitive parameter according to the difference value of the reservoir longitudinal wave velocity and the reservoir transverse wave velocity multiplied by the density.
In further embodiments herein, determining the first reservoir sensitivity parameter based on the ratio of the reservoir density to the shear wave velocity comprises calculating the first reservoir sensitivity parameter using the following equation:
X=ρ/vs
and determining a second reservoir sensitivity parameter according to the difference between the longitudinal wave velocity and the transverse wave velocity and the product of the density, wherein the second reservoir sensitivity parameter is determined by using the following formula:
Figure BDA0002821276910000031
where ρ is the reservoir density, vsIs the reservoir shear wave velocity, vpIs the reservoir compressional velocity.
In further embodiments herein, constructing a petrophysical template from the first reservoir sensitivity parameter and the second reservoir sensitivity parameter comprises:
determining values of porosity and water saturation according to the value of the first reservoir sensitive parameter and the value of the second reservoir sensitive parameter;
setting the first reservoir sensitive parameter as an X axis and the second reservoir sensitive parameter as a Y axis, and constructing a rock physical template according to the value of the first reservoir sensitive parameter, the value of the second reservoir sensitive parameter, the corresponding porosity and the value of the water saturation.
A second aspect herein provides a petrophysical template building apparatus comprising:
the elastic parameter calculation module is used for calculating to obtain elastic parameters of the rock reservoir according to rock physical parameters of the reservoir in the research work area and a two-phase medium rock physical theoretical model, wherein the elastic parameters of the rock reservoir comprise longitudinal wave velocity, transverse wave velocity and density of the reservoir under the conditions of different porosities and different water saturation degrees;
the sensitive parameter calculation module is used for constructing a first reservoir sensitive parameter and a second reservoir sensitive parameter which are used for eliminating the influence of a rock framework and keeping the influence of pores and pore fluid according to the elastic parameters of the rock reservoir;
and the construction module is used for constructing the petrophysical template according to the first reservoir sensitive parameters and the second reservoir sensitive parameters.
A third aspect of the present disclosure provides a computer apparatus comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the petrophysical template construction method of any of the preceding embodiments when executing the computer program.
A fourth aspect herein provides a computer readable storage medium storing an executable computer program which, when executed by a processor, implements the petrophysical template construction method of any of the preceding embodiments.
According to the rock physical template construction method and device, the elastic parameters of a rock reservoir are calculated by combining a two-phase medium rock physical theoretical model with reservoir rock physical parameters of a research work area; constructing a first reservoir sensitive parameter and a second reservoir sensitive parameter which are influenced by rock frameworks and reserved pores and pore fluid according to the elastic parameters of the rock reservoirs; and constructing a petrophysical template according to the first reservoir sensitive parameters and the second reservoir sensitive parameters, so that the water-containing sandstone and the oil-containing sandstone in the constructed petrophysical template have obvious distribution difference, and the water-containing sandstone and the oil-containing sandstone still have strong resolution even under the condition of medium-low porosity, thereby improving the precision of reservoir prediction quantitative interpretation and reducing the exploration risk.
In order to make the aforementioned and other objects, features and advantages of the present invention comprehensible, preferred embodiments accompanied with figures are described in detail below.
Drawings
In order to more clearly illustrate the embodiments or technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 shows a flow diagram of a petrophysical template construction method of embodiments herein;
FIG. 2 illustrates a flow chart of elasticity parameter calculation according to embodiments herein;
FIG. 3 illustrates a flow diagram of a sensitive parameter construction process of embodiments herein;
FIG. 4 illustrates a distribution law graph of a first reservoir sensitivity parameter according to embodiments herein;
FIG. 5 is a graph illustrating a distribution law of a second reservoir sensitivity parameter according to embodiments herein;
FIG. 6 shows a schematic diagram of a petrophysical template of an embodiment herein;
FIG. 7 shows a schematic diagram of a prior art petrophysical template;
FIG. 8 is a block diagram of a petrophysical template building apparatus according to an embodiment herein;
FIG. 9 shows a block diagram of a computer device according to an embodiment of the present disclosure.
Description of the symbols of the drawings:
810. an elasticity parameter calculation module;
820. a sensitive parameter calculation module;
830. building a module;
902. a computer device;
904. a processor;
906. a memory;
908. a drive mechanism;
910. an input/output module;
912. an input device;
914. an output device;
916. a presentation device;
918. a graphical user interface;
920. a network interface;
922. a communication link;
924. a communication bus.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments herein without making any creative effort, shall fall within the scope of protection.
In an embodiment herein, as shown in fig. 1, fig. 1 shows a flowchart of a rock physical template construction method in the embodiment herein, and the rock physical template construction method provided in this embodiment may be executed in an intelligent terminal, including a smart phone, a tablet computer, a desktop computer, and the like, and may be an individual application program, an applet embedded in another program, and the like, or may also be in a web page form, and the specific implementation manner is not limited herein. The embodiment is used for solving the problem that when the petrophysical template used in the prior art is applied to a water-bearing oil layer with medium and low porosity, the distribution of the water-bearing oil layer and the oil-bearing oil layer cannot be effectively distinguished, so that the risk of reservoir prediction is increased. Specifically, the rock physical template construction method comprises the following steps:
110, calculating elastic parameters of the rock reservoir according to rock physical parameters of the reservoir in the research work area and a two-phase medium rock physical theoretical model, wherein the elastic parameters of the rock reservoir comprise longitudinal wave velocity, transverse wave velocity and density of the reservoir under the conditions of different porosities and different water saturation degrees;
step 120, constructing a first reservoir sensitive parameter and a second reservoir sensitive parameter which are influenced by eliminating a rock framework and retaining pores and pore fluid according to the elastic parameters of the rock reservoir;
and step 130, constructing a rock physical template according to the first reservoir sensitive parameters and the second reservoir sensitive parameters.
In detail, the reservoir of the research work area is, for example, a sand shale reservoir, and may also be other two-phase medium reservoirs, which is not limited herein. Reservoir petrophysical parameters of the research work area comprise matrix bulk modulus, matrix shear modulus, matrix density, water bulk modulus, oil bulk modulus, water density, oil density, porosity and water saturation. The two-phase medium rock physical theory model is used for calculating longitudinal wave velocity, transverse wave velocity and density of the reservoir under the conditions of different porosities and different water saturation based on reservoir rock physical parameters in a research work area, and the specific implementation process refers to the following embodiments.
The first reservoir sensitive parameter and the second reservoir sensitive parameter are not particularly limited, and the first reservoir sensitive parameter and the second reservoir sensitive parameter which are not affected by the rock framework and only affected by the pore and pore fluid can be used as the first reservoir sensitive parameter and the second reservoir sensitive parameter.
In the embodiment, the water-containing sandstone and the oil-containing sandstone in the rock physical template constructed based on the first reservoir sensitive parameters and the second reservoir sensitive parameters have obvious distribution difference, and the water-containing sandstone and the oil-containing sandstone still have strong resolution even under the condition of medium-low porosity, so that the precision of reservoir prediction quantitative interpretation is improved, and the exploration risk is reduced.
In an embodiment of this document, as shown in fig. 2, the step 110 of calculating the elastic parameter of the rock reservoir according to the petrophysical parameter of the reservoir in the research work area and the two-phase medium petrophysical theoretical model includes: based on the two-phase medium rock physical theory model, the following calculation is carried out:
step 210, calculating the matrix volume modulus and the matrix shear modulus according to the matrix volume modulus, the matrix shear modulus and the porosity;
step 220, calculating the bulk modulus and the density of the mixed fluid according to the water saturation, the bulk modulus of water, the bulk modulus of oil, the water density and the oil density;
step 230, calculating the rock equivalent volume modulus of the saturated two-phase fluid according to the skeleton volume modulus, the porosity, the matrix volume modulus and the mixed fluid volume modulus;
step 240, enabling the skeleton shear modulus to be equivalent to the rock equivalent shear modulus of saturated two-phase fluid;
step 250, calculating rock equivalent density of saturated two-phase fluid according to the density of the matrix, the porosity and the density of the mixed fluid;
and step 260, calculating to obtain longitudinal wave velocity and transverse wave velocity of the reservoir according to the rock equivalent volume modulus of the saturated two-phase fluid, the rock equivalent shear modulus of the saturated two-phase fluid and the rock equivalent density of the saturated two-phase fluid.
In step 210, calculating the bulk modulus and the shear modulus of the framework according to the bulk modulus, the shear modulus and the porosity of the matrix includes calculating the bulk modulus and the shear modulus of the framework by using the following formulas:
Km=Kg(1-φ)4/(1-φ)m=Kmμg/Kg
wherein, KgAnd mugRespectively, the bulk modulus and shear modulus of the matrix, phi is the porosity, KmAnd mumThe bulk modulus and shear modulus of the framework are shown.
In the step 220, the bulk modulus of the mixed fluid and the density of the mixed fluid are calculated according to the water saturation, the bulk modulus of water, the bulk modulus of oil, the water density and the oil density, and the bulk modulus of the mixed fluid and the density of the mixed fluid are calculated by using the following formulas:
Figure BDA0002821276910000071
ρf=Swρw+(1-Swo
wherein S iswIs the water saturation, KwIs the water bulk modulus, KoIs the oil volume modulus, ρwIs the water density, ρoIs the oil density, KfFor bulk modulus of mixed fluids, pfIs the mixed fluid density.
In step 230, calculating the rock equivalent bulk modulus of the saturated two-phase fluid according to the skeleton bulk modulus, the porosity, the matrix bulk modulus, and the bulk modulus of the mixed fluid, including calculating the rock equivalent bulk modulus of the saturated two-phase fluid by using the following formula:
Figure BDA0002821276910000072
wherein, KmFor a volume modulus of the frame, KgIs the bulk modulus of the matrix, KfIn order to mix the bulk modulus of the fluid,
Figure BDA0002821276910000073
is porosity.
In step 240, the skeleton shear modulus is equivalent to the rock equivalent shear modulus of the saturated two-phase fluid by using the following formula:
μsat=μm
wherein, mumThe skeletal shear modulus.
In the step 250, calculating the rock equivalent density of the saturated two-phase fluid according to the matrix density, the porosity and the mixed fluid density, including calculating the rock equivalent density of the saturated two-phase fluid by using the following formula:
ρsat=(1-φ)ρg+φρf
where ρ isgThe density of the matrix is used as the density of the matrix,
Figure BDA0002821276910000074
is the porosity, pfIs the mixed fluid density.
In step 260, the longitudinal wave velocity and the transverse wave velocity of the reservoir are calculated according to the rock equivalent volume modulus of the saturated two-phase fluid, the rock equivalent shear modulus of the saturated two-phase fluid and the rock equivalent density of the saturated two-phase fluid, and the method includes the following steps:
Figure BDA0002821276910000081
wherein v ispIs the velocity of longitudinal wave, vsIs the transverse wave velocity, KsatRock equivalent bulk modulus, mu, of saturated two-phase fluidsatRock equivalent shear modulus, ρ, for saturated two-phase fluidssatIs the rock equivalent density of a two-phase fluid.
In an embodiment of this document, as shown in fig. 3, the step 120 includes constructing a first reservoir sensitive parameter and a second reservoir sensitive parameter that eliminate the influence of the rock skeleton and retain the influence of the pores and the pore fluids according to the rock reservoir elastic parameters, including:
step 121, determining a first reservoir sensitivity parameter according to the ratio of the reservoir shear wave velocity to the reservoir density;
and step 122, determining a second reservoir sensitive parameter according to the difference value between the reservoir longitudinal wave velocity and the reservoir transverse wave velocity and the product of the density.
Specifically, the step 121 of determining the first reservoir sensitivity parameter according to the ratio of the reservoir density to the reservoir shear wave velocity includes calculating the first reservoir sensitivity parameter by using the following formula:
Xaρ/vs
wherein X is a first reservoir sensitivity parameter, a is a coefficient, ρ is reservoir density, and v issIs the reservoir shear wave velocity. Where a may be any number greater than 0, and in a preferred embodiment, a is 1.
Reservoir density ρ reflects the comprehensive response of the rock skeleton, rock pores and pore fluids; reservoir shear wave velocity vsThe response of the rock framework is reflected, so that the first reservoir sensitivity parameter formed by dividing the reservoir density by the shear wave speed can eliminate the common part of the reservoir density and the shear wave speed, namely the rock framework, and only the action of pores and fluid is reserved.
Determining the second reservoir sensitivity parameter in step 122 according to the difference between the product of the reservoir compressional wave velocity and the reservoir shear wave velocity and the density comprises calculating a second sensitivity coefficient by using the following formula:
Y=vpρ-bvsρ;
where ρ isIs the reservoir density, vsIs the reservoir shear wave velocity, vpIs the reservoir compressional velocity and b is the coefficient.
The reservoir longitudinal wave velocity reflects the comprehensive response of a rock framework, rock pores and pore fluid as the same as the density, and the longitudinal wave velocity x the density-b transverse wave velocity x the density is also used for eliminating the common part of the longitudinal wave velocity x the density and retaining the effects of the pores and the fluid.
Preferably, the first and second liquid crystal materials are,
Figure BDA0002821276910000082
the reason is as follows: the sandstone containing water generally has a ratio of longitudinal wave velocity to transverse wave velocity
Figure BDA0002821276910000083
In a coordinate system, if the abscissa is the transverse wave velocity and the ordinate is the longitudinal wave velocity, when the slope is
Figure BDA0002821276910000084
And then fitting the longitudinal and transverse wave velocities of the water-containing sandstone to obtain a water-containing baseline.
In an embodiment, as shown in fig. 4 and 5, fig. 4 and 4 respectively show a distribution law graph of a first reservoir sensitive parameter and a second reservoir sensitive parameter according to the embodiment of the present disclosure, in which an abscissa represents porosity and an ordinate represents water saturation, and filling colors in the graph represent values of the sensitive parameters, wherein dark colors represent high values of the property, and light colors represent low values of the property. As can be seen in fig. 4, the value of the first reservoir sensitivity parameter monotonically increases as the porosity progressively increases from 0 to 0.35, regardless of the change in water saturation; for a particular porosity, the value of the first reservoir sensitivity parameter does not change as the water saturation changes, indicating that the first reservoir sensitivity parameter is not sensitive to reservoir fluids, i.e., water saturation, and is only very sensitive to porosity. As can be seen in fig. 5, as the porosity and water saturation progressively increase, the value of the second reservoir sensitivity parameter progressively increases, indicating that the second reservoir sensitivity parameter has sensitivity to both reservoir porosity and water saturation.
In one embodiment herein, constructing a petrophysical template from the first reservoir sensitivity parameter and the second reservoir sensitivity parameter comprises:
determining values of porosity and water saturation according to the value of the first reservoir sensitive parameter and the value of the second reservoir sensitive parameter;
setting the first reservoir sensitive parameter as an X axis and the second reservoir sensitive parameter as a Y axis, and constructing a rock physical template according to the value of the first reservoir sensitive parameter, the value of the second reservoir sensitive parameter, the corresponding porosity and the value of the water saturation.
In a specific embodiment, based on the first reservoir sensitive parameter and the second reservoir sensitive parameter obtained by calculation, the constructed petrophysical template is shown in fig. 6, where an X axis is the first reservoir sensitive parameter (X), a Y axis is the second reservoir sensitive parameter (Y), both the porosity and the water saturation in the petrophysical template are variable values, the variation range of the porosity is 0 to 0.35, and the variation range of the water saturation is 0 to 1. In fig. 6, the vertical solid line is an equal porosity trend line, the porosity gradually increases from left to right, the porosity values are respectively 0, 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, and 0.35, the horizontal solid line is an equal water saturation trend line, the water saturation gradually increases from bottom to top, and the water saturation values are respectively 0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.0. The filling color in the graph characterizes the numerical size of the hydrocarbon-containing pores, and the calculation formula is HC ═ ρ (1-S)w) Wherein HC is hydrocarbon-bearing pore space, HC high represents a high porosity, low water saturation premium reservoir, and HC low represents a tight or water layer. From fig. 6, it can be seen that the high-porosity, low-water-saturation and high-quality reservoir is located in the lower right corner region of the new petrophysical template (the petrophysical template constructed herein), and the compact layer and the water layer are located in the upper left part region of the new petrophysical template, which indicates that the new petrophysical template can effectively distinguish the hydrocarbon-containing reservoir, the water layer and the compact layer. When the porosity is less than 0.15, corresponding to a medium-low porosity reservoir, the trend lines of different water saturation degrees in the new rock physical template are distributed uniformly and divergently, and the reservoir fluid properties can still be effectively identified, which indicates that the new rock physical template has higher reservoir porosity and water saturation degreeSensitivity, improved quantitative reservoir prediction accuracy, and reduced risk of exploration and development.
As shown in fig. 7, fig. 7 shows a conventional petrophysical template, the same data as those in fig. 6 are used in the construction process, the abscissa in the figure is longitudinal wave impedance, and the ordinate is a longitudinal-transverse wave velocity ratio, so that it can be seen that the distribution of trend lines of different water saturation in the conventional petrophysical template is relatively concentrated, especially for medium-low porosity reservoirs, the trend lines of different water saturation are basically distributed in an overlapping manner, so that the water-bearing reservoir and the oil-bearing reservoir cannot be effectively identified, and the risk of quantitative prediction of the reservoir is increased.
Based on the same inventive concept, a petrophysical template building apparatus is also provided herein, as described in the following embodiments. Because the principle of solving the problems of the rock physical template construction device is similar to that of the rock physical template construction method, the implementation of the rock physical template construction device can refer to the rock physical template construction method, and repeated parts are not described again. The petrophysical template construction device provided by this embodiment includes a plurality of functional modules, which may be implemented by dedicated or general chips, and may also be implemented by software programs, which is not limited herein.
Specifically, as shown in fig. 8, the petrophysical template building apparatus includes:
the elastic parameter calculation module 810 is used for calculating and obtaining elastic parameters of a rock reservoir according to rock physical parameters of the reservoir in the research work area and a biphase medium rock physical theory model, wherein the elastic parameters of the rock reservoir comprise longitudinal wave velocity, transverse wave velocity and density of the reservoir under the conditions of different porosities and different water saturations, and the rock physical parameters of the reservoir in the research work area comprise matrix bulk modulus, matrix shear modulus, matrix density, water bulk modulus, oil bulk modulus, water density, oil density, porosity and water saturations;
the sensitive parameter calculation module 820 is used for constructing a first reservoir sensitive parameter and a second reservoir sensitive parameter which are used for eliminating the influence of a rock framework and keeping the influence of pores and pore fluid according to the elastic parameters of the rock reservoir;
and the building module 830 is configured to build a petrophysical template according to the first reservoir sensitivity parameter and the second reservoir sensitivity parameter.
The petrophysical template constructed by the embodiment has the advantages that the water-containing sandstone and the oil-containing sandstone in the constructed petrophysical template have obvious distribution difference, and the two have strong distinguishability even under the condition of medium and low porosity, so that the precision of reservoir prediction quantitative interpretation is improved, and the exploration risk is reduced.
In an embodiment of the present invention, the elastic parameter calculating module 810, according to the petrophysical parameters of the reservoir in the research work area and the two-phase medium petrophysical theoretical model, calculates the elastic parameters of the rock reservoir to obtain an elastic parameter of the rock reservoir, including: based on the two-phase medium rock physical theory model, the following calculation is carried out:
calculating the volume modulus and the shear modulus of the framework according to the volume modulus, the shear modulus and the porosity of the matrix;
calculating the bulk modulus and the density of the mixed fluid according to the water saturation, the bulk modulus of water, the bulk modulus of oil, the water density and the oil density;
calculating the rock equivalent volume modulus of the saturated two-phase fluid according to the skeleton volume modulus, the porosity, the matrix volume modulus and the mixed fluid volume modulus;
equating the skeleton shear modulus to the rock equivalent shear modulus of saturated two-phase fluid;
calculating the rock equivalent density of the saturated two-phase fluid according to the matrix density, the porosity and the mixed fluid density;
calculating longitudinal wave velocity and transverse wave velocity of the reservoir according to the rock equivalent volume modulus of saturated two-phase fluid, the rock equivalent shear modulus of saturated two-phase fluid and the rock equivalent density of saturated two-phase fluid;
the rock equivalent density of the saturated two-phase fluid is taken as the reservoir density.
The method comprises the following steps of calculating the longitudinal wave velocity and the transverse wave velocity of a reservoir according to the rock equivalent volume modulus of saturated two-phase fluid, the rock equivalent shear modulus of saturated two-phase fluid and the rock equivalent density of saturated two-phase fluid, wherein the longitudinal wave velocity and the transverse wave velocity of the reservoir are calculated by using the following formulas:
Figure BDA0002821276910000111
wherein v ispIs the velocity of longitudinal wave, vsIs the transverse wave velocity, KsatRock equivalent bulk modulus, μ, for saturated two-phase fluidssatRock equivalent shear modulus, ρ, for saturated two-phase fluidssatIs the rock equivalent density of a two-phase fluid.
In an embodiment of this document, the sensitive parameter calculating module 820 constructs a first reservoir sensitive parameter and a second reservoir sensitive parameter that eliminate the influence of the rock skeleton and retain the influence of pores and pore fluids according to the elastic parameters of the rock reservoir, including:
determining a first reservoir sensitive parameter according to the ratio of the reservoir density to the reservoir shear wave velocity;
and determining a second reservoir sensitive parameter according to the difference value of the product of the reservoir longitudinal wave velocity and the reservoir transverse wave velocity and the density.
In one embodiment, determining the first reservoir sensitivity parameter according to the ratio of the reservoir density to the reservoir shear wave velocity includes calculating the first reservoir sensitivity parameter using the following formula:
X=ρ/vs
determining a second reservoir sensitivity parameter according to the difference between the reservoir compressional wave velocity and the reservoir shear wave velocity multiplied by the reservoir density, wherein the second reservoir sensitivity parameter comprises the following formula:
Figure BDA0002821276910000112
where ρ is the reservoir density, vsIs the reservoir shear wave velocity, vpIs the reservoir compressional velocity.
In an embodiment of this document, the constructing module 830 constructs the petrophysical template according to the first reservoir sensitivity parameter and the second reservoir sensitivity parameter, including:
determining values of porosity and water saturation according to the value of the first reservoir sensitive parameter and the value of the second reservoir sensitive parameter;
setting the first reservoir sensitive parameter as an X axis and the second reservoir sensitive parameter as a Y axis, and constructing a rock physical template according to the value of the first reservoir sensitive parameter, the value of the second reservoir sensitive parameter, the corresponding porosity and the value of the water saturation.
In the prediction of the sandstone reservoir with medium and low porosity, aiming at the problem that the resolution capability of the conventional rock physical template for quantitatively explaining the water-containing sandstone reservoir and the oil-containing sandstone reservoir is limited, the rock physical template construction method and the rock physical template construction device provided by the invention construct a new rock physical template based on forward calculation of a biphase medium rock physical theory model, can improve the precision of quantitative explanation of reservoir prediction, and reduce the exploration risk.
In an embodiment herein, there is also provided a computer device, as shown in fig. 9, the computer device 902 may include one or more processors 904, such as one or more Central Processing Units (CPUs), each of which may implement one or more hardware threads. The computer device 902 may also include any memory 906 for storing any kind of information, such as code, settings, data, etc. For example, and without limitation, memory 906 may include any one or more of the following in combination: any type of RAM, any type of ROM, flash memory devices, hard disks, optical disks, etc. More generally, any memory may use any technology to store information. Further, any memory may provide volatile or non-volatile retention of information. Further, any memory may represent fixed or removable components of computer device 902. In one case, when the processor 904 executes the associated instructions, which are stored in any memory or combination of memories, the computer device 902 can perform any of the operations of the associated instructions. The computer device 902 also includes one or more drive mechanisms 908, such as a hard disk drive mechanism, an optical disk drive mechanism, etc., for interacting with any memory.
Computer device 902 may also include an input/output module 910(I/O) for receiving various inputs (via input device 912) and for providing various outputs (via output device 914)). One particular output mechanism may include a presentation device 916 and an associated graphical user interface 918 (GUI). In other embodiments, input/output module 910(I/O), input device 912, and output device 914 may also be excluded, acting as only one computer device in a network. Computer device 902 may also include one or more network interfaces 920 for exchanging data with other devices via one or more communication links 922. One or more communication buses 924 couple the above-described components together.
Communication link 922 may be implemented in any manner, such as over a local area network, a wide area network (e.g., the Internet), a point-to-point connection, etc., or any combination thereof. Communication link 922 may include any combination of hardwired links, wireless links, routers, gateway functions, name servers, etc., governed by any protocol or combination of protocols.
Embodiments herein also provide a computer-readable storage medium having a computer program stored thereon, where the computer program is executed by a processor to execute the steps of the petrophysical template construction method according to any of the above embodiments.
Embodiments herein also provide computer readable instructions, wherein when executed by a processor, the program causes the processor to perform a petrophysical template construction method as described in any of the above embodiments.
It should be understood that, in various embodiments herein, the sequence numbers of the above-mentioned processes do not mean the execution sequence, and the execution sequence of each process should be determined by its function and inherent logic, and should not constitute any limitation to the implementation process of the embodiments herein.
It should also be understood that, in the embodiments herein, the term "and/or" is only one kind of association relation describing an associated object, meaning that three kinds of relations may exist. For example, a and/or B, may represent: a exists alone, A and B exist simultaneously, and B exists alone. In addition, the character "/" herein generally indicates that the former and latter related objects are in an "or" relationship.
Those of ordinary skill in the art will appreciate that the elements and algorithm steps of the examples described in connection with the embodiments disclosed herein may be embodied in electronic hardware, computer software, or combinations of both, and that the components and steps of the examples have been described in a functional general in the foregoing description for the purpose of illustrating clearly the interchangeability of hardware and software. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present disclosure.
It is clear to those skilled in the art that, for convenience and brevity of description, the specific working processes of the above-described systems, apparatuses and units may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
In the several embodiments provided herein, it should be understood that the disclosed system, apparatus, and method may be implemented in other ways. For example, the above-described apparatus embodiments are merely illustrative, and for example, the division of the units is only one logical division, and other divisions may be realized in practice, for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may also be an electric, mechanical or other form of connection.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purposes of the embodiments herein.
In addition, functional units in the embodiments herein may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit. The integrated unit can be realized in a form of hardware, and can also be realized in a form of a software functional unit.
The integrated unit, if implemented in the form of a software functional unit and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solutions of the present invention may be implemented in a form of a software product, which is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the methods described in the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
The principles and embodiments of this document are explained herein using specific examples, which are presented only to aid in understanding the methods and their core concepts; meanwhile, for the general technical personnel in the field, according to the idea of this document, there may be changes in the concrete implementation and the application scope, in summary, this description should not be understood as the limitation of this document.

Claims (10)

1. A petrophysical template construction method, comprising:
calculating to obtain rock reservoir elastic parameters according to reservoir rock physical parameters of a research work area and a two-phase medium rock physical theoretical model, wherein the rock reservoir elastic parameters comprise reservoir longitudinal wave velocity, transverse wave velocity and density under the conditions of different porosities and different water saturation;
constructing a first reservoir sensitive parameter and a second reservoir sensitive parameter which are influenced by eliminating a rock framework and retaining pores and pore fluid according to the rock reservoir elastic parameters;
and constructing a rock physical template according to the first reservoir sensitive parameter and the second reservoir sensitive parameter.
2. The method of claim 1, wherein the studying the zonal reservoir petrophysical parameters comprises: matrix bulk modulus, matrix shear modulus, matrix density, water bulk modulus, oil bulk modulus, water density, oil density, porosity, and water saturation.
3. The method of claim 2, wherein the calculating the elastic parameters of the rock reservoir according to the petrophysical parameters of the reservoir in the research work area and the rock physical theory model of the dual-phase medium comprises: based on the two-phase medium rock physical theory model, the following calculation is carried out:
calculating the matrix bulk modulus and the matrix shear modulus according to the matrix bulk modulus, the matrix shear modulus and the porosity;
calculating the bulk modulus and the density of the mixed fluid according to the water saturation, the bulk modulus of water, the bulk modulus of oil, the water density and the oil density;
calculating rock equivalent bulk modulus of saturated two-phase fluid according to the skeleton bulk modulus, the porosity, the matrix bulk modulus and the bulk modulus of the mixed fluid;
equating the skeletal shear modulus to the rock equivalent shear modulus of a saturated two-phase fluid;
calculating rock equivalent density of saturated two-phase fluid according to the matrix density, the porosity and the mixed fluid density;
calculating longitudinal wave velocity and transverse wave velocity of the reservoir according to the rock equivalent volume modulus of the saturated two-phase fluid, the rock equivalent shear modulus of the saturated two-phase fluid and the rock equivalent density of the saturated two-phase fluid;
the rock equivalent density of the saturated two-phase fluid is taken as the reservoir density.
4. The method of claim 3, wherein calculating reservoir compressional and shear velocities from the petrophysical equivalent bulk modulus of the saturated two-phase fluid, the petrophysical equivalent shear modulus of the saturated two-phase fluid, and the petrophysical equivalent density of the saturated two-phase fluid comprises calculating reservoir compressional and shear velocities using the following formulas:
Figure FDA0002821276900000021
wherein v ispIs the reservoir longitudinal wave velocity, vsIs the reservoir shear wave velocity, KsatIs the rock equivalent bulk modulus, mu, of the saturated two-phase fluidsatIs the rock equivalent shear modulus, ρ, of the saturated two-phase fluidsatIs the rock equivalent density of the two-phase fluid.
5. The method of claim 1, wherein constructing the first reservoir sensitivity parameter and the second reservoir sensitivity parameter which eliminate rock skeleton influence and reserve pore and pore fluid influence according to the rock reservoir elasticity parameter comprises:
determining a first reservoir sensitive parameter according to the ratio of the reservoir density to the reservoir shear wave velocity;
and determining a second reservoir sensitive parameter according to the difference value of the reservoir longitudinal wave velocity and the reservoir transverse wave velocity multiplied by the density.
6. The method of claim 5, wherein determining the first reservoir sensitivity parameter based on the ratio of the reservoir density to the reservoir shear wave velocity comprises calculating the first reservoir sensitivity parameter using the following equation:
X=ρ/vs
and determining a second reservoir sensitive parameter according to the difference value obtained by multiplying the reservoir longitudinal wave velocity and the reservoir transverse wave velocity with the reservoir density respectively, wherein the second reservoir sensitive parameter comprises the following formula:
Figure FDA0002821276900000022
where ρ is the reservoir density, vsIs the reservoir shear wave velocity, vpIs the reservoir compressional velocity.
7. The method of claim 6, wherein constructing a petrophysical template from the first reservoir sensitivity parameter and the second reservoir sensitivity parameter comprises:
determining values of porosity and water saturation according to the value of the first reservoir sensitive parameter and the value of the second reservoir sensitive parameter;
setting the first reservoir sensitive parameter as an X axis and the second reservoir sensitive parameter as a Y axis, and constructing a rock physical template according to the value of the first reservoir sensitive parameter, the value of the second reservoir sensitive parameter, the corresponding porosity and the value of the water saturation.
8. A petrophysical template building apparatus comprising:
the elastic parameter calculation module is used for calculating to obtain elastic parameters of the rock reservoir according to rock physical parameters of the reservoir in the research work area and a two-phase medium rock physical theoretical model, wherein the elastic parameters of the rock reservoir comprise longitudinal wave velocity, transverse wave velocity and density of the reservoir under the conditions of different porosities and different water saturation degrees;
the sensitive parameter calculation module is used for constructing a first reservoir sensitive parameter and a second reservoir sensitive parameter which are used for eliminating the influence of a rock framework and keeping the influence of pores and pore fluid according to the elastic parameters of the rock reservoir;
and the construction module is used for constructing the petrophysical template according to the first reservoir sensitive parameters and the second reservoir sensitive parameters.
9. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the petrophysical template construction method of any one of claims 1 to 7 when executing the computer program.
10. A computer-readable storage medium, characterized in that the computer-readable storage medium stores an executable computer program, which when executed by a processor implements the petrophysical template construction method of any one of claims 1 to 7.
CN202011418723.5A 2020-12-07 2020-12-07 Rock physical template construction method, device, equipment and storage medium Active CN112698401B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011418723.5A CN112698401B (en) 2020-12-07 2020-12-07 Rock physical template construction method, device, equipment and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011418723.5A CN112698401B (en) 2020-12-07 2020-12-07 Rock physical template construction method, device, equipment and storage medium

Publications (2)

Publication Number Publication Date
CN112698401A true CN112698401A (en) 2021-04-23
CN112698401B CN112698401B (en) 2024-03-01

Family

ID=75506307

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011418723.5A Active CN112698401B (en) 2020-12-07 2020-12-07 Rock physical template construction method, device, equipment and storage medium

Country Status (1)

Country Link
CN (1) CN112698401B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115201910A (en) * 2022-07-07 2022-10-18 中国石油化工股份有限公司 Rock physical dynamic template establishing method, system, storage medium and equipment

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109115987A (en) * 2018-07-20 2019-01-01 中国石油天然气股份有限公司 A kind of evaluation method and device of the fluid factor based on petrophysical model
CN111381292A (en) * 2019-07-31 2020-07-07 中国石油天然气股份有限公司 Logging interpretation method and device for predicting sandstone hydrocarbon-bearing reservoir

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109115987A (en) * 2018-07-20 2019-01-01 中国石油天然气股份有限公司 A kind of evaluation method and device of the fluid factor based on petrophysical model
CN111381292A (en) * 2019-07-31 2020-07-07 中国石油天然气股份有限公司 Logging interpretation method and device for predicting sandstone hydrocarbon-bearing reservoir

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
WANG DA-XING: ""A STUDY ON THE ROCK PHYSICS MODEL OF GAS RESERVOIR IN TIGHT SANDSTONE"", 《CHINESE JOURNAL OF GEOPHYSICS》, vol. 60, no. 1, pages 64 - 83 *
何剑: ""双相介质储层参数反演"", 中国优秀硕士学位论文全文数据库 基础科技辑》, no. 2, pages 011 - 819 *
周水生;宜伟;郝召兵;黄为清;伍向阳;: "基于叠前反演的流体敏感属性实验研究及应用", 地球物理学报, no. 06, pages 1985 - 1992 *
黄捍东;汪佳蓓;郭飞;: "敏感参数分析在叠前反演流体识别中的应用", 物探与化探, no. 06, pages 55 - 60 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115201910A (en) * 2022-07-07 2022-10-18 中国石油化工股份有限公司 Rock physical dynamic template establishing method, system, storage medium and equipment
CN115201910B (en) * 2022-07-07 2023-05-26 中国石油化工股份有限公司 Rock physical dynamic template establishment method, system, storage medium and equipment

Also Published As

Publication number Publication date
CN112698401B (en) 2024-03-01

Similar Documents

Publication Publication Date Title
Bisdom et al. Calibrating discrete fracture-network models with a carbonate three-dimensional outcrop fracture network: Implications for naturally fractured reservoir modeling
Joekar‐Niasar et al. Trapping and hysteresis in two‐phase flow in porous media: A pore‐network study
EP1729152B1 (en) Method for changing the scale of absolute permeabilities to build a flow simulation model
Zhang et al. A new practical method for upscaling in highly heterogeneous reservoir models
Illiassov et al. Field-scale characterization of permeability and saturation distribution using partitioning tracer tests: The Ranger field, Texas
Porta et al. Inverse modeling of geochemical and mechanical compaction in sedimentary basins through P olynomial C haos E xpansion
Fanchi Integrated flow modeling
Correia et al. Multiscale integration for Karst-reservoir flow-simulation models
Nakashima et al. Accurate representation of near-well effects in coarse-scale models of primary oil production
Massart et al. Effective flow properties of heterolithic, cross-bedded tidal sandstones: Part 2. Flow simulation
CN112684503B (en) Reservoir fluid identification method and device based on post-stack seismic attributes
CN112698401A (en) Rock physical template construction method, device, equipment and storage medium
Ringrose et al. Use of geology in the interpretation of core-scale relative permeability data
Kaczmaryk et al. Interference pumping tests in a fractured limestone (Poitiers–France): Inversion of data by means of dual-medium approaches
Stephen et al. Reservoir simulations developed from an outcrop of incised valley fill strata
Garrison Jr et al. The Fractal Nature of Geological Data Sets: Power Law Processes Everywhere!
Habrat et al. The concept of a computer system for interpretation of tight rocks using X-ray computed tomography results
Agarwal et al. A streamline-based method for assisted history matching applied to an Arabian Gulf field
Brun et al. Use of prior information in gradient-based history-matching
Li et al. Implicit space-time domain decomposition approach for solving multiphase miscible flow: accuracy and scalability
CN111381279B (en) Reservoir porosity quantitative prediction method and device
US11682167B2 (en) Seamless scaling geomodeling
O'Neill et al. Diagenesis and burial history modeling of heterogeneous marginal marine to shoreface Paleocene glauconitic sandstones, Taranaki Basin, New Zealand
Irving et al. Optimization of uncertain structural parameters with production and observation well data
Irving et al. Optimisation of uncertain structural parameters using production and observation well data

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant