CN105445794A - A double-phase anisotropic medium staggered-grid finite-difference analogue method and apparatus - Google Patents

Abstract

The invention provides a double-phase anisotropic medium staggered-grid finite-difference analogue method and apparatus. Through the double-phase anisotropic medium staggered-grid finite-difference analogue method and apparatus provided by the embodiment of the invention, in a process of performing simulation of spreading of seismic waves in a double-phase anisotropic medium via a staggered grid, interpolation calculation does not need to be applied to shear forces, and the precision of obtained speed-stress data of medium points of double-phase anisotropic mediums can be raised.

Description

Two-phase Anisotropic Media staggering mesh finite-difference analogy method and device

Technical field

The present invention relates to calculating technical field of data processing, in particular to Two-phase Anisotropic Media staggering mesh finite-difference analogy method and device.

Background technology

At present, in order to study the feature of earth interior Two-phase Anisotropic Media, understanding being deepened to reservoir medium, a series of forward simulation need be done to Two-phase Anisotropic Media.Frequent employing staggered-mesh is simulated the propagation of seismic event in Two-phase Anisotropic Media.

In correlation technique, in the process adopting staggered-mesh to simulate the propagation of seismic event in Two-phase Anisotropic Media, net point is divided into whole net point and half net point, Particle Vibration Velocity and suffered stress are defined in two differences respectively and on adjacent time horizon, usually all normal stress is positioned over whole net point place, and shearing stress is positioned over half net point place.Utilize the field value of adjacent whole net point and half net point to respectively computing velocity component and the components of stress along the difference of horizontal component and vertical component, differential is replaced by difference, physical quantity simultaneously on adjacent two time horizons is lucky staggered half grid in space distribution, to reach the staggered object of Time and place.

In the process propagation of seismic event in Two-phase Anisotropic Media simulated by staggered-mesh, must by carrying out interpolation to shearing stress component, just can simulate the propagation of seismic event in Two-phase Anisotropic Media, interpolation calculation itself is out of true also, so the error of calculation of simulation process can be increased, reduce the precision of the speed-stress data of the medium particle of the Two-phase Anisotropic Media obtained.

Summary of the invention

In view of this, the object of the embodiment of the present invention is to provide Two-phase Anisotropic Media staggering mesh finite-difference analogy method and device, in the process propagation of seismic event in Two-phase Anisotropic Media simulated by staggered-mesh, without the need to using interpolation calculation, to improve the precision of the speed-stress data of the Two-phase Anisotropic Media particle obtained to shearing stress.

First aspect, embodiments provides a kind of Two-phase Anisotropic Media staggering mesh finite-difference analogy method, comprising:

According to the mode of motion of described Two-phase Anisotropic Media, the speed of described Two-phase Anisotropic Media and the components of stress are separately positioned on the net point place of staggered-mesh, wherein, the shearing stress component of the medium particle of Two-phase Anisotropic Media is arranged on the whole net point place in staggered-mesh;

According to the speed of the medium particle determined in described staggered-mesh and components of stress position, the discretize data of computing velocity and the components of stress;

By the described staggered-mesh obtained, the propagation of seismic event in Two-phase Anisotropic Media is simulated, obtain the speed-stress data of the medium particle of described Two-phase Anisotropic Media.

In conjunction with first aspect, embodiments provide the first possible embodiment of first aspect, wherein, the shearing stress component of the medium particle of Two-phase Anisotropic Media be arranged on the whole net point place in staggered-mesh, comprise:

Pass through formula

$\frac{\∂{\mathrm{\σ}}_{xy}}{\∂x}=c_1\frac{{\mathrm{\σ}}_{xy}(i+1,j)-{\mathrm{\σ}}_{xy}(i,j)}{\mathrm{\Δ}x}+c_2\frac{{\mathrm{\σ}}_{xy}(i+2,j)-{\mathrm{\σ}}_{xy}(i-1,j)}{\mathrm{\Δ}x}$ Represent that the shearing stress component of the medium particle of described Two-phase Anisotropic Media is arranged on the whole net point place in staggered-mesh;

Wherein, σ _xyrepresent x direction stress component in the y-direction, c ₁and c ₂representation space difference coefficient respectively, i, j get xyz respectively, and x represents horizontal direction, and y represents vertical direction, and z represents vertical direction.

In conjunction with first aspect, embodiments provide the embodiment that the second of first aspect is possible, wherein, according to the mode of motion of described Two-phase Anisotropic Media, the shearing stress component of the medium particle of Two-phase Anisotropic Media is arranged on the whole net point place in staggered-mesh, comprises:

Pass through formula $\left\{\begin{array}{c}{\mathrm{\ρ}}_{11}\frac{{\∂}^2{w}_i}{\∂t^2}+{\mathrm{\ρ}}_{12}\frac{{\∂}^2{W}_i}{\∂t^2}={\mathrm{\τ}}_{{\mathrm{ij}}^{\′}j}+b_{ij}\left(\frac{\∂W_j}{\∂t}-\frac{\∂w_j}{\∂t}\right)\\ {\mathrm{\ρ}}_{12}\frac{{\∂}^2{w}_i}{\∂t^2}+{\mathrm{\ρ}}_{22}\frac{{\∂}^2{W}_i}{\∂t^2}=s_{{}^{\′}i}-b_{ij}\left(\frac{\∂W_j}{\∂t}-\frac{\∂w_j}{\∂t}\right)\end{array}\right.$ Represent the mode of motion of described Two-phase Anisotropic Media;

Wherein, w _jand W _jrepresent the component of displacement in j direction of solid phase and stream phase respectively, τ _ij'jfor the solid phase components of stress are at the local derviation in j direction, b _ijdissipation factor during solid phase motion relative to stream, s represents the stress acting on stream phase, and i, j get xyz respectively, and x represents horizontal direction, and y represents vertical direction, and z represents vertical direction.

In conjunction with first aspect, embodiments provide the third possible embodiment of first aspect, wherein, according to the speed of the medium particle determined in staggered-mesh and components of stress position, the discretize data of computing velocity and the components of stress, comprising:

Pass through formula $\left\{\begin{array}{c}\frac{\∂{\mathrm{\υ}}_y}{\∂t}=(D_2+D_3)b_{22}({\mathrm{\υ}}_y-V_y)-D_3(\frac{\∂{\mathrm{\τ}}_{xy}}{\∂x}+\frac{\∂{\mathrm{\τ}}_{yz}}{\∂z})\\ \frac{\∂V_y}{\∂t}=-(D_1+D_2)b_{22}({\mathrm{\υ}}_y-V_y)+D_2(\frac{\∂{\mathrm{\τ}}_{xy}}{\∂x}+\frac{\∂{\mathrm{\τ}}_{yz}}{\∂z})\end{array}\right.$ Represent described vertical direction speed component;

Wherein, ν _y, V _yrepresent solid-phase component and stream phase constituent vertically speed component respectively, τ _ijrepresent the component of i direction stress along j direction, D _ifor the multi-term expression about density, b ₂₂represent dissipation factor, i, j all represent x, y and z, x represent horizontal direction, and y represents vertical direction, and z represents vertical direction.

In conjunction with first aspect, embodiments provide the 4th kind of possible embodiment of first aspect, wherein, according to the speed of the medium particle determined in staggered-mesh and components of stress position, the discretize data of computing velocity and the components of stress, comprising:

Pass through formula

$\left\{\begin{array}{c}{\mathrm{\υ}}_{yi,j}^{n+\frac{1}{2}}={\mathrm{\υ}}_{yi,j}^{n-\frac{1}{2}}-D_{3i,j}(Q_x^{-}\left({\mathrm{\τ}}_{xyi,j}^n\right)+Q_z^{-}\left({\mathrm{\τ}}_{yzi,j}^n\right))+\mathrm{\Δ}t(D_{2i,j}+D_{3i,j})b_{22i,j}({\mathrm{\υ}}_{yi,j}^{n-\frac{1}{2}}-V_{yi,j}^{n-\frac{1}{2}})\\ V_{yi,j}^{n+\frac{1}{2}}=V_{yi,j}^{n-\frac{1}{2}}+D_{2i,j}(Q_x^{-}\left({\mathrm{\τ}}_{xyi,j}^n\right)+Q_z^{-}\left({\mathrm{\τ}}_{yzi,j}^n\right))-\mathrm{\Δ}t(D_{1i,j}+D_{2i,j})b_{22i,j}({\mathrm{\υ}}_{yi,j}^{n-\frac{1}{2}}-V_{yi,j}^{n-\frac{1}{2}})\end{array}\right.$ Calculate the discretize data of vertical direction speed component;

Wherein, Q representation space difference operator, Δ t represents sampling time interval, ν _y, V _yrepresent solid-phase component and stream phase constituent vertically speed component respectively, τ _ijrepresent the component of i direction stress along j direction, D _ifor the multi-term expression about density, b ₂₂represent dissipation factor.

Second aspect, embodiments provides a kind of Two-phase Anisotropic Media staggering mesh finite-difference analogue means, comprising:

Module is set, for the mode of motion according to described Two-phase Anisotropic Media, the speed of described Two-phase Anisotropic Media and the components of stress are separately positioned on the net point place of staggered-mesh, wherein, the shearing stress component of the medium particle of Two-phase Anisotropic Media is arranged on the whole net point place in staggered-mesh;

Computing module, for according to the speed of the medium particle determined in described staggered-mesh and components of stress position, the discretize data of computing velocity and the components of stress;

Analog module, simulates the propagation of seismic event in Two-phase Anisotropic Media for the described staggered-mesh by obtaining, and obtains the speed-stress data of the medium particle of described Two-phase Anisotropic Media.

In conjunction with second aspect, embodiments provide the first possible embodiment of second aspect, wherein, described module is set, for:

Pass through formula

$\frac{\∂{\mathrm{\σ}}_{xy}}{\∂x}=c_1\frac{{\mathrm{\σ}}_{xy}(i+1,j)-{\mathrm{\σ}}_{xy}(i,j)}{\mathrm{\Δ}x}+c_2\frac{{\mathrm{\σ}}_{xy}(i+2,j)-{\mathrm{\σ}}_{xy}(i-1,j)}{\mathrm{\Δ}x}$ Represent that the shearing stress component of the medium particle of described Two-phase Anisotropic Media is arranged on the whole net point place in staggered-mesh;

Wherein, σ _xyrepresent x direction stress component in the y-direction, c ₁and c ₂representation space difference coefficient respectively, i, j get xyz respectively, and x represents horizontal direction, and y represents vertical direction, and z represents vertical direction.

In conjunction with second aspect, embodiments provide the embodiment that the second of second aspect is possible, wherein, described module is set, for:

Pass through formula $\left\{\begin{array}{c}{\mathrm{\ρ}}_{11}\frac{{\∂}^2{w}_i}{\∂t^2}+{\mathrm{\ρ}}_{12}\frac{{\∂}^2{W}_i}{\∂t^2}={\mathrm{\τ}}_{{\mathrm{ij}}^{\′}j}+b_{ij}\left(\frac{\∂W_j}{\∂t}-\frac{\∂w_j}{\∂t}\right)\\ {\mathrm{\ρ}}_{12}\frac{{\∂}^2{w}_i}{\∂t^2}+{\mathrm{\ρ}}_{22}\frac{{\∂}^2{W}_i}{\∂t^2}=s_{{}^{\′}i}-b_{ij}\left(\frac{\∂W_j}{\∂t}-\frac{\∂w_j}{\∂t}\right)\end{array}\right.$ Represent the mode of motion of described Two-phase Anisotropic Media;

Wherein, w _jand W _jrepresent the component of displacement in j direction of solid phase and stream phase respectively, τ _ij'jfor the solid phase components of stress are at the local derviation in j direction, b _ijdissipation factor during solid phase motion relative to stream, s represents the stress acting on stream phase, and i, j get xyz respectively, and x represents horizontal direction, and y represents vertical direction, and z represents vertical direction.

In conjunction with second aspect, embodiments provide the third possible embodiment of second aspect, wherein, described computing module, for:

Pass through formula $\left\{\begin{array}{c}\frac{\∂{\mathrm{\υ}}_y}{\∂t}=(D_2+D_3)b_{22}({\mathrm{\υ}}_y-V_y)-D_3(\frac{\∂{\mathrm{\τ}}_{xy}}{\∂x}+\frac{\∂{\mathrm{\τ}}_{yz}}{\∂z})\\ \frac{\∂V_y}{\∂t}=-(D_1+D_2)b_{22}({\mathrm{\υ}}_y-V_y)+D_2(\frac{\∂{\mathrm{\τ}}_{xy}}{\∂x}+\frac{\∂{\mathrm{\τ}}_{yz}}{\∂z})\end{array}\right.$ Represent described vertical direction speed component;

Wherein, ν _y, V _yrepresent solid-phase component and stream phase constituent vertically speed component respectively, τ _ijrepresent the component of i direction stress along j direction, D _ifor the multi-term expression about density, b ₂₂represent dissipation factor, i, j all represent x, y and z, x represent horizontal direction, and y represents vertical direction, and z represents vertical direction.

In conjunction with second aspect, embodiments provide the 4th kind of possible embodiment of second aspect, wherein, described computing module, for:

Pass through formula

$\left\{\begin{array}{c}{\mathrm{\υ}}_{yi,j}^{n+\frac{1}{2}}={\mathrm{\υ}}_{yi,j}^{n-\frac{1}{2}}-D_{3i,j}(Q_x^{-}\left({\mathrm{\τ}}_{xyi,j}^n\right)+Q_z^{-}\left({\mathrm{\τ}}_{yzi,j}^n\right))+\mathrm{\Δ}t(D_{2i,j}+D_{3i,j})b_{22i,j}({\mathrm{\υ}}_{yi,j}^{n-\frac{1}{2}}-V_{yi,j}^{n-\frac{1}{2}})\\ V_{yi,j}^{n+\frac{1}{2}}=V_{yi,j}^{n-\frac{1}{2}}+D_{2i,j}(Q_x^{-}\left({\mathrm{\τ}}_{xyi,j}^n\right)+Q_z^{-}\left({\mathrm{\τ}}_{yzi,j}^n\right))-\mathrm{\Δ}t(D_{1i,j}+D_{2i,j})b_{22i,j}({\mathrm{\υ}}_{yi,j}^{n-\frac{1}{2}}-V_{yi,j}^{n-\frac{1}{2}})\end{array}\right.$

Calculate the discretize data of vertical direction speed component;

Wherein, Q representation space difference operator, Δ t represents sampling time interval, ν _y, V _yrepresent solid-phase component and stream phase constituent vertically speed component respectively, τ _ijrepresent the component of i direction stress along j direction, D _ifor the multi-term expression about density, b ₂₂represent dissipation factor.

The Two-phase Anisotropic Media staggering mesh finite-difference analogy method that the embodiment of the present invention provides and device, the whole net point place in staggered-mesh is arranged on by the shearing stress component of the medium particle by Two-phase Anisotropic Media, when Simulating Seismic Wave, the shearing stress component of medium particle is used to be positioned at the staggered-mesh at whole net point place, the seismic event propagated in Two-phase Anisotropic Media is simulated, thus without the need to using interpolation calculation to shearing stress component in simulation process, just can be simulated the propagation of seismic event in Two-phase Anisotropic Media by staggered-mesh, improve the precision of the speed-stress data of the medium particle of the Two-phase Anisotropic Media obtained, so that the propagation law of Study of Seismic ripple in Two-phase Anisotropic Media, instruct the production in reality.

For making above-mentioned purpose of the present invention, feature and advantage become apparent, preferred embodiment cited below particularly, and coordinate appended accompanying drawing, be described in detail below.

Accompanying drawing explanation

In order to be illustrated more clearly in the technical scheme of the embodiment of the present invention, be briefly described to the accompanying drawing used required in embodiment below, be to be understood that, the following drawings illustrate only some embodiment of the present invention, therefore the restriction to scope should be counted as, for those of ordinary skill in the art, under the prerequisite not paying creative work, other relevant accompanying drawings can also be obtained according to these accompanying drawings.

Fig. 1 shows the process flow diagram of a kind of Two-phase Anisotropic Media staggering mesh finite-difference analogy method that the embodiment of the present invention 1 provides;

Fig. 2 shows the existing staggered-mesh structure described in the embodiment of the present invention 1;

Fig. 3 shows the staggered-mesh structure after the improvement described in the embodiment of the present invention 1;

Fig. 4 shows the structural representation of a kind of Two-phase Anisotropic Media staggering mesh finite-difference analogue means that the embodiment of the present invention 2 provides.

Embodiment

Below in conjunction with accompanying drawing in the embodiment of the present invention, be clearly and completely described the technical scheme in the embodiment of the present invention, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.The assembly of the embodiment of the present invention describing and illustrate in usual accompanying drawing herein can be arranged with various different configuration and design.Therefore, below to the detailed description of the embodiments of the invention provided in the accompanying drawings and the claimed scope of the present invention of not intended to be limiting, but selected embodiment of the present invention is only represented.Based on embodiments of the invention, the every other embodiment that those skilled in the art obtain under the prerequisite not making creative work, all belongs to the scope of protection of the invention.

To consider in correlation technique in the process propagation of seismic event in Two-phase Anisotropic Media simulated by staggered-mesh, must by carrying out interpolation to shearing stress component, just can simulate the propagation of seismic event in Two-phase Anisotropic Media, interpolation calculation itself is out of true also, so the error of calculation of simulation process can be increased, reduce the precision of the speed-stress data of the medium particle of the Two-phase Anisotropic Media obtained.Based on this, embodiments provide Two-phase Anisotropic Media staggering mesh finite-difference analogy method and device, be described below by embodiment.

Embodiment 1

Present embodiments provide a kind of Two-phase Anisotropic Media staggering mesh finite-difference analogy method.The executive agent of the embodiment of the present invention is server, when server receives the dummy instruction of technician's triggering, by obtaining the discretization equation of stress and speed component according to the stress of the particle of the Two-phase Anisotropic Media determined in staggered-mesh and speed component, and when Simulating Seismic Wave, use discretize stress in staggered-mesh of medium particle and speed component, the seismic event propagated in Two-phase Anisotropic Media is simulated.

Server, can use existing any computing equipment or computing machine, simulate, repeat no longer one by one here the seismic event propagated in Two-phase Anisotropic Media.

See Fig. 1, a kind of Two-phase Anisotropic Media staggering mesh finite-difference analogy method, comprises the following steps:

Step 100, mode of motion according to Two-phase Anisotropic Media, the speed of Two-phase Anisotropic Media and the components of stress are separately positioned on the net point place of staggered-mesh, wherein, the shearing stress component of the medium particle of Two-phase Anisotropic Media is arranged on the whole net point place in staggered-mesh.

In correlation technique, in staggered-mesh, usually all normal stress is positioned over whole net point place, and shearing stress is positioned over half net point place.The definition mode of stress and speed as shown in Figure 2.Wherein, position 1 place represents normal stress τ _iiposition, 2 places are shearing stress τ _ijposition, 3 places are speed component position, x and y direction, and 4 places are z (vertically) direction speed component position.

And in the present embodiment, as shown in Figure 3, the shearing stress component of medium particle is positioned at the whole net point place of staggered-mesh, thus in simulation process, without the need to using interpolation calculation to shearing stress component, just can be simulated the propagation of seismic event in Two-phase Anisotropic Media by staggered-mesh.

In step 100, formula is passed through

$\frac{\∂{\mathrm{\σ}}_{xy}}{\∂x}=c_1\frac{{\mathrm{\σ}}_{xy}(i+1,j)-{\mathrm{\σ}}_{xy}(i,j)}{\mathrm{\Δ}x}+c_2\frac{{\mathrm{\σ}}_{xy}(i+2,j)-{\mathrm{\σ}}_{xy}(i-1,j)}{\mathrm{\Δ}x}$ Represent that the shearing stress component of the medium particle of Two-phase Anisotropic Media is arranged on the whole net point place in staggered-mesh;

Wherein, σ _xyrepresent x direction stress component in the y-direction, c ₁and c ₂representation space difference coefficient respectively, i, j get xyz respectively, and x represents horizontal direction, and y represents vertical direction, and z represents vertical direction.

Pass through formula $\left\{\begin{array}{c}{\mathrm{\ρ}}_{11}\frac{{\∂}^2{w}_i}{\∂t^2}+{\mathrm{\ρ}}_{12}\frac{{\∂}^2{W}_i}{\∂t^2}={\mathrm{\τ}}_{{\mathrm{ij}}^{\′}j}+b_{ij}\left(\frac{\∂W_j}{\∂t}-\frac{\∂w_j}{\∂t}\right)\\ {\mathrm{\ρ}}_{12}\frac{{\∂}^2{w}_i}{\∂t^2}+{\mathrm{\ρ}}_{22}\frac{{\∂}^2{W}_i}{\∂t^2}=s_{{}^{\′}i}-b_{ij}\left(\frac{\∂W_j}{\∂t}-\frac{\∂w_j}{\∂t}\right)\end{array}\right.$ Represent the mode of motion of Two-phase Anisotropic Media;

Wherein, w _jand W _jrepresent the component of displacement in j direction of solid phase and stream phase respectively,

τ _ij'jfor the solid phase components of stress are at the local derviation in j direction, b _ijdissipation factor during solid phase motion relative to stream, s represents the stress acting on stream phase, and i, j get xyz respectively.

The mode of motion of Two-phase Anisotropic Media is theoretical based on Biot two-phase media, in conjunction with the characteristic of anisotropic medium, and the formula drawn.

Step 102, according to the speed of the medium particle determined in staggered-mesh and components of stress position, the discretize data of computing velocity and the components of stress.

Wherein, formula is passed through $\left\{\begin{array}{c}\frac{\∂{\mathrm{\υ}}_y}{\∂t}=(D_2+D_3)b_{22}({\mathrm{\υ}}_y-V_y)-D_3(\frac{\∂{\mathrm{\τ}}_{xy}}{\∂x}+\frac{\∂{\mathrm{\τ}}_{yz}}{\∂z})\\ \frac{\∂V_y}{\∂t}=-(D_1+D_2)b_{22}({\mathrm{\υ}}_y-V_y)+D_2(\frac{\∂{\mathrm{\τ}}_{xy}}{\∂x}+\frac{\∂{\mathrm{\τ}}_{yz}}{\∂z})\end{array}\right.$ Represent vertical direction speed component;

Wherein, ν _y, V _yrepresent solid-phase component and stream phase constituent vertically speed component respectively, τ _ijrepresent the component of i direction stress along j direction, D _ifor the multi-term expression about density, b ₂₂represent dissipation factor, i, j all represent x, y and z.

Pass through formula

$\left\{\begin{array}{c}{\mathrm{\υ}}_{yi,j}^{n+\frac{1}{2}}={\mathrm{\υ}}_{yi,j}^{n-\frac{1}{2}}-D_{3i,j}(Q_x^{-}\left({\mathrm{\τ}}_{xyi,j}^n\right)+Q_z^{-}\left({\mathrm{\τ}}_{yzi,j}^n\right))+\mathrm{\Δ}t(D_{2i,j}+D_{3i,j})b_{22i,j}({\mathrm{\υ}}_{yi,j}^{n-\frac{1}{2}}-V_{yi,j}^{n-\frac{1}{2}})\\ V_{yi,j}^{n+\frac{1}{2}}=V_{yi,j}^{n-\frac{1}{2}}+D_{2i,j}(Q_x^{-}\left({\mathrm{\τ}}_{xyi,j}^n\right)+Q_z^{-}\left({\mathrm{\τ}}_{yzi,j}^n\right))-\mathrm{\Δ}t(D_{1i,j}+D_{2i,j})b_{22i,j}({\mathrm{\υ}}_{yi,j}^{n-\frac{1}{2}}-V_{yi,j}^{n-\frac{1}{2}})\end{array}\right.$

Calculate the discretize data of vertical direction speed component;

Wherein, Q representation space difference operator, Δ t represents sampling time interval, ν _y, V _yrepresent solid-phase component and stream phase constituent vertically speed component respectively, τ _ijrepresent the component of i direction stress along j direction, D _ifor the multi-term expression about density, b ₂₂represent dissipation factor.

The components of stress, can by the existing formulae discovery that can be used for calculated stress component out, repeat here no longer one by one.

Step 104, by the staggered-mesh obtained, the propagation of seismic event in Two-phase Anisotropic Media to be simulated, obtain the speed-stress data of the medium particle of Two-phase Anisotropic Media.

In sum, the Two-phase Anisotropic Media staggering mesh finite-difference analogy method that the present embodiment provides, the whole net point place in staggered-mesh is arranged on by the shearing stress component of the medium particle by Two-phase Anisotropic Media, when Simulating Seismic Wave, the shearing stress component of medium particle is used to be positioned at the staggered-mesh at whole net point place, the seismic event propagated in Two-phase Anisotropic Media is simulated, thus without the need to using interpolation calculation to shearing stress component in simulation process, just can be simulated the propagation of seismic event in Two-phase Anisotropic Media by staggered-mesh, improve the precision of the speed-stress data of the medium particle of the Two-phase Anisotropic Media obtained, so that the propagation law of Study of Seismic ripple in Two-phase Anisotropic Media, instruct the production in reality.

Embodiment 2

See Fig. 4, the present embodiment provides a kind of Two-phase Anisotropic Media staggering mesh finite-difference analogue means, for performing above-mentioned Two-phase Anisotropic Media staggering mesh finite-difference simulation, comprising:

Module 400 is set, for the mode of motion according to Two-phase Anisotropic Media, the speed of Two-phase Anisotropic Media and the components of stress are separately positioned on the net point place of staggered-mesh, wherein, the shearing stress component of the medium particle of Two-phase Anisotropic Media is arranged on the whole net point place in staggered-mesh;

Computing module 402, for according to the speed of the medium particle determined in staggered-mesh and components of stress position, the discretize data of computing velocity and the components of stress;

Analog module 404, simulates the propagation of seismic event in Two-phase Anisotropic Media for the staggered-mesh by obtaining, and obtains the speed-stress data of the medium particle of Two-phase Anisotropic Media.

Wherein, module is set, for:

Pass through formula

$\frac{\∂{\mathrm{\σ}}_{xy}}{\∂x}=c_1\frac{{\mathrm{\σ}}_{xy}(i+1,j)-{\mathrm{\σ}}_{xy}(i,j)}{\mathrm{\Δ}x}+c_2\frac{{\mathrm{\σ}}_{xy}(i+2,j)-{\mathrm{\σ}}_{xy}(i-1,j)}{\mathrm{\Δ}x}$ Represent that the shearing stress component of the medium particle of Two-phase Anisotropic Media is arranged on the whole net point place in staggered-mesh;

Wherein, σ _xyrepresent x direction stress component in the y-direction, c ₁and c ₂representation space difference coefficient respectively, i, j get xyz respectively, and x represents horizontal direction, and y represents vertical direction, and z represents vertical direction.

Module is set, for:

Pass through formula $\left\{\begin{array}{c}{\mathrm{\ρ}}_{11}\frac{{\∂}^2{w}_i}{\∂t^2}+{\mathrm{\ρ}}_{12}\frac{{\∂}^2{W}_i}{\∂t^2}={\mathrm{\τ}}_{{\mathrm{ij}}^{\′}j}+b_{ij}\left(\frac{\∂W_j}{\∂t}-\frac{\∂w_j}{\∂t}\right)\\ {\mathrm{\ρ}}_{12}\frac{{\∂}^2{w}_i}{\∂t^2}+{\mathrm{\ρ}}_{22}\frac{{\∂}^2{W}_i}{\∂t^2}=s_{{}^{\′}i}-b_{ij}\left(\frac{\∂W_j}{\∂t}-\frac{\∂w_j}{\∂t}\right)\end{array}\right.$ Represent the mode of motion of Two-phase Anisotropic Media;

Wherein, w _jand W _jrepresent the component of displacement in j direction of solid phase and stream phase respectively, τ _ij'jfor the solid phase components of stress are at the local derviation in j direction, b _ijdissipation factor during solid phase motion relative to stream, s represents the stress acting on stream phase, and i, j get xyz respectively.

Wherein, computing module, for:

Pass through formula $\left\{\begin{array}{c}\frac{\∂{\mathrm{\υ}}_y}{\∂t}=(D_2+D_3)b_{22}({\mathrm{\υ}}_y-V_y)-D_3(\frac{\∂{\mathrm{\τ}}_{xy}}{\∂x}+\frac{\∂{\mathrm{\τ}}_{yz}}{\∂z})\\ \frac{\∂V_y}{\∂t}=-(D_1+D_2)b_{22}({\mathrm{\υ}}_y-V_y)+D_2(\frac{\∂{\mathrm{\τ}}_{xy}}{\∂x}+\frac{\∂{\mathrm{\τ}}_{yz}}{\∂z})\end{array}\right.$ Represent vertical direction speed component;

Wherein, ν _y, V _yrepresent solid-phase component and stream phase constituent vertically speed component respectively, τ _ijrepresent the component of i direction stress along j direction, D _ifor the multi-term expression about density, b ₂₂represent dissipation factor, i, j all represent x, y and z, x represent horizontal direction, and y represents vertical direction, and z represents vertical direction.

Wherein, computing module, for:

Pass through formula

$\left\{\begin{array}{c}{\mathrm{\υ}}_{yi,j}^{n+\frac{1}{2}}={\mathrm{\υ}}_{yi,j}^{n-\frac{1}{2}}-D_{3i,j}(Q_x^{-}\left({\mathrm{\τ}}_{xyi,j}^n\right)+Q_z^{-}\left({\mathrm{\τ}}_{yzi,j}^n\right))+\mathrm{\Δ}t(D_{2i,j}+D_{3i,j})b_{22i,j}({\mathrm{\υ}}_{yi,j}^{n-\frac{1}{2}}-V_{yi,j}^{n-\frac{1}{2}})\\ V_{yi,j}^{n+\frac{1}{2}}=V_{yi,j}^{n-\frac{1}{2}}+D_{2i,j}(Q_x^{-}\left({\mathrm{\τ}}_{xyi,j}^n\right)+Q_z^{-}\left({\mathrm{\τ}}_{yzi,j}^n\right))-\mathrm{\Δ}t(D_{1i,j}+D_{2i,j})b_{22i,j}({\mathrm{\υ}}_{yi,j}^{n-\frac{1}{2}}-V_{yi,j}^{n-\frac{1}{2}})\end{array}\right.$

Calculate the discretize data of vertical direction speed component;

Wherein, Q representation space difference operator, Δ t represents sampling time interval, ν _y, V _yrepresent solid-phase component and stream phase constituent vertically speed component respectively, τ _ijrepresent the component of i direction stress along j direction, D _ifor the multi-term expression about density, b ₂₂represent dissipation factor.

In sum, the Two-phase Anisotropic Media staggering mesh finite-difference analogy method that the present embodiment provides, the whole net point place in staggered-mesh is arranged on by the shearing stress component of the medium particle by Two-phase Anisotropic Media, when Simulating Seismic Wave, the shearing stress component of medium particle is used to be positioned at the staggered-mesh at whole net point place, the seismic event propagated in Two-phase Anisotropic Media is simulated, thus without the need to using interpolation calculation to shearing stress component in simulation process, just can be simulated the propagation of seismic event in Two-phase Anisotropic Media by staggered-mesh, improve the precision of the speed-stress data of the medium particle of the Two-phase Anisotropic Media obtained, so that the propagation law of Study of Seismic ripple in Two-phase Anisotropic Media, instruct the production in reality.

The computer program carrying out Two-phase Anisotropic Media staggering mesh finite-difference analogy method that the embodiment of the present invention provides, comprise the computer-readable recording medium storing program code, the instruction that described program code comprises can be used for performing the method described in previous methods embodiment, specific implementation see embodiment of the method, can not repeat them here.

Those skilled in the art can be well understood to, and for convenience and simplicity of description, the specific works process of the system of foregoing description, device and unit, with reference to the corresponding process in preceding method embodiment, can not repeat them here.

In several embodiments that the application provides, should be understood that disclosed system, apparatus and method can realize by another way.Device embodiment described above is only schematic, such as, the division of described unit, be only a kind of logic function to divide, actual can have other dividing mode when realizing, again such as, multiple unit or assembly can in conjunction with or another system can be integrated into, or some features can be ignored, or do not perform.Another point, shown or discussed coupling each other or direct-coupling or communication connection can be by some communication interfaces, and the indirect coupling of device or unit or communication connection can be electrical, machinery or other form.

The described unit illustrated as separating component or can may not be and physically separates, and the parts as unit display can be or may not be physical location, namely can be positioned at a place, or also can be distributed in multiple network element.Some or all of unit wherein can be selected according to the actual needs to realize the object of the present embodiment scheme.

In addition, each functional unit in each embodiment of the present invention can be integrated in a processing unit, also can be that the independent physics of unit exists, also can two or more unit in a unit integrated.

If described function using the form of SFU software functional unit realize and as independently production marketing or use time, can be stored in a computer read/write memory medium.Based on such understanding, the part of the part that technical scheme of the present invention contributes to prior art in essence in other words or this technical scheme can embody with the form of software product, this computer software product is stored in a storage medium, comprising some instructions in order to make a computer equipment (can be personal computer, server, or the network equipment etc.) perform all or part of step of method described in each embodiment of the present invention.And aforesaid storage medium comprises: USB flash disk, portable hard drive, ROM (read-only memory) (ROM, Read-OnlyMemory), random access memory (RAM, RandomAccessMemory), magnetic disc or CD etc. various can be program code stored medium.

The above; be only the specific embodiment of the present invention, but protection scope of the present invention is not limited thereto, is anyly familiar with those skilled in the art in the technical scope that the present invention discloses; change can be expected easily or replace, all should be encompassed within protection scope of the present invention.Therefore, protection scope of the present invention should described be as the criterion with the protection domain of claim.

Claims (10)

1. a Two-phase Anisotropic Media staggering mesh finite-difference analogy method, is characterized in that, comprising:

According to the mode of motion of described Two-phase Anisotropic Media, the speed of described Two-phase Anisotropic Media and the components of stress are separately positioned on the net point place of staggered-mesh, wherein, the shearing stress component of the medium particle of Two-phase Anisotropic Media is arranged on the whole net point place in staggered-mesh;

According to the speed of the medium particle determined in described staggered-mesh and components of stress position, the discretize data of computing velocity and the components of stress;

By the described staggered-mesh obtained, the propagation of seismic event in Two-phase Anisotropic Media is simulated, obtain the speed-stress data of the medium particle of described Two-phase Anisotropic Media.

2. method according to claim 1, is characterized in that, the shearing stress component of the medium particle of Two-phase Anisotropic Media is arranged on the whole net point place in staggered-mesh, comprises:

Pass through formula

$\frac{\∂{\mathrm{\σ}}_{xy}}{\∂x}=c_1\frac{{\mathrm{\σ}}_{xy}\left(i+1,j\right)-{\mathrm{\σ}}_{xy}\left(i,j\right)}{\mathrm{\Δ}x}+c_2\frac{{\mathrm{\σ}}_{xy}\left(i+2,j\right)-{\mathrm{\σ}}_{xy}\left(i-1,j\right)}{\mathrm{\Δ}x}$ Represent that the shearing stress component of the medium particle of described Two-phase Anisotropic Media is arranged on the whole net point place in staggered-mesh;

Wherein, σ _xyrepresent x direction stress component in the y-direction, c ₁and c ₂representation space difference coefficient respectively, i, j get xyz respectively, and x represents horizontal direction, and y represents vertical direction, and z represents vertical direction.

3. method according to claim 1, is characterized in that, according to the mode of motion of described Two-phase Anisotropic Media, the shearing stress component of the medium particle of Two-phase Anisotropic Media is arranged on the whole net point place in staggered-mesh, comprises:

Pass through formula $\left\{\begin{array}{c}{\mathrm{\ρ}}_{11}\frac{{\∂}^2{w}_i}{\∂t^2}+{\mathrm{\ρ}}_{12}\frac{{\∂}^2{W}_i}{\∂t^2}={\mathrm{\τ}}_{{\mathrm{ij}}^{\′}j}+b_{ij}\left(\frac{\∂W_j}{\∂t}-\frac{\∂w_j}{\∂t}\right)\\ {\mathrm{\ρ}}_{12}\frac{{\∂}^2{w}_i}{\∂t^2}+{\mathrm{\ρ}}_{22}\frac{{\∂}^2{W}_i}{\∂t^2}=s_{\′i}-b_{ij}\left(\frac{\∂W_j}{\∂t}-\frac{\∂w_j}{\∂t}\right)\end{array}\right.$ Represent the mode of motion of described Two-phase Anisotropic Media;

Wherein, w _jand W _jrepresent the component of displacement in j direction of solid phase and stream phase respectively, τ _ij'jfor the solid phase components of stress are at the local derviation in j direction, b _ijdissipation factor during solid phase motion relative to stream, s represents the stress acting on stream phase, and i, j get xyz respectively, and x represents horizontal direction, and y represents vertical direction, and z represents vertical direction.

4. method according to claim 1, is characterized in that, according to the speed of the medium particle determined in staggered-mesh and components of stress position, the discretize data of computing velocity and the components of stress, comprising:

Pass through formula represent vertical direction speed component;

Wherein, ν _y, V _yrepresent solid-phase component and stream phase constituent vertically speed component respectively, τ _ijrepresent the component of i direction stress along j direction, D _ifor the multi-term expression about density, b ₂₂represent dissipation factor, i, j all represent x, y and z, x represent horizontal direction, and y represents vertical direction, and z represents vertical direction.

5. method according to claim 4, is characterized in that, according to the speed of the medium particle determined in staggered-mesh and components of stress position, the discretize data of computing velocity and the components of stress, comprising:

Pass through formula

$\left\{\begin{array}{c}{\mathrm{\υ}}_{yi,j}^{n+\frac{1}{2}}={\mathrm{\υ}}_{yi,j}^{n-\frac{1}{2}}-D_{3i,j}\left(Q_x^{-}\left({\mathrm{\τ}}_{xyi,j}^n\right)+Q_z^{-}\left({\mathrm{\τ}}_{yzi,j}^n\right)\right)+\mathrm{\Δ}t\left(D_{2i,j}+D_{3i,j}\right)b_{22i,j}\left({\mathrm{\υ}}_{yi,j}^{n-\frac{1}{2}}-V_{yi,j}^{n-\frac{1}{2}}\right)\\ V_{yi,j}^{n+\frac{1}{2}}=V_{yi,j}^{n-\frac{1}{2}}+D_{2i,j}\left(Q_x^{-}\left({\mathrm{\τ}}_{xyi,j}^n\right)+Q_z^{-}\left({\mathrm{\τ}}_{yzi,j}^n\right)\right)-\mathrm{\Δ}t\left(D_{1i,j}+D_{2i,j}\right)b_{22i,j}\left({\mathrm{\υ}}_{yi,j}^{n-\frac{1}{2}}-V_{yi,j}^{n-\frac{1}{2}}\right)\end{array}\right.$

Calculate the discretize data of vertical direction speed component;

Wherein, Q representation space difference operator, Δ t represents sampling time interval, ν _y, V _yrepresent solid-phase component and stream phase constituent vertically speed component respectively, τ _ijrepresent the component of i direction stress along j direction, D _ifor the multi-term expression about density, b ₂₂represent dissipation factor.

6. a Two-phase Anisotropic Media staggering mesh finite-difference analogue means, is characterized in that, comprising:

Module is set, for the mode of motion according to described Two-phase Anisotropic Media, the speed of described Two-phase Anisotropic Media and the components of stress are separately positioned on the net point place of staggered-mesh, wherein, the shearing stress component of the medium particle of Two-phase Anisotropic Media is arranged on the whole net point place in staggered-mesh;

Computing module, for according to the speed of the medium particle determined in described staggered-mesh and components of stress position, the discretize data of computing velocity and the components of stress;

Analog module, simulates the propagation of seismic event in Two-phase Anisotropic Media for the described staggered-mesh by obtaining, and obtains the speed-stress data of the medium particle of described Two-phase Anisotropic Media.

7. device according to claim 6, is characterized in that, describedly arranges module, for:

Pass through formula

$\frac{\∂{\mathrm{\σ}}_{xy}}{\∂x}=c_1\frac{{\mathrm{\σ}}_{xy}\left(i+1,j\right)-{\mathrm{\σ}}_{xy}\left(i,j\right)}{\mathrm{\Δ}x}+c_2\frac{{\mathrm{\σ}}_{xy}\left(i+2,j\right)-{\mathrm{\σ}}_{xy}\left(i-1,j\right)}{\mathrm{\Δ}x}$ Represent that the shearing stress component of the medium particle of described Two-phase Anisotropic Media is arranged on the whole net point place in staggered-mesh;

Wherein, σ _xyrepresent x direction stress component in the y-direction, c ₁and c ₂representation space difference coefficient respectively, i, j get xyz respectively, and x represents horizontal direction, and y represents vertical direction, and z represents vertical direction.

8. device according to claim 6, is characterized in that, describedly arranges module, for:

Pass through formula $\left\{\begin{array}{c}{\mathrm{\ρ}}_{11}\frac{{\∂}^2{w}_i}{\∂t^2}+{\mathrm{\ρ}}_{12}\frac{{\∂}^2{W}_i}{\∂t^2}={\mathrm{\τ}}_{{\mathrm{ij}}^{\′}j}+b_{ij}\left(\frac{\∂W_j}{\∂t}-\frac{\∂w_j}{\∂t}\right)\\ {\mathrm{\ρ}}_{12}\frac{{\∂}^2{w}_i}{\∂t^2}+{\mathrm{\ρ}}_{22}\frac{{\∂}^2{W}_i}{\∂t^2}=s_{\′i}-b_{ij}\left(\frac{\∂W_j}{\∂t}-\frac{\∂w_j}{\∂t}\right)\end{array}\right.$ Represent the mode of motion of described Two-phase Anisotropic Media;

Wherein, w _jand W _jrepresent the component of displacement in j direction of solid phase and stream phase respectively, τ _ij'jfor the solid phase components of stress are at the local derviation in j direction, b _ijdissipation factor during solid phase motion relative to stream, s represents the stress acting on stream phase, and i, j get xyz respectively, and x represents horizontal direction, and y represents vertical direction, and z represents vertical direction.

9. device according to claim 6, is characterized in that, described computing module, for:

Pass through formula $\left\{\begin{array}{c}\frac{\∂{\mathrm{\υ}}_y}{\∂t}=\left(D_2+D_3\right)b_{22}\left({\mathrm{\υ}}_y-V_y\right)-D_3\left(\frac{\∂{\mathrm{\τ}}_{xy}}{\∂x}+\frac{\∂{\mathrm{\τ}}_{yz}}{\∂z}\right)\\ \frac{\∂V_y}{\∂t}=-\left(D_1+D_2\right)b_{22}\left({\mathrm{\υ}}_y-V_y\right)+D_2\left(\frac{\∂{\mathrm{\τ}}_{xy}}{\∂x}+\frac{\∂{\mathrm{\τ}}_{yz}}{\∂z}\right)\end{array}\right.$ Represent described vertical direction speed component;

Wherein, ν _y, V _yrepresent solid-phase component and stream phase constituent vertically speed component respectively, τ _ijrepresent the component of i direction stress along j direction, D _ifor the multi-term expression about density, b ₂₂represent dissipation factor, i, j all represent x, y and z, x represent horizontal direction, and y represents vertical direction, and z represents vertical direction.

10. device according to claim 9, is characterized in that, described computing module, for:

Pass through formula

$\left\{\begin{array}{c}{\mathrm{\υ}}_{yi,j}^{n+\frac{1}{2}}={\mathrm{\υ}}_{yi,j}^{n-\frac{1}{2}}-D_{3i,j}\left(Q_x^{-}\left({\mathrm{\τ}}_{xyi,j}^n\right)+Q_z^{-}\left({\mathrm{\τ}}_{yzi,j}^n\right)\right)+\mathrm{\Δ}t\left(D_{2i,j}+D_{3i,j}\right)b_{22i,j}\left({\mathrm{\υ}}_{yi,j}^{n-\frac{1}{2}}-V_{yi,j}^{n-\frac{1}{2}}\right)\\ V_{yi,j}^{n+\frac{1}{2}}=V_{yi,j}^{n-\frac{1}{2}}+D_{2i,j}\left(Q_x^{-}\left({\mathrm{\τ}}_{xyi,j}^n\right)+Q_z^{-}\left({\mathrm{\τ}}_{yzi,j}^n\right)\right)-\mathrm{\Δ}t\left(D_{1i,j}+D_{2i,j}\right)b_{22i,j}\left({\mathrm{\υ}}_{yi,j}^{n-\frac{1}{2}}-V_{yi,j}^{n-\frac{1}{2}}\right)\end{array}\right.$

Calculate the discretize data of vertical direction speed component;

Wherein, Q representation space difference operator, Δ t represents sampling time interval, ν _y, V _yrepresent solid-phase component and stream phase constituent vertically speed component respectively, τ _ijrepresent the component of i direction stress along j direction, D _ifor the multi-term expression about density, b ₂₂represent dissipation factor.