CN114200538B - Ground stress direction determining method, device, storage medium and computer equipment - Google Patents
Ground stress direction determining method, device, storage medium and computer equipment Download PDFInfo
- Publication number
- CN114200538B CN114200538B CN202010980973.1A CN202010980973A CN114200538B CN 114200538 B CN114200538 B CN 114200538B CN 202010980973 A CN202010980973 A CN 202010980973A CN 114200538 B CN114200538 B CN 114200538B
- Authority
- CN
- China
- Prior art keywords
- core sample
- determining
- target depth
- ground stress
- sampling well
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 61
- 238000003860 storage Methods 0.000 title claims abstract description 17
- 238000005070 sampling Methods 0.000 claims abstract description 66
- 238000003384 imaging method Methods 0.000 claims description 20
- 238000005096 rolling process Methods 0.000 claims description 15
- 230000003068 static effect Effects 0.000 claims description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 5
- 230000001902 propagating effect Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 abstract description 6
- 230000009286 beneficial effect Effects 0.000 abstract description 4
- 238000011161 development Methods 0.000 description 7
- 238000009826 distribution Methods 0.000 description 6
- 239000011435 rock Substances 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 238000004088 simulation Methods 0.000 description 4
- 239000003208 petroleum Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 230000002776 aggregation Effects 0.000 description 2
- 238000004220 aggregation Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 239000003550 marker Substances 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 238000010408 sweeping Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000005465 channeling Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011157 data evaluation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000001965 increasing effect Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000005405 multipole Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000012496 stress study Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- 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
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Abstract
The invention provides a ground stress direction determining method, a ground stress direction determining device, a storage medium and computer equipment. The method comprises the following steps: acquiring a core sample of a stratum at a target depth, measuring longitudinal wave speeds of longitudinal waves in a plurality of different diameter directions when the longitudinal waves propagate in the core sample, and determining the direction of ground stress in the core sample according to the longitudinal wave speeds of the different diameter directions; comparing the grain characteristics of the side surface of the core sample with the grain characteristics of the sampling well wall of the core sample at the target depth, and determining the relative orientation of the core sample and the sampling well wall of the core sample at the target depth according to the comparison result; and determining the direction of the ground stress in the stratum at the target depth according to the direction of the ground stress in the core sample and the relative direction of the core sample and the wall of the sampling well at the target depth. The method has strong theoretical basis, can be used for judging the ground stress direction of various reservoirs, particularly the stratum with small plastic change, is beneficial to optimizing the well path and reasonably arranging the hydraulic fracturing direction.
Description
Technical Field
The invention relates to the technical field of petroleum geology and petroleum engineering, in particular to a method and a device for determining a ground stress direction, a storage medium and computer equipment.
Background
In oil exploration and development, ground stress plays an increasingly important role. The ground stress is a key factor influencing the fracturing modification effect. In the reconstruction of oil field reservoirs, the stress field state determines the form, azimuth, height and width of a hydraulic fracturing fracture, and has an influence on the fracturing yield increasing effect; the well pattern arrangement and adjustment in the development need to study the stress field, otherwise, the oil displacement efficiency is possibly reduced, the anhydrous oil extraction period is shortened, and phenomena such as water channeling, flooding and the like appear. Thus, developing earth stress studies for formation characteristics is an essential element in reservoir recovery.
In a certain sense, the oil and gas migration and aggregation in the oil and gas exploration and development process, the stability of a well wall in the drilling process, the design of a horizontal well, the transformation of an oil layer, the arrangement of a well pattern in water injection development and the like are closely related to the ground stress, and how to effectively and accurately predict the ground stress field of a reservoir has important significance for the oil and gas exploration and development. At present, expert students have widely recognized the role of ground stress in oil and gas exploration and development, and how to effectively conduct research and application of ground stress has become a research hotspot in the geophysical field.
Fairhu-rst proposes a hydraulic fracture crustal stress measurement method, which is the most common method for the current crustal deep crustal stress measurement application; mayer et al propose a stress recovery method, namely, grooving the borehole wall, putting a pressure pillow into the groove, and recovering the stress state of the measuring points at two sides of the groove to the original state by pressurizing, namely recovering the distance between the measuring points to the distance before grooving, wherein the pressure value can be regarded as the original stress value. Li Hong et al experimentally obtained the Kaiser points for each azimuth of the subsurface rock, resulting in the original stress state of the subsurface rock. In the 70 s of the 20 th century, the company stonelen Bei Xiece began to study the method of solving the problems related to formation mechanics (ground stress is an intermediate process) by using logging data, and applied the method to explain the problems of formation collapse pressure, fracture pressure, oil layer sand production and the like in the petroleum exploration and development process. Zhao Liangxiao et al, which discloses that after the width and depth of the well wall breakout are obtained by using a double well diameter curve and a well logging imaging diagram, the maximum and minimum horizontal main stresses can be obtained according to the mechanical properties of the rock. Liu Zhi and the like determine a well wall induced fracturing balance point by using imaging logging data, so that a formula is established to obtain the ground stress. Fu Haicheng and the like study the direction of the ground stress in the ancient areas by using electric imaging data; xue Ru the direction of the ground stress in the ancient areas of the wheel was studied using the electrical imaging data. Wang Xiaojie and the like extract the speed and the azimuth of the fast and slow transverse waves from the waveform data of the multipole acoustic logging, and then determine the azimuth of the maximum stress according to the azimuth of the fast transverse waves. The research shows that the rock mass complete coefficient can be calculated by utilizing acoustic logging data, then the reduction coefficient is calculated, and finally the dynamic mechanical parameter is converted into the static mechanical parameter, so that the continuous ground stress can be calculated based on the static mechanical parameter. Ge Hongkui et al propose two sets of empirical formulas for ground stress calculations for hydraulic fracturing vertical and horizontal joints. Zhao Yongjiang propose a method of identifying stress type elliptical wellbores. Beaumont et al set up a set of relatively complete numerical simulation system by analyzing the geologic formation cause, selecting applicable boundary conditions, and simulating the stress field distribution law. Chen Shu equally simulate the stress distribution rule of each period by analyzing the construction motion evolution process of a certain area and utilizing a finite element method, and summarize the relation between the stress field and the oil gas distribution. Yang Ke et al propose a strain function method for the problems of the stress function method, and the method is used for analyzing the distribution situation of the ground stress around the underground cavern group. The influence of ground stress on migration and aggregation of oil and gas is discussed quantitatively based on finite element law. PaRsons analyze boundary constraints of the target zone based on GPS data, and analyze the structural stress field characteristics of the zone by building a finite element model. Zhang Anzhi establishing a stratum rock physical model through priori information, obtaining an underground medium stiffness matrix based on the model, and calculating underground medium stress distribution according to the stiffness matrix; zong Zhaoyun proposes a method of predicting ground stress using fractured rock parameters; ma Ni considers the influence of HTI medium vertical cracks and the effect of VTI medium horizontal layer at the same time, derives the ground stress expression of the Orthotropic (OA) medium based on orthotropic medium theory, and realizes the earthquake prediction of ground stress by using azimuth pre-stack earthquake data.
The prior results of comprehensive analysis are mainly concentrated on four methods aiming at the evaluation of the ground stress direction. 1. The direct measurement method has the advantages that the ground stress obtained by using an instrument is accurate, but the data size is too small, the cost is very high, and the continuity is lacking, so that the method is not suitable for obtaining the ground stress data in a large range; 2. the acoustic-electric imaging logging data evaluation method has the difficulty that for the stratum with small stress intensity difference, the characteristics of reactive stress deformation, such as induced seams and the like, are difficult to obtain from an acoustic-electric image; 3. the numerical simulation method has the advantages that the numerical simulation ground stress distribution is based on measured data, and meanwhile, the characteristics of topography and geology are considered, a model is built for inversion back calculation and simulation, so that the difficulty is that the geological condition of a region is difficult to be cleaned, and single well data is needed to be used as constraint; 4. the earthquake prediction method has the difficulty that the influence of parameters such as structural stress, pore pressure and the like on the ground stress is required to be considered. The related data of the stratum stress is the scientific basis for horizontal well design and drilling, and the related data of the stratum stress is used for evaluating the direction of the stratum stress, and a long path is needed.
Disclosure of Invention
The invention mainly aims to provide a method, a device, a storage medium and computer equipment for determining the direction of ground stress so as to accurately determine the direction of ground stress in a stratum.
In a first aspect, the present application provides a method for determining a direction of ground stress, including the steps of: acquiring a core sample of a stratum at a target depth, measuring longitudinal wave speeds of longitudinal waves in a plurality of different diameter directions when the longitudinal waves propagate in the core sample, and determining the direction of ground stress in the core sample according to the longitudinal wave speeds of the different diameter directions; comparing the grain characteristics of the side surface of the core sample with the grain characteristics of the sampling well wall of the core sample at the target depth, and determining the relative orientation of the core sample and the sampling well wall of the core sample at the target depth according to the comparison result; and determining the direction of the ground stress in the stratum at the target depth according to the direction of the ground stress in the core sample and the relative direction of the core sample and the wall of the sampling well at the target depth.
In one embodiment, the core sample is cylindrical, and after obtaining the core sample of the formation at the target depth and before measuring the longitudinal wave velocity along a plurality of different diametric directions as the longitudinal wave propagates in the core sample, the method further comprises: setting a reference line parallel to the axial direction of the core sample on the side surface of the core sample, wherein the axial direction of the core sample is a straight line direction passing through the centers of any two circular planes of the core sample; after determining the direction of the earth stress in the core sample from the longitudinal wave velocities of the respective different diametric directions, the method further includes: determining a central angle of the direction of the ground stress in the core sample relative to the reference line, and describing the direction of the ground stress in the core sample by using the central angle.
In one embodiment, the method for determining the relative orientation of the core sample and the sampling well wall thereof at the target depth according to the comparison result comprises the following steps: acquiring a rolling scanning image of a core sample and an electric imaging dynamic and static image of a sampling well wall of the core sample at a target depth; and comparing the rolling scanning image with the electric imaging dynamic-static image, determining the posture of the core sample at the target depth in the sampling well when the grain characteristics of the side surface of the core sample in the rolling scanning image are consistent with the grain characteristics of the sampling well wall of the core sample at the target depth in the electric imaging dynamic-static image, and determining the relative direction of the core sample and the sampling well wall at the target depth according to the posture of the core sample at the target depth in the sampling well.
In one embodiment, determining the relative orientation of the core sample to its sampling well wall at the target depth from the attitude of the core sample at the target depth within its sampling well comprises: determining a connecting line of the projection of a reference line of the core sample in a circular plane perpendicular to the reference line and the center of the circular plane, determining an included angle between the connecting line and the north direction under the posture of the core sample at the target depth in a sampling well, and describing the relative orientation of the core sample and the sampling well wall at the target depth by using the included angle.
In one embodiment, determining the direction of the earth stress in the formation at the target depth based on the direction of the earth stress in the core sample and the relative orientation of the core sample and its sampling well wall at the target depth comprises: and determining the azimuth angle of the direction of the ground stress in the stratum at the target depth relative to the north direction according to the central angle of the direction of the ground stress in the core sample relative to the reference line and the included angle.
In one embodiment, the core sample is a full diameter core sample.
In one embodiment, determining the direction of the earth stress in the core sample from the longitudinal wave velocities of the respective different diametric directions includes: taking the diameter direction of the maximum longitudinal wave speed during propagation as the direction of the ground stress in the core sample; or the direction which is in the same circular plane with the diameter direction of the minimum longitudinal wave speed when the minimum longitudinal wave speed is transmitted and is perpendicular to the diameter direction of the minimum longitudinal wave speed when the minimum longitudinal wave speed is transmitted is taken as the direction of the ground stress in the core sample.
In a second aspect, the present application provides a ground stress direction determining apparatus comprising: the sample processing module is used for acquiring a core sample of the stratum at the target depth, measuring longitudinal wave speeds along a plurality of different diameter directions when longitudinal waves propagate in the core sample, and determining the direction of the ground stress in the core sample according to the longitudinal wave speeds in the different diameter directions; the azimuth determining module is used for comparing the grain characteristics of the side face of the core sample with the grain characteristics of the sampling well wall of the core sample at the target depth, and determining the relative azimuth of the core sample and the sampling well wall of the core sample at the target depth according to the comparison result; and the ground stress determining module is used for determining the ground stress direction in the stratum at the target depth according to the direction of the ground stress in the core sample and the relative direction of the core sample and the sampling well wall of the core sample at the target depth.
In a third aspect, the present application provides a storage medium storing a computer program, characterized in that the computer program, when executed by a processor, implements the steps of the method for determining a direction of ground stress as described above.
In a fourth aspect, the present application provides a computer device comprising a processor and a storage medium storing program code which, when executed by the processor, implements the steps of a method for determining a direction of ground stress as described above.
According to the application, the wave velocity anisotropy and the electric imaging logging data are combined to evaluate the ground stress direction in the stratum, so that the problem that an electric imaging image cannot acquire accurate crack positions and a short plate with larger error in the ground stress direction is determined by a paleomagnetic method is avoided, the ground stress direction can be accurately determined, and a reliable single well basis is provided for well position deployment and engineering transformation and regional stress prediction by utilizing seismic data. The technical scheme of the application uses logging data to determine the direction of the earth stress, and the method has strong theoretical basis, simple operation and higher precision, can be used for judging the direction of the earth stress of various reservoirs, particularly the stratum with little plastic change, is beneficial to optimizing the track of the well body, reasonably arranges the hydraulic fracturing direction and has higher practicability.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this specification, illustrate embodiments of the application and together with the description serve to explain the application and do not constitute a undue limitation on the application, wherein:
FIG. 1 is a flowchart of a method of determining a direction of ground stress according to an exemplary embodiment of the present application;
FIG. 2 is a flow chart of a method for determining a direction of ground stress according to an embodiment of the application;
FIG. 3A is a schematic illustration of a core sample with reference lines according to an embodiment of the present application;
FIG. 3B is a rolling scan image of a full diameter core sample with reference lines according to an embodiment of the present application;
FIG. 4 is a graph showing the variation of the velocity of longitudinal waves of different diameters in a full diameter core sample with respect to the average velocity of the longitudinal waves according to an embodiment of the present application;
FIG. 5 is a contrast plot of an electrically imaged dynamic and static image versus a rolling scan image of a full diameter core sample, according to an embodiment of the present application;
fig. 6 is a graph of the result of determining the direction of ground stress according to an embodiment of the present application.
Detailed Description
It should be noted that, without conflict, the embodiments of the present application and features of the embodiments may be combined with each other. The application will be described in detail below with reference to the drawings in connection with embodiments.
Example 1
Fig. 1 is a flowchart of a method of determining a direction of ground stress according to an exemplary embodiment of the present application. As shown in fig. 1, the present embodiment provides a method for determining a direction of ground stress, which includes the following steps:
S100: acquiring a core sample of a stratum at a target depth, measuring longitudinal wave speeds of longitudinal waves in a plurality of different diameter directions when the longitudinal waves propagate in the core sample, and determining the direction of ground stress in the core sample according to the longitudinal wave speeds of the different diameter directions;
S200: comparing the grain characteristics of the side surface of the core sample with the grain characteristics of the sampling well wall of the core sample at the target depth, and determining the relative orientation of the core sample and the sampling well wall of the core sample at the target depth according to the comparison result;
s300: and determining the direction of the ground stress in the stratum at the target depth according to the direction of the ground stress in the core sample and the relative direction of the core sample and the wall of the sampling well at the target depth.
Example two
The embodiment provides a method for determining a ground stress direction, which comprises the following steps:
The first step: and obtaining a core sample of the stratum at the target depth, and setting a reference line parallel to the axial direction of the core sample at the side surface of the core sample, wherein the core sample is cylindrical and can be a full-diameter core sample.
And a second step of: and measuring longitudinal wave speeds along a plurality of different diameter directions when longitudinal waves propagate in the core sample, determining the direction of the ground stress in the core sample according to the longitudinal wave speeds in the different diameter directions, and specifically determining the central angle of the direction of the ground stress in the core sample relative to the reference line, wherein the direction of the ground stress in the core sample is described by the central angle.
The determining the direction of the ground stress in the core sample according to the longitudinal wave speeds of different diameter directions may include: the diameter direction of the maximum longitudinal wave speed during propagation can be used as the direction of the ground stress in the core sample; or a direction which is in the same circular plane as the diameter direction of the minimum longitudinal wave speed when propagating and is perpendicular to the diameter direction of the minimum longitudinal wave speed when propagating can be used as the direction of the ground stress in the core sample.
And a third step of: comparing the grain characteristics of the side face of the core sample with the grain characteristics of the sampling well wall of the core sample at the target depth, and determining the relative orientation of the core sample and the sampling well wall of the core sample at the target depth according to the comparison result, wherein the method specifically comprises the following steps:
Acquiring a rolling scanning image of a core sample and an electric imaging dynamic and static image of a sampling well wall of the core sample at a target depth; and comparing the rolling scanning image with the electric imaging dynamic-static image, determining the posture of the core sample at the target depth in the sampling well when the grain characteristics of the side surface of the core sample in the rolling scanning image are consistent with the grain characteristics of the sampling well wall of the core sample at the target depth in the electric imaging dynamic-static image, and determining the relative direction of the core sample and the sampling well wall at the target depth according to the posture of the core sample at the target depth in the sampling well.
The determining the relative orientation of the core sample and the wall of the sampling well at the target depth according to the posture of the core sample at the target depth in the sampling well can comprise: determining a connecting line of the projection of a reference line of the core sample in a circular plane perpendicular to the reference line and the center of the circular plane, determining an included angle between the connecting line and the north direction under the posture of the core sample at the target depth in a sampling well, and describing the relative orientation of the core sample and the sampling well wall at the target depth by using the included angle.
Fourth step: determining the direction of the ground stress in the stratum at the target depth according to the direction of the ground stress in the core sample and the relative orientation of the core sample and the wall of the sampling well at the target depth can comprise: and determining the azimuth angle of the direction of the ground stress in the stratum at the target depth relative to the north direction according to the central angle of the direction of the ground stress in the core sample relative to the reference line and the included angle.
According to the application, the wave velocity anisotropy and the electric imaging logging data are combined to evaluate the ground stress direction in the stratum, so that the problem that an electric imaging image cannot acquire accurate crack positions and a short plate with larger error in the ground stress direction determined by a paleomagnetic method is avoided, the ground stress direction can be determined more accurately, and a reliable basis is provided for well position deployment and engineering transformation and regional stress prediction by utilizing seismic data. The technical scheme of the application uses logging data to determine the direction of the earth stress, and the method has strong theoretical basis, simple operation and higher precision, can be used for judging the direction of the earth stress of various reservoirs, particularly the stratum with little plastic change, is beneficial to optimizing the track of the well body, reasonably arranges the hydraulic fracturing direction and has higher practicability.
Example III
Fig. 2 is a flowchart of a method for determining a direction of ground stress according to an embodiment of the present application. As shown in fig. 2, the ground stress direction determining method includes the steps of:
S1: and acquiring a full-diameter core sample of the stratum at the target depth, wherein the full-diameter core sample is cylindrical, and a reference line (shown in fig. 3A and 3B) parallel to the axial direction of the full-diameter core sample is arranged on the side surface of the full-diameter core sample. The numerical marks in fig. 3B are marks made by a logging staff after the core sample is taken out from the ground, the numbers in front of the score lines indicate the number of times, the denominator indicates how many core samples are in total in the number of times, and the numerator indicates how many blocks are the current core sample. The core sample in fig. 3B shows this as a 4 th barrel coring with a total of 40 cores, 7 th-9 th cores in order from top to bottom.
Specifically, the core sample with the reference line drawn can be used for measuring the time for propagation of the longitudinal wave along the diameter direction at the reference line by using an ultrasonic instrument, then the time for propagation of the longitudinal wave along the diameter direction is measured every 10 degrees of clockwise rotation on the circumference of the core sample, and the longitudinal wave speed is calculated by combining the diameter size of the core sample.
S2: and measuring longitudinal wave speeds along a plurality of different diameter directions when longitudinal waves propagate in the full-diameter core sample, determining the direction of the ground stress in the full-diameter core sample according to the longitudinal wave speeds in the different diameter directions, and specifically determining the central angle of the direction of the ground stress in the full-diameter core sample relative to the reference line, wherein the central angle is used for describing the direction of the ground stress in the full-diameter core sample.
The determining the direction of the ground stress in the full-diameter core sample according to the longitudinal wave speeds of different diameter directions may include: taking the diameter direction of the maximum longitudinal wave speed during propagation as the direction of the ground stress in the core sample; or the direction which is in the same circular plane with the diameter direction of the minimum longitudinal wave speed when the minimum longitudinal wave speed is transmitted and is perpendicular to the diameter direction of the minimum longitudinal wave speed when the minimum longitudinal wave speed is transmitted is taken as the direction of the ground stress in the core sample.
Specifically, in experiments, the directions of the maximum longitudinal wave velocity and the minimum longitudinal wave velocity can be determined by the wave velocity deviation. The wave velocity deviation for each diametric direction of the full diameter core sample can be determined using the following equation:
wherein θ represents a central angle between the diameter direction of the full-diameter core sample and the marker line, S θ represents a wave velocity deviation, V θ represents a longitudinal wave velocity in the diameter direction at a central angle θ with the marker line, The average value of longitudinal wave velocities in each diameter direction in the full-diameter core sample is shown.
And (3) making a change curve of the wave speed deviation about the central angle theta, as shown in fig. 4, finding out an angle corresponding to the maximum wave speed or the minimum wave speed from the curve, wherein the central angle corresponding to the maximum wave speed is the central angle of the maximum horizontal main stress direction obtained through experiments relative to the reference line, and in fig. 4, the central angle of the maximum horizontal main stress direction relative to the reference line is 120 degrees.
S3: and rolling and sweeping the full-diameter core sample to obtain a rolling and sweeping image of the full-diameter core sample.
S4: and processing the electric imaging data to obtain dynamic and static images of the sampling well wall of the full-diameter core sample, wherein the electric imaging data can be a logging resistivity curve of the sampling well of the full-diameter core sample.
S5: and comparing the rolling sweep image with the dynamic and static images, and when the grain characteristics of the side surface of the core sample in the rolling sweep image are consistent with the grain characteristics of the sampling well wall of the core sample in the dynamic and static images at the target depth, determining the posture of the core sample at the target depth in the sampling well, and determining the relative direction of the core sample and the sampling well wall at the target depth according to the posture of the core sample at the target depth in the sampling well, as shown in fig. 5.
The determining the relative orientation of the core sample and the wall of the sampling well at the target depth according to the posture of the core sample at the target depth in the sampling well specifically may include: determining a connecting line of the projection of a reference line of the core sample in a circular plane perpendicular to the reference line and the center of the circular plane, determining an included angle between the connecting line and the north direction under the posture of the core sample at the target depth in a sampling well, and describing the relative orientation of the core sample and the sampling well wall at the target depth by using the included angle. As can be seen from fig. 5, the angle is 288.2 °.
S6: the method for determining the direction of the ground stress in the stratum at the target depth can determine the direction of the ground stress in the stratum at the target depth according to the direction of the ground stress in the core sample and the relative direction of the core sample and the sampling well wall of the core sample at the target depth, and specifically comprises the following steps: and determining the azimuth angle of the direction of the ground stress in the stratum at the target depth relative to the north direction according to the central angle of the direction of the ground stress in the core sample relative to the reference line and the included angle.
As shown in fig. 6, the central angle of the maximum horizontal principal stress direction of the stratum at the target depth relative to the reference line is 120 °, and the included angle between the connecting line of the reference line of the core sample in the circular plane perpendicular to the reference line and the center of the circular plane and the north direction is 288.2 °, so that the included angle between the ground stress direction of the stratum at the target depth and the north direction is 48.2 ° can be determined.
The technical scheme of the application uses logging data to determine the direction of the earth stress, and the method has strong theoretical basis, simple operation and higher precision, can be used for judging the direction of the earth stress of various reservoirs, particularly the stratum with little plastic change, is beneficial to optimizing the track of the well body, reasonably arranges the hydraulic fracturing direction and has higher practicability.
Example IV
The present embodiment provides a ground stress direction determining device, including: the sample processing module is used for acquiring a core sample of the stratum at the target depth, measuring longitudinal wave speeds along a plurality of different diameter directions when longitudinal waves propagate in the core sample, and determining the direction of the ground stress in the core sample according to the longitudinal wave speeds in the different diameter directions; the azimuth determining module is used for comparing the grain characteristics of the side face of the core sample with the grain characteristics of the sampling well wall of the core sample at the target depth, and determining the relative azimuth of the core sample and the sampling well wall of the core sample at the target depth according to the comparison result; and the ground stress determining module is used for determining the ground stress direction in the stratum at the target depth according to the direction of the ground stress in the core sample and the relative direction of the core sample and the sampling well wall of the core sample at the target depth.
In this embodiment, the ground stress direction determining apparatus may further include: a processor and a memory, wherein the processor is configured to execute the following program modules stored in the memory: the system comprises a sample processing module, an azimuth determining module and a ground stress determining module, so as to accurately measure the ground stress direction in the stratum at the target depth.
Example five
The present embodiment provides a storage medium storing a computer program, characterized in that the computer program, when executed by a processor, implements the steps of the method for determining a direction of ground stress as described above.
Storage media, including both permanent and non-permanent, removable and non-removable media, may be implemented in any method or technology for storage of information. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of storage media include, but are not limited to, phase change memory (PRAM), static Random Access Memory (SRAM), dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), read Only Memory (ROM), electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device.
Example six
The present embodiment provides a computer device comprising a processor and a storage medium storing program code which, when executed by the processor, implements the steps of the method for determining a direction of ground stress as described above.
In one embodiment, a computer device includes one or more processors (CPUs), an input/output interface, a network interface, and memory.
The memory may include volatile memory in a computer-readable medium, random Access Memory (RAM) and/or nonvolatile memory, etc., such as Read Only Memory (ROM) or flash memory (FLASH FLASH RAM). Memory is an example of computer-readable media.
It is noted that the terms used herein are used merely to describe particular embodiments and are not intended to limit exemplary embodiments in accordance with the present application, when the terms "comprising" and/or "including" are used in this specification, they specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof.
It should be noted that the terms "first," "second," and the like in the description and the claims and drawings of the present application are used for distinguishing between similar objects and not for describing a particular sequential or chronological order. It is to be understood that the terms so used are interchangeable under appropriate circumstances such that the embodiments of the application described herein are, for example, capable of operation in sequences other than those illustrated or otherwise described herein.
It should be understood that the exemplary embodiments in this specification may be embodied in many different forms and should not be construed as limited to only the embodiments set forth herein. These embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the concept of these exemplary embodiments to those skilled in the art, and should not be construed as limiting the application.
Claims (7)
1. A method for determining a direction of ground stress, comprising the steps of:
acquiring a core sample of a stratum at a target depth, measuring longitudinal wave speeds of longitudinal waves in a plurality of different diameter directions when the longitudinal waves propagate in the core sample, and determining the direction of ground stress in the core sample according to the longitudinal wave speeds of the different diameter directions;
comparing the grain characteristics of the side surface of the core sample with the grain characteristics of the sampling well wall of the core sample at the target depth, and determining the relative orientation of the core sample and the sampling well wall of the core sample at the target depth according to the comparison result;
Determining the direction of the ground stress in the stratum at the target depth according to the direction of the ground stress in the core sample and the relative direction of the core sample and the wall of the sampling well at the target depth;
Wherein the core sample is cylindrical, after obtaining the core sample of the formation at the target depth and before measuring the longitudinal wave velocity along a plurality of different diameter directions when the longitudinal wave propagates in the core sample, the method further comprises:
Setting a reference line parallel to the axial direction of the core sample on the side surface of the core sample;
After determining the direction of the earth stress in the core sample from the longitudinal wave velocities of the respective different diametric directions, the method further includes:
determining a central angle of the direction of the ground stress in the core sample relative to the reference line, and describing the direction of the ground stress in the core sample by using the central angle;
The method for determining the relative orientation of the core sample and the sampling well wall at the target depth comprises the following steps:
acquiring a rolling scanning image of a core sample and an electric imaging dynamic and static image of a sampling well wall of the core sample at a target depth;
Comparing the rolling scanning image with the electric imaging dynamic-static image, when the grain characteristics of the side surface of the core sample in the rolling scanning image are consistent with the grain characteristics of the sampling well wall of the core sample in the electric imaging dynamic-static image at the target depth, determining the posture of the core sample at the target depth in the sampling well, and determining the relative direction of the core sample and the sampling well wall at the target depth according to the posture of the core sample at the target depth in the sampling well;
The method for determining the relative orientation of the core sample and the sampling well wall at the target depth according to the attitude of the core sample at the target depth in the sampling well comprises the following steps:
determining a connecting line of the projection of a reference line of the core sample in a circular plane perpendicular to the reference line and the center of the circular plane, determining an included angle between the connecting line and the north direction under the posture of the core sample at the target depth in a sampling well, and describing the relative orientation of the core sample and the sampling well wall at the target depth by using the included angle.
2. The method of claim 1, wherein determining the direction of the earth stress in the formation at the target depth based on the direction of the earth stress in the core sample and the relative orientation of the core sample and its sampling wall of the well at the target depth comprises:
and determining the azimuth angle of the direction of the ground stress in the stratum at the target depth relative to the north direction according to the central angle of the direction of the ground stress in the core sample relative to the reference line and the included angle.
3. The method of claim 1, wherein the core sample is a full diameter core sample.
4. The method of determining the direction of ground stress according to claim 1, wherein determining the direction of ground stress in the core sample from the longitudinal wave velocities of the respective different diametric directions comprises:
Taking the diameter direction of the maximum longitudinal wave speed during propagation as the direction of the ground stress in the core sample; or (b)
The direction which is in the same circular plane as the diameter direction of the minimum longitudinal wave speed when propagating and is perpendicular to the diameter direction of the minimum longitudinal wave speed when propagating is taken as the direction of the ground stress in the core sample.
5. A ground stress direction determining apparatus based on the ground stress direction determining method according to any one of claims 1 to 4, characterized by comprising:
The sample processing module is used for acquiring a core sample of the stratum at the target depth, measuring longitudinal wave speeds along a plurality of different diameter directions when longitudinal waves propagate in the core sample, and determining the direction of the ground stress in the core sample according to the longitudinal wave speeds in the different diameter directions;
the azimuth determining module is used for comparing the grain characteristics of the side face of the core sample with the grain characteristics of the sampling well wall of the core sample at the target depth, and determining the relative azimuth of the core sample and the sampling well wall of the core sample at the target depth according to the comparison result;
And the ground stress determining module is used for determining the ground stress direction in the stratum at the target depth according to the direction of the ground stress in the core sample and the relative direction of the core sample and the sampling well wall of the core sample at the target depth.
6. A storage medium storing a computer program, wherein the computer program, when executed by a processor, implements the steps of the ground stress direction determining method according to any one of claims 1-4.
7. A computer device comprising a processor and a storage medium storing program code which, when executed by the processor, implements the steps of the ground stress direction determining method of any one of claims 1-4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010980973.1A CN114200538B (en) | 2020-09-17 | 2020-09-17 | Ground stress direction determining method, device, storage medium and computer equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010980973.1A CN114200538B (en) | 2020-09-17 | 2020-09-17 | Ground stress direction determining method, device, storage medium and computer equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114200538A CN114200538A (en) | 2022-03-18 |
| CN114200538B true CN114200538B (en) | 2024-04-30 |
Family
ID=80644815
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010980973.1A Active CN114200538B (en) | 2020-09-17 | 2020-09-17 | Ground stress direction determining method, device, storage medium and computer equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114200538B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115435936B (en) * | 2022-09-20 | 2023-10-10 | 成都理工大学 | Calculation method for ground stress of deep tight reservoir |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106033127A (en) * | 2016-06-29 | 2016-10-19 | 中国石油化工股份有限公司 | Ground stress orientation seismic prediction method based on shear wave speed variation rate |
| CN107829727A (en) * | 2017-09-12 | 2018-03-23 | 中国石油天然气股份有限公司 | A kind of method and device for determining crustal stress |
| CN110320277A (en) * | 2018-03-30 | 2019-10-11 | 中国石油天然气股份有限公司 | In-situ stress state test method |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100305927A1 (en) * | 2009-05-27 | 2010-12-02 | Schlumberger Technology Corporation | Updating a reservoir model using oriented core measurements |
| WO2015102508A1 (en) * | 2013-12-30 | 2015-07-09 | Schlumberger Holdings Limited | Method for estimating petrophysical properties of a hydrocarbon reservoir |
-
2020
- 2020-09-17 CN CN202010980973.1A patent/CN114200538B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106033127A (en) * | 2016-06-29 | 2016-10-19 | 中国石油化工股份有限公司 | Ground stress orientation seismic prediction method based on shear wave speed variation rate |
| CN107829727A (en) * | 2017-09-12 | 2018-03-23 | 中国石油天然气股份有限公司 | A kind of method and device for determining crustal stress |
| CN110320277A (en) * | 2018-03-30 | 2019-10-11 | 中国石油天然气股份有限公司 | In-situ stress state test method |
Non-Patent Citations (3)
| Title |
|---|
| 基于小波变换图像分割技术的电成像测井资料裂缝、孔洞面孔率提取方法;刘瑞林;谢芳;肖承文;张丽莉;;地球物理学报(第12期);全文 * |
| 深层地应力地理方位确定的新方法;卢运虎;陈勉;金衍;吴文;张家振;李远;;岩石力学与工程学报(第02期);全文 * |
| 黔南地区～3km油气深孔地应力测量与构造应力场分析;曹慧等;《中国地质》;第第47卷卷(第第1期期);第88-98页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114200538A (en) | 2022-03-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103175730B (en) | Method for determining jointed rock mechanical parameters | |
| CN105550448B (en) | Pre- drilling well three-dimensional well eye modeling method based on wellbore trace design parameter and device | |
| WO2016041189A1 (en) | Method for evaluating shale gas reservoir and seeking desert area | |
| US20100204972A1 (en) | Method For Predicting Well Reliability By Computer Simulation | |
| CN113919196A (en) | Reservoir three-dimensional stress field simulation method, simulation system, terminal and storage medium | |
| CN113901681B (en) | Three-dimensional compressibility evaluation method for dual desserts of shale gas reservoir in whole life cycle | |
| CN106014399B (en) | Method for establishing high-precision three-dimensional ground stress model of heterogeneous stratum | |
| Ouenes et al. | Integrated fractured reservoir modeling using both discrete and continuum approaches | |
| CN105549079A (en) | Method and device for establishing full-waveform inversion model for geophysics parameters | |
| CN109799540B (en) | Volcanic rock type uranium deposit magnetic susceptibility inversion method based on geological information constraint | |
| CN106154322A (en) | Log curve correction method and apparatus | |
| CN107045143A (en) | A kind of method and device of Prediction of fracture | |
| CN109033504A (en) | A kind of casing damage in oil-water well prediction technique | |
| CN105114063A (en) | Inclined-well dual laterolog curve correction method | |
| CN106295119A (en) | A kind of shale gas stratum crustal stress computational methods | |
| CN102943669A (en) | High-accuracy depth domain stochastic simulation inversion method suitable for designing track of horizontal well | |
| CN107728205A (en) | A kind of Formation pressure prediction method | |
| CN113779811B (en) | Quantitative analysis method for structural deformation to current ground stress orientation disturbance | |
| CN114066084B (en) | Method and system for predicting phase permeation curve based on machine learning | |
| CN114200538B (en) | Ground stress direction determining method, device, storage medium and computer equipment | |
| CN106021840A (en) | Method for inverting transverse isotropic rock mass ground stress | |
| CN102877828A (en) | CT (Computed Tomography) imaging method of three-dimensional multi-well combined well land | |
| Lofts et al. | Three-axis borehole gravity logging for reservoir surveillance | |
| CN110500089B (en) | Shale gas horizontal well bedding crack evaluation method based on array induction logging data | |
| CN103615239A (en) | Logging formation imaging system and method |
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 |