CN102841374A - Pseudo three-dimensional fast microseism forward modeling method based on scanning surface forward modeling - Google Patents
Pseudo three-dimensional fast microseism forward modeling method based on scanning surface forward modeling Download PDFInfo
- Publication number
- CN102841374A CN102841374A CN2012103078136A CN201210307813A CN102841374A CN 102841374 A CN102841374 A CN 102841374A CN 2012103078136 A CN2012103078136 A CN 2012103078136A CN 201210307813 A CN201210307813 A CN 201210307813A CN 102841374 A CN102841374 A CN 102841374A
- Authority
- CN
- China
- Prior art keywords
- tour
- microearthquake
- whilst
- vertical plane
- net
- 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
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000005553 drilling Methods 0.000 claims description 15
- 230000000644 propagated effect Effects 0.000 claims description 7
- 230000001902 propagating effect Effects 0.000 claims description 6
- 230000009466 transformation Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 12
- 238000012544 monitoring process Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 239000004215 Carbon black (E152) Substances 0.000 description 5
- 229930195733 hydrocarbon Natural products 0.000 description 5
- 125000001183 hydrocarbyl group Chemical group 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Toys (AREA)
Abstract
A pseudo three-dimensional fast microseism forward modeling method based on scanning surface forward modeling includes that a three-dimensional (3D) horizontal isotropic anisotropic model is set up and provided with a coordinate system (X, Y, Z), a microseism origin is disposed inside the anisotropic model, and a plurality of seismic geophones are arranged on the upper surface of the anisotropic model; a two-dimensional (2D) vertical plane is set up for the microseism origin and the seismic geophones so that the microseism origin and the seismic geophones are all arranged on the 2D vertical plane; the 2D vertical plane is vertically projected in a horizontal coordinate system (X, Y) of the anisotropic model, and according to horizontal coordinates of the microseism origin and the seismic geophones in the horizontal coordinate system (X, Y), an azimuth angle theta is obtained; according to the azimuth angle theta and the coordinate transformation formula, the horizontal coordinate system (X, Y) is rotated and then a first new coordinate system X<1> Y<1> is obtained; a Z-axis of the coordinate system (X, Y, Z) and an X<1>-axis of the first new coordinate system X<1> Y<1> are used to form a second new coordinate system (X1, Z); and first arrival travel-times of microseismic waves generated from the microseism origin to every seismic geophone are calculated in the second new coordinate system (X<1>, Z).
Description
Technical field
The present invention relates to the petroleum gas field of seismic exploration, specifically, relate to a kind of pseudo-three-dimensional microearthquake forward modeling method fast of just drilling, be mainly used in the microearthquake monitoring of oil seismic exploration based on scanning plane.
Background technology
Microearthquake pressure break monitoring technology is the important technology of the fast-developing in recent years solution low-permeability oil gas reservoir development that gets up.Although carried out the monitoring of microearthquake pressure break the seventies in 20th century, at that time just as experimental.Along with further attention and the demand of country to energy development, the microearthquake monitoring is played crucial effects with research to the exploitation of unconventional hydrocarbon-bearing pools such as tight sand, shale.China has at least up to ten thousand mouthfuls of wells need do the monitoring of microearthquake pressure break every year.The microearthquake technology is to analyze from the microearthquake signal (at present mainly refer to microearthquake first arrival whilst on tour) of fractured well in fracturing process through what the seismoreceiver in the adjacent well was received, describes the geometric distributions and the fluid migration characteristic (the for example height in crack, length and orientation) of crack growth in the fracturing process.These information can be optimized FRACTURING DESIGN and improve Reservoir Management, thereby improve the production capacity of oil gas field.
Just drilling the technological part that has two keys in the research in microearthquake: microearthquake modelling and microearthquake first arrival whilst on tour are just being drilled calculating.Because descend the media property in hydrocarbon-bearing pool zone (even the whole earth's crust is inner) to show as complicated anisotropic character truly, therefore original microearthquake isotropic model can not accurately be described the physical attribute in hydrocarbon-bearing pool zone.Adopt original isotropic model if can not develop the microearthquake anisotropic model, then be difficult to effectively utilize the first arrival whilst on tour data (these data loggings are from the hydrocarbon-bearing pool with anisotropic characteristic) in the real data.The microearthquake source location and the fracture spaces geometric distributions that may lead to errors based on the result of calculation of isotropic model.
Just drilling the efficient that calculating has determined that microearthquake focus location and fracture distribution calculate based on the microearthquake first arrival whilst on tour of three-dimensional anisotropic medium model.Usually, three-dimensional microearthquake ripple first arrival whilst on tour is just being drilled to calculate and is being had the slow problem of counting yield.How improving the arithmetic speed of microearthquake ripple first arrival whilst on tour and keeping corresponding precision is an important subject in the microearthquake research.
Therefore, the present invention proposes the method that a kind of pseudo-three-dimensional microearthquake of just drilling based on scanning plane fast of anisotropy is just being drilled.
Summary of the invention
Microearthquake ripple first arrival whilst on tour is just being drilled and calculated is current basis and prerequisite of carrying out microearthquake research.Effectively microearthquake first arrival whilst on tour calculates and can guarantee carrying out smoothly of the research of microearthquake later stage (the for example focus of microearthquake location).In addition, just drill calculating based on the microearthquake ripple whilst on tour of three-dimensional anisotropic medium and having the slow problem of counting yield.Therefore, propose the technology that a kind of pseudo-three-dimensional microearthquake of just drilling based on scanning plane is fast just being drilled, underground geologic bodies is assumed to be anisotropic medium.This technology have calculate microearthquake first arrival whilst on tour accurately, advantages such as fast operation and calculation stability, can improve development and utilization to hydrocarbon-bearing pool.The present invention relates to the methods such as foundation, coordinate system transformation and numerical simulation of microearthquake model.
According to an aspect of the present invention, a kind of pseudo-three-dimensional microearthquake forward modeling method fast of just drilling based on scanning plane is provided, said method comprising the steps of: set up three-dimensional transversely isotropic anisotropic model; The coordinate system of anisotropic model is (X; Y, Z), anisotropic model is divided into a plurality of square grids; Set inside at anisotropic model has the microearthquake focus, on the upper surface of anisotropic model, is laid with a plurality of seismoreceivers; Microearthquake focus and each seismoreceiver are set up two-dimentional vertical plane, and two-dimentional vertical plane is vertical with the upper surface of anisotropic model, makes microearthquake focus and said each seismoreceiver be positioned on the said two-dimentional vertical plane; With two-dimentional vertical plane vertical projection to the horizontal coordinates of anisotropic model (X, Y) on, according to microearthquake focus and seismoreceiver horizontal coordinates (X, the horizontal coordinate under Y) obtains azimuth angle theta; (X Y) is rotated, and obtains the first new coordinate system (X to horizontal coordinates according to azimuth angle theta and coordinate transform formula
1, Y
1), make microearthquake focus and seismoreceiver be positioned at the first new coordinate system (X
1, Y
1) X
1On the coordinate axis; Utilize coordinate system (X, Y, Z axle Z) and the first new coordinate system (X
1, Y
1) X
1Axle forms the second new coordinate system (X
1, Z); At the second new coordinate system (X
1, Z) under, calculate the microearthquake ripple that sends from the microearthquake focus first arrival whilst on tour to each seismoreceiver point.
Said anisotropic model can comprise the information about the anisotropic parameters of moulded dimension, size of mesh opening, each net point, microearthquake source location and seismoreceiver position.
Said anisotropic parameters can comprise vertical phase velocity and the Thomsen parameter ε and the δ of microearthquake ripple.
The coordinate transform formula can be:
The microearthquake ripple that can adopt the eikonal equation finite difference method to calculate to send from the microearthquake focus is to the first arrival whilst on tour of each seismoreceiver point.
Adopting the eikonal equation finite difference method to calculate the microearthquake ripple can comprise to the step of the first arrival whilst on tour of each seismoreceiver point: (1) is divided into a plurality of square nets with two-dimentional vertical plane, and microearthquake focus and seismoreceiver are positioned on the net point of two-dimentional vertical plane; (2) with the perpendicular line at the microearthquake focus on two-dimentional vertical plane place as the vertical edges boundary line, with the horizontal line at the place of the microearthquake focus on the two-dimentional vertical plane as the horizontal sides boundary line; (3) calculate the whilst on tour that the microearthquake ripple arrives each net point on horizontal sides boundary line and the vertical edges boundary line according to time of origin, mesh width and the microearthquake wave propagation velocity of microearthquake focus; (4) calculate the slope that the microearthquake ripple is propagated according to the whilst on tour that lays respectively at two net points adjacent on vertical edges boundary line and the horizontal sides boundary line in the time of origin of microearthquake focus and the said a plurality of square net with the microearthquake focus; (5) slope of propagating according to the microearthquake ripple obtains the phase angle that the microearthquake ripple is propagated through arctan function; (6) phase angle of propagating according to the microearthquake ripple, vertical phase velocity and the Thomsen parameter ε and the δ of microearthquake ripple calculate the microearthquake phase velocity of wave; (7) according to the whilst on tour of said two net points adjacent in the time of origin of the width of microearthquake phase velocity of wave, square net, microearthquake focus and the square net with the microearthquake focus; Calculate the whilst on tour of another net point in said a plurality of square nets of microearthquake ripple to the two-dimentional vertical plane, said another net point is adjacent with said two net points and not on horizontal sides boundary line and vertical edges boundary line; (8) according to the identical mode of step (7); On two-dimentional vertical plane, carry out iterative computation from top to bottom and from left to right from the microearthquake focus; For not other net point on horizontal sides boundary line and vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model top on the two-dimentional vertical plane with said other net point; (9) according to the identical mode of step (7); On two-dimentional vertical plane, carry out iterative computation from top to bottom and from right to left from the microearthquake focus; For not other net point on horizontal sides boundary line and vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model top on the two-dimentional vertical plane with said other net point; (10) according to the identical mode of step (7); On two-dimentional vertical plane, carry out iterative computation from top to bottom and from left to right from the microearthquake focus; For not other net point on horizontal sides boundary line and vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model bottom on the two-dimentional vertical plane with said other net point; (11) according to the identical mode of step (7); On two-dimentional vertical plane, carry out iterative computation from top to bottom and from right to left from the microearthquake focus; For not other net point on horizontal sides boundary line and vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated with said other net point; Till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model bottom on the two-dimentional vertical plane, accomplish the scanning first time of two-dimentional vertical plane thus; (12) according to the identical mode of step (7); The vertical edges boundary line from the two-dimentional vertical plane and the net point of anisotropic model roof intersection carry out iterative computation from top to bottom and from left to right; With the first arrival whilst on tour of the anisotropic model top net lattice point of previous calculating as initial value; For not other net point on the vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model bottom on the two-dimentional vertical plane with said other net point; (13) according to the identical mode of step (7); The vertical edges boundary line from the two-dimentional vertical plane and the net point of anisotropic model roof intersection carry out iterative computation from top to bottom and from right to left; With the first arrival whilst on tour of the anisotropic model top net lattice point of previous calculating as initial value; For not other net point on the vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated with said other net point; Till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model bottom on the two-dimentional vertical plane, accomplish the scanning second time of two-dimentional vertical plane thus; (14) according to the identical mode of step (7); The crossing net point in vertical edges boundary line and anisotropic model bottom from two-dimentional vertical plane carries out iterative computation from top to bottom and from left to right; With the first arrival whilst on tour of the anisotropic model bottom web lattice point of previous calculating as initial value; For not other net point on the vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model top on the two-dimentional vertical plane with said other net point; (15) according to the identical mode of step (7); The crossing net point in vertical edges boundary line and anisotropic model bottom from two-dimentional vertical plane carries out iterative computation from top to bottom and from right to left; With the first arrival whilst on tour of the anisotropic model bottom web lattice point of previous calculating as initial value; For not other net point on the vertical edges boundary line on the two-dimentional vertical plane, calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model top on the two-dimentional vertical plane with said other net point; Accomplish the scanning for the third time of two-dimentional vertical plane thus; Wherein, in the scanning process each time of two-dimentional vertical plane, the whilst on tour on each net point that calculates compares with the whilst on tour that previous scanning is calculated on the resulting said net point; Keep minimum whilst on tour, obtain the first arrival whilst on tour that the microearthquake ripple arrives seismoreceiver thus.
In above-mentioned steps, through three scannings of two-dimentional vertical plane, get the minimum whilst on tour of iterative computation at last, can further improve microearthquake and just drilling precision.
Can be with the time of origin of microearthquake focus initial whilst on tour as the net point on the horizontal sides boundary line; Iterative computation microearthquake ripple arrives the whilst on tour of the net point on the horizontal sides boundary line; Wherein, For the ad hoc networks lattice point that is positioned on the horizontal sides boundary line, calculate the whilst on tour that the microearthquake ripple arrives the said ad hoc networks lattice point on the horizontal sides boundary line according to whilst on tour, mesh width and the microearthquake wave propagation velocity of the net point adjacent on the horizontal sides boundary line that has calculated with said ad hoc networks lattice point.
Can be with the time of origin of microearthquake focus initial whilst on tour as the net point on the vertical edges boundary line; Arrive the whilst on tour of the net point on the vertical edges boundary line with iterative computation microearthquake ripple; Wherein, For the ad hoc networks lattice point that is positioned on the vertical edges boundary line, calculate the whilst on tour that the microearthquake ripple arrives the said ad hoc networks lattice point on the vertical edges boundary line according to whilst on tour, mesh width and the microearthquake wave propagation velocity of the net point adjacent on the vertical edges boundary line that has calculated with said ad hoc networks lattice point.
Description of drawings
In conjunction with the drawings, from the description of following embodiment, the present invention these and/or others and advantage will become clear, and are easier to understand, wherein:
Fig. 1 is the pseudo-three-dimensional process flow diagram of microearthquake forward modeling method fast of anisotropy of just drilling based on scanning plane according to of the present invention;
Fig. 2 shows according to the microearthquake focus of the transversely isotropic anisotropic model of three-dimensional of the present invention and the synoptic diagram of seismoreceiver relation;
Fig. 3 shows the (x through microearthquake focus S
0, y
0, z
0) and seismoreceiver R (x
1, y
1, z
1) two-dimentional vertical plane;
Fig. 4 is the subpoint S that connects microearthquake focus and seismoreceiver
1And R
1Azimuthal synoptic diagram;
Fig. 5 is that microearthquake focus and seismoreceiver are at the first new coordinate system (X
1, Y
1) under the signal W figure of position;
Fig. 6 is that microearthquake focus and seismoreceiver are at the second new coordinate system (X
1, the synoptic diagram of the position under Z);
Fig. 7 is the initialized synoptic diagram of whilst on tour;
Fig. 8 illustrates the geometric representation of microearthquake ripple first arrival whilst on tour;
Fig. 9 illustrates the synoptic diagram of three-dimensional five layers anisotropic model;
Figure 10 is the isogram of the first arrival whilst on tour in the anisotropic model;
Figure 11 is the curve of the first arrival whilst on tour that obtains of face of land survey line.
Embodiment
At the pseudo-three-dimensional microearthquake forward modeling method fast of anisotropy of just drilling according to the present invention,, obtain many two-dimentional vertical planes through microearthquake focus and each seismoreceiver reception geometry arrangement (being arranged in the face of land or well) are divided into groups based on scanning plane.To each two-dimentional vertical plane, through vertical projection vertical plane is projected to surface level, at this moment on surface level straight line.Further be rotated coordinate system transformation, utilize anisotropy eikonal equation method of finite difference to calculate microearthquake first arrival whilst on tour then.In computation process, model top and model bottom design three scanning calculative strategies, the whilst on tour of guaranteeing to obtain is for minimum.
Describe the pseudo-three-dimensional microseism forward modeling method fast of anisotropy of just drilling based on scanning plane in detail according to of the present invention below with reference to Fig. 1.Fig. 1 is the pseudo-three-dimensional process flow diagram of microearthquake forward modeling method fast of anisotropy of just drilling based on scanning plane according to of the present invention.
With reference to Fig. 1, in step 101, set up three-dimensional transversely isotropic anisotropic model, coordinate system be (X, Y, Z).Anisotropic model is divided into a plurality of square grids; Set inside at anisotropic model has the microearthquake focus, is being laid with a plurality of seismoreceivers (each seismoreceiver is positioned on the respective grid points on the upper surface of anisotropic model) on the upper surface (the expression face of land) of anisotropic model.Said anisotropic model comprises the information such as anisotropic parameters, microearthquake source location and seismoreceiver position of moulded dimension, size of mesh opening, each net point.Anisotropic parameters comprises the vertical phase velocity v of microearthquake ripple
0And Thomsen parameter (ε and δ).
Fig. 2 shows according to the synoptic diagram of the microearthquake focus of the transversely isotropic anisotropic model of three-dimensional of the present invention (microseism point) with the seismoreceiver relation.In Fig. 2, S (x
0, y
0, z
0) expression microearthquake focus, R (x
1, y
1, z
1) the expression seismoreceiver.The microearthquake ripple that the microearthquake focus sends is received and detects by seismoreceiver.
In step 102; Microearthquake focus and each seismoreceiver (being arranged in the face of land or well) are set up two-dimentional vertical plane; The two dimension vertical plane is vertical with the upper surface of anisotropic model, makes microearthquake focus and seismoreceiver be positioned on this two dimension vertical plane.Fig. 3 shows one through microearthquake focus S (x
0, y
0, z
0) and seismoreceiver R (x
1, y
1, z
1) two-dimentional vertical plane, that is, and microearthquake focus S (x
0, y
0, z
0) and seismoreceiver R (x
1, y
1, z
1) be positioned on this two dimension vertical plane.
In step 103, with two-dimentional vertical plane vertical projection to the horizontal coordinates of anisotropic model (X, Y) on, (X, the subpoint coordinate on Y) is respectively S at horizontal coordinates for microearthquake focus and seismoreceiver
1(x
0, y
0) and R
1(x
1, y
1), according to the horizontal coordinate of microearthquake focus and seismoreceiver, obtain azimuth angle theta then.
The position angle
Fig. 4 is the subpoint S that connects microearthquake focus and seismoreceiver
1And R
1Azimuthal synoptic diagram.
In step 104, (X Y) is rotated, and obtains the first new coordinate system (X to the horizontal coordinates of anisotropic model according to azimuth angle theta and coordinate transform formula
1, Y
1), make microearthquake focus and seismoreceiver be positioned at the first new coordinate system (X
1, Y
1) X
1On the coordinate axis.The coordinate transform formula is following:
Through coordinate transform, microearthquake focus and seismoreceiver are at the first new coordinate system (X
1, Y
1) under coordinate be divided into
With
Fig. 5 is that microearthquake focus and seismoreceiver are at the first new coordinate system (X
1, Y
1) under the synoptic diagram of position.
In step 105, because microearthquake focus and seismoreceiver are positioned at the first new coordinate system (X
1, Y
1) X
1On the axle, so utilize coordinate system (X, Y, Z axle Z) and the first new coordinate system (X
1, Y
1) X
1Axle further forms new vertical coordinate system (X
1, Z), i.e. the second new coordinate system (X
1, Z).At the second new coordinate system (X
1, Z) under, the coordinate of microearthquake focus and seismoreceiver becomes
With
Fig. 6 be for microearthquake focus and seismoreceiver at the second new coordinate system (X
1, the synoptic diagram of the position under Z).
In step 106, at the second new coordinate system (X
1, Z) under, calculate from the microearthquake focus
The microearthquake ripple that sends is to each seismoreceiver point
The first arrival whilst on tour.Eikonal equation can be used to compute the finite difference method from the micro-seismic source
issued to each micro-seismic geophone point
The first break when traveling.It is fast and stable that the eikonal equation finite difference method is compared other method computing velocity.
Described in detail below eikonal equation using finite difference method to calculate from the micro-seismic source
issued to each micro-seismic geophone point
The initial step to travel.
Particularly, the step that adopts the eikonal equation finite difference method to calculate the first arrival whilst on tour can comprise: (1) is divided into a plurality of square nets with two-dimentional vertical plane, and microearthquake focus and seismoreceiver are positioned on the net point of two-dimentional vertical plane; (2) with the perpendicular line at the microearthquake focus on two-dimentional vertical plane place as the vertical edges boundary line, with the horizontal line at the place of the microearthquake focus on the two-dimentional vertical plane as the horizontal sides boundary line; (3) calculate the whilst on tour that the microearthquake ripple arrives each net point on horizontal sides boundary line and the vertical edges boundary line according to time of origin, mesh width and the microearthquake wave propagation velocity of microearthquake focus; (4) calculate the slope that the microearthquake ripple is propagated according to the whilst on tour that lays respectively at two net points adjacent on vertical edges boundary line and the horizontal sides boundary line in the time of origin of microearthquake focus and the said a plurality of square net with the microearthquake focus; (5) slope of propagating according to the microearthquake ripple obtains the phase angle that the microearthquake ripple is propagated through arctan function; (6) phase angle of propagating according to the microearthquake ripple, vertical phase velocity and the Thomsen parameter ε and the δ of microearthquake ripple calculate the microearthquake phase velocity of wave; (7) according to the whilst on tour of said two net points adjacent in the time of origin of the width of microearthquake phase velocity of wave, square net, microearthquake focus and the square net with the microearthquake focus; Calculate the whilst on tour of another net point in said a plurality of square nets of microearthquake ripple to the two-dimentional vertical plane, said another net point is adjacent with said two net points and not on horizontal sides boundary line and vertical edges boundary line; (8) according to the identical mode of step (7); On two-dimentional vertical plane, carry out iterative computation from top to bottom and from left to right from the microearthquake focus; For not other net point on horizontal sides boundary line and vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model top on the two-dimentional vertical plane with said other net point; (9) according to the identical mode of step (7); On two-dimentional vertical plane, carry out iterative computation from top to bottom and from right to left from the microearthquake focus; For not other net point on horizontal sides boundary line and vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model top on the two-dimentional vertical plane with said other net point; (10) according to the identical mode of step (7); On two-dimentional vertical plane, carry out iterative computation from top to bottom and from left to right from the microearthquake focus; For not other net point on horizontal sides boundary line and vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model bottom on the two-dimentional vertical plane with said other net point; (11) according to the identical mode of step (7); On two-dimentional vertical plane, carry out iterative computation from top to bottom and from right to left from the microearthquake focus; For not other net point on horizontal sides boundary line and vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated with said other net point; Till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model bottom on the two-dimentional vertical plane, accomplish the scanning first time of two-dimentional vertical plane thus; (12) according to the identical mode of step (7); The vertical edges boundary line from the two-dimentional vertical plane and the net point of anisotropic model roof intersection carry out iterative computation from top to bottom and from left to right; With the first arrival whilst on tour of the anisotropic model top net lattice point of previous calculating as initial value; For not other net point on the vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model bottom on the two-dimentional vertical plane with said other net point; (13) according to the identical mode of step (7); The vertical edges boundary line from the two-dimentional vertical plane and the net point of anisotropic model roof intersection carry out iterative computation from top to bottom and from right to left; With the first arrival whilst on tour of the anisotropic model top net lattice point of previous calculating as initial value; For not other net point on the vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated with said other net point; Till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model bottom on the two-dimentional vertical plane, accomplish the scanning second time of two-dimentional vertical plane thus; (14) according to the identical mode of step (7); The crossing net point in vertical edges boundary line and anisotropic model bottom from two-dimentional vertical plane carries out iterative computation from top to bottom and from left to right; With the first arrival whilst on tour of the anisotropic model bottom web lattice point of previous calculating as initial value; For not other net point on the vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model top on the two-dimentional vertical plane with said other net point; (15) according to the identical mode of step (7); The crossing net point in vertical edges boundary line and anisotropic model bottom from two-dimentional vertical plane carries out iterative computation from top to bottom and from right to left; With the first arrival whilst on tour of the anisotropic model bottom web lattice point of previous calculating as initial value; For not other net point on the vertical edges boundary line on the two-dimentional vertical plane, calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model top on the two-dimentional vertical plane with said other net point; Accomplish the scanning for the third time of two-dimentional vertical plane thus; Wherein, in the scanning process each time of two-dimentional vertical plane, the whilst on tour on each net point that calculates compares with the whilst on tour that previous scanning is calculated on the resulting said net point; Keep minimum whilst on tour, obtain the first arrival whilst on tour that the microearthquake ripple arrives seismoreceiver thus.
In above-mentioned steps, through three scannings of two-dimentional vertical plane, get the minimum whilst on tour of iterative computation at last, can further improve microearthquake and just drilling precision.
Can the time of origin of the microearthquake focus initial whilst on tour as the net point on the horizontal sides boundary line be arrived the whilst on tour of the net point on the horizontal sides boundary line with iterative computation microearthquake ripple; For be positioned on the horizontal sides boundary line the ad hoc networks lattice point (correspondingly; Microearthquake phase of wave angle is 0 °), calculate the whilst on tour that the microearthquake ripple arrives the said ad hoc networks lattice point on the horizontal sides boundary line according to whilst on tour, mesh width and the microearthquake wave propagation velocity of the net point adjacent on the horizontal sides boundary line that has calculated with said ad hoc networks lattice point.
Likewise, can the time of origin of the microearthquake focus initial whilst on tour as the net point on the vertical edges boundary line be arrived the whilst on tour of the net point on the vertical edges boundary line with iterative computation microearthquake ripple; For be positioned on the vertical edges boundary line the ad hoc networks lattice point (correspondingly; Microearthquake phase of wave angle is 90 °), calculate the whilst on tour that the microearthquake ripple arrives the said ad hoc networks lattice point on the vertical edges boundary line according to whilst on tour, mesh width and the microearthquake wave propagation velocity of the net point adjacent on the vertical edges boundary line that has calculated with said ad hoc networks lattice point.
Fig. 7 is the initialized synoptic diagram of whilst on tour.As shown in Figure 7, suppose that the microearthquake focus is positioned at T
0Net point, T
0Time of origin for the microearthquake focus.The mesh width h and the microearthquake wave propagation velocity of two dimension vertical plane are known, then can utilize equality T
1=T
0+ hs (0 °) calculates the first initial whilst on tour T
1, and the like can calculate the whilst on tour (black triangle is represented) that all L1 lines (horizontal sides boundary line) are gone up net point.S (0 °) refers to the slowness on the horizontal direction.Can utilize T
2=T
0+ hs (90 °) calculates the second initial whilst on tour T
2, and the like can calculate the whilst on tour (black circle is represented) that all C0 lines (vertical edges boundary line) are gone up net point.S (90 °) refers to the slowness on the vertical direction.More than the whilst on tour of each net point just as initial value, can in net point scanning process subsequently, be upgraded by iteration.
Specifically describe the iterative computation of the whilst on tour of each net point below in conjunction with accompanying drawing and formula.The T of front
0And the T that calculates
1And T
2Can be used as the whilst on tour initial value of iterative computation.
Suppose anisotropic parameters (that is vertical phase velocity v of microearthquake ripple, of microearthquake ripple
0And Thomsen parameter ε and δ) in each square net unit, be constant, and hypothesis microearthquake ripple outwards propagates with the plane wave form according to phase velocity (V (α)), and wherein, α is microearthquake phase of wave angle.In case phase angle α is determined, then can calculate phase velocity V (α), further calculate the first arrival whilst on tour that the microearthquake ripple arrives each seismoreceiver then.
Under the weak-moderate anisotropy situation of Method in Transverse Isotropic Medium, the calculation expression of microearthquake phase velocity of wave V (α) is:
V(α)=v
0(1+δsin
2α+(ε-δ)sin
4α) (2)
Can find out that the microearthquake ripple is different along the phase velocity of each direction of propagation.
The eikonal equation of two dimension isotropic medium is following:
Wherein, τ express time.
The approximate form of eikonal equation that can obtain the Two-Dimensional Anisotropic medium from the eikonal equation of two-dimentional isotropic medium is following:
In the eikonal equation finite difference method, arrive the unknown whilst on tour (said another net point is not on horizontal sides boundary line and vertical edge boundary line) of another net point adjacent of two-dimentional vertical plane to the whilst on tour calculating microearthquake ripple of three net points that are positioned at two-dimentional vertical plane according to the microearthquake ripple with said three net points.
As shown in Figure 7, t
0, t
1And t
2Be the known whilst on tour of known microearthquake ripple to three net points of two-dimentional vertical plane, the microearthquake ripple is asked for waiting to another one net point whilst on tour t, and the square net width of two-dimentional vertical plane is h.
Slope through the propagation of computes microearthquake ripple:
Obtain the phase angle α that the microearthquake ripple is propagated through arctan function then.Phase angle α is updated to phase velocity computing formula (2), thereby obtains phase velocity V (α).Further with phase velocity V (α) substitution following formula:
Thus, can draw the whilst on tour t of microearthquake ripple to the another one net point.
Then according to identical mode; On two-dimentional vertical plane, carry out iterative computation from top to bottom and from left to right from the microearthquake focus; For not other net point on horizontal sides boundary line and vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model top on the two-dimentional vertical plane with said other net point; According to identical mode; On two-dimentional vertical plane, carry out iterative computation from top to bottom and from right to left from the microearthquake focus; For not other net point on horizontal sides boundary line and vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model top on the two-dimentional vertical plane with said other net point; According to identical mode; On two-dimentional vertical plane, carry out iterative computation from top to bottom and from left to right from the microearthquake focus; For not other net point on horizontal sides boundary line and vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model bottom on the two-dimentional vertical plane with said other net point; According to identical mode; On two-dimentional vertical plane, carry out iterative computation from top to bottom and from right to left from the microearthquake focus; For not other net point on horizontal sides boundary line and vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated with said other net point; Till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model bottom on the two-dimentional vertical plane, accomplish the scanning first time of two-dimentional vertical plane thus.
Then according to identical mode; The vertical edges boundary line from the two-dimentional vertical plane and the net point of anisotropic model roof intersection carry out iterative computation from top to bottom and from left to right; With the first arrival whilst on tour of the anisotropic model top net lattice point of previous calculating as initial value; For not other net point on the vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model bottom on the two-dimentional vertical plane with said other net point; According to identical mode; The vertical edges boundary line from the two-dimentional vertical plane and the net point of anisotropic model roof intersection carry out iterative computation from top to bottom and from right to left; With the first arrival whilst on tour of the anisotropic model top net lattice point of previous calculating as initial value; For not other net point on the vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model bottom on the two-dimentional vertical plane, accomplish the scanning second time of two-dimentional vertical plane thus with said other net point.
Then according to identical mode; The crossing net point in vertical edges boundary line and anisotropic model bottom from two-dimentional vertical plane carries out iterative computation from top to bottom and from left to right; With the first arrival whilst on tour of the anisotropic model bottom web lattice point of previous calculating as initial value; For not other net point on the vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model top on the two-dimentional vertical plane with said other net point; According to identical mode; The crossing net point in vertical edges boundary line and anisotropic model bottom from two-dimentional vertical plane carries out iterative computation from top to bottom and from right to left; With the first arrival whilst on tour of the anisotropic model bottom web lattice point of previous calculating as initial value; For not other net point on the vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model top on the two-dimentional vertical plane, accomplish the scanning for the third time of two-dimentional vertical plane thus with said other net point.
In the scanning process each time of two-dimentional vertical plane; Whilst on tour on each net point that calculates compares with the whilst on tour that previous scanning is calculated on the resulting said net point; Keep minimum whilst on tour, obtain the first arrival whilst on tour that the microearthquake ripple arrives seismoreceiver thus.
As stated, in scanning (iterative computation) for the first time, T capable of using
0, T
1And T
2Respectively as t
0, t
1And t
2Initial value calculate.In addition; For be positioned on horizontal sides boundary line or the vertical edges boundary line the ad hoc networks lattice point (correspondingly; Microearthquake phase of wave angle is 0 ° or 90 °), can calculate the whilst on tour that the microearthquake ripple arrives the said ad hoc networks lattice point on horizontal sides boundary line or the vertical edges boundary line according to whilst on tour, mesh width and the microearthquake wave propagation velocity of the net point adjacent on horizontal sides boundary line that has calculated or the vertical edges boundary line with said ad hoc networks lattice point; For the net point on horizontal sides boundary line and vertical edge boundary line not, can calculate the whilst on tour that the microearthquake ripple arrive said net point according to iterative formula (6), till the scanning of accomplishing whole two-dimentional vertical plane.
The present invention is applied in three-dimensional five layers of anisotropic medium model, sets up three-dimensional transversely isotropic anisotropic model, each layer model anisotropic parameters is as shown in the table:
Table 1 model anisotropic parameter
| The number of plies | v 0 | ε | |
| 1 | 2000 | 0 | 0 |
| 2 | 2400 | 0.05 | 0.05 |
| 3 | 2800 | 0.1 | 0.1 |
| 4 | 3000 | 0.1 | 0.1 |
| 5 | 3500 | 0.15 | 0.15 |
Obtain the second new coordinate system (X through coordinate transform
1, the moulded dimension under Z) is 3000 * 2000, i.e. long 3000 meters, dark 2000 meters (referring to the five layers of anisotropic model of three-dimensional shown in Figure 9) of model.Size of mesh opening is 2 meters, and the microearthquake source location is net point (1000,980), shown in the asterisk among Fig. 9.Each net point of the face of land is provided with a seismoreceiver, and dotted line is represented survey line.Microearthquake forward modeling method through the present invention proposes can calculate the first arrival whilst on tour quickly and accurately.
Figure 10 is the isogram of the first arrival whilst on tour in the anisotropic model, and Figure 11 is the curve of the first arrival whilst on tour that obtains of face of land survey line.
Model experiment shows; The pseudo-three-dimensional positive artistic skills art of just drilling based on scanning plane of microearthquake fast according to the present invention has not only kept the accurate of microearthquake first arrival whilst on tour calculating; And have advantages such as fast operation and calculation stability, greatly improved the bearing accuracy of microearthquake monitoring.
The present invention is applicable to the independent development of microearthquake monitoring software system.The shale gas reservoir has stronger anisotropic character; And to the bearing accuracy important influence; The present invention can be through based on anisotropic forward modeling method, to the later microearthquake monitoring in real time of microearthquake and fine processing, recording geometry design etc. have stronger actual directive function afterwards.
Though the present invention is specifically described with reference to its exemplary embodiment and is shown; But will be understood by those skilled in the art that; Under the situation that does not break away from the spirit and scope of the present invention that are defined by the claims, can carry out the various changes of form and details to it.
Claims (8)
1. pseudo-three-dimensional microearthquake forward modeling method fast of just drilling based on scanning plane may further comprise the steps:
Set up three-dimensional transversely isotropic anisotropic model; The coordinate system of anisotropic model is (X; Y, Z), anisotropic model is divided into a plurality of square grids; Set inside at anisotropic model has the microearthquake focus, and the upper surface of anisotropic model is provided with a plurality of seismoreceivers;
Microearthquake focus and each seismoreceiver are set up two-dimentional vertical plane, and two-dimentional vertical plane is vertical with the upper surface of anisotropic model, makes microearthquake focus and said each seismoreceiver be positioned on the said two-dimentional vertical plane;
With two-dimentional vertical plane vertical projection to the horizontal coordinates of anisotropic model (X, Y) on, according to microearthquake focus and seismoreceiver horizontal coordinates (X, the horizontal coordinate under Y) obtains azimuth angle theta;
(X Y) is rotated, and obtains the first new coordinate system (X to horizontal coordinates according to azimuth angle theta and coordinate transform formula
1, Y
1), make microearthquake focus and seismoreceiver be positioned at the first new coordinate system (X
1, Y
1) X
1On the coordinate axis;
Utilize coordinate system (X, Y, Z axle Z) and the first new coordinate system (X
1, Y
1) X
1Axle forms the second new coordinate system (X
1, Z);
At the second new coordinate system (X
1, Z) under, calculate the microearthquake ripple that sends from the microearthquake focus first arrival whilst on tour to each seismoreceiver point.
2. pseudo-three-dimensional microearthquake forward modeling method fast according to claim 1; Wherein, said anisotropic model comprises the information of anisotropic parameters, microearthquake source location and seismoreceiver position about moulded dimension, size of mesh opening, each net point.
3. pseudo-three-dimensional microearthquake forward modeling method fast according to claim 2, wherein, said anisotropic parameters comprises vertical phase velocity and the Thomsen parameter ε and the δ of microearthquake ripple.
4. pseudo-three-dimensional microearthquake forward modeling method fast according to claim 1, wherein, the coordinate transform formula is:
5. pseudo-three-dimensional microearthquake forward modeling method fast according to claim 1, wherein, the microearthquake ripple that adopts the eikonal equation finite difference method to calculate to send from the microearthquake focus is to the first arrival whilst on tour of each seismoreceiver point.
6. pseudo-three-dimensional microearthquake forward modeling method fast according to claim 5, wherein, adopt the eikonal equation finite difference method to calculate the microearthquake ripple and comprise to the step of the first arrival whilst on tour of each seismoreceiver point:
(1) two-dimentional vertical plane is divided into a plurality of square nets, microearthquake focus and seismoreceiver are positioned on the net point of two-dimentional vertical plane;
(2) with the perpendicular line at the microearthquake focus on two-dimentional vertical plane place as the vertical edges boundary line, with the horizontal line at the place of the microearthquake focus on the two-dimentional vertical plane as the horizontal sides boundary line;
(3) calculate the whilst on tour that the microearthquake ripple arrives each net point on horizontal sides boundary line and the vertical edges boundary line according to time of origin, mesh width and the microearthquake wave propagation velocity of microearthquake focus;
(4) calculate the slope that the microearthquake ripple is propagated according to the whilst on tour that lays respectively at two net points adjacent on vertical edges boundary line and the horizontal sides boundary line in the time of origin of microearthquake focus and the said a plurality of square net with the microearthquake focus;
(5) slope of propagating according to the microearthquake ripple obtains the phase angle that the microearthquake ripple is propagated through arctan function;
(6) phase angle of propagating according to the microearthquake ripple, vertical phase velocity and the Thomsen parameter ε and the δ of microearthquake ripple calculate the microearthquake phase velocity of wave;
(7) based on the whilst on tour of said two mesh points adjacent in the time of origin of the width of microseism phase velocity of wave, square net, microseism focus and the square net with the microseism focus; Calculate the whilst on tour of another mesh point in said a plurality of square nets of microseism ripple to the two-dimentional vertical plane, said another mesh point is adjacent with said two mesh points and not on horizontal sides boundary line and vertical edges boundary line;
(8) according to the identical mode of step (7); On two-dimentional vertical plane, carry out iterative computation from top to bottom and from left to right from the microearthquake focus; For not other net point on horizontal sides boundary line and vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model top on the two-dimentional vertical plane with said other net point;
(9) according to the identical mode of step (7); On two-dimentional vertical plane, carry out iterative computation from top to bottom and from right to left from the microearthquake focus; For not other net point on horizontal sides boundary line and vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model top on the two-dimentional vertical plane with said other net point;
(10) according to the identical mode of step (7); On two-dimentional vertical plane, carry out iterative computation from top to bottom and from left to right from the microearthquake focus; For not other net point on horizontal sides boundary line and vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model bottom on the two-dimentional vertical plane with said other net point;
(11) according to the identical mode of step (7); On two-dimentional vertical plane, carry out iterative computation from top to bottom and from right to left from the microearthquake focus; For not other net point on horizontal sides boundary line and vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated with said other net point; Till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model bottom on the two-dimentional vertical plane, accomplish the scanning first time of two-dimentional vertical plane thus;
(12) according to the identical mode of step (7); The vertical edges boundary line from the two-dimentional vertical plane and the net point of anisotropic model roof intersection carry out iterative computation from top to bottom and from left to right; With the first arrival whilst on tour of the anisotropic model top net lattice point of previous calculating as initial value; For not other net point on the vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model bottom on the two-dimentional vertical plane with said other net point;
(13) according to the identical mode of step (7); The vertical edges boundary line from the two-dimentional vertical plane and the net point of anisotropic model roof intersection carry out iterative computation from top to bottom and from right to left; With the first arrival whilst on tour of the anisotropic model top net lattice point of previous calculating as initial value; For not other net point on the vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated with said other net point; Till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model bottom on the two-dimentional vertical plane, accomplish the scanning second time of two-dimentional vertical plane thus;
(14) according to the identical mode of step (7); The crossing net point in vertical edges boundary line and anisotropic model bottom from two-dimentional vertical plane carries out iterative computation from top to bottom and from left to right; With the first arrival whilst on tour of the anisotropic model bottom web lattice point of previous calculating as initial value; For not other net point on the vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated, till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model top on the two-dimentional vertical plane with said other net point;
(15) according to the identical mode of step (7); The crossing net point in vertical edges boundary line and anisotropic model bottom from two-dimentional vertical plane carries out iterative computation from top to bottom and from right to left; With the first arrival whilst on tour of the anisotropic model bottom web lattice point of previous calculating as initial value; For not other net point on the vertical edges boundary line on the two-dimentional vertical plane; Calculate the whilst on tour of said other net point according to the whilst on tour of three net points adjacent that calculated with said other net point; Till the first arrival whilst on tour that calculates the net point that is positioned at the anisotropic model top on the two-dimentional vertical plane, accomplish the scanning for the third time of two-dimentional vertical plane thus
Wherein, In the scanning process each time of two-dimentional vertical plane; Whilst on tour on each net point that calculates compares with the whilst on tour that previous scanning is calculated on the resulting said net point, keeps minimum whilst on tour, obtains the first arrival whilst on tour that the microearthquake ripple arrives seismoreceiver thus.
7. pseudo-three-dimensional microearthquake forward modeling method fast according to claim 6, wherein, with the time of origin of the microearthquake focus initial whilst on tour as the net point on the horizontal sides boundary line, iterative computation microearthquake ripple arrives the whilst on tour of the net point on the horizontal sides boundary line,
Wherein, For the ad hoc networks lattice point that is positioned on the horizontal sides boundary line, calculate the whilst on tour that the microearthquake ripple arrives the said ad hoc networks lattice point on the horizontal sides boundary line according to whilst on tour, mesh width and the microearthquake wave propagation velocity of the net point adjacent on the horizontal sides boundary line that has calculated with said ad hoc networks lattice point.
8. pseudo-three-dimensional microearthquake forward modeling method fast according to claim 6; Wherein, With the time of origin of microearthquake focus initial whilst on tour, arrive the whilst on tour of the net point on the vertical edges boundary line with iterative computation microearthquake ripple as the net point on the vertical edges boundary line
Wherein, For the ad hoc networks lattice point that is positioned on the vertical edges boundary line, calculate the whilst on tour that the microearthquake ripple arrives the said ad hoc networks lattice point on the vertical edges boundary line according to whilst on tour, mesh width and the microearthquake wave propagation velocity of the net point adjacent on the vertical edges boundary line that has calculated with said ad hoc networks lattice point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210307813.6A CN102841374B (en) | 2012-08-27 | 2012-08-27 | Pseudo three-dimensional fast microseism forward modeling method based on scanning surface forward modeling |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210307813.6A CN102841374B (en) | 2012-08-27 | 2012-08-27 | Pseudo three-dimensional fast microseism forward modeling method based on scanning surface forward modeling |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102841374A true CN102841374A (en) | 2012-12-26 |
| CN102841374B CN102841374B (en) | 2014-12-17 |
Family
ID=47368907
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210307813.6A Active CN102841374B (en) | 2012-08-27 | 2012-08-27 | Pseudo three-dimensional fast microseism forward modeling method based on scanning surface forward modeling |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102841374B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104122578A (en) * | 2013-04-23 | 2014-10-29 | 中国石油化工股份有限公司 | Earthquake excitation simulation method based on rock microstructure and medium elastic parameters |
| CN105893723A (en) * | 2014-10-15 | 2016-08-24 | 长沙矿山研究院有限责任公司 | Rock mass fault gliding plane occurrence calculation method based on microseism event cluster PCA method |
| CN106569279A (en) * | 2015-10-12 | 2017-04-19 | 中国石油化工股份有限公司 | Method for correcting three-dimensional shortest ray tracking path and device thereof |
| CN109725345A (en) * | 2018-11-15 | 2019-05-07 | 中国石油天然气集团有限公司 | A kind of preliminary wave the Forward Modeling and device |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007044459A2 (en) * | 2005-10-07 | 2007-04-19 | Wireless Seismic | Wireless exploration seismic system |
| CN102213769A (en) * | 2010-04-07 | 2011-10-12 | 中国石油天然气集团公司 | Method for determining anisotropic parameters by utilizing data of three-dimensional VSP (Vertical Seismic Profile) |
-
2012
- 2012-08-27 CN CN201210307813.6A patent/CN102841374B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007044459A2 (en) * | 2005-10-07 | 2007-04-19 | Wireless Seismic | Wireless exploration seismic system |
| CN102213769A (en) * | 2010-04-07 | 2011-10-12 | 中国石油天然气集团公司 | Method for determining anisotropic parameters by utilizing data of three-dimensional VSP (Vertical Seismic Profile) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104122578A (en) * | 2013-04-23 | 2014-10-29 | 中国石油化工股份有限公司 | Earthquake excitation simulation method based on rock microstructure and medium elastic parameters |
| CN105893723A (en) * | 2014-10-15 | 2016-08-24 | 长沙矿山研究院有限责任公司 | Rock mass fault gliding plane occurrence calculation method based on microseism event cluster PCA method |
| CN106569279A (en) * | 2015-10-12 | 2017-04-19 | 中国石油化工股份有限公司 | Method for correcting three-dimensional shortest ray tracking path and device thereof |
| CN106569279B (en) * | 2015-10-12 | 2018-08-07 | 中国石油化工股份有限公司 | Method and apparatus for correcting three-dimensional most short ray tracing path |
| CN109725345A (en) * | 2018-11-15 | 2019-05-07 | 中国石油天然气集团有限公司 | A kind of preliminary wave the Forward Modeling and device |
| CN109725345B (en) * | 2018-11-15 | 2020-08-11 | 中国石油天然气集团有限公司 | First-arrival forward modeling method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102841374B (en) | 2014-12-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10795053B2 (en) | Systems and methods of multi-scale meshing for geologic time modeling | |
| CN104977618B (en) | A kind of method evaluated shale gas reservoir and find dessert area | |
| Liu et al. | Three-dimensional seismic ahead-prospecting method and application in TBM tunneling | |
| CN106772577B (en) | Source inversion method based on microseism data and SPSA optimization algorithm | |
| US11175434B2 (en) | Geologic stratigraphy via implicit and jump functions | |
| CN102914791B (en) | Kirchhoff prestack time migration method for processing seismic data of undulating surface | |
| CN104199090B (en) | A kind of rate pattern of ground monitoring microseism positioning builds and method for solving | |
| US20180322232A1 (en) | Representing structural uncertainty in a mesh representing a geological environment | |
| CN103499835A (en) | Method for inverting near-surface velocity model by utilizing preliminary waveforms | |
| CN105277978A (en) | Method and device for determining near-ground-surface speed model | |
| CN105093274A (en) | Hydrofracture focal mechanism inversion method and system | |
| CN102116869A (en) | High-precision prestack domain least square migration seismic imaging technology | |
| CN102636809B (en) | Method for generating spreading angle domain common image point gathers | |
| CN104678434A (en) | Method for predicting storage layer crack development parameters | |
| CN105549077B (en) | The microseism seismic source location method calculated based on multistage multiple dimensioned grid likeness coefficient | |
| CN103576200A (en) | Low signal-to-noise ratio zone shallow wave impedance interface static correction method | |
| CN104730574A (en) | Method for constructing near-surface structural model | |
| CN107817516A (en) | Near surface modeling method and system based on preliminary wave information | |
| CN102841374B (en) | Pseudo three-dimensional fast microseism forward modeling method based on scanning surface forward modeling | |
| CN102877828A (en) | CT (Computed Tomography) imaging method of three-dimensional multi-well combined well land | |
| CN104749630B (en) | Method for constructing microseism monitoring velocity model | |
| CN104199088B (en) | Incident angle gather extraction method and system | |
| CN113376695B (en) | Full waveform inversion method suitable for complex collapse column of coal seam floor | |
| Zhou et al. | An iterative factored topography-dependent eikonal solver for anisotropic media | |
| CN104111476A (en) | Method and device for building formation velocity field |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20180212 Address after: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189 Patentee after: BGP INC., CHINA NATIONAL PETROLEUM Corp. Address before: 610213 No. 1, No. 1, No. 1, Huayang Avenue, Huayang Town, Shuangliu County, Chengdu, Sichuan Patentee before: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20200922 Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No. Patentee after: CHINA NATIONAL PETROLEUM Corp. Patentee after: BGP Inc., China National Petroleum Corp. Address before: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189 Patentee before: BGP Inc., China National Petroleum Corp. |