CN106321081B - For seeking the method and system of salt constituent stratigraphic anormaly crustal stress - Google Patents
For seeking the method and system of salt constituent stratigraphic anormaly crustal stress Download PDFInfo
- Publication number
- CN106321081B CN106321081B CN201510379614.XA CN201510379614A CN106321081B CN 106321081 B CN106321081 B CN 106321081B CN 201510379614 A CN201510379614 A CN 201510379614A CN 106321081 B CN106321081 B CN 106321081B
- Authority
- CN
- China
- Prior art keywords
- stress
- region
- crustal
- coefficient
- weakened
- 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
Abstract
The invention discloses a kind of for seeking the method and system of salt constituent stratigraphic anormaly crustal stress, this method comprises: determining the stress enhancement region, stress weakened region and unaffected zone on salt constituent stratum according to distribution of faults feature and geometric shape；The crustal stress of unaffected zone, stress enhancement region and stress weakened region is determined according to the well location crustal stress of unaffected zone；Salt constituent stratum ground stress model is established according to the crustal stress of unaffected zone, stress enhancement region and stress weakened region.The present invention considers various factors by the way of combination of qualitative and quantitative analysis, and the real brill data of integrated use numerical analysis and scene are modified existing computation model, improve ground stress model computational accuracy.
Description
Technical field
The invention belongs to petroleum engineering well drilling industries, specifically, being related to a kind of for seeking salt constituent stratigraphic anormaly
The method and system of stress.
Background technique
Crustal stress is internal stress present in stratum, is an important factor for influencing drilling engineering.In drilling engineering and open
In hair, the size and orientation of the tectonic stress in oil gas reservoir region are grasped, drilling engineering design can be optimized, oil gas field is opened
Well pattern arrangement and control wellbore stability are sent out, serious financial consequences and the person caused by accidents such as leakage is reduced or avoided, sprays, collapse, blocking
Accident etc..
Often assume that stratum is under the effect of three axis crustal stress, the principal stress of three of them principal direction is maximum horizontal principal stress
σH, minimum level crustal stress σhWith vertical stress σz, wherein vertical stress σzIt is horizontal mainly as caused by the self gravity of stratum
Crustal stress is changed by rock mass self weight, tectonic movement, formation fluid pressure and formation temperature to be generated.
Currently, having carried out some researchs to the calculation method for being influenced salt constituent stratigraphic anormaly crustal stress by local fracture.Text
Offer: the Finite-Element Study that influence and tomography based on measured data and numerical simulation tomography to crustal stress influence stress field refers to
It is influenced out by being broken, different well locations not far from one another show the biggish crustal stress numerical value of gap.The difference of crustal stress states
It is caused by being influenced due to tomography.It is anti-in conjunction with optimization using three-dimensional finite element simulation method based on FEM-software ANSYS
Artistic skills art carries out inverting to stress field and shows: since tomography influences and the superposition of compund fault disturbs, studying area's crustal stress
State table reveals complicated non-uniform Distribution feature.Document: fault region stress field Study on Forecasting Method is pointed out containing fault zone
Crustal stress distribution research in domain is the challenge with multi-solution, is difficult to adopt numerical method and is calculated, needs to use
Mechanical model, which carrys out sunykatuib analysis, to be attempted to apply discontinuous trend surface analysis, from the angle of geometry deformation, using small deflection
Plate theory is carried out according to the distribution that discontinuous trend surface principal curvatures analyzes the landform trend surface and tectonic stress field of fault region
Stress field Study on Forecasting Method.Document: influence and its engineering significance of the rift structure to stress field are pointed out to determine rift structure
The range that 1. rift structures influence is specifically included that on the parameter that stress field influences；2. local fault construction and regional stress field
Relationship, the influence of the angle α between influence and boundary stress direction including boundary stress ratio (Kb) and fracture；3. fracture
The influence of the physico-mechanical properties (deformation modulus, modulus of shearing, internal friction angle, cohesive force) of the rock of two sides；4. the power of fracture
Learn the influence of property (fracture normal stiffness, shear stiffness, internal friction angle, cohesive force)；5. the geometric shape of tomography and compound feelings
Condition.
Document above, which all shows to be broken, bigger influence to local crustal stress, and impacted factor also compares more.But
It is the crustal stress Numerical Predicting Method all without a set of suitable field application of proposition, it is more not this to saline bed by crustal stress shadow
The realistic problem for ringing biggish stratum carries out exact value analysis.
Summary of the invention
In order to solve the above problem, the method for seeking salt constituent stratigraphic anormaly crustal stress that the present invention provides a kind of and it is
System is carried out quasi- to provide a kind of crustal stress Numerical Predicting Method of suitable field application with the crustal stress to salt constituent stratum
True numerical analysis.
According to an aspect of the invention, there is provided a kind of method for seeking salt constituent stratigraphic anormaly crustal stress, packet
It includes:
According to distribution of faults feature and geometric shape determine the stress enhancement region on salt constituent stratum, stress weakened region and not by
The zone of influence；
The unaffected zone, the stress enhancement region and described are determined according to the well location crustal stress of the unaffected zone
The crustal stress of stress weakened region；
With establishing salt constituent according to the crustal stress of the unaffected zone, the stress enhancement region and the stress weakened region
Layer ground stress model.
According to one embodiment of present invention, determine that the unaffected zone, the stress enhancement region and the stress subtract
The crustal stress in weak area further comprises:
Well location crustal stress based on the unaffected zone determines regional structure stress coefficient, region crustal stress ratio and breaks
The angle of layer trend and region maximum and minimum level stress direction；
Based on the regional structure stress coefficient, the region crustal stress ratio and the fault strike and region maximum and most
The angle of small horizontal stress direction determines the crustal stress of the stress enhancement region and the stress weakened region.
According to one embodiment of present invention, the regional structure stress coefficient, the region crustal stress ratio and institute are based on
The angle for stating fault strike and region maximum and minimum level stress direction, determines that the stress enhancement region and the stress subtract
The crustal stress in weak area further comprises:
Well location crustal stress, the region crustal stress ratio and the fault strike and region based on the unaffected zone are most
Big and minimum level stress direction angle establishes numerical simulator；
Based on the well location crustal stress of the numerical simulator and the unaffected zone obtain the stress enhancement region and
The crustal stress variable quantity by braking effect of the stress weakened region；
Based on the stress enhancement region and the stress weakened region by braking effect crustal stress variable quantity, it is described not by
The well location crustal stress of the zone of influence obtains the tectonic stress coefficient of the stress enhancement region and the tectonic stress of the stress weakened region
Coefficient；
The tectonic stress coefficient or described of tectonic stress coefficient and the stress enhancement region based on the unaffected zone
Braking effect coefficient is calculated in the tectonic stress coefficient of stress weakened region；
The maximum horizontal of the stress enhancement region and the stress weakened region is calculated based on the braking effect coefficient
Crustal stress and minimum level crustal stress.
According to one embodiment of present invention, the numerical simulator includes stress enhancement region numerical simulator and answers
Power weakened region numerical simulator, stress enhancement region numerical simulator include the maximum horizontal crustal stress of stress enhancement region
Increment percentage-region crustal stress is than submodel, maximum horizontal stress-fault strike and maximum horizontal stress direction
Angle submodel, minimum level crustal stress increment percentage-region crustal stress are more disconnected than submodel and minimum level stress-
Layer trend and minimum level stress direction angle submodel；
Stress weakened region numerical simulator includes maximum horizontal crustal stress increment percentage-area of stress weakened region
Domain crustal stress than submodel and maximum horizontal stress-fault strike and maximum horizontal stress direction angle submodel,
Minimum level crustal stress increment percentage-region crustal stress is than submodel and minimum level stress-fault strike and most
Small horizontal stress direction angle submodel.
According to one embodiment of present invention, it answers to well location based on the numerical simulator and the unaffected zone
Power obtains the stress enhancement region and the crustal stress variable quantity by braking effect of the stress weakened region further comprises:
Pass through the stress enhancement region or maximum horizontal crustal stress increment percentage-region of stress weakened region
Crustal stress obtains the maximum horizontal crustal stress by braking effect of the stress enhancement region or the stress weakened region than submodel
Increment percentage；
Pass through maximum horizontal stress-rift direction of the stress enhancement region or the stress weakened region and maximum
Horizontal stress direction angle submodel obtains the maximum by braking effect of the stress enhancement region or the stress weakened region
Horizontal stress；
It is answered to region by the minimum level crustal stress increment percentage-of the stress enhancement region or the stress weakened region
Power obtains the minimum level crustal stress increment by braking effect of the stress enhancement region or the stress weakened region than submodel
Percentage；
Pass through minimum level stress-minimum level stress direction of the stress enhancement region or stress weakened region
The minimum level by braking effect of the stress enhancement region or stress weakened region is obtained with the angle submodel of rift direction
Stress difference；
Pass through maximum horizontal crustal stress increment percentage, the maximum horizontal of the stress enhancement region or the stress weakened region
The maximum horizontal principal stress of stress and the unaffected zone calculates the stress enhancement region or the stress weakened region
The maximum horizontal crustal stress variable quantity by braking effect, pass through the minimum water of the stress enhancement region or the stress weakened region
Described in the minimum level crustal stress of flat principal stress increment percentage, minimum level stress and the unaffected zone calculates
The minimum level crustal stress variable quantity by braking effect of stress enhancement region or the stress weakened region.
According to one embodiment of present invention, the stress enhancement region and institute are calculated based on the braking effect coefficient
The maximum horizontal crustal stress and minimum level crustal stress for stating stress weakened region further comprise introducing the braking effect coefficient
The maximum horizontal crustal stress and minimum level crustal stress of the stress enhancement region and the stress weakened region is calculated in following formula:
Wherein, σ 'HFor the maximum horizontal crustal stress of stress enhancement region or stress weakened region, σ 'hFor stress enhancement region or stress
The minimum level crustal stress of weakened region, ξ1For the maximum horizontal tectonic stress coefficient of unaffected zone, ξ2Most for unaffected zone
Small horizontal structural s tress coefficient, σ 'zFor the burden pressure of stress enhancement region or weakened region, P'pFor stress enhancement region or weakened region
Pore pressure, υ 'sFor the stratum static state Poisson's ratio of stress enhancement region or weakened region, E'sFor stress enhancement region or the weakened region
Elasticity modulus, α ' are the effective stress coefficient of stress enhancement region or weakened region, and T is the braking effect of stress enhancement region or weakened region
Coefficient.
According to one embodiment of present invention, the tectonic stress coefficient of the stress enhancement region and the stress weakened region is logical
Following formula is crossed to be calculated:
Wherein, σ 'HIt is the stress enhancement region that is calculated by numerical simulator or stress weakened region by fracture shadow
Loud maximum horizontal crustal stress, σ 'hFor the stress enhancement region that is calculated by numerical simulator or stress weakened region by
The minimum level crustal stress of braking effect, σ 'zFor the burden pressure of stress enhancement region or weakened region, υs' stress enhancement region or to answer
The stratum static state Poisson's ratio of power weakened region, E'sFor the elasticity modulus of stress enhancement region or weakened region, α ' is stress enhancement region or subtracts
The effective stress coefficient in weak area, P'pFor the pore pressure of stress enhancement region or weakened region.
According to one embodiment of present invention, the braking effect coefficient T is calculate by the following formula to obtain:
Wherein, ξ '1With ξ '2For the tectonic stress coefficient of stress enhancement region or stress weakened region, ξ '1For unaffected zone
Maximum horizontal tectonic stress coefficient, ξ '2For the minimum level tectonic stress coefficient of unaffected zone, ξ1And ξ2For unaffected zone
Tectonic stress coefficient, ξ1For the maximum horizontal tectonic stress coefficient of unaffected zone, ξ2For the minimum level structure of unaffected zone
Make stress coefficient.
According to one embodiment of present invention, the crustal stress for determining the unaffected zone further comprises:
Well location crustal stress based on the unaffected zone obtains formation rock fracture pressure, instantaneous characteristics, stratum
Pressure, Tensile Strength of Rock, breakpoint depth, drilling fluid density, effective stress coefficient, break area；
Based on the instantaneous characteristics, the drilling fluid density, described ground breakpoint depth and described break areal calculation
Obtain the minimum level crustal stress of the unaffected zone；
Minimum level crustal stress, the formation rock fracture pressure, the effective stress based on the unaffected zone
Coefficient, the strata pressure and the Tensile Strength of Rock calculate the maximum horizontal crustal stress of the unaffected zone.
According to another aspect of the present invention, it additionally provides a kind of for seeking by salt constituent stratigraphic anormaly crustal stress system
System, comprising:
Stress area division module determines that the stress on salt constituent stratum enhances according to distribution of faults feature and geometric shape
Area, stress weakened region and unaffected zone；
Crustal stress computing module determines the unaffected zone, institute according to the well location crustal stress in the uninfluenced region
State the crustal stress of stress enhancement region and the stress weakened region；
Ground stress model establishes module, according to the unaffected zone, the stress enhancement region and the stress weakened region
Crustal stress establish salt constituent stratum ground stress model.
Beneficial effects of the present invention:
The present invention proposes that a set of reasonable relatively accurate crustal stress numerical value is asked by the way of combination of qualitative and quantitative analysis
Method is taken, considers various factors, the real brill data of integrated use numerical analysis and scene repair existing computation model
Just, ground stress model computational accuracy is improved.
Other features and advantages of the present invention will be illustrated in the following description, also, partly becomes from specification
It obtains it is clear that understand through the implementation of the invention.The objectives and other advantages of the invention can be by specification, right
Specifically noted structure is achieved and obtained in claim and attached drawing.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is required attached drawing in technical description to do simple introduction:
Fig. 1 is method flow diagram according to an embodiment of the invention；
Fig. 2 is a fracture and well location schematic diagram；
Maximum horizontal crustal stress increment percentage-region of the stress enhancement region according to an embodiment of the invention Fig. 3 a
Crustal stress is than submodel schematic diagram；
Maximum horizontal stress-fault strike of the stress enhancement region according to an embodiment of the invention Fig. 3 b with
Maximum horizontal stress direction angle submodel schematic diagram；
Minimum level crustal stress increment percentage-region of the stress enhancement region according to an embodiment of the invention Fig. 4 a
Crustal stress is than submodel schematic diagram；And
Minimum level stress-fault strike of the stress enhancement region according to an embodiment of the invention Fig. 4 b with
Minimum level stress direction angle submodel schematic diagram.
Specific embodiment
Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings and examples, how to apply to the present invention whereby
Technological means solves technical problem, and the realization process for reaching technical effect can fully understand and implement.It needs to illustrate
As long as not constituting conflict, each feature in each embodiment and each embodiment in the present invention can be combined with each other,
It is within the scope of the present invention to be formed by technical solution.
As shown in Figure 1 for according to the method flow diagram of one embodiment of the present of invention, below with reference to Fig. 1 come to the present invention into
Row is described in detail.
Firstly, in step s 110, determining that the stress on salt constituent stratum enhances according to distribution of faults feature and geometric shape
Area, stress weakened region and unaffected zone.
In this step, by geological exploration techniques, form, shape and the combination side in salt constituent stratum breaking region are obtained
Formula, and find out according to the distribution of faults feature and geometric shape on stratum stress enhancement region, the stress weakened region in stratum breaking region
The unaffected zone and.Wherein, it is stress enhancement region that crustal stress, which is significantly greater than the region of the normal crustal stress in this area,.Apart from tomography compared with
Far, being influenced lesser region by tomography is unaffected zone.The region that stress is significantly less than the normal crustal stress in this area is stress
Weakened region.
Next, in the step s 120, determining unaffected zone according to the well location crustal stress situation in uninfluenced region, answering
The crustal stress of power enhancement region and stress weakened region.
In this step, firstly, the well location crustal stress situation based on unaffected zone, determines that the ground of unaffected zone is answered
Power, region crustal stress are than the angle with fault strike and region maximum horizontal stress direction.Wherein, the well location of unaffected zone
Crustal stress situation include according to well data, survey logging data and by break formation data that experiment and core experiment obtain etc.
The crustal stress size and Orientation obtained, specifically, well data includes, well depth opens secondary, strata division situation；Survey logging data
Including, sound wave, gamma, density, lithology, wellbore shape etc.；It includes formation rock fracture pressure, instantaneous termination of pumping that ground, which breaks experimental data,
Pressure, strata pressure, Tensile Strength of Rock, breakpoint depth, drilling fluid density, effective stress coefficient, break area；Rock core is real
The data for testing acquisition include the crustal stress size for the specific location that coring point rock core Kaiser effect obtains.
Specifically, in this step, in conjunction with experiment and laboratory core experiment is broken lively, determining the maximum of unaffected zone
Horizontal crustal stress and minimum level crustal stress, fault strike and region maximum horizontal stress direction and minimum level crustal stress side
To angle, region crustal stress ratio.
Specifically, fracture pressure and maximum, the relationship of minimum level crustal stress are as follows:
Pf=3 σh-σH-αPp+St (1)
Wherein, fracture pressure PfExperimental calculation is broken by ground to obtain.
Next, based on instantaneous characteristics, drilling fluid density, breakpoint depth and ground break areal calculation obtain it is described not
The minimum level crustal stress of affected area is calculate by the following formula to obtain minimum level crustal stress:
σh=Ps+10-6×ρmgHw (2)
Finally, minimum level crustal stress, formation rock fracture pressure, effective stress coefficient based on unaffected zone,
Stressor layer and Tensile Strength of Rock calculate the maximum horizontal crustal stress of unaffected zone, maximum horizontal is calculated based on following formula
Stress:
σH=3 σh-Pf-α·Pp+|St| (3)
Wherein, StFor Tensile Strength of Rock, acquired by rock mechanics parameters calculating；PsFor instantaneous characteristics, reality is broken by ground
Rupture curve data figure is tested to read；PfFor formation rock fracture pressure, obtained by rupture curve；PpFor strata pressure；α is effective
Stress coefficient is obtained by formation pressure calculation related data；ρm, HwRespectively drilling fluid density and ground breakpoint depth.Thus it obtains
The minimum level crustal stress σ of ground breakpointh, maximum horizontal crustal stress σHNumerical value, bring following formula into and obtain:
Tectonic stress coefficient (regional structure stress coefficient) ξ of unaffected zone can be obtained1And ξ2；σH,σh,σzFor region
Maximum horizontal, minimum level crustal stress and burden pressure, maximum burden pressure according in well-log information density log or sound
The value of wave well logging is calculated；PpFor pore pressure, calculated according to log data；υs、EsFor regional stratum static state Poisson
Than and elasticity modulus, calculated according to well-log information；α is region effective stress coefficient, is calculated according to well-log information.
Zone boundary stress axis is calculate by the following formula to obtain:
Unaffected zone will be obtained in relevant well location crustal stress situation substitution formula (1) of unaffected zone, (2), (3)
Minimum level crustal stress and maximum horizontal crustal stress, and pass through the zone boundary stress ratio that unaffected zone is calculated in formula (5).
The fault strike of unaffected zone and maximum and maximum horizontal stress direction angle are calculated by log data
It obtains.
Next, crustal stress, region crustal stress ratio (zone boundary stress axis), fault strike based on unaffected zone
With the angle of region maximum and maximum horizontal stress direction, the crustal stress of identified sign enhancement region and stress weakened region.Herein
The angle of region maximum and maximum horizontal stress direction is that the fault strike of unaffected zone is answered with maximum and maximum horizontal
The angle in power direction.
In this step, firstly, according to Theory of Fracture Mechanics, well location using discrete element method based on unaffected zone
Stress, region crustal stress are found out than establishing numerical simulator with the angle of fault strike and region maximum stress direction
Under various influence conditions, by the stress situation of braking effect.The numerical simulator includes stress enhancement region simulation mould
Type and stress weakened region simulation model.Wherein, stress enhancement region simulation model includes the maximum horizontal crustal stress of stress enhancement region
Increment percentage-region crustal stress is than submodel and maximum horizontal stress-fault strike and maximum horizontal crustal stress side
To angle submodel, minimum level crustal stress increment percentage-region crustal stress than submodel and minimum level stress-
Fault strike and minimum level stress direction angle submodel.
Stress weakened region simulation model includes maximum horizontal crustal stress increment percentage-region crustal stress of stress weakened region
Than submodel and maximum horizontal stress-fault strike and maximum horizontal stress direction angle submodel, minimum level
Crustal stress increment percentage-region crustal stress is than submodel and minimum level stress-fault strike and minimum level
Stress direction angle submodel.
Pass through maximum horizontal crustal stress increment percentage-region crustal stress ratio of stress enhancement region or the stress weakened region
Submodel is obtained in certain area crustal stress ratio, certain rift direction and maximum horizontal stress direction angle and centainly not by shadow
The maximum horizontal by braking effect of stress enhancement region or stress weakened region under the conditions of the maximum horizontal crustal stress in sound region
Stress increment percentage；
It with maximum horizontal is answered by the maximum horizontal of stress enhancement region or stress weakened region stress-rift direction
Power angular separation submodel is obtained in certain-rift direction and maximum horizontal stress direction angle, in certain area crustal stress
The maximum horizontal stress by braking effect of stress enhancement region or stress weakened region than under the conditions of；
Compare submodule by minimum level crustal stress increment percentage-region crustal stress of stress enhancement region or stress weakened region
Type obtains in certain area crustal stress ratio, certain-rift direction and minimum level stress direction angle and uninfluenced region
Minimum level crustal stress under the conditions of stress enhancement region or stress weakened region by braking effect minimum level crustal stress increase
Measure percentage；
By the minimum level of stress enhancement region or stress weakened region stress-minimum level stress direction and break
The angle submodel for splitting direction, which is obtained, answers in certain-rift direction with minimum level stress direction angle, in certain area
The minimum level stress by braking effect of stress enhancement region or stress weakened region under the conditions of power ratio；
By stress enhancement region or stress weakened region by the maximum horizontal crustal stress increment percentage of braking effect, by disconnected
The maximum horizontal principal stress of the maximum horizontal stress and unaffected zone that split influence calculates stress enhancement region or stress subtracts
The maximum horizontal crustal stress variable quantity by braking effect in weak area, by stress enhancement region or stress weakened region by braking effect
Minimum horizontal principal stress increment percentage, by the minimum level stress of braking effect and the minimum water of unaffected zone
The minimum level crustal stress variable quantity by braking effect of level land Stress calculation stress enhancement region or stress weakened region.
The maximum horizontal crustal stress variable quantity and minimum level crustal stress by braking effect based on stress enhancement region changes
The maximum horizontal principal stress and minimum level crustal stress of stress enhancement region is calculated in amount.By the stress enhancement region being calculated
Maximum horizontal crustal stress and minimum level crustal stress substitute into the Seventh Five-Year Plan model inverse, i.e. substitution formula (4) obtains the structure of stress enhancement region
Make stress coefficient ξ '1With ξ '2:
Wherein, σ 'HIt is the stress enhancement region that is calculated by numerical simulator or stress weakened region by fracture shadow
Loud maximum horizontal crustal stress variable quantity, σ 'hWeaken for the stress enhancement region being calculated by numerical simulator or stress
The minimum level crustal stress variable quantity by braking effect in area, υs' moored for the stratum static state of stress enhancement region or stress weakened region
Loose ratio, E'sFor the elasticity modulus of stress enhancement region or weakened region, α ' is the effective stress coefficient of stress enhancement region or weakened region,
P'pFor the pore pressure of stress enhancement region or weakened region.
Stress coefficient ξ ' through this construction1Or ξ '2, the tectonic stress coefficient ξ of unaffected zone1Or ξ2, counted by formula (6)
Calculation obtains braking effect coefficient:
Wherein, the tectonic stress coefficient ξ ' of stress enhancement region is calculated by formula (4-1)1With ξ '2When, σ ' thereinHTo pass through
The maximum horizontal crustal stress by braking effect for the stress enhancement region that simulation model is calculated；σ'hTo pass through simulation model meter
The minimum level crustal stress by braking effect of obtained stress enhancement region；υ'sFor the stratum static state Poisson of stress enhancement region
Than；E'sFor the elasticity modulus of stress enhancement region；α ' is the effective stress coefficient of stress enhancement region, P'pFor the hole of stress enhancement region
Gap pressure.
When stress enhancement region or when having drilling well of stress weakened region, the well location crustal stress data that can use drilling well are asked
Take the corresponding tectonic stress coefficient of the well location, then the tectonic stress coefficient obtained based on the well location is to the tectonic stress being calculated
Coefficient ξ '1Or ξ '2It is modified.Specifically, the average value of these tectonic stress coefficients can be used, or design weight coefficient
Seek final tectonic stress coefficient ξ '1With ξ '2。
Finally, in step s 130, establishing salt according to the crustal stress of unaffected zone, stress enhancement region and stress weakened region
Cream rock stratum ground stress model.
In this step, using the maximum horizontal crustal stress and minimum level of the stress enhancement region of the Seventh Five-Year Plan model pair of optimization
Crustal stress, the maximum horizontal crustal stress of stress weakened region and minimum level crustal stress optimize.
The Seventh Five-Year Plan model is to assume that stratum is homogeneous isotropic line elastomer, and assume that phase geological structure is transported after deposit
During dynamic, relative displacement does not occur between stratum and stratum, the strain of all two horizontal directions of stratum is constant, then maximum
Horizontal crustal stress σH, minimum level crustal stress σhIt is expressed as follows:
In formula: ξ1, ξ2For the tectonic stress coefficient in the region, ξ1For maximum horizontal tectonic stress coefficient, ξ2For minimum level
Tectonic stress coefficient, region can be unaffected zone, stress enhancement region or stress weakened region；σzFor the burden pressure in the region,
Maximum burden pressure is calculated according to the value of density log or acoustic logging in well-log information；PpFor pore pressure, root
It is calculated according to log data；υs、EsFor stratum static state Poisson's ratio and elasticity modulus, calculated according to well-log information；α is to have
Effect force coefficient, is calculated according to well-log information.Wherein, maximum horizontal crustal stress and minimum level crustal stress pass through above
Formula (1), (2), (3) are calculated.
Next, optimizing to obtain to the Seventh Five-Year Plan model by braking effect coefficient T:
Wherein, T is braking effect coefficient.
When specifically being calculated using formula (8), the corresponding numerical value in stress enhancement region or stress weakened region is substituted into formula (8) and is obtained
It arrives:
Wherein, σ 'HFor the maximum horizontal of the stress enhancement region or stress weakened region that are calculated by numerical simulator
Crustal stress, σ 'hFor the minimum level crustal stress of the stress enhancement region or stress weakened region that are calculated by numerical simulator,
ξ1For the maximum horizontal tectonic stress coefficient of unaffected zone, ξ2For the minimum level tectonic stress coefficient of unaffected zone, σ 'z
For the burden pressure of stress enhancement region or stress weakened region, P'pFor the pore pressure of stress enhancement region or stress weakened region, υ 'sFor
The stratum static state Poisson's ratio of stress enhancement region or stress weakened region, E'sFor the elasticity modulus of stress enhancement region or stress weakened region,
α ' is the effective stress coefficient of stress enhancement region or stress weakened region, and T is the braking effect of stress enhancement region or stress weakened region
Coefficient.
Answer to the maximum horizontal of maximum horizontal crustal stress and minimum level crustal stress, stress weakened region to stress enhancement region
After power and minimum level crustal stress optimize, next, (such as according to the well location crustal stress situation of stress enhancement region and weakened region
The area has well location crustal stress) to the maximum horizontal crustal stress and minimum level crustal stress, stress of the stress enhancement region after optimization
The maximum horizontal crustal stress and minimum level crustal stress of weakened region are tested and are calibrated.Specifically, the rock obtained using scene
The heart carries out Rock Mechanics Test, carries out Kaiser effect and tests to obtain two groups to three groups of parameters and earth stress of coring point, with calculating
Result test and calibrate.
Finally, maximum horizontal crustal stress and minimum level crustal stress based on the stress enhancement region after inspection and calibration, answering
The maximum horizontal crustal stress and minimum level crustal stress of power weakened region establish salt constituent reservoir stress model with interpolation algorithm.
Feasibility of the invention is verified below by way of a specific example.It is illustrated in figure 2 and is locally broken
It splits the specific well location of salt constituent stratum block of influence and is broken distribution mode and geometry schematic diagram.As shown in Fig. 2, influencing wheat
The fracture of Gai Tiyi block is mainly the fracture of Ma south, and No. 5 structural belts, main body is two parallel fractures, and outside is lower wall, is broken it
Interior is upper disk.
It as shown in Fig. 2, jade Bei1Jing, Yu Bei 1-1 well, Yu Bei 1-2 well are between two parallel fractures, and is fracture
Upper disk position.Beautiful 6 well of Bei5Jing, Yu Bei is in the two sides of parallel fracture, and belongs to the lower wall position of fracture.According to fracture
Influence to local crustal stress can judge beautiful Bei5Jing and beautiful northern 6 well stress situations and jade Bei1Jing, Yu Bei 1-1 well, Yu Bei
The stress situation of 1-2 well has a difference, and the size of difference is by the corner dimension of rift direction and region biggest principal stress direction, disconnected
Split the influence of the rock strength between two sides and fracture, the internal friction angle size for being broken itself.
According to the actual conditions of beautiful northern 1 block, there is different degrees of casing collapse in 5 well of Yu Bei and jade Bei6Jing
Situation illustrates that the crustal stress for the two side areas being broken under conditions of local structure influences close to Ma south will be significantly greater than the ground
The normal crustal stress in area illustrates that this two well location enhances region in stress.And locating for jade Bei1Jing, Yu Bei 1-1 well and beautiful north 1-2 well
Two fracture belts inside influenced by local structure stress, the crustal stress in the region should be lower than the normal crustal stress in this area
Situation, this three well location is in stress weakened region.
Beautiful 8 well of Bei4Jing, Yu Bei and jade Bei9Jing, farther out apart from tomography, crustal stress is mainly by overburden pressure and region
Construction determines, is influenced by local structure smaller, belongs to the normal crustal stress situation in this area, the well not influenced by local structure
Position is located at unaffected zone.
Each well is as shown in table 1 with respect to the specific location of tomography geometric shape in Fig. 2.
Beautiful Bei7Jing is broken north side in Ma south, is broken northwest side in a single line, should answer on the whole with beautiful northern 5 wells
The trend of power is consistent, but relatively weak, by local structure influenced the result is that normal crustal stress should be slightly higher than.Beautiful north 3
Well from Ma south is main due to being broken farther out, and nearby offset well is less, and concrete condition also needs to continue to judge.
Table 1
Next, according to the well location crustal stress situation in uninfluenced region, based on formula (1), (2), (3) be calculated not by
The minimum level crustal stress and maximum horizontal crustal stress of the zone of influence.It is then based on the construction that unaffected zone is calculated in formula (4)
Stress coefficient ξ1=0.778 and ξ2=0.368, and then zone boundary stress axis is calculated by formula (5)Can also obtain region maximum horizontal stress direction by well location crustal stress situation is nearly North and South direction
(N4 °/184 °), fracture direction and maximum horizontal crustal stress angle are about 45 °, and the direction of crustal stress refers to uninfluenced region
Crustal stress direction and azimuthal form, according to Image Logging Data read.
Based on numerical simulator according to beautiful backlands region stress ratio KbValue about 1.26, fault strike and maximum horizontal
Crustal stress angle is about 45 °, according to the numerical simulation result of Fig. 3 a, 3b, according to the concrete condition of embodiment block, it is known that fracture
Influence to crustal stress is about 25% or so of region crustal stress, by the maximum principal stress difference under the conditions of the angle multiplied by 25%
The maximum horizontal crustal stress that stress enhancement region can be obtained is added with the maximum horizontal crustal stress of unaffected zone again.Stress enhancing
The minimum level crustal stress numerical model in area obtains the minimum level of stress enhancement region as shown in Fig. 4 a, 4b by the above process
Crustal stress.Bring numerical value into formula (4) to obtain: it is about 0.085 that tomography, which influences coefficient,；Stress enhancement region tectonic stress coefficient: ξ1=
0.844, ξ2=0.399；Stress weakened region tectonic stress coefficient: ξ1=0.712, ξ2=0.337.
Next, based on tomography influence coefficient about T, using the Seventh Five-Year Plan modular form (8-1) of optimization to stress enhancement region most
Horizontally stress and minimum level crustal stress, the maximum horizontal crustal stress of stress weakened region and the progress of minimum level crustal stress are excellent greatly
Change.And then examined and calibrated, salt constituent reservoir stress model is established with interpolation algorithm.
According to another aspect of the present invention, a kind of saline bed for seeking by braking effect is additionally provided to answer singularly
The system of power, including stress area division module, crustal stress computing module and ground stress model establish module.
Wherein, stress area division module determines the stress on salt constituent stratum according to distribution of faults feature and geometric shape
Enhancement region, stress weakened region and unaffected zone.Crustal stress computing module is determined according to the well location crustal stress in uninfluenced region
Unaffected zone, stress enhancement region and stress weakened region crustal stress.Ground stress model establishes module, according to unaffected zone,
The crustal stress of stress enhancement region and stress weakened region establishes salt constituent stratum ground stress model.
The mode that the present invention is combined by the way of qualitative and quantitative proposes a set of reasonable relatively accurate crustal stress
Numerical value acquiring method considers various factors, and integrated use numerical analysis, scene is real to bore data, to existing computation model into
Row amendment, improves the computational accuracy of crustal stress.
While it is disclosed that embodiment content as above but described only to facilitate understanding the present invention and adopting
Embodiment is not intended to limit the invention.Any those skilled in the art to which this invention pertains are not departing from this
Under the premise of the disclosed spirit and scope of invention, any modification and change can be made in the implementing form and in details,
But scope of patent protection of the invention, still should be subject to the scope of the claims as defined in the appended claims.
Claims (7)
1. a kind of method for seeking salt constituent stratigraphic anormaly crustal stress, comprising:
The stress enhancement region on salt constituent stratum, stress weakened region and uninfluenced are determined according to distribution of faults feature and geometric shape
Area；
The unaffected zone, the stress enhancement region and the stress are determined according to the well location crustal stress of the unaffected zone
The crustal stress of weakened region；
With establishing salt constituent stratum according to the crustal stress of the unaffected zone, the stress enhancement region and the stress weakened region
Stress model；
Wherein it is determined that the crustal stress of the unaffected zone, the stress enhancement region and the stress weakened region further comprises:
Well location crustal stress based on the unaffected zone determines that regional structure stress coefficient, region crustal stress ratio and tomography are walked
To the angle with region maximum and minimum level stress direction；
Well location crustal stress, the region crustal stress ratio and the fault strike based on the unaffected zone and region it is maximum and
The angle of minimum level stress direction establishes numerical simulator；
The stress enhancement region and described is obtained based on the well location crustal stress of the numerical simulator and the unaffected zone
The crustal stress variable quantity by braking effect of stress weakened region；
The crustal stress variable quantity, described uninfluenced by braking effect based on the stress enhancement region and the stress weakened region
The well location crustal stress in area obtains the tectonic stress coefficient of the stress enhancement region and the tectonic stress coefficient of the stress weakened region；
The tectonic stress coefficient or the stress of tectonic stress coefficient and the stress enhancement region based on the unaffected zone
Braking effect coefficient is calculated in the tectonic stress coefficient of weakened region；
It is answered the maximum horizontal of the stress enhancement region and the stress weakened region is calculated based on the braking effect coefficient
Power and minimum level crustal stress；
Wherein, the braking effect coefficient T is calculate by the following formula to obtain:
Wherein, ξ '1With ξ '2For the tectonic stress coefficient of stress enhancement region or stress weakened region, ξ1For the most flood of unaffected zone
Flat tectonic stress coefficient, ξ2For the minimum level tectonic stress coefficient of unaffected zone.
2. the method according to claim 1, wherein the numerical simulator includes stress enhancement region Numerical-Mode
Analog model and stress weakened region numerical simulator, stress enhancement region numerical simulator include the maximum of stress enhancement region
Horizontally stress increment percentage-region crustal stress is than submodel, maximum horizontal stress-fault strike and maximum horizontal
Stress direction angle submodel, minimum level crustal stress increment percentage-region crustal stress are than submodel and minimum level
Stress difference-fault strike and minimum level stress direction angle submodel；
Stress weakened region numerical simulator includes maximum horizontal crustal stress increment percentage-region of stress weakened region
Stress ratio submodel and maximum horizontal stress-fault strike and maximum horizontal stress direction angle submodel, minimum
Horizontally stress increment percentage-region crustal stress is than submodel and minimum level stress-fault strike and minimum water
Level land stress direction angle submodel.
3. the method according to claim 1, wherein being based on the numerical simulator and the unaffected zone
Well location crustal stress obtain the crustal stress variable quantity by braking effect of the stress enhancement region and the stress weakened region into one
Step includes:
It is answered to region by the maximum horizontal crustal stress increment percentage-of the stress enhancement region or stress weakened region
Power obtains the maximum horizontal crustal stress increment by braking effect of the stress enhancement region or the stress weakened region than submodel
Percentage；
Pass through the maximum horizontal stress-rift direction and maximum horizontal of the stress enhancement region or the stress weakened region
Stress direction angle submodel obtains the maximum horizontal by braking effect of the stress enhancement region or the stress weakened region
Stress；
Pass through minimum level crustal stress increment percentage-region crustal stress ratio of the stress enhancement region or the stress weakened region
Submodel obtains the minimum level crustal stress increment percentage by braking effect of the stress enhancement region or the stress weakened region
Than；
By minimum level stress-minimum level stress direction of the stress enhancement region or stress weakened region and break
The angle submodel for splitting direction obtains the minimum level crustal stress by braking effect of the stress enhancement region or stress weakened region
Difference；
By the maximum horizontal crustal stress increment percentage of the stress enhancement region or the stress weakened region, answer to maximum horizontal
Power difference and the maximum horizontal principal stress of the unaffected zone calculate the stress enhancement region or the stress weakened region by
The maximum horizontal crustal stress variable quantity of braking effect, passes through the minimum level master of the stress enhancement region or the stress weakened region
The minimum level crustal stress of stress increment percentage, minimum level stress and the unaffected zone calculates the stress
The minimum level crustal stress variable quantity by braking effect of enhancement region or the stress weakened region.
4. the method according to claim 1, wherein the stress is calculated based on the braking effect coefficient
The maximum horizontal crustal stress and minimum level crustal stress of enhancement region and the stress weakened region further comprise, by the fracture shadow
It rings coefficient and introduces the stress enhancement region and the stress weakened region is calculated in following formula maximum horizontal crustal stress and minimum water
Level land stress:
Wherein, σ 'HFor the maximum horizontal crustal stress of stress enhancement region or stress weakened region, σ 'hWeaken for stress enhancement region or stress
The minimum level crustal stress in area, ξ1For the maximum horizontal tectonic stress coefficient of unaffected zone, ξ2For the minimum water of unaffected zone
Flat tectonic stress coefficient, σ 'zFor the burden pressure of stress enhancement region or weakened region, P'pFor the hole of stress enhancement region or weakened region
Pressure, υ 'sFor the stratum static state Poisson's ratio of stress enhancement region or weakened region, E'sFor the elasticity of stress enhancement region or the weakened region
Modulus, α ' are the effective stress coefficient of stress enhancement region or weakened region, and T is the braking effect system of stress enhancement region or weakened region
Number.
5. the method according to claim 1, wherein the construction of the stress enhancement region and the stress weakened region
Stress coefficient is calculate by the following formula to obtain:
Wherein, σ 'HFor the stress enhancement region that is calculated by numerical simulator or stress weakened region by braking effect most
Big horizontally stress, σ 'hIt is the stress enhancement region that is calculated by numerical simulator or stress weakened region by fracture shadow
Loud minimum level crustal stress, σ 'zFor the burden pressure of stress enhancement region or weakened region, υs' weaken for stress enhancement region or stress
The stratum static state Poisson's ratio in area, E'sFor the elasticity modulus of stress enhancement region or weakened region, α ' is stress enhancement region or weakened region
Effective stress coefficient, P'pFor the pore pressure of stress enhancement region or weakened region.
6. the method according to claim 1, wherein determining that the crustal stress of the unaffected zone further wraps
It includes:
Well location crustal stress based on the unaffected zone obtain formation rock fracture pressure, instantaneous characteristics, strata pressure,
Tensile Strength of Rock, breakpoint depth, drilling fluid density, effective stress coefficient, break area；
Based on the instantaneous characteristics, the drilling fluid density, described ground breakpoint depth and described break areal calculation and obtain
The minimum level crustal stress of the unaffected zone；
Minimum level crustal stress, the formation rock fracture pressure, the effective stress coefficient based on the unaffected zone,
The strata pressure and the Tensile Strength of Rock calculate the maximum horizontal crustal stress of the unaffected zone.
7. a kind of system for seeking salt constituent stratigraphic anormaly crustal stress, comprising:
Stress area division module determines the stress enhancement region on salt constituent stratum according to distribution of faults feature and geometric shape, answers
Power weakened region and unaffected zone；
Crustal stress computing module determines the unaffected zone according to the well location crustal stress in the uninfluenced region, described answers
The crustal stress of power enhancement region and the stress weakened region；
Ground stress model establishes module, according to the ground of the unaffected zone, the stress enhancement region and the stress weakened region
Stress establishes salt constituent stratum ground stress model；Wherein, the crustal stress computing module by following steps determine it is described not by
The crustal stress of the zone of influence, the stress enhancement region and the stress weakened region:
Well location crustal stress based on the unaffected zone determines that regional structure stress coefficient, region crustal stress ratio and tomography are walked
To the angle with region maximum and minimum level stress direction；
Well location crustal stress, the region crustal stress ratio and the fault strike based on the unaffected zone and region it is maximum and
The angle of minimum level stress direction establishes numerical simulator；
The stress enhancement region and described is obtained based on the well location crustal stress of the numerical simulator and the unaffected zone
The crustal stress variable quantity by braking effect of stress weakened region；
The crustal stress variable quantity, described uninfluenced by braking effect based on the stress enhancement region and the stress weakened region
The well location crustal stress in area obtains the tectonic stress coefficient of the stress enhancement region and the tectonic stress coefficient of the stress weakened region；
The tectonic stress coefficient or the stress of tectonic stress coefficient and the stress enhancement region based on the unaffected zone
Braking effect coefficient is calculated in the tectonic stress coefficient of weakened region；
It is answered the maximum horizontal of the stress enhancement region and the stress weakened region is calculated based on the braking effect coefficient
Power and minimum level crustal stress；
Wherein, the braking effect coefficient T is calculate by the following formula to obtain:
Wherein, ξ '1With ξ '2For the tectonic stress coefficient of stress enhancement region or stress weakened region, ξ1For the most flood of unaffected zone
Flat tectonic stress coefficient, ξ2For the minimum level tectonic stress coefficient of unaffected zone.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510379614.XA CN106321081B (en) | 2015-07-01 | 2015-07-01 | For seeking the method and system of salt constituent stratigraphic anormaly crustal stress |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510379614.XA CN106321081B (en) | 2015-07-01 | 2015-07-01 | For seeking the method and system of salt constituent stratigraphic anormaly crustal stress |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106321081A CN106321081A (en) | 2017-01-11 |
| CN106321081B true CN106321081B (en) | 2019-03-15 |
Family
ID=57726876
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510379614.XA Active CN106321081B (en) | 2015-07-01 | 2015-07-01 | For seeking the method and system of salt constituent stratigraphic anormaly crustal stress |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106321081B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110107276B (en) * | 2018-01-30 | 2022-06-03 | 中国石油天然气股份有限公司 | Casing design method and apparatus |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4665984A (en) * | 1985-08-29 | 1987-05-19 | Tohoku University | Method of measuring crustal stress by hydraulic fracture based on analysis of crack growth in rock |
| JP2004239640A (en) * | 2003-02-03 | 2004-08-26 | Mutsuto Kawahara | Method of predicting geological properties in front of ground excavation part |
| CN101162177A (en) * | 2007-11-14 | 2008-04-16 | 南京银茂铅锌矿业有限公司 | Method for measuring ground stress |
| CN103091727A (en) * | 2013-01-11 | 2013-05-08 | 上海交通大学 | Geostress field distributed data management (DDM) retrieval method based on actual measurement results and linear programming method |
| CN103838936A (en) * | 2014-03-25 | 2014-06-04 | 中国石油大学(华东) | High-precision tectonic stress field simulation method applicable to turbidite sand low-permeability reservoirs |
| CN103983490A (en) * | 2014-05-26 | 2014-08-13 | 中国科学院武汉岩土力学研究所 | Preparation method of disturbed belt sample for reflecting stress state of primary rock |
| CN104088666A (en) * | 2014-07-09 | 2014-10-08 | 中国电建集团华东勘测设计研究院有限公司 | In-situ test method and structure for deeply-buried soft-rock tunnel |
| CN104360013A (en) * | 2014-09-10 | 2015-02-18 | 山东科技大学 | Analytical method of deep mining floor rock stratum damage and fracturation evolution and failure law |
| CN104463717A (en) * | 2014-12-09 | 2015-03-25 | 中国水利水电第十四工程局有限公司 | Construction expert information system applicable to extra-large underground cavern group |
| CN105804731A (en) * | 2014-12-30 | 2016-07-27 | 中国石油天然气股份有限公司 | Rock ground stress detection method and system |
-
2015
- 2015-07-01 CN CN201510379614.XA patent/CN106321081B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4665984A (en) * | 1985-08-29 | 1987-05-19 | Tohoku University | Method of measuring crustal stress by hydraulic fracture based on analysis of crack growth in rock |
| JP2004239640A (en) * | 2003-02-03 | 2004-08-26 | Mutsuto Kawahara | Method of predicting geological properties in front of ground excavation part |
| CN101162177A (en) * | 2007-11-14 | 2008-04-16 | 南京银茂铅锌矿业有限公司 | Method for measuring ground stress |
| CN103091727A (en) * | 2013-01-11 | 2013-05-08 | 上海交通大学 | Geostress field distributed data management (DDM) retrieval method based on actual measurement results and linear programming method |
| CN103838936A (en) * | 2014-03-25 | 2014-06-04 | 中国石油大学(华东) | High-precision tectonic stress field simulation method applicable to turbidite sand low-permeability reservoirs |
| CN103983490A (en) * | 2014-05-26 | 2014-08-13 | 中国科学院武汉岩土力学研究所 | Preparation method of disturbed belt sample for reflecting stress state of primary rock |
| CN104088666A (en) * | 2014-07-09 | 2014-10-08 | 中国电建集团华东勘测设计研究院有限公司 | In-situ test method and structure for deeply-buried soft-rock tunnel |
| CN104360013A (en) * | 2014-09-10 | 2015-02-18 | 山东科技大学 | Analytical method of deep mining floor rock stratum damage and fracturation evolution and failure law |
| CN104463717A (en) * | 2014-12-09 | 2015-03-25 | 中国水利水电第十四工程局有限公司 | Construction expert information system applicable to extra-large underground cavern group |
| CN105804731A (en) * | 2014-12-30 | 2016-07-27 | 中国石油天然气股份有限公司 | Rock ground stress detection method and system |
Non-Patent Citations (2)
| Title |
|---|
| 深埋大断层构造区三维地应力场反演分析;黄耀光;《煤矿开采》;20140630;第19卷(第3期);第24-28页 |
| 玉北古近系膏泥层地应力分布规律探索;唐世春;《石油钻探技术》;20130131;第41卷(第1期);第31-34页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106321081A (en) | 2017-01-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gong et al. | A new approach for characterization and prediction of natural fracture occurrence in tight oil sandstones with intense anisotropy | |
| Ju et al. | Tectonic fractures in the Lower Cretaceous Xiagou Formation of Qingxi Oilfield, Jiuxi Basin, NW China. Part two: Numerical simulation of tectonic stress field and prediction of tectonic fractures | |
| CN108680952B (en) | A kind of strike-slip fault structural evolution analytic method | |
| Yıldırım | Relative tectonic activity assessment of the Tuz Gölü fault zone; Central Anatolia, Turkey | |
| Negredo et al. | Neotectonic modeling of the Ibero‐Maghrebian region | |
| CN106484958B (en) | A kind of system of the three-dimensional cracking permeability tensor computation model based on pit shaft gap observation | |
| Pinel et al. | A two‐step model for dynamical dike propagation in two dimensions: Application to the July 2001 Etna eruption | |
| CN103914620B (en) | Method for computing distribution of opening spaces of fractures of fault rupture zones | |
| CN106321082B (en) | For seeking the method and system of salt constituent stratum maximum outside casing load | |
| CN110969699B (en) | High-precision three-dimensional ground stress modeling method and application thereof | |
| Feng et al. | Near-surface stress measurements in the Longmenshan fault belt after the 2008 Wenchuan Ms8. 0 earthquake | |
| CN106815412B (en) | Simulation method and device for structural stress field | |
| Neves et al. | The effect of sea level changes on fault reactivation potential in Portugal | |
| Rodriguez-Herrera et al. | Field-scale geomechanical characterization of the Haynesville shale | |
| Xiao‐Guang et al. | Effects of topography and local geological structure on computing co‐seismic deformation—A case study of the 2011 Japan Tohoku earthquake (Mw 9.0) | |
| Mejia et al. | Geomechanical effects of natural fractures on fluid flow in a pre-salt field | |
| CN106321081B (en) | For seeking the method and system of salt constituent stratigraphic anormaly crustal stress | |
| Feng et al. | An improved geomechanical model for the prediction of fracture generation and distribution in brittle reservoirs | |
| Roberts | Fault orientation variations along the strike of active normal fault systems in Italy and Greece: Implications for predicting the orientations of subseismic-resolution faults in hydrocarbon reservoirs | |
| Savage et al. | The effect of non-parallel thrust fault interaction on fold patterns | |
| Zeng et al. | Variation in the orientation of the maximum horizontal stress in thick channel-fill sandstones with low-permeability: A case of the Bonan Oilfield in the Bohai Bay Basin, eastern China | |
| Manchanda et al. | Creation of multi-stranded fractures in homogeneous rocks | |
| Ulutaş | The May 11 Paphos, Cyprus, earthquake: implications for stress regime and tsunami modelling for the Eastern Mediterranean shorelines | |
| Romano et al. | Non-isothermal momentum transfer and ground displacements rate at Campi Flegrei caldera (Southern Italy) | |
| Petrie et al. | Modeling strain across mechanical sedimentary lithologic interfaces: geomechanical models derived from outcrop analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |