CN109492238A - A method of using ground stress analysis tomography and fisstured flow characteristic - Google Patents
A method of using ground stress analysis tomography and fisstured flow characteristic Download PDFInfo
- Publication number
- CN109492238A CN109492238A CN201710818821.XA CN201710818821A CN109492238A CN 109492238 A CN109492238 A CN 109492238A CN 201710818821 A CN201710818821 A CN 201710818821A CN 109492238 A CN109492238 A CN 109492238A
- Authority
- CN
- China
- Prior art keywords
- stress
- formula
- cff
- tomography
- coulomb
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Abstract
The present invention provides a kind of method using ground stress analysis tomography and fisstured flow characteristic, includes the following steps: the dynamic 3 D geomechanics model and discrete fracture network model of 1) establishing target work area respectively；2) by six stress parameters σ in dynamic 3 D geomechanics modelx、σy、σz、τxy、τxzAnd τyzAccordingly it is loaded on the tomography and fracture surface of discrete fracture network model；3) according to the tilt angle gamma in tomography and crack and tendency θ, three stress cosine coefficients l, m and n are obtained；4) according to three stress cosine coefficients and six stress parameters, three components of stress p are obtainedx、py、pz；5) according to three stress cosine coefficients and three components of stress, direct stress σ is obtainedN；6) according to three components of stress and direct stress σN, obtain shear stress τN；7) according to the direct stress σNWith shear stress τN, obtain a mole coulomb function CFF.Method of the invention can accurately analyze the seepage characteristic in tomography and crack.
Description
Technical field
The present invention relates to a kind of Seepage characteristic analysis methods, and in particular to a kind of application ground stress analysis tomography and crack are seeped
The method of properties of flow.
Background technique
Fisstured flow characteristic is a basic research of fractured reservoirs, passes through geotechnic stress field control critical eigenvalue in the past
Formation, distribution and development degree, Present Stress Field not only controls the form and extending direction in artificial fracturing crack, also controls
Percolation ability and validity of the intrinsic fracture in underground.Oil and gas development practice have shown that, trend and orientation of principal stress
Parallel crack is larger with aperture, connectivity is universal preferably, relatively strong etc. advantages of flow conductivity；And trend and maximum horizontal
The aperture in the vertical crack of principal direction of stress is generally smaller or even closure, and connectivity is poor, and flow conductivity is also relatively
It is weak.It is a kind of qualitative method according to fracture strike and the relationship assessment fisstured flow of stress direction, in view of fisstured flow
It is strong and weak also closely related with the factors such as present daygeodynamics size, pore pressure, therefore carry out the seepage flow energy that crack is quantitatively evaluated
The correlative studys such as power are that fractured reservoirs develop a most important urgently to be solved.
During oil and gas development, the variation of strata pressure can significantly change the crustal stress states of oil-gas reservoir, to make ground
The permeability in layer interrupting layer and crack changes.In natural petroleum gas field geological research field, geomechanics modeling and
Hydrocarbon reservoirs modeling is usually independent from each other two parts；In hydrocarbon reservoirs modeling, to the description master of intrinsic fracture
If characterizing intrinsic fracture by crack modeling technique in the distribution in space, however which does not account for crustal stress counterincision
Stitch the influence of seepage flow.
Geomechanics modeling mainly includes one-dimensional geology mechanical modeling and three-dimensional geological mechanical modeling.Three-dimensional geological mechanics is built
For mould based on various initial data, foundation is able to reflect Geological Structural Forms, tectonic relationship and the variation of geologic body built-in attribute
The digital model of rule；According to the different disposal to time dimension, three-dimensional geological mechanical modeling can be divided into three-dimensional static geology power
Model and Three-Dimensional Dynamic geomechanics model are learned, wherein the time is constant in static modelling, and model is only used for description object a certain
The static nature at moment, and dynamic modeling can be used for the consecutive variations of simulated object whithin a period of time.
Currently, three-dimensional geological mechanical modeling is also immature to the quantitative evaluation technology of intrinsic fracture percolation ability；It is one-dimensional
Geomechanics modeling can carry out quantitative evaluation by percolation ability of the data such as imaging logging to the intrinsic fracture in individual well,
However it is difficult to realize the mechanical analysis to three-dimensional space a large amount of crack and tomography, therefore more difficult qualitative assessment three-dimensional space is split
Stitch influence of the seepage field variation to oil and gas development.
Summary of the invention
The present invention provides a kind of method using ground stress analysis tomography and fisstured flow characteristic, can overcome existing infiltration
Properties of flow analysis method because lack crustal stress factor to can not the seepage characteristic accurately to tomography and crack analyze etc.
Defect.
The present invention provides a kind of method using ground stress analysis tomography and fisstured flow characteristic, includes the following steps:
1) the dynamic 3 D geomechanics model and discrete fracture network model (DFN model) in target work area are established respectively；
2) by six stress parameters σ in dynamic 3 D geomechanics modelx、σy、σz、τxy、τxzAnd τyzAccordingly it is loaded into
On the tomography and fracture surface of discrete fracture network model, wherein σx、σy、σzThe respectively direct stress in x, y, z direction, τxy、τxz、τyz
The respectively shear stress in the direction xy, xz, yz；
3) according to the tilt angle gamma in tomography and crack and tendency θ, three stress cosine coefficient l, m and n are obtained, wherein l, m, n
The respectively stress cosine coefficient in x, y, z direction；
4) according to three stress cosine coefficients and six stress parameters, three components of stress p are obtainedx、py、pz,
Wherein px、py、pzThe respectively components of stress in x, y, z direction；
5) according to three stress cosine coefficients and three components of stress, direct stress σ is obtainedN；
6) according to three components of stress and direct stress σN, obtain shear stress τN；
7) according to the direct stress σNWith shear stress τN, obtain a mole coulomb function CFF.
The method of the present invention application ground stress analysis tomography and fisstured flow characteristic, dynamic 3 D geomechanics modeling and
On the basis of the modeling of crack, quickly, in bulk by the crustal stress achievement assignment in dynamic 3 D geomechanics model to tomography and
On fracture surface, the mechanical state of quantitative analysis three-dimensional space interrupting layer different location and the mechanical state of millions fracture surface become
The shear stress of law, especially tomography and fracture surface and the relationship of direct stress are seeped using a mole coulomb mechanical equation with three-dimensional
The form in flow field describes the activity and fisstured flow Spatial Variation of tomography, realizes the seepage flow point that crack is quantitatively evaluated
Boot is sought peace changing rule, and analysis visual result is illustrated, is easily understood, accuracy is high, good reliability, to optimize Reservoir Development
Scheme and measure provide scientific basis and guidance, to be conducive to that oil-gas field development personnel is assisted to fully consider seepage field change
Optimal Development production decision and well stimulation on the basis of law realize the lasting of oil gas field, high yield and scientific development.
In step 1) of the invention, the conventional method that can commonly use this field establishes the dynamic 3 D geology in target work area
Mechanical model and discrete fracture network model.
In step 2) of the invention, by six stress parameters σ in dynamic 3 D geomechanics modelx、σy、σz、τxy、
τxzAnd τyzIt is accordingly loaded on the tomography and fracture surface of discrete fracture network model, it may be assumed that right respectively in corresponding all directions
The tomography and fracture surface of discrete fracture network model apply and six stress parameters σx、σy、σz、τxy、τxzAnd τyzIdentical stress,
To make dynamic 3 D geomechanics model and discrete fracture network Model coupling.
In step 3) of the invention, formula 1, formula 2 and formula 3 can be passed through respectively and obtain three stress cosine coefficients
L, m and n,
Formula 1:l=sin (γ) × cos (θ),
Formula 2:m=sin (γ) × sin (θ),
Formula 3:n=cos (γ),
Wherein, l2+m2+n2=1.
In step 4) of the invention, formula 4, formula 5 and formula 6 can be passed through respectively and obtain three components of stress px、
py、pz,
Formula 4:px=l σx+mτxy+nτxz,
Formula 5:py=l τxy+mσy+nτyz,
Formula 6:pz=l τxz+mτyz+nσz。
In step 5) of the invention, direct stress σ can be obtained by formula 7N,
Formula 7: σN=lpx+mpy+npz。
In step 6) of the invention, shear stress τ can be obtained by formula 8N,
Formula 8: τN 2=px 2+py 2+pz 2-σN 2。
In step 7) of the invention, a mole coulomb function CFF can be obtained by formula 9,
Formula 9:CFF=τN–μ(σN-pp),
Wherein, μ is conversion coefficient, ppFor Pore Pressure force value.
It is understood that mole coulomb function CFF can be used for assessing the seepage characteristic in tomography and crack；Specifically, it rubs
Your coulomb function CFF is smaller, shows that the closure degree in tomography and crack is bigger, i.e. the seepage flow of tomography and crack is poorer；Conversely, rubbing
Your coulomb function CFF is bigger, shows that the opening degree in tomography and crack is bigger, i.e., the seepage flow of tomography and crack is better.
Further, the method for the invention using ground stress analysis tomography and fisstured flow characteristic can also include step
It is rapid 8), according to described mole of coulomb function CFF, obtain a mole coulomb function variable quantity △ CFF.
In step 8) of the invention, a mole coulomb function variable quantity △ CFF can be obtained by formula 10,
10: △ CFF=CFF-CFF of formula0,
Wherein, CFF is current mole of coulomb function, CFF0For initial molar coulomb function.
Further, the method for the invention using ground stress analysis tomography and fisstured flow characteristic can also include step
It is rapid 9), according to a mole coulomb function variable quantity △ CFF, obtain the current permeability K in target work areai。
In step 9) of the invention, the current permeability K in target work area can be obtained by formula 11i,
Formula 11:Ki=K0×(1-a×△CFF)-1.092,
Wherein, KiFor current permeability, K0For original permeability, △ CFF is mole coulomb function variable quantity, and a is correction system
Number；Specifically, core experiment data be can use, the relationship of mechanics parameter CFF and crack core permeability power function is fitted, obtains
To the specific value of correction coefficient a.
With the seepage characteristic of permeability K reflection tomography and crack, analysis visual result is illustrated, is simple easy method of the invention
Understand；It is understood that current permeability KiIt is smaller, show that the closure degree in tomography and crack is bigger, i.e. tomography and crack
Seepage flow is poorer；Conversely, current permeability KiIt is bigger, show that the opening degree in tomography and crack is bigger, i.e. the infiltration of tomography and crack
Stream is better.
Method of the invention can analyze the mechanics shape of the different times of the crack piece of large scale tomography and ten million number of stages
State, to accurately assess the changing rule of the sealing characteristics of tomography and the seepage flow in crack；This method is according in geomechanics
Mole coulomb theory assesses the mechanics of tomography and crack by the relationship of the shear stress and direct stress of analysis tomography and fracture surface
(i.e. closure or unlatching, seepage flow is poor if closure for state；If opening and opening degree being bigger, seepage flow is better), accuracy
High, good reliability.
Method of the invention is suitable for tomography and fractured reservoirs, changes especially for strata pressure in development and production
Big and oil-gas reservoir crustal stress changes greatly that stability so as to cause tomography changes or the seepage flow in crack changes greatly
Oil gas field.
Further, the present invention can by it is above-mentioned using the method for ground stress analysis tomography and fisstured flow characteristic with software mould
The modes such as block are applied.
It specifically, can be according to the dynamic 3 D geomechanics model and discrete fracture network model in target work area, automatically
Generating Jewel Suite software (learning software platform to Baker Hughes Inc) identifiable script file, (script file is usual
For XML format file), it is run after which is imported into Jewel Suite software, tomography and fracture surface can be obtained
CFF value；Working principle is the above-mentioned method using ground stress analysis tomography and fisstured flow characteristic.
Implementation of the invention, at least has the advantage that
1, method of the invention, on the basis of the modeling of dynamic 3 D geomechanics and crack model, by dynamic 3 D
Matter mechanical model is coupled with discrete fracture network model, because lacking is answered to overcome existing Seepage characteristic analysis method
Power factor to can not the seepage characteristic accurately to tomography and crack analyze the defects of.
2, method of the invention can quickly, in bulk by the crustal stress achievement assignment in dynamic 3 D geomechanics model
Onto tomography and fracture surface, the quantitative analysis mechanical state of three-dimensional space interrupting layer different location and the power of millions fracture surface
State change rule is learned, and describes the activity and fisstured flow Spatial Variation of tomography in the form of three-dimensional seepage field,
Realize the seepage flow distribution characteristics and changing rule that crack is quantitatively evaluated.
3, method of the invention is illustrated with the seepage characteristic of permeability K reflection tomography and crack, analysis visual result, is simple
It is understandable, accuracy is high, good reliability, thus for optimization Reservoir Development scheme and measure provide scientific basis and guidance, have
Conducive to assistance oil-gas field development personnel Optimal Development production decision and volume increase on the basis of fully considering seepage field changing rule
Measure realizes the lasting of oil gas field, high yield and scientific development.
4, method of the invention can be applied in a manner of software module etc., according to the dynamic 3 D geology power in target work area
Model and discrete fracture network model are learned, the identifiable script file of Jewel Suite software is automatically generated, by the script file
It is run after importeding into Jewel Suite software, the CFF value of tomography and fracture surface can be obtained, so as to conveniently and efficiently right
Tomography and the seepage characteristic in crack are quantitatively evaluated.
5, method of the invention is suitable for tomography and fractured reservoirs, becomes especially suitable for strata pressure in development and production
Change that the stability that big and oil-gas reservoir crustal stress changes greatly so as to cause tomography changes or the seepage flow in crack changes greatly
Oil gas field.
Detailed description of the invention
The spatial distribution of crack penetration parameter CFF parameter field when Fig. 1 is target work area starting in the embodiment of the present invention 1
Figure；
Fig. 2 is the spatial distribution of the crack penetration parameter CFF parameter field in the embodiment of the present invention 1 at the end of target work area
Figure；
The spatial distribution of permeability of fault parameter CFF parameter field when Fig. 3 is target work area starting in the embodiment of the present invention 1
Figure；
Fig. 4 is the spatial distribution of the permeability of fault parameter CFF parameter field in the embodiment of the present invention 1 at the end of target work area
Figure；
The spatial distribution map of permeability K parameter field when Fig. 5 is target work area starting in the embodiment of the present invention 2；
Fig. 6 is the spatial distribution map of the permeability K parameter field in the embodiment of the present invention 2 at the end of target work area.
Specific embodiment
To make the object, technical solutions and advantages of the present invention clearer, below in conjunction with attached drawing of the invention and implementation
Example, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is the present invention
A part of the embodiment, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art are not having
Every other embodiment obtained under the premise of creative work is made, shall fall within the protection scope of the present invention.
Embodiment 1
The present embodiment is using Dina gas field as target work area, using the side of ground stress analysis tomography and fisstured flow characteristic
Method includes the following steps:
1) the dynamic 3 D geomechanics model and discrete fracture network model in target work area are established respectively.
Specifically, the dynamic 3 D geomechanics model in above-mentioned target work area can be established using the usual manner of this field
With discrete fracture network model；Wherein, the simulated object of dynamic 3 D geomechanics model was above-mentioned target work area in 2009
Consecutive variations in 2026.
2) by six stress parameters σ in dynamic 3 D geomechanics modelx、σy、σz、τxy、τxzAnd τyzAccordingly it is loaded into
On the tomography and fracture surface of discrete fracture network model, wherein σx、σy、σzThe respectively direct stress in x, y, z direction, τxy、τxz、τyz
The respectively shear stress in the direction xy, xz, yz.
Specifically, respectively to the tomography of the discrete fracture network model of above-mentioned foundation on x, y, z, the direction xy, xz and yz
With the corresponding application of fracture surface and above-mentioned six stress parameters σx、σy、σz、τxy、τxzAnd τyzIdentical stress, to make dynamic 3 D
Geomechanics model and discrete fracture network Model coupling.
3) according to the tilt angle gamma in tomography and crack and tendency θ, three stress cosine coefficient l, m and n are obtained, wherein l, m, n
The respectively stress cosine coefficient in x, y, z direction.
Specifically, three stress cosine coefficient l, m and n are obtained by formula 1, formula 2 and formula 3 respectively,
Formula 1:l=sin (γ) × cos (θ),
Formula 2:m=sin (γ) × sin (θ),
Formula 3:n=cos (γ),
Wherein, l2+m2+n2=1.
4) according to three stress cosine coefficients and six stress parameters, three components of stress p are obtainedx、py、pz,
Wherein px、py、pzThe respectively components of stress in x, y, z direction.
Specifically, three components of stress p are obtained by formula 4, formula 5 and formula 6 respectivelyx、py、pz,
Formula 4:px=l σx+mτxy+nτxz,
Formula 5:py=l τxy+mσy+nτyz,
Formula 6:pz=l τxz+mτyz+nσz。
5) according to three stress cosine coefficients and three components of stress, direct stress σ is obtainedN。
Specifically, direct stress σ is obtained by formula 7N,
Formula 7: σN=lpx+mpy+npz。
6) according to three components of stress and direct stress σN, obtain shear stress τN。
Specifically, shear stress τ is obtained by formula 8N,
Formula 8: τN 2=px 2+py 2+pz 2-σN 2。
7) according to the direct stress σNWith shear stress τN, obtain a mole coulomb function CFF.
Specifically, a mole coulomb function CFF is obtained by formula 9,
Formula 9:CFF=τN–μ(σN-pp),
Wherein, μ is conversion coefficient, ppFor Pore Pressure force value.
The spatial distribution of crack penetration parameter CFF parameter field when Fig. 1 and Fig. 2 is respectively target work area starting and ending
Figure；The spatial distribution map of permeability of fault parameter CFF parameter field when Fig. 3 and Fig. 4 is respectively target work area starting and ending.Its
In, mole coulomb function CFF is smaller, shows that the closure degree in tomography and crack is bigger, i.e. the seepage flow of tomography and crack is poorer；Instead
It, mole coulomb function CFF is bigger, shows that the opening degree in tomography and crack is bigger, i.e., the seepage flow of tomography and crack is better.
The result shows that: Dina gas field construction top crack mechanical state change greatly, CFF value gradually becomes smaller, crack by
Gradually in closure, and the variation of the crack piece mechanical state of edge and alar part is smaller.
Embodiment 2
The present embodiment is using Dina gas field as target work area, using the side of ground stress analysis tomography and fisstured flow characteristic
Method includes the following steps:
1) the dynamic 3 D geomechanics model and discrete fracture network model in target work area are established respectively.
Specifically, the dynamic 3 D geomechanics model in above-mentioned target work area can be established using the usual manner of this field
With discrete fracture network model；Wherein, the simulated object of dynamic 3 D geomechanics model was above-mentioned target work area in 2009
Consecutive variations in 2026.
2) by six stress parameters σ in dynamic 3 D geomechanics modelx、σy、σz、τxy、τxzAnd τyzAccordingly it is loaded into
On the tomography and fracture surface of discrete fracture network model, wherein σx、σy、σzThe respectively direct stress in x, y, z direction, τxy、τxz、τyz
The respectively shear stress in the direction xy, xz, yz.
Specifically, respectively to the tomography of the discrete fracture network model of above-mentioned foundation on x, y, z, the direction xy, xz and yz
With the corresponding application of fracture surface and above-mentioned six stress parameters σx、σy、σz、τxy、τxzAnd τyzIdentical stress, to make dynamic 3 D
Geomechanics model and discrete fracture network Model coupling.
3) according to the tilt angle gamma in tomography and crack and tendency θ, three stress cosine coefficient l, m and n are obtained, wherein l, m, n
The respectively stress cosine coefficient in x, y, z direction.
Specifically, three stress cosine coefficient l, m and n are obtained by formula 1, formula 2 and formula 3 respectively,
Formula 1:l=sin (γ) × cos (θ),
Formula 2:m=sin (γ) × sin (θ),
Formula 3:n=cos (γ),
Wherein, l2+m2+n2=1.
4) according to three stress cosine coefficients and six stress parameters, three components of stress p are obtainedx、py、pz,
Wherein px、py、pzThe respectively components of stress in x, y, z direction.
Specifically, three components of stress p are obtained by formula 4, formula 5 and formula 6 respectivelyx、py、pz,
Formula 4:px=l σx+mτxy+nτxz,
Formula 5:py=l τxy+mσy+nτyz,
Formula 6:pz=l τxz+mτyz+nσz。
5) according to three stress cosine coefficients and three components of stress, direct stress σ is obtainedN。
Specifically, direct stress σ is obtained by formula 7N,
Formula 7: σN=lpx+mpy+npz。
6) according to three components of stress and direct stress σN, obtain shear stress τN。
Specifically, shear stress τ is obtained by formula 8N,
Formula 8: τN 2=px 2+py 2+pz 2-σN 2。
7) according to the direct stress σNWith shear stress τN, obtain a mole coulomb function CFF.
Specifically, a mole coulomb function CFF is obtained by formula 9,
Formula 9:CFF=τN–μ(σN-pp),
Wherein, μ is conversion coefficient, ppFor Pore Pressure force value.
8) according to mole coulomb function CFF, a mole coulomb function variable quantity △ CFF is obtained.
Specifically, a mole coulomb function variable quantity △ CFF is obtained by formula 10,
10: △ CFF=CFF-CFF of formula0,
Wherein, CFF is current mole of coulomb function, CFF0For initial molar coulomb function.
9) according to mole coulomb function variable quantity △ CFF, the current permeability K in target work area is obtainedi。
Specifically, the current permeability K in target work area is obtained by formula 11i,
Formula 11:Ki=K0×(1-a×△CFF)-1.092,
Wherein, KiFor current permeability, K0For original permeability, △ CFF is mole coulomb function variable quantity, and a is correction system
Number.
The spatial distribution map of permeability K parameter field when Fig. 5 and Fig. 6 is respectively target work area starting and ending；Wherein, when
Preceding permeability KiIt is smaller, show that the closure degree in tomography and crack is bigger, i.e. the seepage flow of tomography and crack is poorer；Conversely, current
Permeability KiIt is bigger, show that the opening degree in tomography and crack is bigger, i.e., the seepage flow of tomography and crack is better.
The result shows that: the simulation permeability variation at Dina gas field construction top is larger, and edge and alar part variation are smaller.
Finally, it should be noted that the above embodiments are only used to illustrate the technical solution of the present invention., rather than its limitations；To the greatest extent
Pipe present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: its according to
So be possible to modify the technical solutions described in the foregoing embodiments, or to some or all of the technical features into
Row equivalent replacement；And these are modified or replaceed, various embodiments of the present invention technology that it does not separate the essence of the corresponding technical solution
The range of scheme.
Claims (10)
1. a kind of method using ground stress analysis tomography and fisstured flow characteristic, which comprises the steps of:
1) the dynamic 3 D geomechanics model and discrete fracture network model in target work area are established respectively；
2) by six stress parameters σ in dynamic 3 D geomechanics modelx、σy、σz、τxy、τxzAnd τyzIt is accordingly loaded into discrete
On the tomography and fracture surface of fracture network model of analogous outcrop, wherein σx、σy、σzThe respectively direct stress in x, y, z direction, τxy、τxz、τyzRespectively
For the shear stress in the direction xy, xz, yz；
3) according to the tilt angle gamma in tomography and crack and tendency θ, three stress cosine coefficient l, m and n are obtained, wherein l, m, n distinguish
For the stress cosine coefficient in x, y, z direction；
4) according to three stress cosine coefficients and six stress parameters, three components of stress p are obtainedx、py、pz, wherein
px、py、pzThe respectively components of stress in x, y, z direction；
5) according to three stress cosine coefficients and three components of stress, direct stress σ is obtainedN；
6) according to three components of stress and direct stress σN, obtain shear stress τN；
7) according to the direct stress σNWith shear stress τN, obtain a mole coulomb function CFF.
2. the method according to claim 1, wherein further including step 8), according to described mole of coulomb function
CFF obtains a mole coulomb function variable quantity △ CFF.
3. according to the method described in claim 2, it is characterized in that, further including step 9), according to a mole coulomb function variable quantity
△ CFF obtains the current permeability K in target work areai。
4. the method according to claim 1, wherein passing through formula 1, formula 2 and formula 3 respectively in step 3)
Three stress cosine coefficient l, m and n are obtained,
Formula 1:l=sin (γ) × cos (θ),
Formula 2:m=sin (γ) × sin (θ),
Formula 3:n=cos (γ),
Wherein, l2+m2+n2=1.
5. the method according to claim 1, wherein passing through formula 4, formula 5 and formula 6 respectively in step 4)
Obtain three components of stress px、py、pz,
Formula 4:px=l σx+mτxy+nτxz,
Formula 5:py=l τxy+mσy+nτyz,
Formula 6:pz=l τxz+mτyz+nσz。
6. the method according to claim 1, wherein obtaining direct stress σ by formula 7 in step 5)N,
Formula 7: σN=lpx+mpy+npz。
7. the method according to claim 1, wherein obtaining shear stress τ by formula 8 in step 6)N,
Formula 8: τN 2=px 2+py 2+pz 2-σN 2。
8. the method according to claim 1, wherein obtaining a mole coulomb function by formula 9 in step 7)
CFF,
Formula 9:CFF=τN–μ(σN-pp),
Wherein, μ is conversion coefficient, ppFor Pore Pressure force value.
9. according to the method described in claim 2, it is characterized in that, obtaining a mole coulomb function by formula 10 in step 8)
Variable quantity △ CFF,
10: △ CFF=CFF-CFF of formula0,
Wherein, CFF is current mole of coulomb function, CFF0For initial molar coulomb function.
10. according to the method described in claim 3, it is characterized in that, obtaining working as target work area by formula 11 in step 9)
Preceding permeability Ki,
Formula 11:Ki=K0×(1-a×△CFF)-1.092,
Wherein, KiFor current permeability, K0For original permeability, △ CFF is mole coulomb function variable quantity, and a is correction coefficient.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710818821.XA CN109492238B (en) | 2017-09-12 | 2017-09-12 | Method for analyzing fault and crack seepage characteristics by applying ground stress |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710818821.XA CN109492238B (en) | 2017-09-12 | 2017-09-12 | Method for analyzing fault and crack seepage characteristics by applying ground stress |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109492238A true CN109492238A (en) | 2019-03-19 |
| CN109492238B CN109492238B (en) | 2022-03-01 |
Family
ID=65688683
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710818821.XA Active CN109492238B (en) | 2017-09-12 | 2017-09-12 | Method for analyzing fault and crack seepage characteristics by applying ground stress |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109492238B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114611339A (en) * | 2022-05-12 | 2022-06-10 | 中国石油大学(华东) | Crack connectivity identification method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104200039A (en) * | 2014-09-17 | 2014-12-10 | 中国石油大学(华东) | Quantitative forecasting method of tectonic fissure occurrence |
| CN105319603A (en) * | 2015-11-06 | 2016-02-10 | 中国石油大学(华东) | Compact sandstone reservoir complex netted fracture prediction method |
| CN105334536A (en) * | 2015-12-01 | 2016-02-17 | 中国石油大学(华东) | Effectiveness evaluation method for compact sandstone reservoir map cracking system |
| US20160245939A1 (en) * | 2013-10-21 | 2016-08-25 | Westerngeco Llc | Seismic data analysis |
-
2017
- 2017-09-12 CN CN201710818821.XA patent/CN109492238B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160245939A1 (en) * | 2013-10-21 | 2016-08-25 | Westerngeco Llc | Seismic data analysis |
| CN104200039A (en) * | 2014-09-17 | 2014-12-10 | 中国石油大学(华东) | Quantitative forecasting method of tectonic fissure occurrence |
| CN105319603A (en) * | 2015-11-06 | 2016-02-10 | 中国石油大学(华东) | Compact sandstone reservoir complex netted fracture prediction method |
| CN105334536A (en) * | 2015-12-01 | 2016-02-17 | 中国石油大学(华东) | Effectiveness evaluation method for compact sandstone reservoir map cracking system |
Non-Patent Citations (1)
| Title |
|---|
| 高英: ""薄互层低渗透油藏压裂开发裂缝扩展规律及产能预测研究"", 《中国博士学位论文全文数据库》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114611339A (en) * | 2022-05-12 | 2022-06-10 | 中国石油大学(华东) | Crack connectivity identification method |
| CN114611339B (en) * | 2022-05-12 | 2022-07-15 | 中国石油大学(华东) | Crack connectivity identification method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109492238B (en) | 2022-03-01 |
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 | |
| CA2895549C (en) | Fracturing and reactivated fracture volumes | |
| US9835746B2 (en) | Formation stability modeling | |
| Wang et al. | Numerical analysis on scale effect of elasticity, strength and failure patterns of jointed rock masses | |
| Close et al. | Integrated workflows for shale gas and case study results for the Horn River Basin, British Columbia, Canada | |
| US10386531B2 (en) | Geological model analysis incorporating cross-well electromagnetic measurements | |
| US10048403B2 (en) | Method and system for generation of upscaled mechanical stratigraphy from petrophysical measurements | |
| US9910173B2 (en) | Saturation end-point adjustment | |
| CA2888249A1 (en) | System and method for characterizing uncertainty in subterranean reservoir fracture networks | |
| NO342764B1 (en) | Model compatible structure reconstruction for geomechanical and petroleum system modeling | |
| US11789170B2 (en) | Induced seismicity | |
| WO2021048597A1 (en) | An integrated geomechanics model for predicting hydrocarbon and migration pathways | |
| Liu et al. | Asymmetric propagation mechanism of hydraulic fracture networks in continental reservoirs | |
| EP3303763B1 (en) | Oilfield reservoir saturation and permeability modeling | |
| CN105089615A (en) | Log data historical retrogression treatment method based on oil reservoir model | |
| Chang et al. | Data assimilation of coupled fluid flow and geomechanics using the ensemble Kalman filter | |
| WO2021029875A1 (en) | Integrated rock mechanics laboratory for predicting stress-strain behavior | |
| CN106355571A (en) | Determination method and device for dolostone reservoir quality | |
| SA120420296B1 (en) | Cluster Efficiency Operation Control | |
| CN113296166A (en) | Method for constructing crack model | |
| US20220082719A1 (en) | Reservoir performance system | |
| US11578564B2 (en) | Systems and methods for predicting shear failure of a rock formation | |
| Williams-Stroud | Using microseismic events to constrain fracture network models and implications for generating fracture flow properties for reservoir simulation | |
| CN109492238A (en) | A method of using ground stress analysis tomography and fisstured flow characteristic | |
| CN107942381B (en) | Quantitative prediction method for tight oil reservoir bedding joints |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |