CN105093319A - Ground micro-seismic static correction method based on three-dimensional seismic data - Google Patents

Ground micro-seismic static correction method based on three-dimensional seismic data Download PDF

Info

Publication number
CN105093319A
CN105093319A CN201410203127.3A CN201410203127A CN105093319A CN 105093319 A CN105093319 A CN 105093319A CN 201410203127 A CN201410203127 A CN 201410203127A CN 105093319 A CN105093319 A CN 105093319A
Authority
CN
China
Prior art keywords
seismic
static correction
correction value
geophone station
ground micro
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
Application number
CN201410203127.3A
Other languages
Chinese (zh)
Other versions
CN105093319B (en
Inventor
胡晓婷
尚新民
王延光
王兴谋
马晓义
赵胜天
崔庆辉
揭景荣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Geophysical Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Geophysical Research Institute of Sinopec Shengli Oilfield Co
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Geophysical Research Institute of Sinopec Shengli Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN201410203127.3A priority Critical patent/CN105093319B/en
Publication of CN105093319A publication Critical patent/CN105093319A/en
Application granted granted Critical
Publication of CN105093319B publication Critical patent/CN105093319B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention discloses a ground micro-seismic static correction method based on three-dimensional seismic data. The method comprises the steps of selecting three-dimensional seismic data capable of including a work area; using a tomographic static correction method for solving a static correction value; as a demodulation point of the three-dimensional seismic data can be connected with a receiving point of a ground micro-seism through a geodetic coordinate, in this way, for a demodulation point of each ground micro-seism, taking the demodulation point as the center, giving a maximum radius and a minimum radius, and taking the seismic demodulation point static correction value in the range of the minimum radius as the static correction value of the micro-seismic demodulation point; and when the range is greater than the minimum radius and less than the maximum radius, using the static correction value of the three-dimensional seismic data for interpolating the static correction value of the micro-seism by the way of inverse ratio distance interpolation. According to the invention, the three-dimensional seismic data is used for solving the static correction value of the ground micro-seism, manual pick-up first breaks are avoided, no priori information is required, only the static correction of the seismic data is required to be well done, and a simple and efficient static correction processing technique is obtained.

Description

Based on the ground micro-seismic static correcting method of 3D seismic data
Technical field
The present invention relates to earthquake and the process of ground micro-seismic data, belong to the static corrections processing of ground micro-seismic data, specifically a kind of ground micro-seismic static correcting method based on 3D seismic data.
Background technology
Because the wave detector of ground micro-seismic monitoring is laid in earth's surface, the first arrival of micro-seismic event can be subject to the impact of surface relief and low velocity layer thickness, velocity variations.This is unfavorable for that microearthquake data improves the treatment effect of signal to noise ratio (S/N ratio), and affects the positioning precision of microearthquake validity event.Therefore, need to carry out static corrections processing before carrying out microearthquake data processing.
Conventional ground micro-seismic static correction is generally utilize the first arrival of pickup perforation validity event to ask for static correction value, and therefore, the levels of precision of first break pickup directly affects the quality of static correction quality.Because ground micro-seismic data SNR is usually very low, even perforation data, conventional automatic first break pickup method is utilized also to be difficult to the needs meeting precision, according to artificial pickup, then there is repeatedly the process of a large amount of amendment, waste time and energy, and each event can the record number of channels of accurate first break picking be limited, the record number of channels making it possible to ask for static correction value is just less.Further, in the process that microearthquake wave is propagated, because the travel path of different event is different, often its static correction value is also different, when the static correction value making to utilize perforation and the strong event of some microearthquake to calculate is applied to and goes in other event, there is the phenomenon of localised jitter through overcorrect funeral affairs part lineups.For these problems, proposing a kind of method utilizing 3D seismic data to ask for ground micro-seismic static correction value is have higher research and apply to be worth.
Summary of the invention
The object of the invention is to provide a kind of ground micro-seismic static correcting method based on 3D seismic data.The method does not need the first break information of a large amount of artificial pickup perforation data, avoids the error in the artificial pick process caused because data signal to noise ratio (S/N ratio) is on the low side, is a kind of static correcting method utilizing 3D seismic data to ask for ground micro-seismic static correction value.
The present invention realizes by following technological means.
Based on the ground micro-seismic static correcting method of 3D seismic data, it is characterized in that: comprise the following steps.
Step 1: the three dimensional seismic data that can comprise this work area is chosen in the ground micro-seismic work area according to gathering.
Step 2: Tomography Velocity inverting is carried out to the three dimensional seismic data chosen.
Step 3: utilize the rate pattern of inverting to carry out tomographic statics.
Step 4: judge whether static correction value tallies with the actual situation, if qualified, performs step 5); Otherwise, repeat step 2-3, until find qualified static correction value.
Step 5: static correction value qualified in selecting step 4.
Step 6: the geophone station choosing a ground micro-seismic, centered by this geophone station, setting least radius R min, maximum radius R max;
If the distance d between a certain 3-D seismics geophone station and the geophone station chosen meets
d≤R min
Then think that the geophone station of this 3-D seismics overlaps with the geophone station of this ground micro-seismic, the static correction value of the two is consistent;
If the distance d between a certain 3-D seismics geophone station and the geophone station chosen meets
R min<d≤R max
Then at R minwith R maxutilize in scope, utilize the method for inverse distance weighting to calculate the static correction value of this ground micro-seismic geophone station;
If the distance d between a certain 3-D seismics geophone station and the geophone station chosen meets
d>R max
Then think this geophone station not in the scope calculating consideration;
According to mentioned above principle, all 3-D seismics geophone stations are chosen one by one and processed;
Step 7: the geophone station choosing other ground micro-seismic one by one, repeats step 6, until all the geophone station of ground micro-seismic completes the calculating of static correction value.
Step 8: the static correction value of all ground micro-seismic geophone stations is applied to ground micro-seismic data.
Further, the concrete methods of realizing of the Tomography Velocity inverting of described step 2 is as follows.
Near surface is separated into grid, and in each grid, speed is equal but unknown.Given initial velocity model carries out first arrival ray tracing, then the primary travel time and real first break time that calculate acquisition is compared.By separating the inversion equation group A △ S=△ b shown in formula 1, calculating the speed correction amount in each grid, and then initial model is modified;
a 11 a 12 . . . a 1 j . . . a 1 J a 21 a 22 . . . a 2 j . . . a 2 J . . . . . . . . . . . . . . . . . . a i 1 a i 2 . . . a ij . . . a iJ . . . . . . . . . . . . . . . . . . a I 1 a I 2 . . . a Ij . . . a IJ &Delta;s 1 &Delta;s 2 . . . &Delta;s j . . . &Delta;s J = &Delta;b 1 &Delta;b 2 . . . &Delta;b i . . . &Delta;b I (formula 1)
Wherein, in inversion equation group, A represents ray spread length matrix within a grid, and △ S is underground slowness correction, and △ b is poor between primary travel time.
Further, described step 2 needs the process of repetition FORWARD AND INVERSE PROBLEMS, approaches the minimum rate pattern of whilst on tour error gradually and stop iteration by successive ignition, and between primary travel time, difference △ b is the difference of true first break time and forward simulation primary travel time.
Further, utilize the rate pattern of inverting to carry out tomographic statics described in step 3 and refer to: after the Tomography Velocity inverting of three dimensional seismic data, carry out the tomographic statics process based on first break time.
Its concrete implementation method is as follows.
Obtain the static correction value that near-surface velocity structure just can calculate shot point and geophone station position later easily.Under earth's surface-consistent prerequisite, think that the ray through earth's surface low velocity layer is perpendicular ray, the correcting value of each physical points is unique.Concrete formula is as follows:
t = t A + t B = - ( &Sigma; i = 1 I &Delta;h i &Delta;v i - h b - h gA V r ) - ( &Sigma; j = 1 J &Delta;h j &Delta;v j - h b - h gB V r )
(formula 2)
In formula, △ h ifor the vertical thickness of i-th grid that ray passes, h bfor uniform reference surface elevation, h gA, h gBfor the elevation on shot point, high speed top that geophone station is corresponding, V rfor replacement velocity, △ v ifor the speed of the grid of i-th after tomographic inversion, I is the grid number of shot point place from earth's surface vertically downward at a high speed top, and J is the grid number of geophone station place from earth's surface vertically downward to top at a high speed; t afor the static correction value of shot point, t bfor the static correction value of geophone station, t is that geophone station B is relative to this static correction value of one of shot point A.
Further, after three dimensional seismic data carries out tomographic statics, generally there are shot statics and receiver static correction amount.When judging whether static correction value tallies with the actual situation in step 4, because the vibration source of ground micro-seismic is in underground, the receiver static correction amount only choosing three dimensional seismic data calculates, and judges whether it tallies with the actual situation.
Further, the concrete methods of realizing of the weighting of inverse distance described in step 6 is:
If interpolation ground micro-seismic geophone station is P (x p, y p, z p), P point minimum with maximum radius within the scope of have known 3-D seismics geophone station Q i(x i, y i, z i), i=1,2 ..., n, utilizes distance weighted inverse ratio method to the static correction value Z of P point pcarry out interpolation.The weighted mean of its Interpolation Principle to be the static correction value of interpolation point be 3-D seismics static correction value in the given radius of point to be inserted, in the size of power and point to be inserted and given radius, the distance dependent of known point, is the inverse of distance k power; That is:
Z p = &Sigma; i = 1 n Z i / d i 2 &Sigma; i = 1 n 1 / d i 2 (formula 3)
Wherein, d ifor the distance in interpolation point and given radius between i-th point.
Further, R in described step 6 minand R maxthe track pitch that arranges according to ground micro-seismic of selection and 3-D seismics arrangement track pitch and determine, namely according to the condition of construction in actual work area, specifically given.
Further, the traversal scope of parameter k is [0,2], parameter d itraversal scope be [R min, R max].
The present invention utilizes the static correction value of 3D seismic data to ask for the static correction value of ground micro-seismic, the first break information of a large amount of artificial pickup perforation data repeatedly can be avoided, more avoid because different event travel path is different, often corrected value different caused with the static correction value that perforation and the strong event of some microearthquake calculate be applied to go in other event time, the problem of the shake existed in local after event lineups correct.The method has that stability is strong, efficiency is high and calculate easy advantage.
Accompanying drawing explanation
Fig. 1 is the ground micro-seismic static correcting method process flow diagram based on 3D seismic data of the present invention.
Fig. 2 is the schematic diagram in the 3-D seismics work area comprising ground micro-seismic work area.
Fig. 3, for utilizing 3D seismic data, carries out the velocity field of tomographic inversion.
Fig. 4 is the receiver static correction amount that model after utilizing Tomography Velocity inverting carries out the three dimensional seismic data that tomographic statics obtains.
Fig. 5 is the elevation curve in ground micro-seismic work area and the static correction value of the ground micro-seismic utilizing the static correction value of 3-D seismics to be tried to achieve by interpolation.
Fig. 6 is the ground micro-seismic section before application ground micro-seismic static correction value.
Fig. 7 is the ground micro-seismic section after applying ground micro-seismic static correction value.
Embodiment
For making above and other object of the present invention, feature and advantage can become apparent, cited below particularly go out preferred embodiment, and coordinate institute's accompanying drawings, elaborate.
Embodiment 1.Based on the ground micro-seismic static correcting method of 3D seismic data, should comprise the following steps based on the ground micro-seismic static correcting method of 3D seismic data.
1) three dimensional seismic data that can comprise this work area is chosen according to the ground micro-seismic work area gathered;
2) Tomography Velocity inverting is carried out to three dimensional seismic data;
The concrete methods of realizing of Tomography Velocity inverting: near surface is separated into grid, in each grid, speed is equal but unknown.Given initial velocity model carries out first arrival ray tracing, then the primary travel time and real first break time that calculate acquisition are compared, calculated the speed correction amount in each grid by solution inversion equation group A △ S=△ b (see formula 1), and then initial model is modified.
a 11 a 12 . . . a 1 j . . . a 1 J a 21 a 22 . . . a 2 j . . . a 2 J . . . . . . . . . . . . . . . . . . a i 1 a i 2 . . . a ij . . . a iJ . . . . . . . . . . . . . . . . . . a I 1 a I 2 . . . a Ij . . . a IJ &Delta;s 1 &Delta;s 2 . . . &Delta;s j . . . &Delta;s J = &Delta;b 1 &Delta;b 2 . . . &Delta;b i . . . &Delta;b I (formula 1)
Wherein, in inversion equation group, A represents ray spread length matrix within a grid, and △ S is underground slowness correction, and △ b is poor between primary travel time.
3), after the Tomography Velocity inverting of three dimensional seismic data, the tomographic statics process based on first break time is carried out;
Concrete implementation method: obtain the static correction value that near-surface velocity structure just can calculate shot point and geophone station position later easily.Under earth's surface-consistent prerequisite, think that the ray through earth's surface low velocity layer is perpendicular ray, the correcting value of each physical points is unique.Concrete formula is as follows:
t = t A + t B = - ( &Sigma; i = 1 I &Delta;h i &Delta;v i - h b - h gA V r ) - ( &Sigma; j = 1 J &Delta;h j &Delta;v j - h b - h gB V r ) (formula 2)
In formula, △ h ifor the vertical thickness of i-th grid that ray passes, h bfor uniform reference surface elevation, h gA, h gBfor the elevation on shot point, high speed top that geophone station is corresponding, V rfor replacement velocity, △ v ifor the speed of the grid of i-th after tomographic inversion, I is the grid number of shot point place from earth's surface vertically downward at a high speed top, and J is the grid number of geophone station place from earth's surface vertically downward to top at a high speed.
4) after three dimensional seismic data carries out tomographic statics, shot statics and receiver static correction amount is generally had, because the vibration source of ground micro-seismic is in underground, so the receiver static correction amount only choosing three dimensional seismic data calculates;
5) repeat step 2 to 4, until the static correction value of the three dimensional seismic data calculated tallies with the actual situation, and guarantee accuracy;
6) because three dimensional seismic data work area scope comprises ground micro-seismic work area, so the geophone station in 3-D seismics work area can be connected by terrestrial coordinate and ground micro-seismic geophone station coordinate, like this, for the geophone station of each ground micro-seismic, centered by this geophone station, set minimum, maximum radius R minand R maxif a certain 3-D seismics geophone station is at R minin scope, then think that the geophone station of this 3-D seismics overlaps with the geophone station of this ground micro-seismic, the static correction value of the two is then consistent; If the geophone station of the three dimensional seismic data do not overlapped with this ground micro-seismic geophone station, then at R minwith R maxthe method of inverse distance weighting is utilized to calculate the static correction value of this ground micro-seismic geophone station in scope.Wherein, R is greater than maxthe geophone station of 3-D seismics then not calculating in the scope considered;
The concrete methods of realizing of inverse distance weighting is: set interpolation ground micro-seismic geophone station as P (x p, y p, z p), P point minimum with maximum radius within the scope of have known 3-D seismics geophone station Q i(x i, y i, z i), i=1,2 ..., n, utilizes distance weighted inverse ratio method to the static correction value Z of P point pcarry out interpolation.The weighted mean of its Interpolation Principle to be the static correction value of interpolation point be 3-D seismics static correction value in the given radius of point to be inserted, in the size of power and point to be inserted and given radius, the distance dependent of known point, is the inverse of distance k power.That is:
Z p = &Sigma; i = 1 n Z i / d i 2 &Sigma; i = 1 n 1 / d i 2 (formula 3)
Wherein, d ifor the distance in interpolation point and given radius between i-th point.
7) step 6 is repeated, until all the geophone station of ground micro-seismic completes the calculating of static correction value;
8) static correction value of all ground micro-seismic geophone stations is applied to ground micro-seismic data.
Embodiment 2.As flow process as shown in Figure 1.Based on the ground micro-seismic static correcting method of 3D seismic data, should comprise the following steps based on the ground micro-seismic static correcting method of 3D seismic data.
1) choose the data in certain ground micro-seismic block of Shengli Oil Field construction, and choose the 3D seismic data comprising this work area, as shown in Figure 2.
2) utilize the 3D seismic data chosen to carry out Tomography Velocity inverting, the degree of depth of inverting is 700m, and velocity range is 1200 ~ 2500m/s, as shown in Figure 3.
3) tomographic statics that the rate pattern after Tomography Velocity inverting carries out based on first break time is utilized.
4) after three dimensional seismic data carries out tomographic statics, shot statics and receiver static correction amount is generally had, because the vibration source of ground micro-seismic is in underground, so the receiver static correction amount only choosing three dimensional seismic data calculates.
5) repeat step 2 to 4, until the static correction value of the three dimensional seismic data calculated tallies with the actual situation, and guarantee accuracy, as shown in Figure 4.
6) for the geophone station of each ground micro-seismic, centered by this geophone station, according to the number of the 3-D seismics geophone station that the respective track pitch in two work areas and computer capacity can comprise, given least radius R minscope is [5,10] and maximum radius R maxscope is [100,300], and unit is m, and this scope should according to the condition of construction in actual work area, specifically given.
7) if first 3-D seismics geophone station is at R minin scope, then think that the geophone station of this 3-D seismics overlaps with the geophone station of this ground micro-seismic, the static correction value of the two is identical, then stop search, and carries out the calculating of next ground micro-seismic geophone station; If search for the geophone station of whole 3-D seismics, and situation about overlapping with this ground micro-seismic geophone station, then at R minwith R maxthe method of inverse distance weighting is utilized to calculate the static correction value of this ground micro-seismic geophone station in scope, till the geophone station searching for whole 3-D seismics.Wherein, R is greater than maxthe geophone station of 3-D seismics then not calculating in the scope considered, and when utilizing inverse distance algorithm to calculate, weighting factor value is 2.
8) step 7 is repeated, until complete the static correction value calculating whole microearthquake geophone station, as shown in Figure 5, Fig. 5 is the elevation of each geophone station of ground micro-seismic and the microearthquake static correction value curve comparison figure of calculating, can find out, calculate static correction value and the elevation goodness of fit higher, only velocity variations district and high velocity slightly different.
9) Fig. 6 is the ground micro-seismic section before application static correction value, the static correction value of all ground micro-seismic geophone stations is applied to ground micro-seismic data, as shown in Figure 7.Comparison diagram 6 and Fig. 7, after application static correction, improve the twisted phenomena of ground micro-seismic section lineups, and do not occur the situation of localised jitter, for the process such as follow-up seismic source location are laid a good foundation.
The present invention, after acquisition seismic data, according to the scope in ground micro-seismic work area, chooses the three dimensional seismic data that can comprise this work area scope.Tomographic statics method is utilized to ask for the static correction value of three dimensional seismic data.Because the geophone station of three dimensional seismic data can be connected by the acceptance point of terrestrial coordinate and ground micro-seismic, like this for the geophone station of each ground micro-seismic, centered by this geophone station, given maximum a, least radius, the seismic detection point static correction value within the scope of least radius is as the static correction value of this microearthquake geophone station; Be greater than least radius, when being less than maximum radius, utilize the static correction value interpolation of three dimensional seismic data to go out the static correction value of microearthquake by the method for inverse distance interpolation.The present invention utilizes three dimensional seismic data to ask for the static correction value of ground micro-seismic, avoids a large amount of artificial first break pickings, and without any need for prior imformation, only need carry out the static correction of seismic data, is a kind of simple static corrections processing technology efficiently.

Claims (8)

1., based on the ground micro-seismic static correcting method of 3D seismic data, it is characterized in that: comprise the following steps:
Step 1: the three dimensional seismic data that can comprise this work area is chosen in the ground micro-seismic work area according to gathering;
Step 2: Tomography Velocity inverting is carried out to the three dimensional seismic data chosen;
Step 3: utilize the rate pattern of inverting to carry out tomographic statics;
Step 4: judge whether static correction value tallies with the actual situation, if qualified, performs step 5; Otherwise, repeat step 2-3, until find qualified static correction value;
Step 5: static correction value qualified in selecting step 4;
Step 6: the geophone station choosing a ground micro-seismic, centered by this geophone station, setting least radius R min, maximum radius R max;
If the distance d between a certain 3-D seismics geophone station and the geophone station chosen meets
d≤R min
Then think that the geophone station of this 3-D seismics overlaps with the geophone station of this ground micro-seismic, the static correction value of the two is consistent;
If the distance d between a certain 3-D seismics geophone station and the geophone station chosen meets
R min<d≤R max
Then at R minwith R maxutilize in scope, utilize the method for inverse distance weighting to calculate the static correction value of this ground micro-seismic geophone station;
If the distance d between a certain 3-D seismics geophone station and the geophone station chosen meets
d>R max
Then think this geophone station not in the scope calculating consideration;
According to mentioned above principle, all 3-D seismics geophone stations are chosen one by one and processed;
Step 7: the geophone station choosing other ground micro-seismic one by one, repeats step 6, until all the geophone station of ground micro-seismic completes the calculating of static correction value;
Step 8: the static correction value of all ground micro-seismic geophone stations is applied to ground micro-seismic data.
2., as claimed in claim 1 based on the ground micro-seismic static correcting method of 3D seismic data, it is characterized in that:
The concrete methods of realizing of the Tomography Velocity inverting of described step 2 is:
Near surface is separated into grid, and in each grid, speed is equal but unknown; Given initial velocity model carries out first arrival ray tracing, then the primary travel time and real first break time that calculate acquisition are compared, by separating the inversion equation group A △ S=△ b shown in formula 1, calculating the speed correction amount in each grid, and then initial model is modified;
a 11 a 12 . . . a 1 j . . . a 1 J a 21 a 22 . . . a 2 j . . . a 2 J . . . . . . . . . . . . . . . . . . a i 1 a i 2 . . . a ij . . . a iJ . . . . . . . . . . . . . . . . . . a I 1 a I 2 . . . a Ij . . . a IJ &Delta;s 1 &Delta;s 2 . . . &Delta;s j . . . &Delta;s J = &Delta;b 1 &Delta;b 2 . . . &Delta;b i . . . &Delta;b I (formula 1)
Wherein, in inversion equation group, A represents ray spread length matrix within a grid, and △ S is underground slowness correction, and △ b is poor between primary travel time.
3. the ground micro-seismic static correcting method based on 3D seismic data according to claim 1, is characterized in that:
Described step 2 needs the process of repetition FORWARD AND INVERSE PROBLEMS, approaches the minimum rate pattern of whilst on tour error gradually and stop iteration by successive ignition, and between primary travel time, difference △ b is the difference of true first break time and forward simulation primary travel time.
4., as claimed in claim 1 based on the ground micro-seismic static correcting method of 3D seismic data, it is characterized in that:
Utilize the rate pattern of inverting to carry out tomographic statics described in step 3 to refer to: after the Tomography Velocity inverting of three dimensional seismic data, carry out the tomographic statics process based on first break time;
Its concrete implementation method is: obtain the static correction value that near-surface velocity structure just can calculate shot point and geophone station position later easily; Under earth's surface-consistent prerequisite, think that the ray through earth's surface low velocity layer is perpendicular ray, the correcting value of each physical points is unique; Concrete formula is as follows:
t = t A + t B = - ( &Sigma; i = 1 I &Delta;h i &Delta;v i - h b - h gA V r ) - ( &Sigma; j = 1 J &Delta;h j &Delta;v j - h b - h gB V r ) (formula 2)
In formula, △ h ifor the vertical thickness of i-th grid that ray passes, h bfor uniform reference surface elevation, h gA, h gBfor the elevation on shot point, high speed top that geophone station is corresponding, V rfor replacement velocity, △ v ifor the speed of the grid of i-th after tomographic inversion, I is the grid number of shot point place from earth's surface vertically downward at a high speed top, and J is the grid number of geophone station place from earth's surface vertically downward to top at a high speed; t afor the static correction value of shot point, t bfor the static correction value of geophone station, t is that geophone station B is relative to this static correction value of one of shot point A.
5., as claimed in claim 1 based on the ground micro-seismic static correcting method of 3D seismic data, it is characterized in that:
After three dimensional seismic data carries out tomographic statics, generally there are shot statics and receiver static correction amount, when judging whether static correction value tallies with the actual situation in step 4, because the vibration source of ground micro-seismic is in underground, the receiver static correction amount only choosing three dimensional seismic data calculates, and judges whether it tallies with the actual situation.
6., as claimed in claim 1 based on the ground micro-seismic static correcting method of 3D seismic data, it is characterized in that:
The concrete methods of realizing of the weighting of inverse distance described in step 6 is:
If interpolation ground micro-seismic geophone station is P (x p, y p, z p), P point minimum with maximum radius within the scope of have known 3-D seismics geophone station Q i(x i, y i, z i), i=1,2 ..., n, utilizes distance weighted inverse ratio method to the static correction value Z of P point pcarry out interpolation; The weighted mean of its Interpolation Principle to be the static correction value of interpolation point be 3-D seismics static correction value in the given radius of point to be inserted, in the size of power and point to be inserted and given radius, the distance dependent of known point, is the inverse of distance k power; That is:
Z p = &Sigma; i = 1 n Z i / d i 2 &Sigma; i = 1 n 1 / d i 2 (formula 3)
Wherein, d ifor the distance in interpolation point and given radius between i-th point.
7., as claimed in claim 6 based on the ground micro-seismic static correcting method of 3D seismic data, it is characterized in that: R in described step 6 minand R maxthe track pitch that arranges according to ground micro-seismic of selection and 3-D seismics arrangement track pitch and determine.
8., as claimed in claim 6 based on the ground micro-seismic static correcting method of 3D seismic data, it is characterized in that: the traversal scope of parameter k is [0,2], parameter d itraversal scope be [R min, R max].
CN201410203127.3A 2014-05-14 2014-05-14 Ground micro-seismic static correcting method based on 3D seismic data Active CN105093319B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410203127.3A CN105093319B (en) 2014-05-14 2014-05-14 Ground micro-seismic static correcting method based on 3D seismic data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410203127.3A CN105093319B (en) 2014-05-14 2014-05-14 Ground micro-seismic static correcting method based on 3D seismic data

Publications (2)

Publication Number Publication Date
CN105093319A true CN105093319A (en) 2015-11-25
CN105093319B CN105093319B (en) 2017-11-28

Family

ID=54574173

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410203127.3A Active CN105093319B (en) 2014-05-14 2014-05-14 Ground micro-seismic static correcting method based on 3D seismic data

Country Status (1)

Country Link
CN (1) CN105093319B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107783185A (en) * 2017-09-14 2018-03-09 中国石油天然气股份有限公司 A kind of processing method and processing device of tomographic statics
CN107918154A (en) * 2016-10-10 2018-04-17 中国石油化工股份有限公司 A kind of combined method of the datum static correction amount of earth's surface-consistent
CN108828667A (en) * 2018-03-23 2018-11-16 中国矿业大学(北京) A kind of microseism complicated earth surface elevation correction method
CN109188527A (en) * 2018-10-31 2019-01-11 中国石油化工股份有限公司 The method that sea shoal quickly establishes three-dimensional near Sea Bottom rate pattern
CN109188526A (en) * 2018-11-02 2019-01-11 中国石油天然气股份有限公司大港油田分公司 A kind of modified earth's surface-consistent residual static corrections
CN109239781A (en) * 2018-10-26 2019-01-18 中国石油天然气股份有限公司 A kind of seismic data corrections method and device

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1773310A (en) * 2004-11-09 2006-05-17 中国石油天然气集团公司 Static correcting method during converted wave delaying
US20120163121A1 (en) * 2010-08-27 2012-06-28 Board Of Regents Of The University Of Texas System Extracting sv shear data from p-wave marine data
CN102819042A (en) * 2011-06-10 2012-12-12 中国石油天然气集团公司 Near-surface chromatography inversion static rectification method of rarefied seismic data

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1773310A (en) * 2004-11-09 2006-05-17 中国石油天然气集团公司 Static correcting method during converted wave delaying
US20120163121A1 (en) * 2010-08-27 2012-06-28 Board Of Regents Of The University Of Texas System Extracting sv shear data from p-wave marine data
CN102819042A (en) * 2011-06-10 2012-12-12 中国石油天然气集团公司 Near-surface chromatography inversion static rectification method of rarefied seismic data

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
LIU LIPING ET AL.: "Near Surface Characteristics and Static Correction in Desert Area of Central Junggar Basin", 《NSGAPC BEIJING 2013: FIRST NEAR SURFACE GEOPHYSICS ASIA PACIFIC CONFERENCE》 *
姜宇东等: "地面微地震监测资料静校正方法研究", 《石油物探》 *
王孝等: "先验信息约束层析成像静校正技术在柴西南地区的应用", 《石油物探》 *
王振华等: "复杂地表条件下的静校正方法", 《石油地球物理勘探》 *
魏亦文等: "混合插值法重构近地表模型", 《计算机辅助设计与图形学学报》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107918154A (en) * 2016-10-10 2018-04-17 中国石油化工股份有限公司 A kind of combined method of the datum static correction amount of earth's surface-consistent
CN107918154B (en) * 2016-10-10 2019-06-04 中国石油化工股份有限公司 A kind of combined method of the datum static correction amount of earth's surface-consistent
CN107783185A (en) * 2017-09-14 2018-03-09 中国石油天然气股份有限公司 A kind of processing method and processing device of tomographic statics
CN107783185B (en) * 2017-09-14 2019-05-07 中国石油天然气股份有限公司 A kind of processing method and processing device of tomographic statics
CN108828667A (en) * 2018-03-23 2018-11-16 中国矿业大学(北京) A kind of microseism complicated earth surface elevation correction method
CN108828667B (en) * 2018-03-23 2019-09-03 中国矿业大学(北京) A kind of microseism complicated earth surface elevation correction method
CN109239781A (en) * 2018-10-26 2019-01-18 中国石油天然气股份有限公司 A kind of seismic data corrections method and device
CN109188527A (en) * 2018-10-31 2019-01-11 中国石油化工股份有限公司 The method that sea shoal quickly establishes three-dimensional near Sea Bottom rate pattern
CN109188527B (en) * 2018-10-31 2021-02-26 中国石油化工股份有限公司 Method for rapidly establishing three-dimensional offshore bottom speed model in beach and shallow sea area
CN109188526A (en) * 2018-11-02 2019-01-11 中国石油天然气股份有限公司大港油田分公司 A kind of modified earth's surface-consistent residual static corrections

Also Published As

Publication number Publication date
CN105093319B (en) 2017-11-28

Similar Documents

Publication Publication Date Title
CN105093319A (en) Ground micro-seismic static correction method based on three-dimensional seismic data
CN102841379B (en) Method for analyzing pre-stack time migration and speed based on common scatter point channel set
CN104133245B (en) The static correcting method and system of a kind of seismic data
CN103399300B (en) Wave packet superposition microseism ground location method
CN105589100A (en) Micro-seismic source location and velocity model simultaneous inversion method
CN100549731C (en) Demodulator probe secondary localization method
CN102879820B (en) Based on the three-dimensional table layer model construction method of triangle gridding
CN104730579B (en) A kind of joint static correcting method of ripple in length and breadth based on calculation of near surface shear velocity inverting
CN104181593B (en) Three-dimensional ray-free tracing contraflexure wave tomographic imaging method and device
CN103995288A (en) Gauss beam prestack depth migration method and device
CN104360396B (en) A kind of three kinds of preliminary wave Zoumaling tunnel methods of TTI medium between offshore well
CN106291725A (en) A kind of method of fast inversion underground geologic bodies locus
CN105277978A (en) Method and device for determining near-ground-surface speed model
CN105807316A (en) Surface observation microseism speed model correcting method based on amplitude stack
CN102901985A (en) Depth domain layer speed correcting method suitable for undulating surface
CN107817516A (en) Near surface modeling method and system based on preliminary wave information
CN105549077A (en) Micro-earthquake epicenter positioning method calculated based on multilevel multi-scale grid similarity coefficient
CN104090301A (en) Three-dimensional high-frequency static correction value obtaining method
CN104199103A (en) Static correction method and static correction device for converted wave
CN106199704B (en) A kind of Three-dimendimal fusion submarine cable seismic data velocity modeling method
CN103513279B (en) A kind of illumination analysis computing method based on seismic wave equation and calculation element
CN105607119B (en) Near-surface model construction method and static correction value acquiring method
CN103576197A (en) Method for extracting converted wave angle channel set
CN105093298B (en) A kind of quick calculation method of the four-dimensional search inverse time superposition of microseism data
CN102819042A (en) Near-surface chromatography inversion static rectification method of rarefied seismic data

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
GR01 Patent grant