CN118068399B - A method, equipment, medium, and product for quantitatively evaluating the lateral sealing capacity of faults. - Google Patents

A method, equipment, medium, and product for quantitatively evaluating the lateral sealing capacity of faults.

Info

Publication number
CN118068399B
CN118068399B CN202410208589.8A CN202410208589A CN118068399B CN 118068399 B CN118068399 B CN 118068399B CN 202410208589 A CN202410208589 A CN 202410208589A CN 118068399 B CN118068399 B CN 118068399B
Authority
CN
China
Prior art keywords
fault
stress
pressure
plane
lateral
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.)
Active
Application number
CN202410208589.8A
Other languages
Chinese (zh)
Other versions
CN118068399A (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.)
Northeast Petroleum University
Original Assignee
Northeast Petroleum University
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 Northeast Petroleum University filed Critical Northeast Petroleum University
Priority to CN202410208589.8A priority Critical patent/CN118068399B/en
Publication of CN118068399A publication Critical patent/CN118068399A/en
Application granted granted Critical
Publication of CN118068399B publication Critical patent/CN118068399B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/282—Application of seismic models, synthetic seismograms
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

The invention discloses a quantitative evaluation method, equipment, medium and product for lateral sealing capacity of faults, and relates to the field of geological exploration and development of oil and gas resources. According to the method, the fault mud ratio and the effective normal stress of the fault surface at each point on the fault surface are calculated, the fault mud ratio and the effective normal stress of the fault surface are combined, the fault mud ratio with normalized stress is determined, the fault mud ratio with normalized stress and floating pressure are utilized for casting points, a fault lateral sealing capacity evaluation relational expression is fitted, and the lateral sealing capacity of the non-drilling trap control fault is predicted according to the fault lateral sealing capacity evaluation relational expression. The method can comprehensively reflect the control effect of the fault mud content and the effective normal stress of the fault surface on the fault lateral sealing performance, more comprehensively consider the main control factors of the fault lateral sealing performance, improve the accuracy of fault lateral sealing performance evaluation, and greatly reduce the risk of fault trap drilling.

Description

Quantitative evaluation method, equipment, medium and product for lateral sealing capability of fault
Technical Field
The invention relates to the field of oil and gas resource geological exploration and development, in particular to a quantitative evaluation method, equipment, medium and product for lateral sealing capacity of faults.
Background
A number of exploration practices have shown that the lateral sealing capability of faults in hydrocarbon-bearing reservoirs plays a critical role in migration, accumulation and distribution of hydrocarbons. Since rocks in a stratum generally have hydrophilicity, when the pore throat radius of rocks (fault rocks) in a fault is smaller than the pore throat radius of a reservoir, the direction of capillary pressure difference points to the reservoir, and the capillary pressure difference prevents oil gas from migrating across the fault, so that the fault forms a seal, and the seal mechanism of the fault is capillary seal. Thus, from the fault blocking mechanism, the fault lateral blocking capability depends on the displacement pressure of the fractured rock. Numerous studies have demonstrated that one of the determining factors in the displacement pressure of fractured rock is the shale content in the fault zone (i.e., the fault mud content).
Based on the above, at present, two types of fault lateral sealing capability evaluation methods are mainly formed. A direct method for evaluating lateral sealing ability of fault features that the underground fault rock sample is sampled, the displacement pressure of fault rock is measured, and a statistical relation is established between the sampled sample and the parameters (SGR, SSF, CSP) reflecting the content of fault mud, and the experimental method is used to determine the pressure difference of fault by dissecting the drilled fault oil and gas reservoirs. Although the direct method can accurately reflect the lateral sealing capability of the fault, the method is controlled by the limit of the number of underground through-fault coring samples of the oil-gas bearing basin, and has relatively poor applicability, so that the most commonly used fault lateral sealing capability evaluation method in the oil area exploration and development field is an empirical method. However, the two methods mainly consider the influence of the fault mud content on the fault lateral sealing capacity, and in fact, the fault lateral sealing capacity is a result under the comprehensive action of various geological factors, and particularly a large amount of field and experimental data prove that the effective normal stress of a fault plane also has a critical control effect on the fault lateral sealing capacity. Therefore, the lateral sealing capacity of the fault is only evaluated in terms of the content of the fault mud, and the influence of the effective normal stress of the fault surface on the sealing capacity is ignored, so that the accuracy of the lateral sealing capacity evaluation of the fault is often greatly reduced.
Disclosure of Invention
The invention aims to provide a quantitative evaluation method, equipment, medium and product for lateral sealing capability of a fault, which comprehensively reflect the control effect of the content of fault mud and the effective normal stress of a fault surface on the lateral sealing capability of the fault, and improve the accuracy of the evaluation of the lateral sealing capability of the fault surface.
In order to achieve the above object, the present invention provides the following solutions:
a method for quantitatively evaluating lateral closure ability of a fault, the method comprising:
The method comprises the steps of establishing a three-dimensional geological model comprising faults and strata, calculating fault mud ratio at each point on a fault plane according to fault distance and drilling data of faults in the three-dimensional geological model, determining effective normal stress of the fault plane at each point on the fault plane, calculating stress normalized fault mud ratio according to the fault mud ratio at each point on the fault plane and the effective normal stress of the fault plane at each point on the fault plane, determining floating pressure of an actually closed fault of a drilled hydrocarbon reservoir, carrying out casting on the floating pressure by using the stress normalized fault mud ratio, fitting a fault lateral closing capability evaluation relation, and predicting the lateral closing capability of an un-drilled closed loop fault according to the fault lateral closing capability evaluation relation.
Optionally, the calculation formula of the fault mud ratio at each point on the fault plane is:
Wherein SGR is a fault mud ratio, V sh is the clay content of the stratum, deltaZ is the stratum thickness, and D is the fault distance of the fault.
Optionally, determining the effective normal stress of the fault plane at each point on the fault plane specifically includes:
the trend of the stress relief joint in the imaging logging information is indicated to the direction of the horizontal maximum principal stress;
reversely calculating a ground stress construction coefficient by using the horizontal maximum principal stress and the horizontal minimum principal stress corresponding to the test points given by the stratum fracture test data;
Calculating the horizontal maximum principal stress and the horizontal minimum principal stress of the underground rock mass by using a Huang model according to the ground stress construction coefficient;
calculating vertical main stress;
according to the direction of the horizontal maximum principal stress, the horizontal maximum principal stress of the underground rock mass, the horizontal minimum principal stress of the underground rock mass and the vertical principal stress, calculating the effective normal stress of the fault plane at each point on the fault plane according to a formula SN＝(sinα·sinθ)2SH+(cosα·sinθ)2Sh+cosθ2Sv-Pp,, wherein S N is the effective normal stress of the fault plane, S H is the horizontal maximum principal stress, S h is the horizontal minimum principal stress, S v is the vertical principal stress, alpha is the included angle between the trend of the fault plane and the direction of the maximum horizontal principal stress, theta is the inclination angle of the fault plane, and P p is the pore fluid pressure.
Optionally, the calculation formula of the stress normalized fault mud ratio is:
Wherein SSGR is a stress normalized fault mud ratio, SGR is a fault mud ratio, S N is a fault plane effective normal stress, S max is a largest fault plane effective normal stress of all fault plane effective normal stresses extracted on a fault, and S min is a smallest fault plane effective normal stress of all fault plane effective normal stresses extracted on a fault.
Optionally, determining the float pressure of the actual closure of the drilled hydrocarbon reservoir control ring fault specifically includes:
longitudinally dividing the drilled hydrocarbon reservoir profile into a plurality of oil-water units;
Establishing formation pressure-depth profile of each oil-water unit by using formation pressure test data of the well drilling;
Determining a pressure trend line of hydrocarbons and a pressure trend line of a water layer in a stratum pressure-depth profile of each oil-water unit;
the difference between the hydrocarbon pressure on the pressure trend line of the hydrocarbon at the same depth and the hydrostatic pressure on the pressure trend line of the water layer is determined as the float pressure generated by the oil gas at the same depth.
Optionally, the method comprises the steps of carrying out point casting on the fault mud ratio normalized by stress and the floating pressure and fitting out a fault lateral sealing capacity evaluation relational expression, and specifically comprises the following steps:
The functional relation between the fault mud ratio normalized by the stress and the floating pressure for carrying out point casting and fitting the change of the floating pressure along with the fault mud ratio normalized by the stress is that P Sealing device = a x ln (SSGR) -b, wherein P Sealing device is floating pressure, SSGR is the fault mud ratio normalized by the stress, and a and b are coefficients;
According to a functional relation of the float pressure along with the change of the fault mud ratio of the stress normalization, and combining the relation of the float pressure, the height of the hydrocarbon column and the density of the hydrocarbon water, determining a fault lateral sealing capacity evaluation relation as follows: Wherein H Hydrocarbons is the height of a hydrocarbon column with a closeable fault, ρ w is the density of water under the formation condition, ρ o is the density of hydrocarbons under the formation condition, and g is the gravity acceleration.
Optionally, the height of the fault closable hydrocarbon column obtained according to the fault lateral closure capability evaluation relation is used for representing the lateral closure capability of the non-drilling trap control fault.
A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor executing the computer program to perform the steps of the above method for quantitatively evaluating the lateral closure ability of a fault.
A computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the above-described quantitative assessment method of fault lateral closure ability.
A computer program product comprising a computer program which, when executed by a processor, implements the steps of the above-described quantitative assessment method of fault lateral closure capability.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
According to the quantitative evaluation method, the quantitative evaluation equipment, the quantitative evaluation medium and the quantitative evaluation product for the lateral sealing capacity of the fault, the fault mud ratio at each point on the fault surface and the effective normal stress of the fault surface are calculated, the fault mud ratio and the effective normal stress of the fault surface are combined, and the fault mud ratio normalized by the stress is calculated, so that a fault lateral sealing capacity evaluation relational expression is determined, and the lateral sealing capacity of the non-drilling trap fault is predicted. The invention comprehensively reflects the control effect of the fault mud content and the fault surface effective normal stress on the fault lateral sealing performance, and improves the accuracy of the fault surface lateral sealing capability evaluation.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are needed in the embodiments 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 other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic flow chart of a quantitative evaluation method for lateral sealing capability of faults provided in embodiment 1 of the present invention;
Fig. 2 is a distribution diagram of SGR values of a fault plane in a Bohai middling and Xiy region F1 provided in embodiment 1 of the present invention;
Fig. 3 is a cross-sectional view of structural ground stress and hydrostatic pressure in the mid-west depression of the Bohai provided in embodiment 1 of the present invention;
fig. 4 is a distribution diagram of effective normal stress of a fault section of a Bohai Zhongxili area F1 provided in embodiment 1 of the present invention;
FIG. 5 is a cross-sectional view of a Bohai mid-west well region B8-4 reservoir provided in example 1 of the present invention;
FIG. 6 is a stratum pressure-depth profile of the Bohai midfoot-to-Xiy-to-4 reservoir provided in example 1 of the present invention;
fig. 7 is a diagram showing a relationship between the closeable float pressure of a fault in a middle and west depression of Bohai as SSGR according to embodiment 1 of the present invention;
FIG. 8 is a plot of the height of a closeable hydrocarbon column for a trap of the Bohai midfoot and western depression C6-2 in accordance with example 1 of the present invention;
fig. 9 is an internal structural view of the computer device.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
In view of the problems in the background art and the defects and shortcomings in the prior art, the invention provides a novel quantitative evaluation method for the lateral sealing capability of a fault, which is suitable for the geological exploration and development field of a hydrocarbon reservoir, so that the accuracy of the evaluation of the lateral sealing capability of the fault is improved, and the success rate of drilling related trap of the fault is improved.
In order that the above-recited objects, features and advantages of the present invention will become more readily apparent, a more particular description of the invention will be rendered by reference to the appended drawings and appended detailed description.
Example 1
As shown in fig. 1, a quantitative evaluation method for lateral sealing capability of a fault in the embodiment includes the following steps:
and 1, establishing a three-dimensional geological model containing faults and stratum.
And establishing a three-dimensional geological model containing faults and strata by utilizing the seismic interpretation data, inputting drilling data into the three-dimensional geological model, and calculating the breaking distance of the faults by utilizing the projection depth of the two strata of the faults on the fault plane.
And 2, calculating the fault mud ratio at each point on the fault surface according to the fault distance and the drilling data of the fault in the three-dimensional geological model.
According to fault distance and drilling data of faults in the three-dimensional geological model, calculating fault mud ratio distribution at each point on a fault plane by using the following formula:
Wherein SGR is fault mud ratio,% >, V sh is the mud content of the stratum,% >, deltaZ is the stratum thickness, m, and D is the fault interval, m.
And 3, determining effective normal stress of the fault plane at each point on the fault plane.
The determination process of the effective positive stress of each point on the fault plane comprises the following steps:
The method comprises the steps of indicating the direction of horizontal maximum principal stress by the trend of stress relief joints in imaging logging data, calculating the ground stress construction coefficient by means of the horizontal maximum principal stress and the horizontal minimum principal stress corresponding to test points given by stratum fracture test data, calculating the horizontal maximum principal stress and the horizontal minimum principal stress of an underground rock body by means of a Huang model according to the ground stress construction coefficient, calculating the vertical principal stress, calculating the effective normal stress of each point on a fault plane according to the direction of the horizontal maximum principal stress, the horizontal maximum principal stress of the underground rock body, the horizontal minimum principal stress of the underground rock body and the vertical principal stress and a formula SN＝(sinα·sinθ)2SH+(cosα·sinθ)2Sh+cosθ2Sv-Pp,, wherein S N is the effective normal stress of the fault plane, S H is the horizontal maximum principal stress, S h is the horizontal minimum principal stress, S v is the vertical principal stress, mpa, θ is the dip angle of the fault plane, α is the included angle between the trend of the fault plane and the maximum horizontal principal stress, and P p is pore fluid pressure and Mpa.
Illustratively, the vertical principal stress S v experienced by the formation is primarily caused by overburden weight, and is defined as the sum of all overburden weights at a given formation burial depth h. Thus, at depth h, it can be expressed as an integral of overburden density over buried depth:
Sv=∫ρgdh
Wherein S v is vertical main stress, MPa, ρ is stratum density, g/cm 3, g is gravitational acceleration, N/m, and h is stratum burial depth, m. Formation rock density is derived from density log data.
And 4, calculating the fault mud ratio of stress normalization according to the fault mud ratio at each point on the fault plane and the effective normal stress of the fault plane at each point on the fault plane.
On the basis of calculating the fault mud ratio and the fault plane effective normal stress distribution, calculating the numerical value of the fault mud ratio of stress normalization according to the following formula:
Wherein SSGR is the fault mud ratio of stress normalization,%, S max is the effective normal stress of the fault plane with the largest data concentration, and S min is the effective normal stress of the fault plane with the largest data concentration, and Mpa.
And 5, determining the float pressure of the actually closed fracture of the drilled hydrocarbon reservoir.
By carrying out fine dissection on a fault-related oil-gas reservoir drilled in a research area, determining the buoyancy pressure (floating pressure) in a fault-related trap, wherein the pressure generated by buoyancy after oil gas gathers in the trap is the floating pressure, and when the floating pressure generated by the oil gas reaches the limit of the lateral sealing capacity of a fault, leakage occurs, the larger the floating pressure sealed in the trap is, the stronger the fault sealing capacity is, so that the floating pressure in the trap can reflect the lateral sealing capacity of the fault. The specific method for determining the float pressure is as follows:
And (3) establishing a stratum pressure-depth profile by using stratum pressure test data, wherein the depth is the elevation depth after the inclined shaft is straightened and subjected to heart tonifying height correction, the stratum pressure is the abscissa by taking the depth as the ordinate, the stratum pressure-depth profile is obtained by casting points, and the pressure trend lines of hydrocarbon and water layers are determined, and the difference between the hydrocarbon pressure and hydrostatic pressure at the same depth is the floating pressure generated by oil gas.
And 6, carrying out point casting by using the fault mud ratio normalized by the stress and the floating pressure, and fitting out a fault lateral sealing capacity evaluation relation.
Dropping points by SSGR and fault-sealing floating pressure (P Sealing device ) and fitting a fault-sealing failure outer envelope curve representing the fault-sealing maximum floating pressure under a certain SSGR value condition and an expression representing the fault-sealing failure outer envelope curve functional relation, wherein the functional relation of the fault-sealing floating pressure along with SSGR value is as follows:
P Sealing device =a*ln(SSGR)-b
Wherein P Sealing device is the size of the fault closeable float pressure, MPa, a and b are constants related to the region, and can be obtained through fitting the actual SSGR value and float pressure data.
In the process of continuously filling oil gas in the fault trap, the float pressure generated by oil gas aggregation is gradually increased, when the float pressure is lower than the fault rock displacement pressure, the limit of fault sealing capability is not reached, at the moment, the float pressure generated by oil gas aggregation in the fault trap can not completely reflect the size of the fault rock displacement pressure, when the oil gas float pressure is equal to the fault rock displacement pressure along with continuous filling of the oil gas, the limit of fault sealing capability is reached, at the moment, the height of oil gas aggregation in the fault trap is the height of a hydrocarbon column which can be sealed by a fault, according to the principle, the relation between the oil gas float pressure and the height of the hydrocarbon column and the density of the hydrocarbon water are utilized, and the height of the hydrocarbon column which can be sealed by the fault can be calculated according to the following formula by combining a function relation of the float pressure changing along with SSGR values:
Wherein H Hydrocarbons is the height of a hydrocarbon column with a closeable fault, m, ρ w is the density of water under the stratum condition, kg/m 3;ρo is the density of hydrocarbons under the stratum condition, kg/m 3, g is the gravity acceleration, and m/s 2.
And 7, predicting the lateral closure capacity of the non-drilling trap control fault according to the fault lateral closure capacity evaluation relation.
According to the technical scheme adopted by the invention, a new parameter-stress normalized fault mud ratio (SSGR) capable of comprehensively reflecting the fault mud content and the fault surface effective normal stress on the fault lateral sealing control function is established according to a fault lateral sealing mechanism and a sealing capacity main control factor, then the floating pressure size of fault sealing is determined by finely dissecting a large number of drilled fault-related hydrocarbon reservoirs in a research area, finally a fault lateral sealing capacity evaluation relation which represents the fault sealing floating pressure size under a certain improved fault mud ratio is established by using a statistical method, and the height of a fault sealing hydrocarbon column can be quantitatively evaluated according to the fault lateral sealing capacity evaluation relation.
Compared with the prior art, the method has the beneficial effects that the influence factors of the fault lateral sealing capability are complex, a large amount of field and experimental data prove that the effective normal stress of the fault plane has a decisive effect on the fault lateral sealing capability, and the evaluation of the fault lateral sealing capability only from the aspect of the fault mud content often leads to the great reduction of the evaluation accuracy. The invention establishes a new evaluation method capable of comprehensively reflecting the fault mud content and the fault surface effective normal stress on the fault lateral sealing control effect, more comprehensively considers the main control factors of the fault lateral sealing capability, and has higher accuracy, higher evaluation result and actual anastomosis rate compared with the existing SGR method, and greatly reduces the risk of fault trap drilling.
The present invention will be described in further detail with reference to the specific examples of the lateral blocking capability of the fault in Bohai Zhongxiqing region of Bohai Bay basin, but the present invention is not limited to the examples.
A quantitative evaluation method for lateral blocking capacity of faults comprises the following specific steps:
firstly, utilizing seismic interpretation data of a research area, using TRAPTESTER software developed by Badleys company to establish a three-dimensional geological model containing faults and strata, simultaneously inputting drilling data including well coordinates, well tracks, geological layering and stratum argillaceous content data into the three-dimensional geological model, and utilizing projection depth of two strata of the faults on the fault plane to calculate the breaking distance of the faults, wherein the relative distance between the two equivalent strata of the faults is the breaking distance.
And secondly, calculating the fault mud ratio distribution at each point on the fault surface according to a formula (1) by using the fault distance of the fault and the drilling data of the fault adjacent drilling in the three-dimensional geological model established by TRAPTESTER software.
Wherein SGR is fault mud ratio,% >, V sh is the mud content of the stratum,% >, deltaZ is the stratum thickness, m, and D is the fault interval, m.
The TRAPTESTER software inputs the formula 1 into the software, and the SGR distribution of fault planes can be directly calculated in the TRAPTESTER software three-dimensional geological model, for example, the SGR distribution of fault planes in F1 in the mid-west depression area of Bohai Bay basin is shown in fig. 2.
And thirdly, calculating the effective positive stress of the fault plane according to the following steps.
The effective normal stress of the fault surface is calculated by determining the magnitude and the direction of the ground stress, and the stress state of the underground rock mass can be represented by three main stresses, namely, the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress.
Firstly, determining that the horizontal maximum main stress direction of Bohai middle and West depression areas is NE 70-80 degrees by using imaging logging data. The orientations of the well wall breakout, the fracturing induced seam and the stress relief seam are given in the imaging logging data, the orientations of the stress relief seam indicate the maximum horizontal main stress direction, and the elliptical well long axis caused by the well wall breakout indicates the minimum horizontal main stress direction.
The direction of the vertical earth stress is consistent with the direction of the gravity, and the directions of the horizontal earth stress are mutually perpendicular. The vertical main stress of the stratum is mainly caused by the gravity of the overburden stratum, and the direction of the gravity is consistent, and the vertical main stress is downward. The direction of the horizontal minimum main stress and the direction of the horizontal maximum main stress are mutually perpendicular on the same horizontal plane, and the direction of the horizontal maximum main stress is determined, namely the direction of the horizontal minimum main stress is increased or decreased by 90 degrees.
And secondly, calculating the vertical ground stress of the research area, wherein the vertical ground stress is mainly caused by the gravity of the overlying strata, and integrating the function representing the density through the burial depth of the strata can obtain the vertical main stress. And integrating a function representing the density through the burial depth of the rock stratum by utilizing the density logging data of a plurality of wells in Bohai Chinese and Western depression areas, and establishing the change relation (such as formula (2)) of the rock density along with the burial depth of the rock stratum.
Sv=0.021854*h-2.288398 (2)
Wherein S v is vertical main stress, MPa, h is the burial depth of the rock stratum, and m.
And then, calculating the horizontal principal stress by using the small fracturing test data and the array acoustic logging data of the Bohai Chinese and Western depression area. The level of the horizontal main stress can be generally obtained by interpretation of small fracturing test data, and the method is simple to operate, strong in adaptability and high in accuracy. However, the Bohai mid-western depression area has relatively less small fracturing test data, and the horizontal main stress can be obtained by combining an array acoustic logging calculation method. According to the method, firstly, static rock mechanical parameters are calculated, experimental researches show that the static rock mechanical parameters and dynamic rock mechanical parameters have a certain linear relation, the dynamic rock mechanical parameters can be obtained through dynamic and static conversion, and the dynamic poisson ratio is converted through reference of experimental results of scholars such as Lin Yingsong. And then comprehensively considering and selecting a Huang model to calculate the horizontal main stress, calculating a ground stress structural coefficient, and obtaining the change relation of the horizontal main stress along with the depth, thereby establishing a ground stress section (figure 3) of the Bohai Chinese and Western depression area, and determining the magnitude of the three-dimensional main stress according to the ground stress section.
The dynamic rock mechanical parameters can be obtained through the conversion of the propagation speed of sound waves in the rock sample, and the dynamic rock mechanical parameters comprise dynamic poisson ratio and dynamic Young modulus.
The calculation formulas of the dynamic Young's modulus and the dynamic Poisson's ratio are respectively:
Wherein E d is dynamic Young modulus, GPa, ρ b is rock density, g/cm 3;Δts is transverse wave time difference, μs/m, Δt p is longitudinal wave time difference, μs/m, and v d is dynamic Poisson's ratio.
Dynamic rock mechanical parameter dynamic and dynamic conversion refers to conversion of dynamic poisson ratio into static poisson ratio, and the conversion formula is as follows:
μs＝0.1268+0.250μd
Where μ s and μ d are static poisson's ratio and dynamic poisson's ratio (dimensionless), respectively.
Comprehensively considering and selecting a Huang model to calculate the horizontal main stress, calculating the structural coefficient of the ground stress, and obtaining the concrete flow of the horizontal main stress changing relation along with the depth:
The formula of the Huang's model calculation horizontal principal stress is as follows:
σv＝∫GzdH
wherein, sigma v、σH、σh -overburden pressure, maximum and minimum horizontal main ground stress, MPa, beta and gamma-construction stress coefficient.
Before the ground stress is calculated by using the horizontal principal stress formula, the ground stress construction coefficient is reversely calculated by using the horizontal maximum principal stress and the horizontal minimum principal stress corresponding to the test points given by the stratum fracture test data. The calculation formula of the ground stress construction coefficient is as follows:
substituting the horizontal maximum principal stress and the horizontal minimum principal stress corresponding to the test points given by the stratum fracture test data into a calculation formula of the ground stress construction coefficient, and calculating to obtain the ground stress construction coefficient of the research area.
Finally, on the basis of determining the three-dimensional main stress and the azimuth and the formation fluid pressure of the Bohai Chinese and Western depression area, inputting a three-dimensional main stress, fluid pressure and fault plane effective normal stress calculation formula (3)) into the established Bohai Chinese and Western depression area three-dimensional fault and formation model, and calculating the effective normal stress of any point on the three-dimensional fault plane. Fig. 4 shows the calculated effective normal stress distribution of the fault fracture surface of the Bohai midfoot and western depression area F1 in the three-dimensional fault and stratum model.
SN＝(sinα·sinθ)2SH+(cosα·sinθ)2Sh+cosθ2Sv-Pp (3)
Wherein S N is the effective normal stress of the fault plane, S H is the horizontal maximum main stress, mpa, S h is the horizontal minimum main stress, mpa, S v is the vertical main stress, mpa, θ is the dip angle of the fault plane, α is the included angle between the trend of the fault plane and the maximum horizontal main stress, and P p is the pore fluid pressure, MPa.
And fourthly, calculating SSGR values of the oil deposit control ring faults drilled in the Bohai mid-west depression area according to a formula (4) on the basis of calculating the fault mud ratio and the effective positive stress distribution of the fault surface.
Wherein SSGR is the stress normalized fault mud ratio,%, "S max" is the largest fault plane effective normal stress among all fault plane effective normal stresses extracted on the fault, and "S min" is the smallest fault plane effective normal stress among all fault plane effective normal stresses extracted on the fault.
And fifthly, carrying out fine dissection on a fault-related oil-gas reservoir drilled in the Bohai Chinese and western depression area, determining the floating pressure in the fault-related trap, wherein the pressure generated by buoyancy after the oil gas is gathered in the trap is the floating pressure, and when the floating pressure generated by the oil gas reaches the limit of the lateral sealing capacity of the fault, leakage occurs, the larger the floating pressure sealed in the trap is, the stronger the fault sealing capacity is, so that the floating pressure in the trap can reflect the lateral sealing capacity of the fault. Taking a B8-4 oil reservoir of a liberal pottery group in Bohai China and Xiyou as an example, determining the actual closed floating pressure of faults, wherein the specific process is as follows:
FIG. 5 is a cross-sectional view of a B8-4 reservoir, which may be divided longitudinally into a plurality of oil-water units. And (3) establishing a stratum pressure-depth profile of each oil-water unit by using stratum pressure test data of the B8-4-B well and the B8-4-G well (figure 6), wherein the depth in the profile is the elevation depth after straightening of the inclined well and after heart tonifying height correction, the stratum pressure is taken as an ordinate, the stratum pressure is taken as an abscissa, the stratum pressure-depth profile is obtained through casting points, the pressure trend lines of hydrocarbon and water layers are determined, and the difference between the hydrocarbon pressure and hydrostatic pressure at the same depth is the floating pressure generated by oil gas. According to the calculation method of the actual closed float pressure of the faults, fine dissection is carried out on fault oil reservoirs drilled in Bohai Chinese and Western depression areas, the actual closed float pressure of each control circle fault is counted, and data support is provided for subsequent fault lateral closed capacity evaluation.
And sixthly, after determining the value and the actual closed float size of each control loop fault SSGR, carrying out point casting by utilizing SSGR and fault closed float (P Sealing device ) and fitting a fault closed failure outer envelope curve representing the fault closed maximum float under the condition of a certain SSGR value, and an expression (5) for representing the fault closed failure outer envelope curve functional relation, wherein the fault closed float is a functional relation of the fault closed float changing along with the SSGR value. The relationship graph of the closed floating pressure of the faults in Bohai mid-west depression along with SSGR is shown in fig. 7, and the points with different gray scales in fig. 7 represent different effective positive stresses of the fault planes.
P Sealing device =0.2228*ln(SSGR)-0.4358 (5)
Wherein P Sealing device is the size of the fault closeable float pressure and MPa.
In the process of continuously filling oil gas in the fault trap, the float pressure generated by oil gas aggregation is gradually increased, when the float pressure is lower than the fault rock displacement pressure, the limit of fault sealing capability is not reached, at the moment, the float pressure generated by oil gas aggregation in the fault trap can not completely reflect the size of the fault rock displacement pressure, when the oil gas float pressure is equal to the fault rock displacement pressure along with continuous filling of the oil gas, the limit of fault sealing capability is reached, at the moment, the height of oil gas aggregation in the fault trap is the height of a hydrocarbon column which can be sealed by faults, and according to the principle, the relation between the oil gas float pressure, the height of the hydrocarbon column and the density of hydrocarbon water is utilized, and a relation (5) is combined, so that a fault lateral sealing capability evaluation relation (formula 6) suitable for a Bohai mid-west depression area can be deduced. The height of the fault-closeable hydrocarbon column can be evaluated by using the relation.
Wherein H Hydrocarbons is the height of a hydrocarbon column with a closeable fault, m, ρ w is the density of water under the stratum condition, kg/m 3;ρo is the density of hydrocarbons under the stratum condition, kg/m 3, g is the gravity acceleration, and m/s 2.
Seventh, in order to verify the reliability and accuracy of the established fault lateral sealing capability evaluation method, C6-2 trap control faults in Bohai Chinese and Western depression areas are selected to evaluate the lateral sealing capability. The C6-2 trap contains two oil-containing sand layers, ng- ④ and Ng- ⑤. The result of the fault surface attribute calculation shows that the SGR values of the Ng- ④ sand layer and the Ng- ⑤ sand layer are respectively 32.3-40.4% and 24.0-31.4%, and the effective normal stress of the fault surface is respectively 12.2-12.9MPa and 13.3-13.9MPa. On the basis, SSGR value distribution is calculated according to a formula (4), and further the heights of hydrocarbon columns which can be closed by C6-2 trap faults are evaluated according to a formula (6), wherein the heights of hydrocarbon columns which can be closed by Ng- ④ sand layers and Ng- ⑤ sand layers are respectively 37.8m and 26.6m, and the corresponding oil-water interfaces are respectively 1662.8m and 1756.6m (figure 8). From the evaluation results, it can be seen that the hydrocarbon column height calculated by the present fault blocking capability evaluation method substantially coincides with the actual hydrocarbon column height. Therefore, the evaluation method is reflected to have strong applicability and reliability in practical application.
Example 2
A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor executing the computer program to perform the steps of the method of quantitatively evaluating lateral closure ability of a fault in embodiment 1.
A computer device may be internally structured as shown in fig. 9. The computer device includes a processor, a memory, an Input/Output interface (I/O) and a communication interface. The processor, the memory and the input/output interface are connected through a system bus, and the communication interface is connected to the system bus through the input/output interface. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, computer programs, and a database. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage media. The database of the computer device is used to store the pending transactions. The input/output interface of the computer device is used to exchange information between the processor and the external device. The communication interface of the computer device is used for communicating with an external terminal through a network connection. The computer program is executed by a processor to implement the method for quantitatively evaluating the lateral closure ability of a fault in embodiment 1.
It should be noted that, the object information (including, but not limited to, object device information, object personal information, etc.) and the data (including, but not limited to, data for analysis, stored data, presented data, etc.) related to the present invention are both information and data authorized by the object or sufficiently authorized by each party, and the collection, use and processing of the related data need to comply with the related laws and regulations and standards of the related countries and regions.
Those skilled in the art will appreciate that implementing all or part of the above described methods may be accomplished by way of a computer program stored on a non-transitory computer readable storage medium, which when executed, may comprise the steps of the embodiments of the methods described above. Any reference to memory, database, or other medium used in embodiments provided herein may include at least one of non-volatile and volatile memory. The nonvolatile Memory may include Read-only Memory (ROM), magnetic tape, floppy disk, flash Memory, optical Memory, high density embedded nonvolatile Memory, resistive random access Memory (ReRAM), magneto-resistive random access Memory (Magnetoresistive Random Access Memory, MRAM), ferroelectric Memory (Ferroelectric Random Access Memory, FRAM), phase change Memory (PHASE CHANGE Memory, PCM), graphene Memory, and the like. Volatile memory can include random access memory (Random Access Memory, RAM) or external cache memory, and the like. By way of illustration, and not limitation, RAM can be in the form of a variety of forms, such as Static Random access memory (Static Random access memory AccessMemory, SRAM) or dynamic Random access memory (Dynamic Random Access Memory, DRAM), and the like. The databases referred to in the embodiments provided herein may include at least one of a relational database and a non-relational database. The non-relational database may include, but is not limited to, a blockchain-based distributed database, and the like. The processor referred to in the embodiments provided in the present invention may be a general-purpose processor, a central processing unit, a graphics processor, a digital signal processor, a programmable logic unit, a data processing logic unit based on quantum computing, or the like, but is not limited thereto.
Example 3
A computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the fault lateral closure ability quantitative evaluation method in embodiment 1.
Example 4
A computer program product comprising a computer program which, when executed by a processor, implements the steps of the quantitative assessment method of fault lateral closure ability in embodiment 1.
The technical features of the above embodiments may be arbitrarily combined, and all possible combinations of the technical features in the above embodiments are not described for brevity of description, however, as long as there is no contradiction between the combinations of the technical features, they should be considered as the scope of the description.
The principles and embodiments of the present invention have been described herein with reference to specific examples, which are intended to facilitate an understanding of the principles and concepts of the invention and are to be varied in scope and detail by persons of ordinary skill in the art based on the teachings herein. In view of the foregoing, this description should not be construed as limiting the invention.

Claims (8)

1. A method for quantitatively evaluating lateral blocking capability of a fault, the method comprising:
establishing a three-dimensional geological model containing faults and stratum;
Calculating the fault mud ratio at each point on the fault surface according to the fault distance and the drilling data of the fault in the three-dimensional geological model, wherein the calculation formula of the fault mud ratio at each point on the fault surface is as follows: x 100%, wherein SGR is fault mud ratio, For the shale content of the formation,Z is the stratum thickness, D is the fault distance of the fault;
determining the effective normal stress of the fault plane at each point on the fault plane;
Calculating the fault mud ratio of stress normalization according to the fault mud ratio at each point on the fault plane and the effective normal stress of the fault plane at each point on the fault plane;
determining the float pressure of the actually closed fault of the drilled oil and gas reservoir control ring;
performing point casting by using the fault mud ratio of stress normalization and the floating pressure and fitting out a fault lateral sealing capacity evaluation relation;
predicting the lateral closure capacity of the non-drilling trap control fault according to the fault lateral closure capacity evaluation relation;
determining the effective normal stress of each point on the fault plane, which specifically comprises the following steps:
the trend of the stress relief joint in the imaging logging information is indicated to the direction of the horizontal maximum principal stress;
reversely calculating a ground stress construction coefficient by using the horizontal maximum principal stress and the horizontal minimum principal stress corresponding to the test points given by the stratum fracture test data;
Calculating the horizontal maximum principal stress and the horizontal minimum principal stress of the underground rock mass by using a Huang model according to the ground stress construction coefficient;
calculating vertical main stress;
according to the direction of the horizontal maximum principal stress, the horizontal maximum principal stress of the underground rock mass, the horizontal minimum principal stress of the underground rock mass and the vertical principal stress, calculating the effective normal stress of the fault plane at each point on the fault plane according to a formula SN=(sinα·sinθ)2SH+(cosα·sinθ)2Sh+cosθ2Sv-Pp,, wherein S N is the effective normal stress of the fault plane, S H is the horizontal maximum principal stress, S h is the horizontal minimum principal stress, S v is the vertical principal stress, alpha is the included angle between the trend of the fault plane and the direction of the maximum horizontal principal stress, theta is the inclination angle of the fault plane, and P p is the pore fluid pressure.
2. The quantitative evaluation method for the lateral sealing capability of a fault according to claim 1, wherein a calculation formula of the stress normalized fault mud ratio is:
;
Wherein SSGR is a stress normalized fault mud ratio, SGR is a fault mud ratio, S N is a fault plane effective normal stress, S max is a largest fault plane effective normal stress of all fault plane effective normal stresses extracted on a fault, and S min is a smallest fault plane effective normal stress of all fault plane effective normal stresses extracted on a fault.
3. The quantitative evaluation method for the lateral sealing capability of a fault according to claim 1, wherein the determination of the float pressure of the actual sealing of the drilled hydrocarbon reservoir control ring fault specifically comprises the following steps:
longitudinally dividing the drilled hydrocarbon reservoir profile into a plurality of oil-water units;
Establishing formation pressure-depth profile of each oil-water unit by using formation pressure test data of the well drilling;
Determining a pressure trend line of hydrocarbons and a pressure trend line of a water layer in a stratum pressure-depth profile of each oil-water unit;
the difference between the hydrocarbon pressure on the pressure trend line of the hydrocarbon at the same depth and the hydrostatic pressure on the pressure trend line of the water layer is determined as the float pressure generated by the oil gas at the same depth.
4. The quantitative evaluation method for the lateral sealing capability of a fault according to claim 1, wherein the method for evaluating the lateral sealing capability of the fault by using the stress normalized fault mud ratio and the floating to perform point casting and fitting comprises the following specific steps:
The functional relation of the fault mud ratio normalized by the stress and the floating to perform point casting and fitting the floating pressure to the fault mud ratio variation normalized by the stress is that P Sealing device =a (SSGR) -b, wherein P Sealing device is float pressure, SSGR is fault mud ratio of stress normalization, and a and b are coefficients;
According to a functional relation of the float pressure along with the change of the fault mud ratio of the stress normalization, and combining the relation of the float pressure, the height of the hydrocarbon column and the density of the hydrocarbon water, determining a fault lateral sealing capacity evaluation relation as follows: X 10 6, wherein H hydrocarbon is the height of a hydrocarbon column with a closeable fault, rho w is the density of water under the formation condition, rho o is the density of hydrocarbons under the formation condition, and g is gravity acceleration.
5. The quantitative evaluation method for the lateral closure ability of a fault according to claim 4, wherein the height of a hydrocarbon column which can be closed by the fault and is obtained according to the evaluation relation of the lateral closure ability of the fault is used for representing the lateral closure ability of an undrilled trap fault.
6. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor executes the computer program to carry out the steps of the method for quantitative evaluation of lateral closure ability of a fault as claimed in any one of claims 1-5.
7. A computer-readable storage medium, on which a computer program is stored, characterized in that the computer program, when being executed by a processor, carries out the steps of the method for quantitative evaluation of lateral closure ability of a fault as claimed in any one of claims 1 to 5.
8. A computer program product comprising a computer program, characterized in that the computer program, when executed by a processor, implements the steps of the method for quantitative evaluation of lateral closure ability of a fault as claimed in any one of claims 1-5.
CN202410208589.8A 2024-02-26 2024-02-26 A method, equipment, medium, and product for quantitatively evaluating the lateral sealing capacity of faults. Active CN118068399B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202410208589.8A CN118068399B (en) 2024-02-26 2024-02-26 A method, equipment, medium, and product for quantitatively evaluating the lateral sealing capacity of faults.

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202410208589.8A CN118068399B (en) 2024-02-26 2024-02-26 A method, equipment, medium, and product for quantitatively evaluating the lateral sealing capacity of faults.

Publications (2)

Publication Number Publication Date
CN118068399A CN118068399A (en) 2024-05-24
CN118068399B true CN118068399B (en) 2026-03-24

Family

ID=91108676

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202410208589.8A Active CN118068399B (en) 2024-02-26 2024-02-26 A method, equipment, medium, and product for quantitatively evaluating the lateral sealing capacity of faults.

Country Status (1)

Country Link
CN (1) CN118068399B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN119129906A (en) * 2024-08-26 2024-12-13 中国石油大学(北京) A method, device and equipment for evaluating oil and gas transport capacity
CN120633260B (en) * 2025-08-13 2026-01-02 中国地质大学（北京） Fault trap capability evaluation method, device, equipment and storage medium

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105760668A (en) * 2016-02-08 2016-07-13 东北石油大学 Quantitative evaluation method for lateral sealing of fault
CN113466931A (en) * 2021-05-31 2021-10-01 中国石油天然气股份有限公司 Geological period fault sealing evaluation method and system

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105866835B (en) * 2016-03-28 2018-11-02 中国石油大学(华东) A kind of tomography three dimensional closure quantitative evaluation method based on crustal stress distribution
CN110321595B (en) * 2019-06-06 2020-03-06 中国石油大学(华东) Fault sealing evaluation method for extracting static quality coefficient by logging
CN114185083B (en) * 2021-12-07 2022-11-22 成都北方石油勘探开发技术有限公司 Quantitative evaluation method for fault sealing in clastic rock stratum
CN117233856A (en) * 2023-06-21 2023-12-15 中海石油(中国)有限公司 Sandstone reservoir fault sealing evaluation method and system based on three-dimensional geological modeling

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105760668A (en) * 2016-02-08 2016-07-13 东北石油大学 Quantitative evaluation method for lateral sealing of fault
CN113466931A (en) * 2021-05-31 2021-10-01 中国石油天然气股份有限公司 Geological period fault sealing evaluation method and system

Also Published As

Publication number Publication date
CN118068399A (en) 2024-05-24

Similar Documents

Publication Publication Date Title
US8548782B2 (en) Method for modeling deformation in subsurface strata
CN118068399B (en) A method, equipment, medium, and product for quantitatively evaluating the lateral sealing capacity of faults.
US20230097859A1 (en) Method and system for determining coarsened grid models using machine-learning models and fracture models
US10167704B2 (en) Reservoir pressure prediction and methods of execution
CN107784150B (en) A method for predicting fracture closure
WO2010047859A1 (en) Method for modeling deformation in subsurface strata
CN120654502B (en) A Computer Simulation-Based Design Optimization Method for Multi-Level Grouting in Overburden
CN114185083B (en) Quantitative evaluation method for fault sealing in clastic rock stratum
US20240093593A1 (en) Automated cell-to-cell calibration of subsidence information map in forward geological models
CN108301825A (en) A kind of high temperature and pressure reservoir pore pressure analysis method
CN117390881B (en) Fracturing design method and device, electronic equipment and storage medium
Haldorsen et al. Review of the stochastic nature of reservoirs
CN115324560A (en) A method for determining the position of oil and gas casing deformation induced by fracturing using in-situ stress field simulation
Ebrahimian et al. Evaluation of undrained shear strength of marine clay using cone penetration resistance at South Pars field in Iran
Gillespie et al. Introduction: geology of fractured reservoirs
Ringrose et al. Reservoir model types
CN114114411B (en) Quantitative determination method and device for three-dimensional fault permeability
Agung et al. Expansive Residual Soil Stability Behavior During Wetting and Drying Process.
CN114154285A (en) Method and device for predicting pressure coefficient of organic-rich shale formation fluid
Nolting et al. The preservation of water-table caves at depth: Observations from subsurface data and numerical modeling
CN111650644B (en) Quantitative prediction method for salt cemented sandstone
CN120633259B (en) Multilayer stacking-dumping type reservoir landslide prediction method and terminal equipment
CN114417587B (en) Stratum pore pressure prediction method suitable for compact clastic rock stratum in low-pressure state
CN119559351B (en) Method, system, medium and equipment for three-dimensional modeling of porosity of carbonate reservoir
Bentley et al. Reservoir Model Types

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