CN103605158B - A kind of maximum offset defining method and device - Google Patents

A kind of maximum offset defining method and device Download PDF

Info

Publication number
CN103605158B
CN103605158B CN201310526104.1A CN201310526104A CN103605158B CN 103605158 B CN103605158 B CN 103605158B CN 201310526104 A CN201310526104 A CN 201310526104A CN 103605158 B CN103605158 B CN 103605158B
Authority
CN
China
Prior art keywords
road
amplitude energy
geological data
rms amplitude
maximum offset
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
CN201310526104.1A
Other languages
Chinese (zh)
Other versions
CN103605158A (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 National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
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 National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201310526104.1A priority Critical patent/CN103605158B/en
Publication of CN103605158A publication Critical patent/CN103605158A/en
Application granted granted Critical
Publication of CN103605158B publication Critical patent/CN103605158B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention provides a kind of maximum offset defining method and device, method comprises: gather each single big gun geological data in work area; Calculate each road geological data selecting the RMS amplitude energy in timing window; RMS amplitude energy is processed, generates new RMS amplitude energy: the total number of channels arranging reception after exponential curve fitting, continuation determination matching continuation are carried out to the new RMS amplitude energy in N-t+1 road to N road; Determine the av eff of the RMS amplitude energy variation of neighboring track geological data, according to av eff, new RMS amplitude energy determination amplitude energy, by the position in road maximum for amplitude energy as single big gun maximum offset, will there is the maximum offset of the maximum position of single big gun maximum offset number as described work area in work area.The present invention does not need to set up seismic model, considers the impact on maximum offset such as field interference wave, can determine suitable maximum offset, obtain high-quality field earthquake prospecting collection data.

Description

A kind of maximum offset defining method and device
Technical field
The present invention relates to geophysical prospecting for oil acquisition technique, particularly based on the method for real data zone of interest energy determination maximum offset, is a kind of maximum offset defining method and device concretely.
Background technology
In seismic data acquisition, the selection of maximum offset is the vital parameter of recording geometry design, because it affects the Disposal quality of latter earthquake data and explain effect in a lot.AndreasCordsen is the factor that will consider of detailed overview maximum offset design once: 1, maximum offset is approximately equal to the buried depth of the darkest zone of interest; 2, the reflection of fundamental purpose layer should avoid the interference of direct wave; 3, the reflection of fundamental purpose layer should avoid the interference of first break refraction; 4, maximum offset is less than the critical refraction geophone offset of the darkest zone of interest, and meets reflection coefficient stability requirement; 5, maximum offset should meet the geophone offset of the darkest low velocity layer (LVL) requirement; 6, velocity analysis permissible accuracy is reached; 7, NMO stretching distortion requirement is met; 8, the requirement of Multiple attenuation is met; 9, the requirement that AVO analyzes is met; 10, the ability (length as existing cable) etc. of data acquisition equipment is considered.As can be seen from these factors above-mentioned, the maximum offset that some factors dictate is larger, some then wishes that maximum offset is the smaller the better, and therefore designer selects a best maximum offset parameter with should taking into account the compromise of each factor.In recent years, a lot of scholar is also had to carry out detailed analysis discussion to above-mentioned factor.The selection of bear in 1992 gold is good when discussing existence " reversal point " from AVO curvilinear characteristic angle maximum offset; Fu in 2005 is excited and spectrum analysis towards Kui etc. by the ray tracing of geologic model, simulation, and driven school stretch distortion, velocity analysis accuracy requirement, interference wave excision and reflection coefficient the aspect such as to be stablized and discussed the impact of maximum offset on seismic imaging respectively; Luo Bin in 2005 etc. discuss the impact of maximum offset and lateral resolution vertical on earthquake according to NMO stretching and Fresnel zone size.Except the above-mentioned comprehensive qualitative analysis to maximum offset selection, the quantitative calculating of a lot of scholar to maximum offset also conducts in-depth research.Gu Peicheng in 1998 etc. based on the relation between maximum offset and nmo stretching distortion, and consider that excision has extrapolated the practical formula of an estimation maximum offset to the restriction of geophone offset; Chen Hao woodss in 2009 etc., based on the concept of generalized seismic spatial resolution, establish the quantitative expression formula between a maximum offset and the zone of interest degree of depth.Be not difficult to find, above scholar is the research carrying out maximum offset method for designing based on HORIZONTAL LAYERED MEDIUM WITH HIGH ACCURACY assumed condition mostly.When the vertical, horizontal change of earth's surface or subsurface formations is violent, the recording geometry parameter designed according to a conventional method is generally difficult to the requirement meeting seismic resolution and imaging.
Based on the deficiency that said method design maximum offset exists, people start to utilize seismic wave illumination energy distribution to carry out design and study seismic acquisition parameter.2006 Nian Dongliang states etc. by continuing to increase spread length to observe the contribution to zone of interest general lighting energy, thus determine an effective spread length scope; Zhao Hu in 2011 etc., on the basis that seismic illumination is simulated, by probing into the relation between wave detector received energy and receiving efficiency, define the maximum offset design receiving illumination energy maximizing efficiency based on wave detector; The people such as Xu Feng in 2012, on the basis of seismic illumination, build geologic model and carry out numerical simulation, receive the maximum and most effective principle of zone of interest energy, accumulation energy and energy derivative are combined according to recording geometry, design maximum offset.Above-mentioned several method based on purpose of model layer energy determination maximum offset, needs to set up geologic model, is difficult to the requirement meeting surveying accuracy in practical application, and simultaneously these methods do not consider the impact of noise jamming ripple in actual excitation process etc.
In sum, the method of design maximum offset has a lot, but it is severe for earth's surface geographical environment in current seismic prospecting, geology of deep part complicated condition, the underground structure horizontal change complex area such as acutely, the defect of existing method all various degrees, is usually difficult to obtain suitable maximum offset.The choice relation of maximum offset, to the effect of earthquake-capturing, designs suitable maximum offset and seems and be even more important.
Summary of the invention
The invention provides a kind of maximum offset defining method, comprising:
Gather single big gun geological data in work area;
The each road geological data calculating described single big gun geological data is selecting the RMS amplitude energy in timing window;
Process according to the RMS amplitude energy of formula (1) to described each road geological data, generate new RMS amplitude energy, formula (1) is:
b i=A i/2+(A i+1+A i-1)/4i=2,…,N-1（1）
Wherein, i is the road sequence number of described single big gun geological data arrangement reception channel, and N is the reception number of channels of single big gun geological data arrangement, A iwei Zhong i-th road geological data is selecting the RMS amplitude energy in timing window;
Total number of channels of reception is arranged after exponential curve fitting, continuation determination matching continuation are carried out to the new RMS amplitude energy in described N-t+1 road to N road;
Determine the av eff of the RMS amplitude energy variation of the neighboring track geological data that each road geological data is corresponding according to formula (2), formula (2) is:
Q i = 1 i Σ m = 1 i | b m - b m + 1 | , i = 2 , . . . , M - 1 - - - ( 2 )
Wherein, M is the total number of channels arranging reception after matching continuation, and m is sequence number, 1≤m≤i;
Determine the amplitude energy of each road geological data according to formula (3) and the av eff determined, new RMS amplitude energy, formula (3) is:
C i = ( 1 i + E i E i - 1 ) E i + Q i , i = 2 , . . . , M - - - ( 3 )
Wherein, E ifor the b in the step 3) that the 1st Dao Dao i-th road receiving arrangement is corresponding mand, namely m is sequence number, 1≤m≤i;
By the position in the road of one maximum for the described amplitude energy determined geological data as single big gun maximum offset;
To there is the maximum offset of the maximum position of single big gun maximum offset number as described work area in work area.
Meanwhile, present invention also offers a kind of maximum offset determining device, comprising:
Acquisition module, for gathering single big gun geological data in work area;
RMS amplitude energy computation module, for each road geological data of calculating described single big gun geological data selecting the RMS amplitude energy in timing window;
New RMS amplitude energy computation module, for processing according to the RMS amplitude energy of formula (1) to described each road geological data, generate new RMS amplitude energy, formula (1) is:
b i=A i/2+(A i+1+A i-1)/4i=2,…,N-1（1）
Wherein, i is the road sequence number of described single big gun geological data arrangement reception channel, and N is the reception number of channels of single big gun geological data arrangement, A iwei Zhong i-th road geological data is selecting the RMS amplitude energy in timing window;
Matching continuation module, for arranging total number of channels of reception after carrying out exponential curve fitting, continuation determination matching continuation to the new RMS amplitude energy in described N-t+1 road to N road;
Av eff computing module, for determining the av eff of the RMS amplitude energy variation of the neighboring track geological data that each road geological data is corresponding according to formula (2), formula (2) is:
Q i = 1 i Σ m = 1 i | b m - b m + 1 | , i = 2 , . . . , M - 1 - - - ( 2 )
Wherein, M is the total number of channels arranging reception after matching continuation, and m is sequence number, 1≤m≤i;
Amplitude energy computation module, for determining the amplitude energy of each road geological data according to formula (3) and the av eff determined, new RMS amplitude energy, formula (3) is:
C i = ( 1 i + E i E i - 1 ) E i + Q i , i = 2 , . . . , M - - - ( 3 )
Wherein, E ifor the b in the step 3) that the 1st Dao Dao i-th road receiving arrangement is corresponding mand, namely m is sequence number, 1≤m≤i;
Single big gun maximum offset determination module, for the position in the road by one maximum for the described amplitude energy determined geological data as single big gun maximum offset;
Maximum offset determination module, will have the maximum offset of the maximum position of single big gun maximum offset number as described work area in work area.
The present invention does not need to set up seismic model, consider the impact on maximum offset such as field interference wave, the complex area severe to earth's surface geographical environment in seismic prospecting, geology of deep part complicated condition, underground structure horizontal change are violent etc., suitable maximum offset can be determined, thus obtain high-quality field earthquake prospecting collection data.
For above and other object of the present invention, feature and advantage can be become apparent, preferred embodiment cited below particularly, and coordinate institute's accompanying drawings, be described in detail below.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.
Fig. 1 is the process flow diagram of a kind of maximum offset defining method disclosed by the invention;
Fig. 2 is the structured flowchart of a kind of maximum offset determining device disclosed by the invention;
Fig. 3 is the original single big gun seismologic record in land somewhere;
The RMS amplitude energy trace figure of right arrangement reception channel in window when Fig. 4 is Fig. 3;
Fig. 5 is the amplitude energy curve map after the continuation of Fig. 4 RMS amplitude energy trace
Fig. 6 is the amplitude energy curve map of the av eff of the neighboring track RMS amplitude energy variation in all roads in Fig. 4;
Fig. 7 is amplitude energy curve map in the embodiment of the present invention;
Fig. 8 is the distribution plan of the maximum offset of all big guns.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, be clearly and completely described the technical scheme in the embodiment of the present invention, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to the scope of protection of the invention.
Core of the present invention first calculates to receive the every RMS amplitude together of arrangement selecting in timing window, then the RMS amplitude energy receiving the road of specifying in arrangement is utilized to carry out exponential curve fitting, according to the average root-mean-square amplitude energy determination matched curve continuation length of background noise, calculate the av eff of every neighboring track RMS amplitude energy variation corresponding together again, finally utilize each trace integration amplitude energy of new formulae discovery, determine the road sequence number that maximal value is corresponding, utilize the distribution of the maximum offset of all big guns to determine the maximum offset in work area.
As shown in Figure 1, the present invention is disclosed a kind of maximum offset defining method, comprising:
Step S101, gathers single big gun geological data in work area;
Step S102, each road geological data calculating described single big gun geological data is selecting the RMS amplitude energy in timing window;
Step S103, process according to the RMS amplitude energy of formula (1) to described each road geological data, generate new RMS amplitude energy, formula (1) is:
b i=A i/2+(A i+1+A i-1)/4i=2,…,N-1（1）
Wherein, i is the road sequence number of described single big gun geological data arrangement reception channel, and N is the reception number of channels of single big gun geological data arrangement, A iwei Zhong i-th road geological data is selecting the RMS amplitude energy in timing window;
Step S104, arranges total number of channels of reception to the new RMS amplitude energy in described N-t+1 road to N road after carrying out exponential curve fitting, continuation determination matching continuation;
Step S105, determine the av eff of the RMS amplitude energy variation of the neighboring track geological data that each road geological data is corresponding according to formula (2), formula (2) is:
Q i = 1 i Σ m = 1 i | b m - b m + 1 | , i = 2 , . . . , M - 1 - - - ( 2 )
Wherein, M is the total number of channels arranging reception after matching continuation, and m is sequence number, 1≤m≤i;
Step S106, determine the amplitude energy of each road geological data according to formula (3) and the av eff determined, new RMS amplitude energy, formula (3) is:
C i = ( 1 i + E i E i - 1 ) E i + Q i , i = 2 , . . . , M - - - ( 3 )
Wherein, E ifor the b in the step 3) that the 1st Dao Dao i-th road receiving arrangement is corresponding mand, namely m is sequence number, 1≤m≤i;
Step S107, by the position in road maximum for amplitude energy in the amplitude energy of described each the road geological data determined as single big gun maximum offset;
Step S108, will have the maximum offset of the maximum position of single big gun maximum offset number as described work area in work area.
Meanwhile, present invention also offers a kind of maximum offset determining device as shown in Figure 2, comprising:
Acquisition module 201, for gathering single big gun geological data in work area;
RMS amplitude energy computation module 202, for each road geological data of calculating described single big gun geological data selecting the RMS amplitude energy in timing window;
New RMS amplitude energy computation module 203, for processing according to the RMS amplitude energy of formula (1) to described each road geological data, generate new RMS amplitude energy, formula (1) is:
b i=A i/2+(A i+1+A i-1)/4i=2,…,N-1（1）
Wherein, i is the road sequence number of described single big gun geological data arrangement reception channel, and N is the reception number of channels of single big gun geological data arrangement, A iwei Zhong i-th road geological data is selecting the RMS amplitude energy in timing window;
Matching continuation module 204, for arranging total number of channels of reception after carrying out exponential curve fitting, continuation determination matching continuation to the new RMS amplitude energy in described N-t+1 road to N road;
Av eff computing module 205, for determining the av eff of the RMS amplitude energy variation of the neighboring track geological data that each road geological data is corresponding according to formula (2), formula (2) is:
Q i = 1 i Σ m = 1 i | b m - b m + 1 | , i = 2 , . . . , M - 1 - - - ( 2 )
Wherein, M is the total number of channels arranging reception after matching continuation, and m is sequence number, 1≤m≤i;
Amplitude energy computation module 206, for determining the amplitude energy of each road geological data according to formula (3) and the av eff determined, new RMS amplitude energy, formula (3) is:
C i = ( 1 i + E i E i - 1 ) E i + Q i , i = 2 , . . . , M - - - ( 3 )
Wherein, E ifor the b in the step 3) that the 1st Dao Dao i-th road receiving arrangement is corresponding mand, namely m is sequence number, 1≤m≤i;
Single big gun maximum offset determination module 207, for by the position in road maximum for amplitude energy in the amplitude energy of described each the road geological data determined as single big gun maximum offset;
Maximum offset determination module 208, will have the maximum offset of the maximum position of single big gun maximum offset number as described work area in work area.
Core of the present invention first calculates to receive the every RMS amplitude together of arrangement selecting in timing window, then the RMS amplitude energy receiving the road of specifying in arrangement is utilized to carry out exponential curve fitting, according to the average root-mean-square amplitude energy determination matched curve continuation length of background noise, calculate the av eff of every neighboring track RMS amplitude energy variation corresponding together again, finally utilize each trace integration amplitude energy of new formulae discovery, determine the road sequence number that maximal value is corresponding, utilize the distribution of the maximum offset of all big guns to determine the maximum offset in work area.
Technology provided by the invention mainly calculates to be specified the RMS amplitude energy in road to carry out exponential curve fitting continuation in the RMS amplitude energy selecting every one in timing window, utilization, calculate the av eff of the neighboring track RMS amplitude energy variation of every one correspondence, calculates each trace integration amplitude energy, determines maximum offset.
In embodiment, the present invention is achieved through the following technical solutions:
1) also acquisition and recording list big gun geological data is excited;
2) on seismologic record, used time window selectes zone of interest region, calculates every together selecting the RMS amplitude energy A in timing window i, i is more than or equal to 1 and is less than or equal to N, and i is the road sequence number of single big gun geological data arrangement reception channel, and N is the number of channels that the arrangement of single big gun geological data receives.
Step 2) in, for unilateral shot, the road receiving arrangement middle distance shot point nearest is the 1st road, sort successively according to offset distance, receiving road farthest in arrangement is N road, for split shooting, centered by shot point, whole arrangement is divided into two arrangements in left and right, and the arrangement of two, left and right calculates according to unilateral shot mode respectively;
Calculate background interference noise files every together selecting the RMS amplitude energy in timing window, then to the RMS amplitude energy in all roads of background interference noise files and divided by record number of channels, draw the average root-mean-square amplitude energy n of ground unrest 0;
3) utilize following formula to every one RMS amplitude energy A iprocess, the corresponding RMS amplitude energy b that must make new advances i:
b i=A i/2+(A i+1+A i-1)/4i=2,…,N-1（1）
In formula: i is the road sequence number of single big gun geological data arrangement reception channel, and i is more than or equal to 2 and is less than or equal to N-1, N is the number of channels that the arrangement of single big gun geological data receives, A ifor step 2) in the i-th road selecting the RMS amplitude energy in timing window;
In step 3), acquiring seismic data is received to the 1st road b of arrangement 1value equal b 2value, N road b nvalue equal b n-1value;
4) to the last t road b of the new RMS amplitude energy in step 3) n-t+1, b n-t+2..., b n-1, b ncarry out exponential curve fitting, carry out continuation to the curve after matching according to offset distance, continuation is not less than the average root-mean-square amplitude energy n of ground unrest record to the RMS amplitude energy arranging reception channel 0, the total number of channels arranging reception after continuation is M;
Last t road described in step 4) carries out exponential curve fitting, and t defaults to 20;
5) the av eff Q of every neighboring track RMS amplitude energy variation corresponding is together calculated with following formula i:
Q i = 1 i Σ m = 1 i | b m - b m + 1 | , i = 2 , . . . , M - 1 - - - ( 2 )
In formula: i is the road sequence number of single big gun geological data arrangement reception channel, and M is the total number of channels arranging reception after matching continuation; The 1st road Q that arrangement described in step 5) receives 1value equal Q 2value, M road Q mvalue equal Q m-1value;
6) with the Q in following formula and step 5) icalculate every amplitude energy C together i:
C i = ( 1 i + E i E i - 1 ) E i + Q i , i = 2 , . . . , M - - - ( 3 )
In formula: i is the road sequence number of single big gun geological data arrangement reception channel, and M is the total number of channels arranging reception after matching continuation, E ifor the b in the step 3) that the 1st Dao Dao i-th road receiving arrangement is corresponding mand, namely m is sequence number, and m is more than or equal to 1 and is less than or equal to i;
Step 6) in, the C that the 1st road is corresponding 1value equal the 2nd road C 2value;
7) to ordered series of numbers C in step 6) i(i=1,2 ..., M) and sequence, find out maximum C icorresponding road sequence number, the position of the maximum offset of the big gun will looked for exactly;
8) according to step 2) method of ~ step 7) calculates the distribution of the maximum offset of all big guns, and the maximum position of corresponding geophone offset number is the position that maximum offset is corresponding;
Below in conjunction with specific embodiment, technical solution of the present invention is described in further details.
1, excite also acquisition and recording list big gun geological data, mode of excitation is split shooting, and as shown in Figure 3, left-right situs Jun You 100 road receives, track pitch 30m;
2, select zone of interest region shaking used time window on record as shown in Figure 3, time window be the hyperbolic curve shown in Fig. 3, calculate the maximum offset receiving arrangement on the right of the 5163rd shot point below, the right receives arrangement Zhong You 100 road, selecting the RMS amplitude A in timing window 1=9.59305, A 2=9.37908 ..., A 99=1.96693, A 100=1.74396.The average RMS amplitude n of ground unrest record 0=0.541892;
3, utilize formula (1) to every one RMS amplitude A iprocess, the corresponding RMS amplitude b that must make new advances 1=9.37078, b 2=9.37078 ..., b 98=2.10295, b 99=1.94672, b 100=1.94672.New RMS amplitude energy trace as shown in Figure 4;
4, exponential curve fitting is carried out to the RMS amplitude energy in last 20 roads that step 3 is obtained, carry out continuation to the curve after matching according to offset distance, continuation is not less than the average root-mean-square amplitude energy n of ground unrest record to the RMS amplitude energy arranging reception channel 0=0.541892, arranging total reception number of channels after continuation is 119 roads, as shown in Figure 5;
5, the av eff of every neighboring track RMS amplitude energy variation corresponding is together calculated with formula (2): Q 1=0.252114, Q 2=0.252114 ..., Q 117=0.175304, Q 118=0.17407, Q 119=0.17407, amplitude energy curve as shown in Figure 6.
6, with formula 3) and step 5 in Q icalculate each C i, i is more than or equal to 1 and is less than or equal to 119, from C iin find out maximal value, each C idivided by the maximal value found out, show that the value of every amplitude energy corresponding is together respectively: 0.842715,0.842715,0.646665 ..., 0.996188,0.995841, the curve that in 0.995438, Fig. 7, black circle is formed is the curve of each trace integration amplitude energy divided by amplitude energy income value maximum in all roads, the curve that real triangle form point is formed be in Fig. 5 every amplitude energy together divided by the curve of amplitude energy income value maximum in all roads.
7, to the value sequence of amplitude energy corresponding together every in step 6, maximal value is the 107th road, in addition, owing to being split shooting, shot point is between two geophone stations, therefore need to add a 15m again, therefore corresponding maximum offset is 107 × 30+15=3225m, position as shown in vertical curve in Fig. 7;
8, calculate the distribution of the maximum offset of all big guns according to the method for step 2 ~ step 7, the geophone offset number that the 107th road is corresponding is maximum, and maximum offset is 3225m, as shown in Figure 8.
The present invention mainly calculates and specifies the RMS amplitude energy in road to carry out exponential curve fitting continuation in the RMS amplitude energy selecting every one in timing window, utilization, calculate the av eff of the neighboring track RMS amplitude energy variation of every one correspondence, calculates each trace integration amplitude energy, determines maximum offset.The present invention does not need to set up seismic model, consider the impact on maximum offset such as field interference wave, the complex area severe to earth's surface geographical environment in seismic prospecting, geology of deep part complicated condition, underground structure horizontal change are violent etc., suitable maximum offset can be determined, thus obtain high-quality field earthquake prospecting collection data.
Apply specific embodiment in the present invention to set forth principle of the present invention and embodiment, the explanation of above embodiment just understands method of the present invention and core concept thereof for helping; Meanwhile, for one of ordinary skill in the art, according to thought of the present invention, all will change in specific embodiments and applications, in sum, this description should not be construed as limitation of the present invention.

Claims (12)

1. a maximum offset defining method, is characterized in that, described method comprises:
Gather each single big gun geological data in work area;
The each road geological data calculating described single big gun geological data is selecting the RMS amplitude energy in timing window;
Process according to the RMS amplitude energy of formula (1) to described each road geological data, generate new RMS amplitude energy, formula (1) is:
b i＝A i/2+(A i+1+A i-1)/4i＝2,…,N-1(1)
Wherein, i is the road sequence number of described single big gun geological data arrangement reception channel, and N is the reception number of channels of single big gun geological data arrangement, A iwei Zhong i-th road geological data is selecting the RMS amplitude energy in timing window;
Total number of channels of reception is arranged after exponential curve fitting, continuation determination matching continuation are carried out to the new RMS amplitude energy in described N-t+1 road to N road;
Determine the av eff of the RMS amplitude energy variation of the neighboring track geological data that each road geological data is corresponding according to formula (2), formula (2) is:
Q i = 1 i Σ m = 1 i | b m - b m + 1 | , i = 2 , ... , M - 1 - - - ( 2 )
Wherein, M is the total number of channels arranging reception after matching continuation, and m is sequence number, 1≤m≤i;
Determine the amplitude energy of each road geological data according to formula (3) and the av eff determined, new RMS amplitude energy, formula (3) is:
C i = ( 1 i + E i E i - 1 ) E i + Q i , i = 2 , ... , M - - - ( 3 )
Wherein, E ifor receiving the step 3 corresponding to the 1st Dao Dao i-th road of arrangement) in b mand, namely m is sequence number, 1≤m≤i;
By the position in the road of one maximum for the described amplitude energy determined geological data as single big gun maximum offset;
To there is the maximum offset of the maximum position of single big gun maximum offset number as described work area in work area;
Wherein, for unilateral shot, receiving road nearest apart from shot point in arrangement is the 1st road, and sort successively according to offset distance, road is farthest N road; For split shooting, centered by shot point, whole arrangement is divided into two arrangements in left and right, and the arrangement of two, left and right calculates according to unilateral shot mode respectively.
2. maximum offset defining method as claimed in claim 1, is characterized in that, described exponential curve fitting, continuation determination matching continuation are carried out to the new RMS amplitude energy in N-t+1 road to N road after arrange reception total number of channels comprise:
Exponential curve fitting is carried out to the new RMS amplitude energy in N-t+1 road to N road and generates the curve after matching;
Carry out continuation to the curve after matching according to offset distance, continuation is not less than the average root-mean-square amplitude energy of ground unrest excitation to the RMS amplitude energy arranging reception channel, determines total number of channels that the arrangement after continuation receives.
3. maximum offset defining method as claimed in claim 2, it is characterized in that, t defaults to 20.
4. maximum offset defining method as claimed in claim 1, it is characterized in that, described step processes according to the RMS amplitude energy of formula (1) to described each road geological data, generates in new RMS amplitude energy,
Acquiring seismic data receives the 1st road b of arrangement 1value equal b 2value, N road b nvalue equal b n-1value.
5. maximum offset defining method as claimed in claim 1, it is characterized in that, described determines in the av eff of RMS amplitude energy variation of the neighboring track geological data that each road geological data is corresponding according to formula (2),
The 1st road Q that arrangement receives 1value equal Q 2value, M road Q mvalue equal Q m-1value.
6. maximum offset defining method as claimed in claim 1, it is characterized in that, described determines in the amplitude energy of each road geological data according to formula (3) and the av eff determined, new RMS amplitude energy,
The C that 1st road is corresponding 1value equal the 2nd road C 2value.
7. a maximum offset determining device, is characterized in that, described device comprises:
Acquisition module, for gathering single big gun geological data in work area;
RMS amplitude energy computation module, for each road geological data of calculating described single big gun geological data selecting the RMS amplitude energy in timing window;
New RMS amplitude energy computation module, for processing according to the RMS amplitude energy of formula (1) to described each road geological data, generate new RMS amplitude energy, formula (1) is:
b i＝A i/2+(A i+1+A i-1)/4i＝2,…,N-1(1)
Wherein, i is the road sequence number of described single big gun geological data arrangement reception channel, and N is the reception number of channels of single big gun geological data arrangement, A iwei Zhong i-th road geological data is selecting the RMS amplitude energy in timing window;
Matching continuation module, for arranging total number of channels of reception after carrying out exponential curve fitting, continuation determination matching continuation to the new RMS amplitude energy in described N-t+1 road to N road;
Av eff computing module, for determining the av eff of the RMS amplitude energy variation of the neighboring track geological data that each road geological data is corresponding according to formula (2), formula (2) is:
Q i = 1 i Σ m = 1 i | b m - b m + 1 | , i = 2 , ... , M - 1 - - - ( 2 )
Wherein, M is the total number of channels arranging reception after matching continuation, and m is sequence number, 1≤m≤i;
Amplitude energy computation module, for determining the amplitude energy of each road geological data according to formula (3) and the av eff determined, new RMS amplitude energy, formula (3) is:
C i = ( 1 i + E i E i - 1 ) E i + Q i , i = 2 , ... , M - - - ( 3 )
Wherein, E ifor receiving the step 3 corresponding to the 1st Dao Dao i-th road of arrangement) in b mand, namely m is sequence number, 1≤m≤i;
Single big gun maximum offset determination module, for the position in the road by one maximum for the described amplitude energy determined geological data as single big gun maximum offset;
Maximum offset determination module, will have the maximum offset of the maximum position of single big gun maximum offset number as described work area in work area;
Road divides module: for unilateral shot, and receiving road nearest apart from shot point in arrangement is the 1st road, and sort successively according to offset distance, road is farthest N road; For split shooting, centered by shot point, whole arrangement is divided into two arrangements in left and right, and the arrangement of two, left and right calculates according to unilateral shot mode respectively.
8. maximum offset determining device as claimed in claim 7, it is characterized in that, described matching continuation module comprises:
Fitting unit, generates the curve after matching for carrying out exponential curve fitting to the new RMS amplitude energy in N-t+1 road to N road;
Continuation unit, for carrying out continuation to the curve after matching according to offset distance, continuation is not less than the average root-mean-square amplitude energy of ground unrest excitation to the RMS amplitude energy arranging reception channel, determines total number of channels that the arrangement after continuation receives.
9. maximum offset determining device as claimed in claim 8, it is characterized in that, t defaults to 20.
10. maximum offset determining device as claimed in claim 7, it is characterized in that, described new RMS amplitude energy determination module processes according to the RMS amplitude energy of formula (1) to described each road geological data, generates in new RMS amplitude energy
Acquiring seismic data receives the 1st road b of arrangement 1value equal b 2value, N road b nvalue equal b n-1value.
11. maximum offset determining devices as claimed in claim 7, it is characterized in that, described av eff computing module is determined in the av eff of RMS amplitude energy variation of the neighboring track geological data that each road geological data is corresponding according to formula (2)
The 1st road Q that arrangement receives 1value equal Q 2value, M road Q mvalue equal Q m-1value.
12. maximum offset determining devices as claimed in claim 7, it is characterized in that, described amplitude energy computation module is determined in the amplitude energy of each road geological data according to formula (3) and the av eff determined, new RMS amplitude energy
The C that 1st road is corresponding 1value equal the 2nd road C 2value.
CN201310526104.1A 2013-10-30 2013-10-30 A kind of maximum offset defining method and device Active CN103605158B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310526104.1A CN103605158B (en) 2013-10-30 2013-10-30 A kind of maximum offset defining method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310526104.1A CN103605158B (en) 2013-10-30 2013-10-30 A kind of maximum offset defining method and device

Publications (2)

Publication Number Publication Date
CN103605158A CN103605158A (en) 2014-02-26
CN103605158B true CN103605158B (en) 2016-04-06

Family

ID=50123398

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310526104.1A Active CN103605158B (en) 2013-10-30 2013-10-30 A kind of maximum offset defining method and device

Country Status (1)

Country Link
CN (1) CN103605158B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104133239A (en) * 2014-07-29 2014-11-05 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Method for defining earthquake single-shot record target stratum time window
CN105093296A (en) * 2015-06-30 2015-11-25 中国石油天然气集团公司 Method and device for optimizing observation system
CN105259568B (en) * 2015-10-19 2017-09-15 中国石油天然气集团公司 A kind of method and apparatus of determination exploratory area shot point maximum offset
CN107085236B (en) * 2017-03-23 2019-02-15 中国石油天然气集团公司 The determination method and apparatus of maximum offset
CN109738944B (en) * 2019-03-05 2020-05-08 中国石油大学(北京) Wide-angle reflection-based seismic acquisition parameter determination method and device
CN114814931A (en) * 2022-04-23 2022-07-29 中国海洋大学 Observation system arrangement length selection method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6243650B1 (en) * 1998-09-11 2001-06-05 Diamond Geoscience Research Corporation Method of determining net reservoir thickness
US6263284B1 (en) * 1999-04-22 2001-07-17 Bp Corporation North America Inc. Selection of seismic modes through amplitude characteristics
FR2858064A1 (en) * 2003-07-21 2005-01-28 Geophysique Cie Gle Determining parameters of speed and anellipticity, for treating seismic traces from collection at common midpoint including anelliptic obliquity correction

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6757217B2 (en) * 2001-08-31 2004-06-29 Exxonmobil Upstream Research Company Method for time-aligning multiple offset seismic data volumes

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6243650B1 (en) * 1998-09-11 2001-06-05 Diamond Geoscience Research Corporation Method of determining net reservoir thickness
US6263284B1 (en) * 1999-04-22 2001-07-17 Bp Corporation North America Inc. Selection of seismic modes through amplitude characteristics
FR2858064A1 (en) * 2003-07-21 2005-01-28 Geophysique Cie Gle Determining parameters of speed and anellipticity, for treating seismic traces from collection at common midpoint including anelliptic obliquity correction

Also Published As

Publication number Publication date
CN103605158A (en) 2014-02-26

Similar Documents

Publication Publication Date Title
CN103605158B (en) A kind of maximum offset defining method and device
CN108037526B (en) Reverse-time migration method based on all-wave wave field VSP/RVSP seismic data
CN108957549B (en) Braided river sediment heterogeneous compact sandstone gas reservoir geological modeling method
CN106094029B (en) Utilize the method for offset distance vector piece geological data Predicating Reservoir Fractures
CN103091710B (en) Reverse time migration imaging method and device
CN104614765B (en) Design method for enhancing seismic waves to stimulate illumination
CN102305941B (en) Method for determining stratum stack quality factor by direct scanning of prestack time migration
CN102645670B (en) Observation system optimization design method based on stack response analysis
CN103645503B (en) A kind of three-dimensional time territory illumination analysis and vibration amplitude compensation method
CN103592698B (en) A kind of recording geometry evaluation method based on seismic properties and device
CN104570125A (en) Method for improving imaging velocity model precision by utilizing well data
CN102419453A (en) Long-lead source transient electromagnetic ground-air detecting method
CN112883564A (en) Water body temperature prediction method and prediction system based on random forest
CN104345336B (en) Observation system optimizing method based on target area illumination level
CN104570116A (en) Geological marker bed-based time difference analyzing and correcting method
CN101852864B (en) Method for processing mass seismic data by using surface consistent statistical spectrum analysis technology
CN107678057B (en) Wave detector or equipment investment amount and the method for arrangement are determined or optimized in 3-D seismics
CN101937102A (en) Focusing performance analysis method of three-dimensional observation system
Guan et al. Selection of noise sources and short-time passive surface wave imaging——A case study on fault investigation
Jun et al. Repeatability enhancement of time-lapse seismic data via a convolutional autoencoder
CN104133244B (en) A kind of near surface static corrections processing method of focus shear wave exploration
CN102636813A (en) Normal moveout stretch cutting method for processing geophysical exploitation seismic data
CN105785441A (en) Signal-to-noise ratio analysis method for seismic data
CN103645501A (en) A method for determining a largest offset
CN103513279B (en) A kind of illumination analysis computing method based on seismic wave equation and calculation element

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant