CN105549075B - Ask for the method and device of shallow-layer velocity field - Google Patents

Ask for the method and device of shallow-layer velocity field Download PDF

Info

Publication number
CN105549075B
CN105549075B CN201410589005.2A CN201410589005A CN105549075B CN 105549075 B CN105549075 B CN 105549075B CN 201410589005 A CN201410589005 A CN 201410589005A CN 105549075 B CN105549075 B CN 105549075B
Authority
CN
China
Prior art keywords
value
floating surface
static correction
geophone station
correction value
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.)
Active
Application number
CN201410589005.2A
Other languages
Chinese (zh)
Other versions
CN105549075A (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 Natural Gas Co Ltd
Original Assignee
China Petroleum and Natural Gas Co Ltd
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 Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN201410589005.2A priority Critical patent/CN105549075B/en
Publication of CN105549075A publication Critical patent/CN105549075A/en
Application granted granted Critical
Publication of CN105549075B publication Critical patent/CN105549075B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The embodiment of the present invention provides a kind of method and device for asking for shallow-layer velocity field.This method includes:The datum static correction amount of geological data is subjected to low-and high-frequency separation;Obtain the value in shot point Relative Floating face and the value in geophone station Relative Floating face;Correcting value of the shot point from earth's surface to floating surface is obtained according to floating surface static correction value, shot point high frequency static correction value and the value in shot point Relative Floating face;Correcting value of the geophone station from earth's surface to floating surface is obtained according to floating surface static correction value, geophone station high frequency static correction value and the value in geophone station Relative Floating face;High-precision shallow-layer velocity field is obtained according to correcting value correction geological data.Each shot point, each geophone station are corrected on corresponding floating surface by the embodiment of the present invention by the correcting value of correcting value and geophone station from earth's surface to floating surface of the shot point from earth's surface to floating surface, recycle first break information to realize shallow-layer velocity modeling by floating surface tomographic inversion technology, improve the precision of shallow-layer velocity field.

Description

Ask for the method and device of shallow-layer velocity field
Technical field
The present embodiments relate to petroleum exploration seismic data processing technology field, more particularly to one kind to ask for shallow-layer speed The method and device of field.
Background technology
It is a committed step in process of seismic data processing to ask for seismic velocity, multiple medium points in normal-moveout spectrum be present, Each medium point is corresponding with seismic velocity, and by All Media point in normal-moveout spectrum, corresponding seismic velocity forms seismic velocity respectively , the precision of seismic velocity field is the key factor of seismic imaging.
There is clear area in usual shallow-layer normal-moveout spectrum, for the phenomenon, shallow-layer speed is upwardly extended using according to speed trend Mode make up clear area corresponding to shallow-layer speed.
But larger error be present by the shallow-layer speed that prior art obtains and real shallow-layer speed, cause shallow-layer Velocity field precision is relatively low.
The content of the invention
The embodiment of the present invention provides a kind of method and device for asking for shallow-layer velocity field, to obtain high-precision shallow-layer speed .
The one side of the embodiment of the present invention is to provide a kind of method for asking for shallow-layer velocity field, including:
The datum static correction amount of geological data is based on floating surface and carries out low-and high-frequency separation, obtains floating surface static correction Amount, shot point high frequency static correction value and geophone station high frequency static correction value;
Obtain the value of the relatively described floating surface of shot point and the value of the relatively described floating surface of geophone station;
According to the floating surface static correction value, the shot point high frequency static correction value and the relatively described floating surface of the shot point Value obtains first correcting value of the shot point from earth's surface to the floating surface;
According to the floating surface static correction value, the geophone station high frequency static correction value and the relatively described floating of the geophone station The value in face obtains second correcting value of the geophone station from earth's surface to the floating surface;
The earthquake number after the geological data is corrected is corrected according to first correcting value and second correcting value According to the geological data first break picking after the correction, according to the first arrival progress tomographic inversion acquisition shallow-layer velocity field.
The other side of the embodiment of the present invention is to provide a kind of device for asking for shallow-layer velocity field, including:
Separation module, low-and high-frequency separation is carried out for the datum static correction amount of geological data to be based on into floating surface, is obtained Floating surface static correction value, shot point high frequency static correction value and geophone station high frequency static correction value;
Acquisition module, for obtaining the value of the relatively described floating surface of shot point and the value of the relatively described floating surface of geophone station;
Processing module, for according to the floating surface static correction value, the shot point high frequency static correction value and the shot point phase First correcting value of the shot point from earth's surface to the floating surface is obtained to the value of the floating surface;According to the quiet school of the floating surface The value of positive quantity, the geophone station high frequency static correction value and the relatively described floating surface of the geophone station obtains the geophone station from earth's surface To the second correcting value of the floating surface;The geological data is corrected according to first correcting value and second correcting value to obtain Geological data after must correcting, to the geological data first break picking after the correction, carry out tomographic inversion according to the first arrival and obtain Obtain the shallow-layer velocity field.
The method and device provided in an embodiment of the present invention for asking for shallow-layer velocity field, by shot point from earth's surface to floating surface Each shot point, each geophone station are corrected on corresponding floating surface by the correcting value of correcting value and geophone station from earth's surface to floating surface, then Shallow-layer velocity modeling is realized by floating surface tomographic inversion technology using first break information, the first arrival tomographic inversion technology of routine is created New opplication improves the precision of shallow-layer velocity field in the first arrival tomographic inversion of floating surface.
Brief description of the drawings
Fig. 1 is the method flow diagram provided in an embodiment of the present invention for asking for shallow-layer velocity field;
Fig. 2 is floating surface provided in an embodiment of the present invention and various correcting value schematic diagrames;
Fig. 3 A are the floating surface static correction value schematic diagram for the geological data that another embodiment of the present invention provides;
Fig. 3 B are the shot point high frequency static correction value schematic diagram that another embodiment of the present invention provides;
Fig. 3 C are the geophone station high frequency static correction value schematic diagram that another embodiment of the present invention provides;
Fig. 4 A are the shot position floating surface plan that another embodiment of the present invention provides;
Fig. 4 B are the geophone station position floating surface plan that another embodiment of the present invention provides;
Fig. 5 is that the seismic data big gun that another embodiment of the present invention provides is cautious from earth's surface to floating surface correcting value plan;
Fig. 6 be another embodiment of the present invention provide using big gun, it is cautious from earth's surface be corrected to floating surface after single shot record;
Fig. 7 A are the radiographic density graph model depth inspection and the reliability of the adjustment model analysis chart that another embodiment of the present invention provides;
Fig. 7 B are the shallow-layer velocity field schematic diagram that another embodiment of the present invention provides;
Fig. 8 A are the high-precision shallow-layer velocity field schematic diagram that another embodiment of the present invention provides;
Fig. 8 B are the shallow-layer velocity field schematic diagram that prior art is asked for;
Fig. 9 is the shallow-layer speed that the high-precision shallow-layer velocity field that another embodiment of the present invention provides is asked for prior art The comparison diagram of field;
Figure 10 A are the CRP gather that the shallow-layer velocity field asked for using prior art carries out migration before stack;
Figure 10 B are that the high-precision shallow-layer velocity field provided using another embodiment of the present invention carries out the anti-altogether of migration before stack Exit point trace gather.
Figure 11 is the structure drawing of device provided in an embodiment of the present invention for asking for shallow-layer velocity field.
Embodiment
Fig. 1 is the method flow diagram provided in an embodiment of the present invention for asking for shallow-layer velocity field.Fig. 2 carries for the embodiment of the present invention The floating surface of confession and various correcting value schematic diagrames.It is provided in an embodiment of the present invention ask for shallow-layer velocity field method and step it is as follows:
Step S101, the datum static correction amount of geological data is based on floating surface and carries out low-and high-frequency separation, floated Face static correction value, shot point high frequency static correction value and geophone station high frequency static correction value;
Determine fixed datum, and using shot point and geophone station elevation information generation floating surface, specially described shot point and The midpoint of the geophone station line is floating point, and multiple points that float are formed described floating corresponding to multiple shot points and multiple geophone stations Dynamic face.As shown in Fig. 2 the midpoint C of shot point S and geophone station D lines is floating point, when shot point is multiple, geophone station is also to be multiple When, the midpoint of any shot point and any geophone station line will form a floating point, and multiple points that float form as shown in Figure 2 float Dynamic face 30.
The floating surface static correction value is that point is floated corresponding to the shot point and the geophone station relative to fixed datum Value.As shown in Fig. 2 floating surface static correction value TcmpIt is that point C is floated corresponding to shot point S and geophone station D relative to fixed datum 10 value;The horizontal line for crossing C points is 40, and it is S ' to cross shot point S vertical line with the intersection point of horizontal line 40, and S to S ' is shot point high frequency Static correction value Tstat-source;The intersection point of the vertical line and horizontal line 40 of crossing geophone station D is D ', and D to D ' is the quiet school of geophone station high frequency Positive quantity Tstat-detect
Step S102, the value of the relatively described floating surface of shot point and the value of the relatively described floating surface of geophone station are obtained;
The value of the relatively described floating surface of the shot point was that the vertical line of the shot point is relative with the intersection point of the floating surface In the half of the value of the fixed datum;The value of the relatively described floating surface of the geophone station was the vertical line of the geophone station Half with the intersection point of the floating surface relative to the value of the fixed datum.
As shown in Fig. 2 shot point S vertical line and the intersection point of floating surface 30 are S 〞, S 〞 relative to the fixed datum excessively The half of 10 value is the value T of the relatively described floating surface of shot pointcmp-source.Cross geophone station D vertical line and the friendship of floating surface 30 Point is D 〞, and D 〞 are relative to the value T that the half of the value of the fixed datum 10 is the relatively described floating surface of geophone stationcmp-detect
Step S103, according to the floating surface static correction value, the shot point high frequency static correction value and the shot point with respect to institute The value for stating floating surface obtains first correcting value of the shot point from earth's surface to the floating surface;
S to S 〞 are first correcting value T of the shot point from earth's surface to the floating surfaces, the first correcting value Ts=Tcmp/2+ Tstat-source-Tcmp-source, wherein, TcmpRepresent the floating surface static correction value, Tstat-sourceRepresent the quiet school of shot point high frequency Positive quantity, Tcmp-sourceRepresent the value of the relatively described floating surface of the shot point.
Step S104, according to the floating surface static correction value, the geophone station high frequency static correction value and the geophone station phase Second correcting value of the geophone station from earth's surface to the floating surface is obtained to the value of the floating surface;
D to D 〞 are second correcting value T of the geophone station from earth's surface to the floating surfaced, the second correcting value Td=Tcmp/2+ Tstat-detect-Tcmp-detect, wherein, TcmpRepresent the floating surface static correction value, Tstat-detectRepresent that the geophone station high frequency is quiet Correcting value, Tcmp-detectRepresent the value of the relatively described floating surface of the geophone station.
Step S105, after being corrected according to first correcting value and second correcting value correction geological data Geological data, to the geological data first break picking after the correction, according to the first arrival carry out tomographic inversion obtain it is described shallow Interval velocity field.
The the first correcting value T obtained according to step S103sThe the second correcting value T obtained with step S104dCorrect the earthquake Data corrected after geological data, to the geological data first break picking after the correction, chromatographed according to the first arrival Inverting obtains the high-precision shallow-layer velocity field.
On the basis of the embodiment of the present invention, the shot point high frequency static correction value is the shot point relative to the fixed base The value in quasi- face subtracts the floating surface static correction value;The geophone station high frequency static correction value is the geophone station relative to described solid The value for determining reference plane subtracts the floating surface static correction value.
As shown in Fig. 2 shot point high frequency static correction value Tstat-sourceSubtracted for shot point S-phase for the value of fixed datum 10 floating Dynamic face static correction value Tcmp;Geophone station high frequency static correction value Tstat-detectValue for geophone station D relative to the fixed datum 10 Subtract the floating surface static correction value Tcmp
Correcting value and geophone station school from earth's surface to floating surface of the embodiment of the present invention by shot point from earth's surface to floating surface Each shot point, each geophone station are corrected on corresponding floating surface by positive quantity, recycle first break information to pass through floating surface tomographic inversion skill Art realizes shallow-layer velocity modeling, and the first arrival tomographic inversion technological innovation of routine is applied into the first arrival tomographic inversion in floating surface, carried The high precision of shallow-layer velocity field.
The present embodiment is carried out using certain actual three dimensional seismic data, and it shares 349 big guns, 2592 geophone stations, obtained The floating surface static correction value T of the real data as shown in Figure 3Acmp, shot point high frequency static correction value as shown in Figure 3 B Tstat-source, geophone station high frequency static correction value T as shown in Figure 3 Cstat-detect
Ask for the value T of the relatively described floating surface of shot pointcmp-sourceWith the value T of the relatively described floating surface of geophone stationcmp-detect, Obtain the value T of the relatively described floating surface of shot point as shown in Figure 4 Acmp-sourceGeophone station as shown in Figure 4 B is relatively described floating The value T in dynamic facecmp-detect, and Tcmp-sourceIt was the intersection point of vertical line and the floating surface of shot point relative to the fixed base The half of the value in quasi- face;Tcmp-detectIt was the intersection point of vertical line and the floating surface of geophone station relative to the fixed reference The half of the value in face.
Ask for first correcting value T of the shot point from earth's surface to floating surfaces, Ts=Tcmp/2+Tstat-source-Tcmp-source, ask for Second correcting value T of the geophone station from earth's surface to the floating surfaced, Td=Tcmp/2+Tstat-detect-Tcmp-detect, Fig. 5, which is shown, to be asked It is the big gun that takes, cautious from earth's surface to floating surface correcting value TsAnd Td
As shown in fig. 6, apply shot correction amount Ts and detection point calibration based on floating surface in original single-shot geological data Measure Td, geological data is corrected on floating surface from collection earth's surface.
To the geological data first break picking after correction, the tomographic inversion based on first arrival is carried out, obtains the interval velocity mould of shallow-layer Type, model depth inspection and the reliability of the adjustment model analysis are carried out according to inverting radiographic density Fig. 7 A to determine the shallow-layer layer of final inverting Rate pattern such as Fig. 7 B.
As shown in Figure 8 A and 8 B, high-precision shallow-layer velocity field provided in an embodiment of the present invention is asked for than prior art Shallow-layer velocity field is more rationally accurate.
Fig. 9 is that the comparison diagram that the conventional speeds field determined asks for velocity field with this method is composed based on velocity analysis, deep in figure Color curve is the shallow-layer speed that this method obtains, and White curves are the speed routinely asked for, it can be seen that the speed that this method is asked for The blank of earthquake shallow-layer normal-moveout spectrum has been filled up in degree field, and more true and reliable.
As shown in figs. 10 a and 10b, the shallow-layer velocity field asked for compared to prior art carries out the common reflection of migration before stack Point trace gather, the CRP gather middle-shallow layer that high-precision shallow-layer velocity field provided in an embodiment of the present invention carries out migration before stack are same Phase axle more focuses on, continuity is more preferable.
Figure 11 is the structure drawing of device provided in an embodiment of the present invention for asking for shallow-layer velocity field.It is provided in an embodiment of the present invention The handling process for the embodiment of the method offer for asking for shallow-layer velocity field can be performed by asking for the device of shallow-layer velocity field.Such as Figure 11 institutes Show, asking for the device 110 of shallow-layer velocity field includes separation module 111, acquisition module 112 and processing module 113, wherein, separation Module 111 is used to the datum static correction amount of geological data being based on floating surface progress low-and high-frequency separation, obtains the quiet school of floating surface Positive quantity, shot point high frequency static correction value and geophone station high frequency static correction value;Acquisition module 112 is used to obtain the relatively described floating of shot point The value of the relatively described floating surface of value and geophone station in face;Processing module 113 is used for according to the floating surface static correction value, the big gun The value of point high frequency static correction value and the relatively described floating surface of the shot point obtains the of the shot point from earth's surface to the floating surface One correcting value;It is relatively described floating according to the floating surface static correction value, the geophone station high frequency static correction value and the geophone station The value in dynamic face obtains second correcting value of the geophone station from earth's surface to the floating surface;According to first correcting value and described Second correcting value corrects the geological data after the geological data is corrected, and the geological data after the correction is picked up just Extremely, carry out tomographic inversion according to the first arrival and obtain the shallow-layer velocity field.
The midpoint of the shot point and the geophone station line is floating point, more as corresponding to multiple shot points and multiple geophone stations The individual point that floats forms the floating surface.
The processing module 113 is specifically used for calculating first correcting value and second correcting value, first school Positive quantity Ts=Tcmp/2+Tstat-source-Tcmp-source, wherein, TcmpRepresent the floating surface static correction value, Tstat-sourceRepresent institute State shot point high frequency static correction value, Tcmp-sourceRepresent the value of the relatively described floating surface of the shot point;The second correcting value Td= Tcmp/2+Tstat-detect-Tcmp-detect, wherein, TcmpRepresent the floating surface static correction value, Tstat-detectRepresent the geophone station High frequency static correction value, Tcmp-detectRepresent the value of the relatively described floating surface of the geophone station.
The floating surface static correction value is that point is floated corresponding to the shot point and the geophone station relative to fixed datum Value;
The value of the relatively described floating surface of the shot point was that the vertical line of the shot point is relative with the intersection point of the floating surface In the half of the value of the fixed datum;
The value of the relatively described floating surface of the geophone station was intersection point of the vertical line with the floating surface of the geophone station Relative to the half of the value of the fixed datum.
The shot point high frequency static correction value is that the shot point subtracts the floating surface relative to the value of the fixed datum Static correction value;
The geophone station high frequency static correction value is that the geophone station subtracts described float relative to the value of the fixed datum Dynamic face static correction value.
Correcting value and geophone station school from earth's surface to floating surface of the embodiment of the present invention by shot point from earth's surface to floating surface Each shot point, each geophone station are corrected on corresponding floating surface by positive quantity, recycle first break information to pass through floating surface tomographic inversion skill Art realizes shallow-layer velocity modeling, and the first arrival tomographic inversion technological innovation of routine is applied into the first arrival tomographic inversion in floating surface, carried The high precision of shallow-layer velocity field.
In summary, the embodiment of the present invention by correcting value of the shot point from earth's surface to floating surface and geophone station from earth's surface to floating Each shot point, each geophone station are corrected on corresponding floating surface by the correcting value in dynamic face, recycle first break information to pass through floating surface layer Analysis inversion technique realizes shallow-layer velocity modeling, and the first arrival tomographic inversion technological innovation of routine is applied into the first arrival chromatography in floating surface Inverting, improve the precision of shallow-layer velocity field;In addition, the embodiment of the present invention combines shallow-layer tomographic inversion speed and deep reflex Both wave velocities advantage, realize and field built based on floating face velocity, filled up SEISMIC VELOCTTY ANALYSIS spectrum blank, obtain finally from Shallow to deep high-precision velocity field, it is very beneficial for migration before stack processing.
In several embodiments provided by the present invention, it should be understood that disclosed apparatus and method, it can be passed through Its mode is realized.For example, device embodiment described above is only schematical, for example, the division of the unit, only Only a kind of division of logic function, there can be other dividing mode when actually realizing, such as multiple units or component can be tied Another system is closed or is desirably integrated into, or some features can be ignored, or do not perform.It is another, it is shown or discussed Mutual coupling or direct-coupling or communication connection can be the INDIRECT COUPLINGs or logical by some interfaces, device or unit Letter connection, can be electrical, mechanical or other forms.
The unit illustrated as separating component can be or may not be physically separate, show as unit The part shown can be or may not be physical location, you can with positioned at a place, or can also be distributed to multiple On NE.Some or all of unit therein can be selected to realize the mesh of this embodiment scheme according to the actual needs 's.
In addition, each functional unit in each embodiment of the present invention can be integrated in a processing unit, can also That unit is individually physically present, can also two or more units it is integrated in a unit.Above-mentioned integrated list Member can both be realized in the form of hardware, can also be realized in the form of hardware adds SFU software functional unit.
The above-mentioned integrated unit realized in the form of SFU software functional unit, can be stored in one and computer-readable deposit In storage media.Above-mentioned SFU software functional unit is stored in a storage medium, including some instructions are causing a computer It is each that equipment (can be personal computer, server, or network equipment etc.) or processor (processor) perform the present invention The part steps of embodiment methods described.And foregoing storage medium includes:USB flash disk, mobile hard disk, read-only storage (Read- Only Memory, ROM), random access memory (Random Access Memory, RAM), magnetic disc or CD etc. it is various Can be with the medium of store program codes.
Those skilled in the art can be understood that, for convenience and simplicity of description, only with above-mentioned each functional module Division progress for example, in practical application, can be complete by different functional modules by above-mentioned function distribution as needed Into the internal structure of device being divided into different functional modules, to complete all or part of function described above.On The specific work process of the device of description is stated, the corresponding process in preceding method embodiment is may be referred to, will not be repeated here.
Finally it should be noted that:Various embodiments above is merely illustrative of the technical solution of the present invention, rather than its limitations;To the greatest extent The present invention is described in detail with reference to foregoing embodiments for pipe, it will be understood by those within the art that:Its according to The technical scheme described in foregoing embodiments can so be modified, either which part or all technical characteristic are entered Row equivalent substitution;And these modifications or replacement, the essence of appropriate technical solution is departed from various embodiments of the present invention technology The scope of scheme.

Claims (8)

  1. A kind of 1. method for asking for shallow-layer velocity field, it is characterised in that including:
    The datum static correction amount of geological data is based on floating surface and carries out low-and high-frequency separation, obtains floating surface static correction value, big gun Point high frequency static correction value and geophone station high frequency static correction value;
    Obtain the value of the relatively described floating surface of shot point and the value of the relatively described floating surface of geophone station;
    Value according to the floating surface static correction value, the shot point high frequency static correction value and the relatively described floating surface of the shot point obtains Obtain first correcting value of the shot point from earth's surface to the floating surface;
    According to the floating surface static correction value, the geophone station high frequency static correction value and the relatively described floating surface of the geophone station Value obtains second correcting value of the geophone station from earth's surface to the floating surface;
    The geological data after the geological data is corrected is corrected according to first correcting value and second correcting value, it is right Geological data first break picking after the correction, carry out tomographic inversion according to the first arrival and obtain the shallow-layer velocity field;
    The floating surface static correction value is value of the point relative to fixed datum of being floated corresponding to the shot point and the geophone station;
    The value of the relatively described floating surface of the shot point was the intersection point of vertical line and the floating surface of the shot point relative to institute State the half of the value of fixed datum;
    The value of the relatively described floating surface of the geophone station was that the vertical line of the geophone station is relative with the intersection point of the floating surface In the half of the value of the fixed datum.
  2. 2. according to the method for claim 1, it is characterised in that the midpoint of the shot point and the geophone station line is floating Point, multiple points that float form the floating surface corresponding to multiple shot points and multiple geophone stations.
  3. 3. according to the method for claim 2, it is characterised in that described according to the floating surface static correction value, the shot point The value of high frequency static correction value and the relatively described floating surface of the shot point obtains first of the shot point from earth's surface to the floating surface Correcting value includes:
    The first correcting value Ts=Tcmp/2+Tstat-source-Tcmp-source, wherein, TcmpThe floating surface static correction value is represented, Tstat-sourceRepresent the shot point high frequency static correction value, Tcmp-sourceRepresent the value of the relatively described floating surface of the shot point;
    It is described according to the floating surface static correction value, the geophone station high frequency static correction value and the relatively described floating of the geophone station The value in face, which obtains second correcting value of the geophone station from earth's surface to the floating surface, to be included:
    The second correcting value Td=Tcmp/2+Tstat-detect-Tcmp-detect, wherein, TcmpThe floating surface static correction value is represented, Tstat-detectRepresent the geophone station high frequency static correction value, Tcmp-detectRepresent the value of the relatively described floating surface of the geophone station.
  4. 4. according to the method for claim 3, it is characterised in that the shot point high frequency static correction value be the shot point relative to The value of the fixed datum subtracts the floating surface static correction value;
    The geophone station high frequency static correction value is that the geophone station subtracts the floating surface relative to the value of the fixed datum Static correction value.
  5. A kind of 5. device for asking for shallow-layer velocity field, it is characterised in that including:
    Separation module, low-and high-frequency separation is carried out for the datum static correction amount of geological data to be based on into floating surface, is floated Face static correction value, shot point high frequency static correction value and geophone station high frequency static correction value;
    Acquisition module, for obtaining the value of the relatively described floating surface of shot point and the value of the relatively described floating surface of geophone station;
    Processing module, for according to the relative institute of the floating surface static correction value, the shot point high frequency static correction value and the shot point The value for stating floating surface obtains first correcting value of the shot point from earth's surface to the floating surface;According to the floating surface static correction Amount, the value of the geophone station high frequency static correction value and the relatively described floating surface of the geophone station obtain the geophone station from earth's surface to Second correcting value of the floating surface;The geological data is corrected according to first correcting value and second correcting value to obtain Geological data after correction, to the geological data first break picking after the correction, tomographic inversion acquisition is carried out according to the first arrival The shallow-layer velocity field;
    The floating surface static correction value is value of the point relative to fixed datum of being floated corresponding to the shot point and the geophone station;
    The value of the relatively described floating surface of the shot point was the intersection point of vertical line and the floating surface of the shot point relative to institute State the half of the value of fixed datum;
    The value of the relatively described floating surface of the geophone station was that the vertical line of the geophone station is relative with the intersection point of the floating surface In the half of the value of the fixed datum.
  6. 6. the device according to claim 5 for asking for shallow-layer velocity field, it is characterised in that the shot point and the geophone station The midpoint of line is floating point, and multiple points that float form the floating surface corresponding to multiple shot points and multiple geophone stations.
  7. 7. the device according to claim 6 for asking for shallow-layer velocity field, it is characterised in that the processing module is specifically used for Calculate first correcting value and second correcting value, the first correcting value Ts=Tcmp/2+Tstat-source-Tcmp-source, Wherein, TcmpRepresent the floating surface static correction value, Tstat-sourceRepresent the shot point high frequency static correction value, Tcmp-sourceRepresent The value of the relatively described floating surface of the shot point;The second correcting value Td=Tcmp/2+Tstat-detect-Tcmp-detect, wherein, Tcmp Represent the floating surface static correction value, Tstat-detectRepresent the geophone station high frequency static correction value, Tcmp-detectRepresent the inspection The value of the relatively described floating surface of wave point.
  8. 8. the device according to claim 7 for asking for shallow-layer velocity field, it is characterised in that the shot point high frequency static correction value It is that the shot point subtracts the floating surface static correction value relative to the value of the fixed datum;
    The geophone station high frequency static correction value is that the geophone station subtracts the floating surface relative to the value of the fixed datum Static correction value.
CN201410589005.2A 2014-10-28 2014-10-28 Ask for the method and device of shallow-layer velocity field Active CN105549075B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410589005.2A CN105549075B (en) 2014-10-28 2014-10-28 Ask for the method and device of shallow-layer velocity field

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410589005.2A CN105549075B (en) 2014-10-28 2014-10-28 Ask for the method and device of shallow-layer velocity field

Publications (2)

Publication Number Publication Date
CN105549075A CN105549075A (en) 2016-05-04
CN105549075B true CN105549075B (en) 2017-11-10

Family

ID=55828375

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410589005.2A Active CN105549075B (en) 2014-10-28 2014-10-28 Ask for the method and device of shallow-layer velocity field

Country Status (1)

Country Link
CN (1) CN105549075B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112394406A (en) * 2019-08-13 2021-02-23 中国石油化工股份有限公司 Method for establishing simulated earth surface depth domain velocity model

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2691811A1 (en) * 1992-05-29 1993-12-03 Total Sa Seismic prospecting model determn. - gives optimum underground seismic wave propagation velocity distribution model
US6501703B2 (en) * 2001-01-21 2002-12-31 Xixiang Zhou Static correction method for exploration seismic data using first arrivals of seismic waves
GB2378507B (en) * 2001-06-08 2005-08-17 Pgs Americas Inc Extraction of p-wave and s-wave velocities from multi-component seismic data by joint velocity inversion
CN1948999A (en) * 2005-10-10 2007-04-18 中国石油天然气股份有限公司 Method of approximating layer displacing static correct
CN103472483A (en) * 2013-09-27 2013-12-25 中国石油化工股份有限公司 Speed modeling method based on true earth surface and floating datum plane
CN103576200A (en) * 2012-07-19 2014-02-12 中国石油天然气集团公司 Low signal-to-noise ratio zone shallow wave impedance interface static correction method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2691811A1 (en) * 1992-05-29 1993-12-03 Total Sa Seismic prospecting model determn. - gives optimum underground seismic wave propagation velocity distribution model
US6501703B2 (en) * 2001-01-21 2002-12-31 Xixiang Zhou Static correction method for exploration seismic data using first arrivals of seismic waves
GB2378507B (en) * 2001-06-08 2005-08-17 Pgs Americas Inc Extraction of p-wave and s-wave velocities from multi-component seismic data by joint velocity inversion
CN1948999A (en) * 2005-10-10 2007-04-18 中国石油天然气股份有限公司 Method of approximating layer displacing static correct
CN103576200A (en) * 2012-07-19 2014-02-12 中国石油天然气集团公司 Low signal-to-noise ratio zone shallow wave impedance interface static correction method
CN103472483A (en) * 2013-09-27 2013-12-25 中国石油化工股份有限公司 Speed modeling method based on true earth surface and floating datum plane

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
复杂地表静校正方法研究及应用;熊定钰;《中国博士学位论文全文数据库.基础科学辑》;20130115(第1期);第A011-25页 *

Also Published As

Publication number Publication date
CN105549075A (en) 2016-05-04

Similar Documents

Publication Publication Date Title
CN105974479B (en) The chromatography 2D/3D anisotropy Depth Domain velocity modeling method of GPU space lattice body
CN105301636B (en) The method for building up and device of rate pattern
CN109917454B (en) True earth surface prestack depth migration imaging method and device based on double reference surfaces
WO2006054181A1 (en) Method for processing at least two sets of seismic data
CN104483705A (en) Three-dimensional residual static correction method
CN113552625B (en) Multi-scale full waveform inversion method for conventional land-domain seismic data
CN105093319B (en) Ground micro-seismic static correcting method based on 3D seismic data
CN101984366A (en) Fluctuating surface pre-stack time migration method and device based on common aperture surface
CN108828668B (en) Pre-stack time migration data processing method and device
CN107229071B (en) A kind of subsurface structure inversion imaging method
CN109239781A (en) A kind of seismic data corrections method and device
CN103217710A (en) Method and system for acquiring shared offset distance and shared azimuth angle domain track sets
CN107817525A (en) A kind of near surface constant speed fill method and device
CN111638551A (en) Seismic first-motion wave travel time chromatography method and device
CN112099082B (en) Seismic folding wave travel time inversion method for coplanar element common azimuth gather
CN105717538B (en) Relief surface earthquake data offset reference plane conversion method and device
CN105549075B (en) Ask for the method and device of shallow-layer velocity field
CN113960532A (en) Microseism positioning method based on secondary positioning calculation of imaginary source
CN105572741B (en) A method of calculating 3D high frequency static correction value
CN102798893B (en) Three-dimensional superimposed tomography inversion static correction method
CN106443829A (en) Method and apparatus for constructing near-surface model
CN106842314A (en) The determination method of formation thickness
CN105445795B (en) The velocity analysis method and its system of tilted interface
US20170176616A1 (en) System and method for dip-guided seismic image stacking
CN103064113A (en) Two-parameter expansion method and device confirming underground velocity structure

Legal Events

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