CN103913772B - Micro-seismic event the Forward Modeling based on reservoir geology mechanics parameter - Google Patents

Abstract

The present invention provides a kind of micro-seismic event the Forward Modeling based on reservoir geology mechanics parameter, should micro-seismic event the Forward Modeling based on reservoir geology mechanics parameter include: for the assigned pressure of initial fracturing, use finite element method (FEM) to carry out modeling；The rock unit Stress calculation of each particular moment is carried out according to time change；According to stress calculating results, in conjunction with criterion in rock, it is judged that whether current rock unit ruptures；Calculate porosity parameter；If rock unit ruptures, calculate focus earthquake magnitude in conjunction with microseism focal mechanism；And complete the simulation in all moment.It is somebody's turn to do micro-seismic event the Forward Modeling result of calculation based on reservoir geology mechanics parameter and can reflect the actual change situation of subsurface reservoir really, the prediction of microseismic event in subterranean fracture and fracturing process can be carried out by the method.

Description

Micro-seismic event the Forward Modeling based on reservoir geology mechanics parameter

Technical field

The microseism simulation that the present invention relates in petroleum works in fracturing simulation and method of seismic prospecting, Especially relate to a kind of micro-seismic event the Forward Modeling based on reservoir geology mechanics parameter.

Background technology

Fracturing is the important means improving the fine and close unconventional petroleum resources yield such as oil gas, shale gas. During pressure break, monitor fracturing effect, Fracturing Project is proposed perfect, for improving yield, fall Low cost has great significance.Microseismic is the most approved for monitoring pressure break in real time Unique effective means of effect.During low-permeability reservoir water filling or gas injection, fluid can be caused The movement of pressure leading edge and the change of pore fluid pressure, thus cause micro-seismic event.Pass through microseism Monitoring, can adjust, optimizing design scheme in real time, is effectively improved the rate of oil and gas recovery, it is achieved oil field The scientific management of whole development.By the micro-seismic monitoring under fracturing, can to pressure break scope, Fracture azimuth and size are effectively followed the trail of, location, the effect of objective appraisal fracturing engineering, Particularly by series crack attributive analysis, the structure trend of location major fracture, and secondary fracture point Cloth, approximate calculation goes out the length in crack, width, highly；By micro-seismic monitoring, can be in rock Portion's fluid front carries out imaging.Reservoir engineer is by fracture imaging and drives leading edge to involve dividing of situation Analysis, provides effective guidance to next step production development.For primary development, crack accurately Orientation and distribution will assist in develops mobile oil gas from well most effectively.The secondary of reservoir is opened For Faing, aid in determining whether oil recovery well location, improve tar productivity, cost-effective, increase benefit.

The mechanism produced for microseism at present is understood not so that microseism data is imitated in actual applications Fruit needs to be further improved.Conventional microseism Forward technology assumes that focus and propagation medium it is known that grind Study carefully when it is walked and the method for the propagation law such as amplitude and frequency, be not bound with reservoir geology mechanics parameter. A kind of new micro-seismic event forward simulation side based on reservoir geology mechanics parameter is we have invented for this Method, solves above technical problem.

Summary of the invention

It is an object of the invention to provide and a kind of can carry out microseismic event in subterranean fracture and fracturing process The micro-seismic event the Forward Modeling based on reservoir geology mechanics parameter of prediction.

The purpose of the present invention can be achieved by the following technical measures: based on reservoir geology mechanics parameter micro- Seismic events the Forward Modeling, should micro-seismic event forward simulation side based on reservoir geology mechanics parameter Method includes: step 1, for the assigned pressure of initial fracturing, uses finite element method (FEM) to carry out pattern die Intend；Step 2, carries out the rock unit Stress calculation of each particular moment according to time change；Step 3, According to stress calculating results, in conjunction with criterion in rock, it is judged that whether current rock unit ruptures； Step 4, calculates porosity parameter；Step 5, if rock unit ruptures, in conjunction with microseism focal mechanism Calculate focus earthquake magnitude；And step 6, repeat step 2, step 3, step 4 and step 5, complete to own The simulation in moment.

The purpose of the present invention realizes also by following technical measures:

Should micro-seismic event the Forward Modeling based on reservoir geology mechanics parameter also include, in step 1 Before, set up the seepage-stress-damage coupling model of subsurface reservoir, compose reservoir geology mechanics ginseng for model Number.

In the step of seepage-stress-damage coupling model setting up subsurface reservoir, this seepage-stress-damage The reservoir parameter governing equation hindering coupling model employing is:

${\stackrel{\·}{\mathrm{\σ}}}_{ij}=2G{\stackrel{\·}{\mathrm{\ϵ}}}_{ij}+(K-\frac{2G}{3}){\stackrel{\·}{\mathrm{\ϵ}}}_{kk}{\mathrm{\δ}}_{ij}+\mathrm{\α}\stackrel{\·}{p}{\mathrm{\δ}}_{ij}+{\mathrm{\γ}}_1\stackrel{\·}{T}{\mathrm{\δ}}_{ij}---\left(1\right)$

$\stackrel{\·}{\mathrm{\ζ}}=\mathrm{\α}{\stackrel{\·}{\mathrm{\ϵ}}}_{ii}+\mathrm{\β}\stackrel{\·}{p}-{\mathrm{\γ}}_2\stackrel{\·}{T}---\left(2\right)$

σ in formula_ijAnd ε_ijRepresenting total stress and overall strain respectively, p and T represents pore pressure and temperature, α respectively Representing Biot coefficient, ζ represents saturation of pore fluid, δ_ijRepresenting Kronecker symbol, K represents volume mould Amount, G represents modulus of shearing, γ₁, γ₂, the value of β is given by formula below:

$\mathrm{\β}=\frac{\mathrm{\α}-\mathrm{\φ}}{K_s}+\frac{\mathrm{\φ}}{K_f}---\left(3\right)$

γ₁=K α_m (4)

γ₂=α α_m+(α_f-α_m)φ (5)

φ represents porosity, α_mAnd α_fRepresent rock and the thermal coefficient of expansion of fluid, the volume of rock respectively Modulus K_sRepresenting, the bulk modulus of fluid uses K_fRepresent；Assume that fluid flows in blowhole Following Darcy's law, conduction of heat follows Fourier law；Then there is a following formula:

$J^f=-{\mathrm{\ρ}}_f\frac{k}{\mathrm{\η}}\▿p---\left(6\right)$

$J^T=-k^T\▿T---\left(7\right)$

ρ_fRepresenting fluid density, k represents that permeability, η represent the coefficient of viscosity, k^TRepresent the coefficient of heat conduction, J^f And J^TRepresent fluid flow and heat respectively；Following formula (8) (9) is that Fluid in Pore flowing is put down with stress Weighing apparatus equation；

σ_ij,j=0 (8)

$\frac{\∂\mathrm{\ζ}}{\∂t}=-\frac{1}{{\mathrm{\ρ}}_f}\▿J^f---\left(9\right)$

Bring governing equation (1) (2) above into equilibrium equation (8) (9), obtain about rock deformation Field equation (10) (11) with fluid flow state；

$(K+\frac{G}{3})\▿(\▿.u)+G{\▿}^{2}u+m(\mathrm{\α}\▿p+{\mathrm{\γ}}_1\▿T)=0---\left(10\right)$

$\mathrm{\α}(\▿.\stackrel{\·}{u})+\mathrm{\β}\stackrel{\·}{p}-\frac{k}{\mathrm{\η}}{\▿}^{2}p-{\mathrm{\γ}}_2\stackrel{\·}{T}=0---\left(11\right)$

$\stackrel{\·}{T}+v(\▿T)-c^T{\▿}^{2}T=0---\left(12\right)$

U is displacement, the m=[1,1,0] when model is two dimension^T, m=[1,1,1,0,0,0] when model is three-dimensional^T, c^TFor thermal diffusion coefficient, speed and the pore pressure of fluid follow Darcy's law

In step 3, in a model, reservoir rock can be described as elastic shape by the process after stress effect Become and rupture two stages；During by less stress, rock generation elastic deformation, this stage rock interior Do not rupture；And when rock stress increases the failure criteria reaching rock, rock starts to rupture, with The development ruptured, ultimately form crack.

In step 3, in elastic deformation-rupture mechanism, the reduction of elastic modelling quantity represents rock and occurs broken The process split,

E=(1-d) E₀ (13)

D represents the degree that rock interior ruptures, and is properly termed as rock fracture coefficient, E₀Represent rock Initial lithology modulus, E is current lithology modulus；If rock is stressed and there occurs and rupture, rock breaks Split coefficient to be stated by the residual compression of rock:

$d=1-\frac{f_{cr}}{E_0\stackrel{\‾}{\mathrm{\ϵ}}}---\left(14\right)$

f_crFor the residual pressure of rock,For under corresponding pressure occur strain, thenIllustrate initial pressure Stress.

In step 4, it is thus achieved that the actual loading situation of a certain each underground unit of moment, and contrast reality Stress size and rock unit keep not rupturing born maximum stress, are made whether the judgement ruptured, If ruptured, then rock unit deforms upon, and porosity parameter changes according to deformation；As Really rock rupture, then porosity is undergone mutation, and is converted into a bigger numerical value, and this numerical value is constant.

In step 5, if it occur that rupture, simulate the intensity of microseismic event according to the situation of change of stress, Microseismic event hypocentral location is current computing unit position, and the vibrations moment is current time.

Should micro-seismic event the Forward Modeling based on reservoir geology mechanics parameter also include, in step 5 Afterwards, all analog results are exported, as the foundation analyzing and explaining reservoir pressure.

Micro-seismic event the Forward Modeling based on reservoir geology mechanics parameter in the present invention, with routine Microseism simulation algorithm compare, it is contemplated that the change of reservoir geology mechanics parameter in fracturing process, therefore Result of calculation can reflect the actual change situation of subsurface reservoir really, can carry out ground by the method The prediction of microseismic event in lower crack and fracturing process.Just drilling is in simulation fracturing implementation process, storage Layer is by the situation after STRESS VARIATION.When using this technology can simulate practice of construction operation, subsurface reservoir The information such as the hypocentral location of microseismic event, source level in the generation in crack, growth and fracturing process. These information have great significance for instructing fracturing work.

Accompanying drawing explanation

Fig. 1 is the one of the micro-seismic event the Forward Modeling based on reservoir geology mechanics parameter of the present invention The flow chart of specific embodiment；

Fig. 2 be the present invention a specific embodiment in model elastic modelling quantity scattergram；

Fig. 3 be the present invention a specific embodiment in model stress envelope primitively；

Fig. 4 be the present invention a specific embodiment in forward simulation obtain model microseismic event schematic diagram.

Detailed description of the invention

For making the above and other objects, features and advantages of the present invention to become apparent, cited below particularly go out Preferred embodiment, and coordinate institute's accompanying drawings, it is described in detail below.

As it is shown in figure 1, the micro-seismic event based on reservoir geology mechanics parameter that Fig. 1 is the present invention is just drilling mould The flow chart of one specific embodiment of plan method.

In step 101, set up the seepage-stress-damage coupling model of subsurface reservoir, compose reservoir for model Geomechanics parameter.

Must be in conjunction with fluid pressure to fracturing simulation, rock is stressed deformation and rock failure process Describe.The generation of subterranean fracture is in the nature rock rupture, and it produces and development is mainly by opening that rock is subject to Power affects.The permeability of reservoir rock and stress are the functions of a coupling, when, after rock rupture, permeating Rate sharply increases；

The reservoir parameter governing equation such as following formula that model of the present invention uses:

${\stackrel{\·}{\mathrm{\σ}}}_{ij}=2G{\stackrel{\·}{\mathrm{\ϵ}}}_{ij}+(K-\frac{2G}{3}){\stackrel{\·}{\mathrm{\ϵ}}}_{kk}{\mathrm{\δ}}_{ij}+\mathrm{\α}\stackrel{\·}{p}{\mathrm{\δ}}_{ij}+{\mathrm{\γ}}_1\stackrel{\·}{T}{\mathrm{\δ}}_{ij}---\left(1\right)$

$\stackrel{\·}{\mathrm{\ζ}}=\mathrm{\α}{\stackrel{\·}{\mathrm{\ϵ}}}_{ii}+\mathrm{\β}\stackrel{\·}{p}-{\mathrm{\γ}}_2\stackrel{\·}{T}---\left(2\right)$

σ in formula_ijAnd ε_ijRepresenting total stress and overall strain respectively, p and T represents pore pressure and temperature, α respectively Representing Biot coefficient, ζ represents saturation of pore fluid, δ_ijRepresenting Kronecker symbol, K represents volume mould Amount, G represents modulus of shearing, γ₁, γ₂, the value of β is given by formula below:

$\mathrm{\β}=\frac{\mathrm{\α}-\mathrm{\φ}}{K_s}+\frac{\mathrm{\φ}}{K_f}---\left(3\right)$

γ₁=K α_m (4)

γ₂=α α_m+(α_f-α_m)φ (5)

φ represents porosity, α_mAnd α_fRepresent rock and the thermal coefficient of expansion of fluid, the volume of rock respectively Modulus K_sRepresenting, the bulk modulus of fluid uses K_fRepresent.Assume that fluid flows in blowhole Following Darcy's law, conduction of heat follows Fourier law.Then there is a following formula:

$J^f=-{\mathrm{\ρ}}_f\frac{k}{\mathrm{\η}}\▿p---\left(6\right)$

$J^T=-k^T\▿T---\left(7\right)$

ρ_fRepresenting fluid density, k represents that permeability, η represent the coefficient of viscosity, k^TRepresent the coefficient of heat conduction, J^f And J^TRepresent fluid flow and heat respectively.Following formula (8) (9) is that Fluid in Pore flowing is put down with stress Weighing apparatus equation.

σ_ij,j=0 (8)

$\frac{\∂\mathrm{\ζ}}{\∂t}=-\frac{1}{{\mathrm{\ρ}}_f}\▿J^f---\left(9\right)$

Bring governing equation (1) (2) above into equilibrium equation (8) (9), obtain about rock deformation Field equation (10) (11) with fluid flow state.

$(K+\frac{G}{3})\▿(\▿.u)+G{\▿}^{2}u+m(\mathrm{\α}\▿p+{\mathrm{\γ}}_1\▿T)=0---\left(10\right)$

$\mathrm{\α}(\▿.\stackrel{\·}{u})+\mathrm{\β}\stackrel{\·}{p}-\frac{k}{\mathrm{\η}}{\▿}^{2}p-{\mathrm{\γ}}_2\stackrel{\·}{T}=0---\left(11\right)$

$\stackrel{\·}{T}+v(\▿T)-c^T{\▿}^{2}T=0---\left(12\right)$

U is displacement, the m=[1,1,0] when model is two dimension^T, m=[1,1,1,0,0,0] when model is three-dimensional^T, c^TFor thermal diffusion coefficient, speed and the pore pressure of fluid follow Darcy's law

Flow process enters into step 102.

In step 102, for the assigned pressure of initial fracturing, reservoir geology mechanics parameter is composed initial value. Use equation (10)-(12) that individual grid cell is described, use finite element method (FEM) to carry out model Simulation.Flow process enters into step 103.

In step 103, carry out the rock unit Stress calculation of each particular moment according to time change.Stream Journey enters into step 104.

In step 104, according to stress calculating results, in conjunction with criterion in rock, it is judged that current rock list Whether unit ruptures.In a model, reservoir rock can be described as elasticity by the process after stress effect Deformation and rupture two stages.During by less stress, rock generation elastic deformation, in this stage rock Portion does not ruptures.And when rock stress increases the failure criteria reaching rock, rock starts to rupture, Along with the development ruptured, ultimately form crack.In elastic deformation-rupture mechanism, elastic modelling quantity Reduction represent the process that rock ruptures.

E=(1-d) E₀ (13)

D represents the degree that rock interior ruptures, and is properly termed as rock fracture coefficient, E₀Represent rock Initial lithology modulus, E is current lithology modulus.If rock is stressed and there occurs and rupture, rock breaks Split coefficient to be stated by the residual compression of rock:

$d=1-\frac{f_{cr}}{E_0\stackrel{\‾}{\mathrm{\ϵ}}}---\left(14\right)$

f_crFor the residual pressure of rock,For under corresponding pressure occur strain, thenIllustrate initial Compressive stress.

When carrying out reservoir geology mechanics parameter and describing, the initial pressure of fluid, and the initial geology of reservoir Mechanics parameter such as Poisson's ratio, elastic modelling quantity, initial field stress etc. are all as known conditions.Needs solve After being fluid pressure change, reservoir internal stress and the situation of change of strain.Use finite element method (FEM) solves Equation (10) (11), and during solving, add mole-coulomb rock rupture judgment criterion, just The response of reservoir under the conditions of certain pressure break can be obtained.Flow process enters into step 105.

In step 105, calculate, it is thus achieved that the actual loading situation of a certain each underground unit of moment, And contrast actual loading size and rock unit keeps not rupturing born maximum stress, be made whether The judgement ruptured, if ruptured, then rock unit deforms upon, and porosity parameter is according to deformation Change.If rock rupture, then porosity is undergone mutation, and being converted into a bigger numerical value (should Numerical value is constant).Flow process enters into step 106.

In step 106, if rock unit ruptures, finally obtain microseism thing in conjunction with focal mechanism Part information.If it occur that rupture, simulate the intensity of microseismic event, microseism according to the situation of change of stress Event hypocentral location is current computing unit position, and the vibrations moment is current time.Flow process enters into step 107。

In step 107, repeat step 103 and arrive step 106, complete the simulation in all moment.Flow process enters To step 108.

In step 108, all analog results are exported, as the foundation analyzing and explaining reservoir pressure. Flow process terminates.

For verification algorithm, establish one 500 meters wide, 500 meters of long two dimensional models, divide Minimum grid unit is the square of the length of side 5 meters.Geomechanics parameter employs certain work area known parameters. Fig. 2 is the elastic modelling quantity scattergram of model, and Fig. 3 is the scattergram of underground stress.This model is used and sends out Bright algorithm is simulated.Obtain the microseismic event distribution results figure shown in Fig. 4.

During the present invention provides a kind of oil gas field fracturing, microseism thing based on reservoir geology mechanics parameter Part analogy method, be mainly used in produce during fracturing the analysis of micro-seismic event feature and predict into The change of one step pressure break subsurface reservoir and the micro-seismic event information etc. that will generate.The method main flow It it is the shape carrying out describing reservoir by setting up the reservoir geology mechanics parameter model of seepage-stress-damage coupling State, by solving model formulation equation under boundary condition, thus realizes hydraulic fracturing process The numerical simulation of subsurface reservoir change, simulates the change by obtaining subsurface reservoir parameter, including geomechanics Parameter change, subterranean fracture are distributed, and combine microseism focal mechanism principle, it is thus achieved that in fracturing process The information such as position, time and the earthquake magnitude that micro-seismic event produces, are just drilling the micro-seismic event information of acquisition, For processing and explaining the reality produced and gather to monitor subterranean fracture to grow in hydraulic fracturing process Microseism data is significant.In order to use the stress of math equation describing reservoir and rupture feelings Condition, method do it is assumed hereinafter that: reservoir has elasticity and fragility simultaneously, and its state uses elastic damage Mechanics is described；Direct stress and shearing stress that reservoir rock is subject to follow Mohr-Coulomb fracture criteria； The permeability of reservoir is the elastically-deformable function that stress causes, and ruptures or tomography generation at rock Time, permeability will sharply increase.Under these assumed conditions, it is established that the mathematics meeting research requirement is retouched State equation.Under given boundary condition, it is possible to according to equation solution reservoir fracture situation, understand and Grasp the distribution of subterranean fracture.

Claims (6)

1. micro-seismic event the Forward Modeling based on reservoir geology mechanics parameter, it is characterised in that should Micro-seismic event the Forward Modeling based on reservoir geology mechanics parameter includes:

Set up the seepage-stress-damage coupling model of subsurface reservoir, compose reservoir geology mechanics ginseng for model Number；The reservoir geology mechanics parameter governing equation that this seepage-stress-damage coupling model uses is:

${\stackrel{\·}{\mathrm{\σ}}}_{ij}=2G{\stackrel{\·}{\mathrm{\ϵ}}}_{ij}+(K-\frac{2G}{3}){\stackrel{\·}{\mathrm{\ϵ}}}_{kk}{\mathrm{\δ}}_{ij}+\mathrm{\α}\stackrel{\·}{p}{\mathrm{\δ}}_{ij}+{\mathrm{\γ}}_1\stackrel{\·}{T}{\mathrm{\δ}}_{ij}---\left(1\right)$

$\stackrel{\·}{\mathrm{\ξ}}=\mathrm{\α}{\stackrel{\·}{\mathrm{\ϵ}}}_{ii}+\mathrm{\β}\stackrel{\·}{p}-{\mathrm{\γ}}_2\stackrel{\·}{T}---\left(2\right)$

σ in formula_ijRepresent total stress, ε_ijRepresent overall strain, ε_kkAnd ε_iiRepresent the strain of underground grid, p and T Representing pore pressure and temperature respectively, α represents Biot coefficient, and ζ represents saturation of pore fluid, δ_ijRepresent Kronecker symbol, K represents that bulk modulus, G represent modulus of shearing, γ₁, γ₂, under the value of β is passed through Face formula is given:

$\mathrm{\β}=\frac{\mathrm{\α}-\mathrm{\φ}}{K_s}+\frac{\mathrm{\φ}}{K_f}---\left(3\right)$

γ₁=K α_m (4)

γ₂=α α_m+(α_f-α_m)φ (5)

φ represents porosity, α_mAnd α_fRepresent rock and the thermal coefficient of expansion of fluid, the volume of rock respectively Modulus K_sRepresenting, the bulk modulus of fluid uses K_fRepresent；Assume that fluid flows in blowhole Following Darcy's law, conduction of heat follows Fourier law；Then there is a following formula:

$J^f=-{\mathrm{\ρ}}_f\frac{k}{\mathrm{\η}}\▿p---\left(6\right)$

$J^T=-k^T\▿T---\left(7\right)$

ρ_fRepresenting fluid density, k represents that permeability, η represent the coefficient of viscosity, k^TRepresent the coefficient of heat conduction, J^f And J^TRepresent fluid flow and fluid heat respectively；Following formula (8) (9) is Fluid in Pore flowing and answers Equilibrium equation；

σ_ij,j=0 (8)

$\frac{\∂\mathrm{\ζ}}{\∂t}=-\frac{1}{{\mathrm{\ρ}}_f}\▿J^f---\left(9\right)$

Bring governing equation (1) (2) above into equilibrium equation (8) (9), obtain about rock deformation Field equation (10) (11) with fluid flow state；

$(K+\frac{G}{3})\▿(\▿.u)+G{\▿}^{2}u+m(\mathrm{\α}\▿p+{\mathrm{\γ}}_1\▿T)=0---\left(10\right)$

$\mathrm{\α}(\▿.\stackrel{\·}{u})+\mathrm{\β}\stackrel{\·}{p}-\frac{k}{\mathrm{\η}}{\▿}^{2}p-{\mathrm{\γ}}_2\stackrel{\·}{T}=0---\left(11\right)$

$\stackrel{\·}{T}+v(\▿T)-c^T{\▿}^{2}T=0---\left(12\right)$

U is displacement, the m=[1,1,0] when model is two dimension^T, m=[1,1,1,0,0,0] when model is three-dimensional^T, c^TFor thermal diffusion coefficient, speed and the pore pressure of fluid follow Darcy's law

Step 1, for the assigned pressure of initial fracturing, uses finite element method (FEM) to carry out modeling；

Step 2, carries out the rock unit Stress calculation of each particular moment according to time change；

Step 3, according to stress calculating results, in conjunction with criterion in rock, it is judged that current rock unit is No rupture；

Step 4, calculates porosity parameter；

Step 5, if rock unit ruptures, calculates focus earthquake magnitude in conjunction with microseism focal mechanism；With And

Step 6, repeats step 2, step 3, step 4 and step 5, completes the forward simulation in all moment.

Micro-seismic event forward simulation based on reservoir geology mechanics parameter the most according to claim 1 Method, it is characterised in that in step 3, in a model, reservoir rock can by the process after stress effect To be described as elastic deformation and to rupture two stages；During by less stress, rock generation elastic deformation, this One stage rock interior does not ruptures；And when rock stress increases the fracture criteria reaching rock, rock Stone starts to rupture, and along with the development ruptured, ultimately forms crack.

Micro-seismic event forward simulation based on reservoir geology mechanics parameter the most according to claim 2 Method, it is characterised in that in step 3, in elastic deformation-rupture mechanism, the reduction generation of elastic modelling quantity The process that table rock ruptures,

E=(1-d) E₀ (13)

D represents the degree that rock interior ruptures, and is properly termed as rock fracture coefficient, E₀Represent rock Initial elastic modulus, E is current elastic modelling quantity；If rock is stressed and there occurs and rupture, rock breaks Split coefficient to be stated by the residual compression of rock:

$d=1-\frac{f_{cr}}{E_0\stackrel{\‾}{\mathrm{\ϵ}}}---\left(14\right)$

f_crFor the residual compression of rock,For under corresponding pressure occur strain, thenIllustrate initial pressure Stress.

Micro-seismic event forward simulation based on reservoir geology mechanics parameter the most according to claim 1 Method, it is characterised in that in step 4, it is thus achieved that the actual loading situation of a certain each underground unit of moment, And contrast actual loading size and rock unit keeps not rupturing born maximum stress, be made whether The judgement ruptured, if ruptured, then rock unit deforms upon, and porosity parameter is according to deformation Change；If rock rupture, then porosity is undergone mutation, and is converted into a bigger numerical value, should Numerical value is constant.

Micro-seismic event forward simulation based on reservoir geology mechanics parameter the most according to claim 1 Method, it is characterised in that in step 5, if it occur that rupture, simulates according to the situation of change of stress The intensity of micro-seismic event, micro-seismic event hypocentral location is current computing unit position, vibrations the moment be Current time.

Micro-seismic event forward simulation based on reservoir geology mechanics parameter the most according to claim 1 Method, it is characterised in that be somebody's turn to do micro-seismic event the Forward Modeling based on reservoir geology mechanics parameter also Including, after step 5, all analog results are exported, as analyzing and explaining depending on of reservoir pressure According to.