CN112943198A - Method for calculating inhomogeneous stress field of deep shale complex formation - Google Patents

Method for calculating inhomogeneous stress field of deep shale complex formation Download PDF

Info

Publication number
CN112943198A
CN112943198A CN202110236279.3A CN202110236279A CN112943198A CN 112943198 A CN112943198 A CN 112943198A CN 202110236279 A CN202110236279 A CN 202110236279A CN 112943198 A CN112943198 A CN 112943198A
Authority
CN
China
Prior art keywords
stress
stratum
reservoir
under
target area
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
CN202110236279.3A
Other languages
Chinese (zh)
Other versions
CN112943198B (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.)
China Petroleum and Chemical Corp
Southwest Petroleum University
Original Assignee
Southwest 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 Southwest Petroleum University filed Critical Southwest Petroleum University
Priority to CN202110236279.3A priority Critical patent/CN112943198B/en
Publication of CN112943198A publication Critical patent/CN112943198A/en
Application granted granted Critical
Publication of CN112943198B publication Critical patent/CN112943198B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/25—Methods for stimulating production
    • E21B43/26—Methods for stimulating production by forming crevices or fractures
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00—Survey of boreholes or wells
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00—Survey of boreholes or wells
    • E21B47/06—Measuring temperature or pressure
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Forestry; Mining

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Geology (AREA)
  • Physics & Mathematics (AREA)
  • Geochemistry & Mineralogy (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Fluid Mechanics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geophysics (AREA)
  • Business, Economics & Management (AREA)
  • Primary Health Care (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Human Resources & Organizations (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Marketing (AREA)
  • General Health & Medical Sciences (AREA)
  • Economics (AREA)
  • Health & Medical Sciences (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Animal Husbandry (AREA)
  • Agronomy & Crop Science (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a method for calculating a heterogeneous stress field of a deep shale complex formation, which comprises the steps of calculating a triaxial stress gradient of the formation under local coordinates of each well position; calculating a stratum triaxial stress value under the middle local coordinate of the reservoir at each well position; calculating the stratum stress tensor component under the central global coordinate of the reservoir at each well position; calculating the stratum stress tensor component under the central global coordinate of the reservoir at any position in the target area; and calculating the stratum triaxial stress value and the maximum principal stress direction under the middle local coordinate of the reservoir at any position in the target area. According to the method, large-scale geological structure mechanics modeling is not needed, the non-uniform distribution condition of the stratum triaxial stress field at any position in the target area can be calculated only by means of limited data such as stratum triaxial stress values and the maximum horizontal main stress direction of drilled wells in the target area, and the problems that the stress field of the complex stratum of the deep shale changes rapidly in a small area and a rapid and accurate non-uniform stratum stress field calculation method is lacked are solved.

Description

Method for calculating inhomogeneous stress field of deep shale complex formation
Technical Field
The invention relates to a method for calculating a heterogeneous stress field of a deep shale complex formation, and belongs to the technical field of shale gas development.
Background
Hydraulic fracturing is one of the major stimulation practices for oil and gas reservoirs. In recent years, the method is widely applied to development of unconventional oil and gas reservoirs represented by shale gas reservoirs. In the shale fracturing process, the minimum level main stress value of the stratum is closely related to the fracturing construction pressure; determining an extension path of the hydraulic fracture according to the maximum horizontal principal stress direction of the stratum; the stress difference affects the construction problem of the complex fracture network, and further affects the production effect of the fractured oil-gas well. Therefore, the hydraulic fracturing construction pressure range and the hydraulic fracture extension path can be accurately predicted through stratum stress field calculation, and guidance is provided for deep shale fracturing optimization design and site construction.
Because of the different degrees of severity of the movement of the geological structure, the distribution of the stratum stress is different. For areas where formation motion is gradual, such as successive sedimentary zones of shallower formations, formation motion has less of an effect on formation stress. The traditional stratum stress calculation method is mainly characterized in that the minimum horizontal principal stress value of a stratum is accurately measured through a hydraulic fracturing test method, the maximum horizontal principal stress value of the stratum is calculated according to a construction coefficient, the orientation of a well wall induced crack is observed through imaging logging, and the direction of the maximum horizontal principal stress is determined. In 2008, the poplar and the Liugong Biao establish a stratum level principal stress explanation chart by utilizing a stratum stress mechanical model under a complex construction condition on the basis of inversion of stratum rock parameters by logging data. In 2016, Yangxiangqiang et al determined the functional relationship between the stratum stress and the rock parameters and the lateral pressure coefficient by geological three-dimensional modeling and adopting orthogonal numerical analysis and genetic programming algorithm, and inversely calculated the stratum stress field. In 2020, the Dong et al propose a three-dimensional modeling-based stratum stress field inversion method, which inverts the stratum stress distribution characteristics by establishing a large-scale three-dimensional geological structure numerical model and by parallel calculation. In 2021, Lijing et al characterized the geologic structure characteristics based on a box-dimension method, and further calculated the formation stress by optimal inversion of boundary constraints using a neural network method.
In summary, most of the existing stratum stress calculation methods mainly aim at the situations that the movement of a geological structure is smooth and the geological structure is simple; the geological structure of the deep shale stratum is usually complex, the stratum principal stress value and the principal stress direction change rapidly in a small area, a large-scale three-dimensional geological structure numerical model needs to be established, model calibration is carried out by depending on a large amount of field measurement data, the workload is large, and the time consumption is long. Therefore, a calculation method for the inhomogeneous stress field of the stratum with the deep shale complex structure needs to be established urgently, large-scale geological structure mechanical modeling is not needed, the inhomogeneous distribution condition of the triaxial stress field of the stratum at any position in a target area can be calculated only by means of a small amount of field measurement data, a foundation is provided for predicting the hydraulic fracturing construction pressure range, the hydraulic fracture extension path and the fracture network complexity, the deep shale fracturing optimization design and the field construction are guided, and the fracturing yield-increasing effect is effectively improved.
Disclosure of Invention
The invention provides a method for calculating a heterogeneous stress field of a deep shale complex formation, aiming at solving the problems that the stress field of the deep shale complex formation in the prior art changes rapidly in a small area and a rapid and accurate calculation method of the heterogeneous stress field is lacked.
The technical scheme provided by the invention for solving the technical problems is as follows: a method for calculating the inhomogeneous stress field of deep shale complex stratum includes
Determining a stratum triaxial stress gradient under local coordinates of each well position according to stratum stress data and well position coordinate data of drilled wells in a target area;
determining a stratum stress tensor under the central local coordinate of the reservoir at each well position according to the stratum triaxial stress gradient under the local coordinate of each well position and the vertical depth data of the central part of the reservoir in the target area;
determining the formation stress tensor under the central global coordinate of the reservoir at each well position according to the formation stress tensor under the central local coordinate of the reservoir at each well position and a tensor coordinate system conversion algorithm;
determining a stratum stress tensor under the central global coordinate of the reservoir at any position in a target area and a stratum stress tensor component distribution cloud picture under the central global coordinate in the target area according to the stratum stress tensor under the central global coordinate of the reservoir at each well position and a natural adjacent point three-dimensional interpolation method;
determining a stratum triaxial main stress value and a direction component thereof under the middle local coordinate of the reservoir at any position in the target area according to the stratum stress tensor under the middle global coordinate of the reservoir at any position in the target area;
according to stratum triaxial main stress values and direction components of the stratum triaxial main stress values under the middle part local coordinates of the reservoir at any position in the target region, drawing a stratum maximum horizontal main stress value, a stratum minimum horizontal main stress value, a stratum vertical main stress value distribution cloud chart and maximum and minimum horizontal stress direction vector charts under the local coordinates in the target region;
and determining the non-uniform distribution condition of the triaxial stress field of the stratum at any position in the target region according to the distribution cloud graph of the maximum horizontal main stress value, the minimum horizontal main stress value and the vertical main stress value of the stratum under the local coordinate in the target region and the direction vector graphs of the maximum horizontal stress and the minimum horizontal stress.
The further technical scheme is that the method for determining the triaxial stress gradient of the stratum under the local coordinate of each well position according to the stratum stress data and the well position coordinate data of the drilled well in the target area comprises the following steps:
collecting stratum stress data obtained by measuring drilled wells in a target area, wherein the stratum stress data comprises a maximum horizontal main stress value, a minimum horizontal main stress value and a vertical main stress value;
collecting coordinate data of drilled well positions in a target area and the vertical depth of each well stratum stress measuring point, and drawing a well position distribution diagram;
respectively calculating three-axis stress gradients of the stratum under the local coordinates of each well position by using the following formula, wherein the three-axis stress gradients of the stratum comprise a maximum horizontal main stress gradient, a minimum horizontal main stress gradient and a vertical main stress gradient;
Figure BDA0002960311050000031
in the formula: gi H、Gi h、Gi vThe maximum horizontal main stress gradient, the minimum horizontal main stress gradient and the vertical main stress gradient of the stratum at the No. i well position are Pa/m; sigmaiH、σi h、σi vThe maximum horizontal main stress value, the minimum horizontal main stress value and the vertical main stress value of the stratum at the I-type well position are Pa; diMeasuring the vertical depth m of a formation stress data measuring point of a No. i well; and i is the number of the formation stress value measuring well.
The further technical scheme is that the step of determining the stratum stress tensor of each well position under the central local coordinate of the reservoir according to the triaxial stress gradient of the reservoir under the local coordinate of each well position and the vertical depth data of the central part of the reservoir in the target area comprises the following steps:
collecting vertical depth data of the middle part of a reservoir in a target area;
respectively calculating the stratum triaxial stress values under the middle local coordinates of the reservoir at each well position by using the following formula;
Figure BDA0002960311050000041
in the formula: gi H、Gi h、Gi vThe maximum horizontal main stress gradient, the minimum horizontal main stress gradient and the vertical main stress gradient of the stratum at the No. i well position are Pa/m; sigmai Hmid、σi hmid、σi vmidThe maximum horizontal main stress value and the minimum horizontal main stress value of the stratum under the middle local coordinate of the reservoir at the position of the No. i wellA vertical main stress value, Pa; dmidThe vertical depth m of the middle part of the reservoir in the target area;
assembling stratum triaxial stress values under the central local coordinates of the reservoir at each well position into a tensor form by using the following formula;
Figure BDA0002960311050000042
in the formula: sigmai Hmid、σi hmid、σi vmidThe maximum horizontal main stress value, the minimum horizontal main stress value and the vertical main stress value of the stratum under the middle local coordinate of the reservoir at the position of the No. i well are Pa; sigmai mid|localIs the stratum stress tensor Pa under the central local coordinate of the reservoir at the position of the No. i well.
The further technical scheme is that the step of determining the formation stress tensor under the central global coordinate of the reservoir at each well position according to the formation stress tensor under the central local coordinate of the reservoir at each well position and a tensor coordinate system conversion algorithm comprises the following steps:
collecting the maximum horizontal main stress direction of each well position obtained by the measurement of the drilled well in the target area, and calculating the included angle between the maximum horizontal main stress direction and the x axis of the global coordinate system;
respectively calculating stratum stress tensor components under the central global coordinates of the reservoir at each well position by using a tensor coordinate system conversion algorithm according to the following formula;
Figure BDA0002960311050000051
in the formula: sigmai xx、σi yy、σi xy、σi xz、σi yz、σi zzThe method comprises the following steps of (1) forming a formation stress tensor component Pa under the central global coordinate of a reservoir at the position of a well I; thetaiThe included angle between the maximum horizontal main stress of the No. i well and the x axis of the global coordinate system is degree; sigmai Hmid、σi hmid、σi vmidThe maximum horizontal main stress value, the minimum horizontal main stress value and the vertical main stress value of the stratum under the middle local coordinate of the reservoir at the position of the No. i well are Pa;
assembling stratum stress tensor components under the central global coordinate of the reservoir at each well position into a tensor form by using the following formula;
Figure BDA0002960311050000052
in the formula: sigmai xx、σi yy、σi xy、σi xz、σi yz、σi zzThe method comprises the following steps of (1) forming a formation stress tensor component Pa under the central global coordinate of a reservoir at the position of a well I; pa; sigmai mid|globalAnd the formation stress tensor is Pa under the central global coordinate of the reservoir at the position of the No. i well.
The further technical scheme is that the method for determining the stratum stress tensor of the reservoir at any position in the target area under the global coordinate in the middle of the reservoir at each well position and the component distribution cloud pictures of the stratum stress tensor under the global coordinate in the target area according to the stratum stress tensor under the global coordinate in the middle of the reservoir at each well position and the three-dimensional interpolation method of the natural adjacent points comprises the following steps:
according to the stratum stress tensor component under the central global coordinate of the reservoir at each well position, calculating the stratum stress tensor component under the central global coordinate of the reservoir at any position in the target area by a natural adjacent point three-dimensional interpolation method by using the following formula;
Figure BDA0002960311050000061
in the formula: sigmai mid|globalThe method comprises the steps of (1) forming a formation stress tensor Pa under the global coordinate of the middle part of a reservoir at any position in a target area; sigmaxx、σyy、σxy、σxz、σyz、σzzThe method comprises the following steps of (1) forming a formation stress tensor component Pa under the global coordinate of the middle part of a reservoir at any position in a target area; interp3N is a natural neighbor three-dimensional interpolation operator；
Assembling the formation stress tensor components under the central global coordinate of the reservoir at any position in the target area into a tensor form by using the following formula;
Figure BDA0002960311050000062
in the formula: sigmai mid|globalThe method comprises the steps of (1) forming a formation stress tensor Pa under the global coordinate of the middle part of a reservoir at any position in a target area; sigmaxx、σyy、σxy、σxz、σyz、σzzThe method comprises the following steps of (1) forming a formation stress tensor component Pa under the global coordinate of the middle part of a reservoir at any position in a target area;
and drawing a distribution cloud picture of the formation stress tensor components under the global coordinate in the target area.
The further technical scheme is that the stratum triaxial principal stress value under the central local coordinate of the reservoir at any position in the target area is determined according to the stratum stress tensor under the central global coordinate of the reservoir at any position in the target area, and the directional components of the stratum triaxial principal stress value comprise:
according to the stratum stress tensor under the central global coordinate of the reservoir at any position in the target area, calculating the stratum stress tensor eigenvalue and eigenvector matrix under the central local coordinate of the reservoir at any position in the target area by utilizing the following formula through the matrix eigenvalue and eigenvector operation;
[V,D]＝eig(σmid|global)
in the formula: eig is a matrix eigenvalue and eigenvector operator; v is a stratum stress tensor eigenvalue matrix Pa under the middle local coordinate of the reservoir at any position in the target area; d is a formation stress tensor eigenvector matrix under the middle local coordinate of the reservoir at any position in the target area, and is dimensionless; sigmamid|globalThe method comprises the steps of (1) forming a formation stress tensor Pa under a global coordinate in the middle of a reservoir;
extracting the maximum horizontal main stress value, the minimum horizontal main stress value, the vertical main stress value and each main stress direction component of the stratum under the middle local coordinate of the reservoir at any position in the target area from the stratum stress tensor eigenvalue and eigenvector matrix under the middle local coordinate of the reservoir at any position in the target area by using the following formula;
V＝[σhmid,σHmid,σvmid]
Figure BDA0002960311050000071
in the formula: v is a stratum stress tensor eigenvalue matrix Pa under the middle local coordinate of the reservoir at any position in the target area; d is a formation stress tensor eigenvector matrix under the middle local coordinate of the reservoir at any position in the target area, and is dimensionless; sigmaHmid、σhmid、σvmidThe maximum horizontal main stress value, the minimum horizontal main stress value and the vertical main stress value Pa of the stratum under the middle local coordinate of the reservoir at any position in the target area are obtained; x is the number ofhmid、yhmid、zhmidThe minimum horizontal principal stress direction component of the stratum under the middle local coordinate of the reservoir at any position in the target area is dimensionless; x is the number ofHmid、yHmid、zHmidThe method is characterized in that the method is a formation maximum horizontal principal stress direction component under the middle local coordinate of a reservoir at any position in a target region and is dimensionless; x is the number ofvmid、yvmid、zvmidThe vertical principal stress direction component of the stratum under the middle local coordinate of the reservoir at any position in the target area is dimensionless.
The invention has the following beneficial effects: compared with the prior art, the method is specially used for solving the problem that the acquisition means of stress field data of the deep shale complex formation is insufficient, and a tensor coordinate system conversion algorithm and a natural adjacent point three-dimensional interpolation method are combined, so that a calculation method of the non-uniform stress field of the deep shale complex formation is provided; because the method can calculate the inhomogeneous distribution condition of the triaxial stress field of the stratum at any position in the target area only by depending on limited data such as the triaxial stress value, the maximum horizontal main stress direction and the like of the stratum drilled with the well in the target area without large-scale geological structure mechanical modeling, the method solves the problems that the stratum stress field of the complex structure of the deep shale changes faster in a smaller area and a fast and accurate inhomogeneous stratum stress field calculation method is lacked.
Drawings
FIG. 1 is a block diagram of the computational flow of the method of the present invention;
FIG. 2 is a diagram of the location of a drilled well within a target area of an embodiment;
FIG. 3 is a cloud of the distribution of the components of the formation stress tensor at the global coordinates of the middle of the reservoir in the target region of the embodiment;
FIG. 4 is a cloud chart of distribution of maximum and minimum horizontal principal stress values of the formation at the local coordinates of the middle portion of the reservoir in the target region according to the embodiment;
FIG. 5 is a diagram of the direction vector of the maximum and minimum horizontal principal stresses of the formation at the local coordinates of the middle portion of the reservoir in the target region according to an embodiment.
Detailed Description
The technical solutions of the present invention will be described clearly and completely with reference to the accompanying drawings, and it should be understood that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As shown in FIG. 1, the method for calculating the inhomogeneous stress field of the deep shale complex formation stratum comprises the following steps:
step A, calculating a stratum triaxial stress gradient under local coordinates of each well position according to drilled well data in a target area;
step A, calculating a flow:
collecting formation stress data obtained by measuring drilled wells in a target area, wherein the formation stress data comprises a maximum horizontal main stress value, a minimum horizontal main stress value and a vertical main stress value;
collecting coordinate data of drilled well positions in the target area and the vertical depth of the stratum stress measuring points of each well, and drawing a well position distribution diagram;
thirdly, calculating three-axis principal stress gradients of the stratum at each well position respectively by using an equation (1), wherein the three-axis principal stress gradients of the stratum comprise a maximum horizontal principal stress gradient, a minimum horizontal principal stress gradient and a vertical principal stress gradient;
Figure BDA0002960311050000091
in the formula: gi H、Gi h、Gi vThe maximum horizontal main stress gradient, the minimum horizontal main stress gradient and the vertical main stress gradient of the stratum at the No. i well position are Pa/m; sigmai H、σi h、σi vThe maximum horizontal main stress value, the minimum horizontal main stress value and the vertical main stress value of the stratum at the I-type well position are Pa; diMeasuring the vertical depth m of a formation stress data measuring point of a No. i well; and i is the number of the formation stress value measuring well.
B, calculating the maximum horizontal principal stress value, the minimum horizontal principal stress value and the vertical principal stress value of the stratum under the local coordinate of the middle part of the reservoir at each well position based on the triaxial principal stress gradient of the stratum and the vertical depth data of the middle part of the reservoir in the target area in the step A, and assembling the maximum horizontal principal stress value, the minimum horizontal principal stress value and the vertical principal stress value into a stratum stress tensor form under the local coordinate;
step B, calculating a flow:
collecting vertical depth data of the middle part of a reservoir in a target area;
secondly, calculating the stratum triaxial stress value under the middle local coordinate of the reservoir at each well position by using an equation (2), wherein the equation comprises the following steps: the maximum horizontal main stress value, the minimum horizontal main stress value and the vertical main stress value of the stratum;
thirdly, assembling stratum triaxial stress values under the central local coordinate of the reservoir at the position of the well I into a tensor form by using an equation (3);
Figure BDA0002960311050000092
Figure BDA0002960311050000093
in the formula: sigmaiHmid、σi hmid、σi vmidThe maximum horizontal main stress value, the minimum horizontal main stress value and the vertical main stress value of the stratum under the middle local coordinate of the reservoir at the position of the No. i well are Pa; dmidThe vertical depth m of the middle part of the reservoir in the target area; sigmai mid|localThe method comprises the following steps of (1) forming a formation stress tensor Pa under a central local coordinate of a reservoir at a well position I;
step C, calculating a stratum stress tensor component under the central global coordinate of the reservoir at each well position through a tensor coordinate system conversion algorithm based on the stratum stress tensor under the central local coordinate of the reservoir at each well position obtained in the step B;
step C, calculating a flow:
collecting the maximum horizontal main stress direction of each well position obtained by drilling and setting well measurement in a target area, and calculating the included angle between the maximum horizontal main stress direction and the x axis of the global coordinate system;
respectively calculating the formation stress tensor components under the central global coordinate of the reservoir at each well position by using an equation (5) through a tensor coordinate system conversion algorithm;
thirdly, assembling the formation stress tensor components under the central global coordinate of the reservoir at each well position into a tensor form by using an equation (4);
Figure BDA0002960311050000101
Figure BDA0002960311050000102
in the formula: sigmai xx、σi yy、σi xy、σixz、σi yz、σi zzThe method comprises the following steps of (1) forming a formation stress tensor component Pa under the central global coordinate of a reservoir at the position of a well I; thetaiThe included angle between the maximum horizontal main stress of the No. i well and the x axis of the global coordinate system is degree;
step D, calculating the stratum stress tensor component under the central global coordinate of the reservoir at any position in the target area through a natural adjacent point three-dimensional interpolation method based on the stratum stress tensor under the central global coordinate of the reservoir at each well position obtained in the step C, assembling the stratum stress tensor component into a stratum stress tensor form under the global coordinate, and drawing a stratum stress tensor component distribution cloud picture under the global coordinate in the target area;
step D, calculating a flow:
calculating the stratum stress tensor component under the central global coordinate of the reservoir at any position in a target area by using an equation (7) through a natural adjacent point three-dimensional interpolation method according to the stratum stress tensor component under the central global coordinate of the reservoir at each well position;
assembling the formation stress tensor components under the global coordinate of the middle part of the reservoir at any position in the target area into a tensor form by using an equation (6);
thirdly, drawing a distribution cloud picture of stratum stress tensor components under the global coordinate in the target area;
Figure BDA0002960311050000111
Figure BDA0002960311050000112
in the formula: sigmai mid|globalThe method comprises the steps of (1) forming a formation stress tensor Pa under the global coordinate of the middle part of a reservoir at any position in a target area; sigmaxx、σyy、σxy、σxz、σyz、σzzThe method comprises the following steps of (1) forming a formation stress tensor component Pa under the global coordinate of the middle part of a reservoir at any position in a target area; interp3N is a natural neighboring point three-dimensional interpolation operator;
step E, calculating a stratum triaxial main stress value under the middle local coordinate of the reservoir at any position in the target area and a direction component thereof based on the stratum stress tensor under the middle global coordinate of the reservoir at any position in the target area obtained in the step D;
step E, calculating a flow:
according to the stratum stress tensor under the central global coordinate of the reservoir at any position in the target area, calculating the stratum stress tensor eigenvalue and eigenvector matrix under the central local coordinate of the reservoir at any position in the target area by using an equation (8) through the matrix eigenvalue and eigenvector operation;
extracting the maximum horizontal principal stress value, the minimum horizontal principal stress value, the vertical principal stress value and each principal stress direction component of the stratum under the middle local coordinate of the reservoir at any position in the target area from the stratum stress tensor eigenvalue and the eigenvector matrix under the middle local coordinate of the reservoir at any position in the target area by using equations (9) and (10);
[V,D]＝eig(σmid|global) (8)
V＝[σhmid,σHmid,σvmid] (9)
Figure BDA0002960311050000121
in the formula: eig is a matrix eigenvalue and eigenvector operator; v is a stratum stress tensor eigenvalue matrix Pa under the middle local coordinate of the reservoir at any position in the target area; d is a formation stress tensor eigenvector matrix under the middle local coordinate of the reservoir at any position in the target area, and is dimensionless; sigmaHmid、σhmid、σvmidThe maximum horizontal main stress value, the minimum horizontal main stress value and the vertical main stress value Pa of the stratum under the middle local coordinate of the reservoir at any position in the target area are obtained; x is the number ofhmid、yhmid、zhmidThe minimum horizontal principal stress direction component of the stratum under the middle local coordinate of the reservoir at any position in the target area is dimensionless; x is the number ofHmid、yHmid、zHmidThe method is characterized in that the method is a formation maximum horizontal principal stress direction component under the middle local coordinate of a reservoir at any position in a target region and is dimensionless; x is the number ofvmid、yvmid、zvmidThe vertical principal stress direction component of the stratum under the middle local coordinate of the reservoir at any position in the target area is dimensionless;
f, drawing a stratum maximum horizontal main stress value, a minimum horizontal main stress value, a vertical main stress value distribution cloud chart and maximum and minimum horizontal stress direction vector charts under local coordinates in the target region;
and G, determining the non-uniform distribution condition of the triaxial stress field of the stratum at any position in the target region according to the stratum stress tensor component distribution cloud picture under the global coordinate in the target region, the stratum maximum horizontal main stress value, the minimum horizontal main stress value and the vertical main stress value distribution cloud picture under the local coordinate in the target region, and the maximum and minimum horizontal stress direction vector pictures.
Example 1
A method for calculating a nonuniform stress field of a deep shale complex formation comprises the following steps:
actual data of 6 drilled wells in a deep shale gas reservoir area with a certain complex structure are shown in table 1.
TABLE 1 actual data of 6 drilled wells in a deep shale gas reservoir area of a complex structure
Figure BDA0002960311050000131
The calculation flow diagram of the method according to the invention (as shown in fig. 1) develops an example calculation:
firstly, the triaxial stress gradient of the stratum under the local coordinate of each well position is calculated by using the drilled well data in the target area.
The method comprises the following specific steps: inputting formation stress data obtained by measuring drilled wells in a target area, wherein the formation stress data comprises a maximum horizontal main stress value, a minimum horizontal main stress value and a vertical main stress value; secondly, inputting the coordinate data of the drilled well position in the target area and the vertical depth of the stratum stress measuring point of each well, and drawing a well position distribution diagram as shown in figure 2; thirdly, calculating the three-axis principal stress gradient of the stratum at each well position by using the equation (1), comprising the following steps: maximum horizontal principal stress gradient, minimum horizontal principal stress gradient, vertical principal stress gradient, as shown in table 2.
TABLE 2 three-axis principal stress gradient of formation at each well location
Figure BDA0002960311050000132
Figure BDA0002960311050000141
And then, calculating the triaxial stress value of the stratum under the middle local coordinate of the reservoir at each well position by using the middle vertical depth data of the reservoir in the target area.
The method comprises the following specific steps: inputting vertical depth data of the middle part of a reservoir in a target area: 2200 m; secondly, calculating the stratum triaxial stress value under the middle local coordinate of the reservoir at each well position by using an equation (2), wherein the equation comprises the following steps: the maximum horizontal principal stress value, the minimum horizontal principal stress value, and the vertical principal stress value are shown in table 3; thirdly, assembling the triaxial stress values of the stratum under the local coordinates of the middle part of the reservoir at the position of the i-type well into a tensor form by using an equation (3).
TABLE 3 triaxial stress values of the formation at the middle local coordinates of the reservoir at each well location
Figure BDA0002960311050000142
And then, calculating the stratum stress tensor component under the central global coordinate of the reservoir at each well position by using a tensor coordinate system conversion algorithm.
The method comprises the following specific steps: firstly, inputting the maximum horizontal main stress direction of each well position obtained by drilling and setting well measurement in a target area, namely the included angle between the maximum horizontal main stress direction and the x axis of the global coordinate system; respectively calculating the formation stress tensor components under the central global coordinate of the reservoir at each well position by using an equation (5) through a tensor coordinate system conversion algorithm, as shown in a table 4; and thirdly, assembling the formation stress tensor components under the central global coordinate of the reservoir at each well position into a tensor form by using an equation (4).
TABLE 4 formation stress tensor components at the central global coordinates of the reservoir at each well location
Figure BDA0002960311050000151
And then, calculating the formation stress tensor component under the central global coordinate of the reservoir at any position in the target area by using a natural adjacent point three-dimensional interpolation algorithm.
The method comprises the following specific steps: calculating the stratum stress tensor component under the central global coordinate of the reservoir at any position in a target area by using an equation (7) through a natural adjacent point three-dimensional interpolation method according to the stratum stress tensor component under the central global coordinate of the reservoir at each well position, and drawing a distribution cloud picture; assembling the formation stress tensor components under the global coordinate of the middle part of the reservoir at any position in the target area into a tensor form by using an equation (6); and thirdly, drawing a distribution cloud picture of each component of the formation stress tensor under the global coordinate in the target area, as shown in fig. 3.
And then, calculating the stratum triaxial stress value and the maximum principal stress direction under the middle local coordinate of the reservoir at any position in the target area by using a tensor coordinate system conversion algorithm.
The method comprises the following specific steps: calculating a stratum stress tensor eigenvalue and an eigenvector matrix under the central local coordinate of the reservoir at any position in the target area by using an equation (8) through the matrix eigenvalue and eigenvector operation according to the stratum stress tensor under the central global coordinate of the reservoir at any position in the target area; extracting the maximum horizontal principal stress value, the minimum horizontal principal stress value, the vertical principal stress value and each principal stress direction component of the stratum under the middle local coordinate of the reservoir at any position in the target area from the stratum stress tensor eigenvalue and the eigenvector matrix under the middle local coordinate of the reservoir at any position in the target area by using equations (9) and (10); and thirdly, drawing a distribution cloud picture of the maximum horizontal principal stress value, the minimum horizontal principal stress value and the vertical principal stress value of the stratum under the local coordinate in the target area, and a direction vector picture of the maximum horizontal stress and the minimum horizontal stress respectively as shown in the figures 4 and 5.
The method is specially used for solving the problem of insufficient data acquisition means of the stress field of the deep shale complex formation, and provides a calculation method of the inhomogeneous stress field of the deep shale complex formation by combining a tensor coordinate system conversion algorithm and a natural adjacent point three-dimensional interpolation method. Because the method can calculate the inhomogeneous distribution condition of the triaxial stress field of the stratum at any position in the target area only by depending on limited data such as the triaxial stress value, the maximum horizontal main stress direction and the like of the stratum drilled with the well in the target area without large-scale geological structure mechanical modeling, the method solves the problems that the stratum stress field of the complex structure of the deep shale changes faster in a smaller area and a fast and accurate inhomogeneous stratum stress field calculation method is lacked.
Although the present invention has been described with reference to the above embodiments, it should be understood that the present invention is not limited to the above embodiments, and those skilled in the art can make various changes and modifications without departing from the scope of the present invention.

Claims (6)

1. A method for calculating a nonuniform stress field of a deep shale complex formation is characterized by comprising the following steps:
determining a stratum triaxial stress gradient under local coordinates of each well position according to stratum stress data and well position coordinate data of drilled wells in a target area;
determining a stratum stress tensor under the central local coordinate of the reservoir at each well position according to the stratum triaxial stress gradient under the local coordinate of each well position and the vertical depth data of the central part of the reservoir in the target area;
determining the formation stress tensor under the central global coordinate of the reservoir at each well position according to the formation stress tensor under the central local coordinate of the reservoir at each well position and a tensor coordinate system conversion algorithm;
determining a stratum stress tensor under the central global coordinate of the reservoir at any position in a target area and a stratum stress tensor component distribution cloud picture under the central global coordinate in the target area according to the stratum stress tensor under the central global coordinate of the reservoir at each well position and a natural adjacent point three-dimensional interpolation method;
determining a stratum triaxial main stress value and a direction component thereof under the middle local coordinate of the reservoir at any position in the target area according to the stratum stress tensor under the middle global coordinate of the reservoir at any position in the target area;
according to stratum triaxial main stress values and direction components of the stratum triaxial main stress values under the middle part local coordinates of the reservoir at any position in the target region, drawing a stratum maximum horizontal main stress value, a stratum minimum horizontal main stress value, a stratum vertical main stress value distribution cloud chart and maximum and minimum horizontal stress direction vector charts under the local coordinates in the target region;
and determining the non-uniform distribution condition of the triaxial stress field of the stratum at any position in the target region according to the distribution cloud graph of the maximum horizontal main stress value, the minimum horizontal main stress value and the vertical main stress value of the stratum under the local coordinate in the target region and the direction vector graphs of the maximum horizontal stress and the minimum horizontal stress.
2. The method for calculating the inhomogeneous stress field of the deep shale complex formation stratum according to claim 1, wherein the determining the triaxial stress gradient of the stratum under the local coordinate of each well position according to the stratum stress data and the well position coordinate data of the drilled well in the target area comprises:
collecting stratum stress data obtained by measuring drilled wells in a target area, wherein the stratum stress data comprises a maximum horizontal main stress value, a minimum horizontal main stress value and a vertical main stress value;
collecting coordinate data of drilled well positions in a target area and the vertical depth of each well stratum stress measuring point, and drawing a well position distribution diagram;
respectively calculating three-axis stress gradients of the stratum under the local coordinates of each well position by using the following formula, wherein the three-axis stress gradients of the stratum comprise a maximum horizontal main stress gradient, a minimum horizontal main stress gradient and a vertical main stress gradient;
Figure FDA0002960311040000021
in the formula: gi H、Gi h、Gi vThe maximum horizontal main stress gradient, the minimum horizontal main stress gradient and the vertical main stress gradient of the stratum at the No. i well position are Pa/m; sigmai H、σi h、σi vThe maximum horizontal main stress value, the minimum horizontal main stress value and the vertical main stress value of the stratum at the I-type well position are Pa; diMeasuring the vertical depth m of a formation stress data measuring point of a No. i well; and i is the number of the formation stress value measuring well.
3. The method for calculating the nonuniform stress field of the deep shale complex formation stratum as claimed in claim 2, wherein the determining the formation stress tensor at the central reservoir local coordinate of each well position according to the three-axis stress gradient of the stratum at the local coordinate of each well position and the vertical depth data at the central reservoir in the target area comprises:
collecting vertical depth data of the middle part of a reservoir in a target area;
respectively calculating the stratum triaxial stress values under the middle local coordinates of the reservoir at each well position by using the following formula;
Figure FDA0002960311040000022
in the formula: gi H、Gi h、Gi vThe maximum horizontal main stress gradient, the minimum horizontal main stress gradient and the vertical main stress gradient of the stratum at the No. i well position are Pa/m; sigmai Hmid、σi hmid、σi vmidThe maximum horizontal main stress value, the minimum horizontal main stress value and the vertical main stress value of the stratum under the middle local coordinate of the reservoir at the position of the No. i well are Pa; dmidThe vertical depth m of the middle part of the reservoir in the target area;
assembling stratum triaxial stress values under the central local coordinates of the reservoir at each well position into a tensor form by using the following formula;
Figure FDA0002960311040000031
in the formula: sigmai Hmid、σi hmid、σi vmidThe maximum horizontal main stress value, the minimum horizontal main stress value and the vertical main stress value of the stratum under the middle local coordinate of the reservoir at the position of the No. i well are Pa; sigmai mid|localIs the stratum stress tensor Pa under the central local coordinate of the reservoir at the position of the No. i well.
4. The method for calculating the nonuniform stress field of the deep shale complex formation stratum according to the claim 3, wherein the step of determining the stratum stress tensor under the central global coordinate of the reservoir at each well position according to the stratum stress tensor under the central local coordinate of the reservoir at each well position and a tensor coordinate system conversion algorithm comprises the following steps:
collecting the maximum horizontal main stress direction of each well position obtained by the measurement of the drilled well in the target area, and calculating the included angle between the maximum horizontal main stress direction and the x axis of the global coordinate system;
respectively calculating stratum stress tensor components under the central global coordinates of the reservoir at each well position by using a tensor coordinate system conversion algorithm according to the following formula;
Figure FDA0002960311040000032
in the formula: sigmai xx、σi yy、σi xy、σi xz、σi yz、σi zzThe method comprises the following steps of (1) forming a formation stress tensor component Pa under the central global coordinate of a reservoir at the position of a well I; thetaiThe included angle between the maximum horizontal main stress of the No. i well and the x axis of the global coordinate system is degree; sigmai Hmid、σi hmid、σi vmidFor the local coordinate of the middle part of the reservoir at the position of the No. i wellThe maximum horizontal main stress value, the minimum horizontal main stress value, the vertical main stress value and Pa of the lower stratum;
assembling stratum stress tensor components under the central global coordinate of the reservoir at each well position into a tensor form by using the following formula;
Figure FDA0002960311040000041
in the formula: sigmai xx、σi yy、σi xy、σi xz、σi yz、σi zzThe method comprises the following steps of (1) forming a formation stress tensor component Pa under the central global coordinate of a reservoir at the position of a well I; pa; sigmai mid|globalAnd the formation stress tensor is Pa under the central global coordinate of the reservoir at the position of the No. i well.
5. The method for calculating the inhomogeneous stress field of the deep shale complex formation stratum according to claim 4, wherein the step of determining the formation stress tensor under the central global coordinate of the reservoir at any position in the target area and the component distribution cloud map of the formation stress tensor under the central global coordinate of the reservoir in the target area according to the formation stress tensor under the central global coordinate of the reservoir at each well position and a natural adjacent point three-dimensional interpolation method comprises the following steps:
according to the stratum stress tensor component under the central global coordinate of the reservoir at each well position, calculating the stratum stress tensor component under the central global coordinate of the reservoir at any position in the target area by a natural adjacent point three-dimensional interpolation method by using the following formula;
Figure FDA0002960311040000042
in the formula: sigmai mid|globalThe method comprises the steps of (1) forming a formation stress tensor Pa under the global coordinate of the middle part of a reservoir at any position in a target area; sigmaxx、σyy、σxy、σxz、σyz、σzzThe method comprises the following steps of (1) forming a formation stress tensor component Pa under the global coordinate of the middle part of a reservoir at any position in a target area; interp3N is a natural neighboring point three-dimensional interpolation operator;
assembling the formation stress tensor components under the central global coordinate of the reservoir at any position in the target area into a tensor form by using the following formula;
Figure FDA0002960311040000043
in the formula: sigmai mid|globalThe method comprises the steps of (1) forming a formation stress tensor Pa under the global coordinate of the middle part of a reservoir at any position in a target area; sigmaxx、σyy、σxy、σxz、σyz、σzzThe method comprises the following steps of (1) forming a formation stress tensor component Pa under the global coordinate of the middle part of a reservoir at any position in a target area;
and drawing a distribution cloud picture of the formation stress tensor components under the global coordinate in the target area.
6. The method for calculating the inhomogeneous stress field of the deep shale complex formation stratum according to claim 5, wherein the determining of the triaxial principal stress value of the stratum under the central local coordinate of the reservoir at any position in the target region according to the stratum stress tensor under the central global coordinate of the reservoir at any position in the target region and the directional component thereof comprises:
according to the stratum stress tensor under the central global coordinate of the reservoir at any position in the target area, calculating the stratum stress tensor eigenvalue and eigenvector matrix under the central local coordinate of the reservoir at any position in the target area by utilizing the following formula through the matrix eigenvalue and eigenvector operation;
[V,D]＝eig(σmidglobal)
in the formula: eig is a matrix eigenvalue and eigenvector operator; v is a stratum stress tensor eigenvalue matrix Pa under the middle local coordinate of the reservoir at any position in the target area; d is a formation stress tensor eigenvector matrix under the middle local coordinate of the reservoir at any position in the target area, and is dimensionless;σmid|globalthe method comprises the steps of (1) forming a formation stress tensor Pa under a global coordinate in the middle of a reservoir;
extracting the maximum horizontal main stress value, the minimum horizontal main stress value, the vertical main stress value and each main stress direction component of the stratum under the middle local coordinate of the reservoir at any position in the target area from the stratum stress tensor eigenvalue and eigenvector matrix under the middle local coordinate of the reservoir at any position in the target area by using the following formula;
V＝[σhmid,σHmid,σvmid]
Figure FDA0002960311040000051
in the formula: v is a stratum stress tensor eigenvalue matrix Pa under the middle local coordinate of the reservoir at any position in the target area; d is a formation stress tensor eigenvector matrix under the middle local coordinate of the reservoir at any position in the target area, and is dimensionless; sigmaHmid、σhmid、σvmidThe maximum horizontal main stress value, the minimum horizontal main stress value and the vertical main stress value Pa of the stratum under the middle local coordinate of the reservoir at any position in the target area are obtained; x is the number ofhmid、yhmid、zhmidThe minimum horizontal principal stress direction component of the stratum under the middle local coordinate of the reservoir at any position in the target area is dimensionless; x is the number ofHmid、yHmid、zHmidThe method is characterized in that the method is a formation maximum horizontal principal stress direction component under the middle local coordinate of a reservoir at any position in a target region and is dimensionless; x is the number ofvmid、yvmid、zvmidThe vertical principal stress direction component of the stratum under the middle local coordinate of the reservoir at any position in the target area is dimensionless.
CN202110236279.3A 2021-03-03 2021-03-03 Method for calculating inhomogeneous stress field of deep shale complex formation Active CN112943198B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110236279.3A CN112943198B (en) 2021-03-03 2021-03-03 Method for calculating inhomogeneous stress field of deep shale complex formation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110236279.3A CN112943198B (en) 2021-03-03 2021-03-03 Method for calculating inhomogeneous stress field of deep shale complex formation

Publications (2)

Publication Number Publication Date
CN112943198A true CN112943198A (en) 2021-06-11
CN112943198B CN112943198B (en) 2022-03-01

Family

ID=76247455

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110236279.3A Active CN112943198B (en) 2021-03-03 2021-03-03 Method for calculating inhomogeneous stress field of deep shale complex formation

Country Status (1)

Country Link
CN (1) CN112943198B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114734640A (en) * 2022-04-22 2022-07-12 上海大学 Printing path generation method, computer equipment and storage medium

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110166843A1 (en) * 2007-08-24 2011-07-07 Sheng-Yuan Hsu Method For Modeling Deformation In Subsurface Strata
US20130138410A1 (en) * 2008-11-12 2013-05-30 Geoscape Analytics, Inc. Methods and systems for constructing and using a subterranean geomechanics model spanning local to zonal scale in complex geological environments
US20160018542A1 (en) * 2014-07-15 2016-01-21 Schlumberger Technology Corporation Method to invert for fault activity and tectonic stress
US20180172857A1 (en) * 2015-05-20 2018-06-21 Schlumberger Technology Corporation Inversion for Tectonic Stress
CN110390175A (en) * 2019-07-30 2019-10-29 西安幔源油气勘探开发研究有限公司 Curtain source Accumulation zone analyzing and predicting method based on crustal stress finite element
CN110705168A (en) * 2019-10-11 2020-01-17 内蒙古科技大学 Simulation method of structural stress field

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110166843A1 (en) * 2007-08-24 2011-07-07 Sheng-Yuan Hsu Method For Modeling Deformation In Subsurface Strata
US20130138410A1 (en) * 2008-11-12 2013-05-30 Geoscape Analytics, Inc. Methods and systems for constructing and using a subterranean geomechanics model spanning local to zonal scale in complex geological environments
US20160018542A1 (en) * 2014-07-15 2016-01-21 Schlumberger Technology Corporation Method to invert for fault activity and tectonic stress
US20180172857A1 (en) * 2015-05-20 2018-06-21 Schlumberger Technology Corporation Inversion for Tectonic Stress
CN110390175A (en) * 2019-07-30 2019-10-29 西安幔源油气勘探开发研究有限公司 Curtain source Accumulation zone analyzing and predicting method based on crustal stress finite element
CN110705168A (en) * 2019-10-11 2020-01-17 内蒙古科技大学 Simulation method of structural stress field

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
GUERRA,C.等: "Stress prediction using 1D and 3D geomechanical models of a tight gas reservoir—A case study from the Lower Magdalena Valley Basin, Colombia", 《GEOMECHANICS FOR ENERGY AND THE ENVIRONMENT》 *
朱传华等: "非均质储层三维构造应力场模拟方法", 《吉林大学学报( 地球科学版)》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114734640A (en) * 2022-04-22 2022-07-12 上海大学 Printing path generation method, computer equipment and storage medium
CN114734640B (en) * 2022-04-22 2024-01-16 上海大学 Printing path generation method, computer equipment and storage medium

Also Published As

Publication number Publication date
CN112943198B (en) 2022-03-01

Similar Documents

Publication Publication Date Title
CN109241588B (en) Simulation method for single crack propagation based on quasi-continuous geomechanical model
US8931580B2 (en) Method for using dynamic target region for well path/drill center optimization
CN113821953B (en) Temporary plugging repeated fracturing design method based on dynamic ground stress and residual oil characteristics
WO2016192077A1 (en) Method for establishing and solving numerical well-testing model of horizontal well for tight gas hydraulic fracturing
CN113919196A (en) Reservoir three-dimensional stress field simulation method, simulation system, terminal and storage medium
CN106484958B (en) A kind of system of the three-dimensional cracking permeability tensor computation model based on pit shaft gap observation
CN106014399B (en) Method for establishing high-precision three-dimensional ground stress model of heterogeneous stratum
CN111456709B (en) Horizontal well multistage fracturing segmentation clustering method based on logging curve
CN110119994A (en) A kind of GIS supports the quick-fried heap displacement extraction of lower metallic ore and prediction technique
CN111695228A (en) Multi-scale fracture modeling method for fracture-cave carbonate reservoir
CN112943198B (en) Method for calculating inhomogeneous stress field of deep shale complex formation
CN113821956B (en) Evaluation method for disturbance quantity of current geostress structure of deep shale reservoir
CN113779811B (en) Quantitative analysis method for structural deformation to current ground stress orientation disturbance
CN111950112A (en) Dynamic analysis method for carbonate reservoir suitable for bottom sealing
CHEN et al. 3D geological modeling of dual porosity carbonate reservoirs: A case from the Kenkiyak pre-salt oilfield, Kazakhstan
CN113356843A (en) Borehole wall stability analysis method, device, medium and equipment for stratum
CN111768104B (en) Quantitative evaluation method for lithological change degree of monoclinic sand shale interbedded dam foundation
CN115964962A (en) Three-dimensional seam network modeling digital-analog integration method based on geological engineering integration idea
CN112696181A (en) Fracturing prediction method and device based on coupling relation between ground stress and natural fractures
CN112558179A (en) Intelligent inversion method of far field stress state
CN111897029A (en) Method for determining spatial distribution of underground anticline cracks through core-logging interactive comparison
CN113236345B (en) Design method of drilling fracture visualization system
CN115828636B (en) Anti-channeling construction parameter optimization method for shale gas well group fracturing
CN109541683B (en) Construction method and system of underground point dam inclined stratum model grid
CN117852416B (en) Multimode grouting precontrolled analysis method and system based on digital geological model

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230411

Address after: 610500, Xindu Avenue, Xindu District, Sichuan, Chengdu, 8

Patentee after: SOUTHWEST PETROLEUM University

Patentee after: CHINA PETROLEUM & CHEMICAL Corp.

Address before: 610500, Xindu Avenue, Xindu District, Sichuan, Chengdu, 8

Patentee before: SOUTHWEST PETROLEUM University