CN103439738B - Seismic prospecting single shot record exception road recognition methods - Google Patents

Seismic prospecting single shot record exception road recognition methods Download PDF

Info

Publication number
CN103439738B
CN103439738B CN201310379292.XA CN201310379292A CN103439738B CN 103439738 B CN103439738 B CN 103439738B CN 201310379292 A CN201310379292 A CN 201310379292A CN 103439738 B CN103439738 B CN 103439738B
Authority
CN
China
Prior art keywords
seismic channel
seismic
energy
road
channel
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
CN201310379292.XA
Other languages
Chinese (zh)
Other versions
CN103439738A (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
Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering 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 Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd filed Critical Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd
Priority to CN201310379292.XA priority Critical patent/CN103439738B/en
Publication of CN103439738A publication Critical patent/CN103439738A/en
Application granted granted Critical
Publication of CN103439738B publication Critical patent/CN103439738B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

Disclose a kind of seismic prospecting single shot record exception road recognition methods.Described method includes: receiving geological data and the positional information of each seismic channel of the multiple seismic channels produced by single-shot, the geological data of each of which seismic channel is by the seismic amplitude of predetermined sampling frequency;Calculate the energy value of each seismic channel;Determine the not working path in the plurality of seismic channel;Determine the abnormal big seismic channel of the energy in the seismic channel in addition to described not working path in the plurality of seismic channel;For each seismic channel in the seismic channel in addition to described not working path and the abnormal big seismic channel of described energy in the plurality of seismic channel, determine energy scale value and determine whether seismic channel is the abnormal road of son based on energy scale value;Export not working path and sub-exception road as final abnormal road.

Description

Seismic prospecting single shot record exception road recognition methods
Technical field
The present invention relates to oil seismic exploration and gather field, more particularly, relate to seismic recognition earthquake and survey The method visiting single shot record exception road.
Background technology
Along with the development of geophysical prospecting technology, when earthquake data acquisition has been enter into high accuracy, high density, high efficiency Phase, reception number of channels is many, job scheduling fast, and geological data amount presents exponential increase year by year, and three-dimensional is daily adopted Collection workload has reached more than 1000 big guns, and two dimension reaches more than 240 big guns.Road working condition is single big gun product The key factor of matter, therefore in traditional artificial evaluation methodology, the judgement to abnormal road is extremely important, The most time-consumingly, 240 records are arranged preliminary assessment from playback and reexamine supervision evaluation altogether Time-consuming 20 hours.Due to time-consuming long, efficiency is low, causes before second day gathers, it is impossible to the previous day The data gathered comprehensively is analyzed, it is difficult to meet the needs monitoring in real time and instructing field to produce.In order to Improve production efficiency, improve the accuracy that single shot record is evaluated, it is necessary to research and development single shot record exception road is certainly Dynamic identification technology.
Summary of the invention
In view of demand present in prior art, this application provides high efficiency and identify that earthquake is surveyed exactly The method visiting single shot record exception road.
In order to solve the above-mentioned problems in the prior art, the invention provides a kind of seismic prospecting list big gun The recognition methods of recording exceptional road, comprises the steps: to receive the multiple earthquakes produced by single-shot The geological data in road and the positional information of each seismic channel, the geological data of each of which seismic channel For the seismic amplitude by pre-fixed sample interval;Calculate the energy value of each seismic channel;Determine the plurality of Not working path in seismic channel;Determine the earthquake in addition to described not working path in the plurality of seismic channel The abnormal big seismic channel of energy in road;For in the plurality of seismic channel in addition to described not working path Seismic channel in each seismic channel, determine energy scale value and determine seismic channel based on energy scale value Whether it is the abnormal road of son;Export not working path and sub-exception road as final abnormal road.
Preferably, the step of the energy value calculating each seismic channel includes: the seismic amplitude to seismic channel The RMS amplitude that carries out calculate the energy value of seismic channel.
Preferably, determine that the step of the not working path in the plurality of seismic channel includes: be 0 by energy value And the seismic channel do not lost in actual production is defined as not working path.
Preferably, the energy in the seismic channel in addition to described not working path in the plurality of seismic channel is determined The step of amount the biggest abnormal seismic channel includes: when in the plurality of seismic channel in addition to described not working path Seismic channel in the energy value of each seismic channel more than at least in the energy value of neighbouring seismic channel During individual energy value prearranged multiple, this seismic channel is defined as the abnormal big seismic channel of energy.
Preferably, determine energy scale value and determine whether seismic channel is the abnormal road of son based on energy scale value Step include: in (a) seismic channel in addition to described not working path from the plurality of seismic channel, Select a seismic channel;B () calculates the seismic channel selected and excites the distance between position;(c) from institute State the seismic channel in addition to described not working path and the abnormal big seismic channel of described energy in multiple seismic channel In, determine and excite the distance of position at the seismic channel selected and to excite the first pre-of distance between position Seismic channel in the range of Ding;D () calculates the meansigma methods of the energy value of the seismic channel determined in step (c); (e) the energy value of the seismic channel selected at the second preset range of meansigma methods calculated time, by select Seismic channel is defined as the abnormal road of son;F () repeats step (a) to (e), until to the plurality of seismic channel In all seismic channels in addition to described not working path all carried out the judgement in sub abnormal road.
Seismic prospecting single shot record exception road recognition methods according to embodiments of the invention described above, can avoid The impact of the subjective factors of people.Additionally, when by software performing said method, labor can be alleviated Power also improves production efficiency.
Accompanying drawing explanation
By the description below in conjunction with the accompanying drawings embodiment carried out, these and/or other aspect of the present invention and Advantage will be made apparent from and it is more readily appreciated that wherein:
Fig. 1 is the seismic prospecting single shot record exception road identification illustrating the exemplary embodiment according to the present invention The flow chart of method;
Fig. 2 is to illustrate that in single big gun section, energy value is the diagram of the seismic channel of 0;
Fig. 3 is to illustrate the diagram of the abnormal big seismic channel of energy in single big gun section;
Fig. 4 is the diagram illustrating and judging the abnormal road of son.
Detailed description of the invention
Being described in detail the embodiment of the present invention now, its example represents in the accompanying drawings, wherein, and phase Same label represents same parts all the time.It is described embodiment below with reference to the accompanying drawings to explain this Invention.
Fig. 1 is the seismic prospecting single shot record exception road identification illustrating the exemplary embodiment according to the present invention The flow chart of method.
In step 110, receive the geological data of the multiple seismic channels produced by single-shot and every The positional information of one seismic channel, the geological data of each of which seismic channel is by pre-fixed sample interval Seismic amplitude.
As it is well known to the skilled in the art, single shot record is for excite in a precalculated position, and Multiple positions receive this and excite the vibration caused.In other words, each receiving position corresponds to a seismic channel, And the geological data of each seismic channel is the seismic amplitude sampled by pre-fixed sample interval.
In step 120, calculate the energy value of each seismic channel.Employing to the seismic amplitude of seismic channel Root-mean-square formula calculates the energy value of seismic channel.
It is assumed that geological data includes n seismic amplitude, respectively a1、a2、……、an, then The energy value in shake road is
σ = a 1 2 + a 2 2 + ...... + a n 2 n 2 .
In step 130, determine the not working path in the plurality of seismic channel.Here, can be 0 by energy value And the seismic channel do not lost in actual production is defined as not working path.As in figure 2 it is shown, by A in Fig. 2 The part represented is that in single big gun section, energy value is the seismic channel of 0.
In step 140, determine in the seismic channel in addition to described not working path in the plurality of seismic channel Energy the biggest abnormal seismic channel.Here, if in the plurality of seismic channel except described not working path it The energy value of a certain seismic channel in outer seismic channel is more than in the energy value of its neighbouring seismic channel at least The prearranged multiple (such as 10 times) of one energy value, then be defined as the abnormal big ground of energy by this seismic channel Shake road.As it is shown on figure 3, the part represented by B in Fig. 3 is the abnormal big ground of energy in single big gun section Shake road.
In step 150, in the seismic channel in addition to described not working path in the plurality of seismic channel Each seismic channel, determine energy scale value and determine that whether seismic channel is that son is different based on energy scale value Often road.
Here, the abnormal road of son can be determined as follows.
In (a) seismic channel in addition to described not working path from the plurality of seismic channel, select one Seismic channel.
B () calculates the seismic channel selected and excites the distance between position.
(c) from the plurality of seismic channel except described not working path and described energy the biggest abnormal seismic channel Outside seismic channel in, determine and excite the distance of position at the seismic channel selected and to excite between position Seismic channel in first preset range of distance.
D () calculates the meansigma methods of the energy value of the seismic channel determined in step (c) as energy scale Value.
(e) the energy value of seismic channel selected at energy scale value the second preset range time, by select Seismic channel is defined as the abnormal road of son.
F () repeats step (a) to (e), until to not working except described in the plurality of seismic channel All seismic channels outside road have all carried out the judgement in sub abnormal road.
As shown in Figure 4, if as seismic channel A1 is judged, calculating seismic channel A1 and excite position The distance of O is 2km, the distance seismic channel at about 2km ± 200m determining with exciting position, it is assumed that Determine that result is that seismic channel B1, C1, D1 meet above-mentioned distance condition.By seismic channel B1, C1, D1 The meansigma methods of energy value as the energy scale value of seismic channel A1, by this energy scale value and seismic channel The energy value of A1 compares, when the energy value of seismic channel A1 exceedes the second predetermined model of energy scale value When enclosing (such as, 0.1-10 times), seismic channel A1 will be confirmed as the abnormal road of son.
In step 160, export not working path and sub-exception road as final abnormal road.Selectively, may be used Output eliminates not working path and the seismic channel in sub-exception road.
Seismic prospecting single shot record exception road recognition methods according to embodiments of the invention described above, can avoid The impact of the subjective factors of people.Additionally, when by software performing said method, labor can be alleviated Power also improves production efficiency.
Although having show and described some embodiments of the present invention, it will be understood by those skilled in the art that Without departing from the principle of the present invention and the situation of spirit being limited its scope by claim and equivalent thereof Under, these embodiments can be modified.

Claims (3)

1. a seismic prospecting single shot record exception road recognition methods, it is characterised in that comprise the steps:
Receive the geological data of multiple seismic channels produced by single-shot and each seismic channel Positional information, the geological data of each of which seismic channel is the seismic amplitude by pre-fixed sample interval;
Calculate the energy value of each seismic channel;
Determine the not working path in the plurality of seismic channel;
Determine that the energy in the seismic channel in addition to described not working path in the plurality of seismic channel is abnormal big Seismic channel;
For each earthquake in the seismic channel in addition to described not working path in the plurality of seismic channel Road, determines energy scale value and determines whether seismic channel is the abnormal road of son based on energy scale value;
Export not working path and sub-exception road as final abnormal road,
Wherein, the step of the energy value calculating each seismic channel includes: enter the seismic amplitude of seismic channel Row RMS amplitude calculates, and calculates the energy value of seismic channel,
Wherein it is determined that energy scale value determine that whether seismic channel is the abnormal road of son based on energy scale value Step includes:
In (a) seismic channel in addition to described not working path from the plurality of seismic channel, select one Seismic channel;
B () calculates the seismic channel selected and excites the distance between position;
(c) from the plurality of seismic channel except described not working path and described energy the biggest abnormal seismic channel Outside seismic channel in, determine and excite the distance of position at the seismic channel selected and to excite between position Seismic channel in first preset range of distance;
D () calculates the meansigma methods of the energy value of the seismic channel determined in step (c);
E () is less than the 0.1 of the meansigma methods calculated or more than calculating at the energy value of the seismic channel selected During 10 times of meansigma methods, the seismic channel of selection is defined as the abnormal road of son;
F () repeats step (a) to (e), until to not working except described in the plurality of seismic channel All seismic channels outside road have all carried out the judgement in sub abnormal road.
Seismic prospecting single shot record exception road recognition methods the most according to claim 1, wherein, really The step of the not working path in fixed the plurality of seismic channel includes: is 0 by energy value and does not has in actual production The seismic channel in You Diu road is defined as not working path.
Seismic prospecting single shot record exception road recognition methods the most according to claim 1, wherein, really The abnormal big earthquake of the energy in the seismic channel in addition to described not working path in fixed the plurality of seismic channel The step in road includes:
Some seismic channel in the seismic channel in addition to described not working path in the plurality of seismic channel Energy value more than seismic channel adjacent thereto energy value at least one energy value prearranged multiple time, This seismic channel is defined as the abnormal big seismic channel of energy.
CN201310379292.XA 2013-08-27 2013-08-27 Seismic prospecting single shot record exception road recognition methods Active CN103439738B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310379292.XA CN103439738B (en) 2013-08-27 2013-08-27 Seismic prospecting single shot record exception road recognition methods

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310379292.XA CN103439738B (en) 2013-08-27 2013-08-27 Seismic prospecting single shot record exception road recognition methods

Publications (2)

Publication Number Publication Date
CN103439738A CN103439738A (en) 2013-12-11
CN103439738B true CN103439738B (en) 2016-09-28

Family

ID=49693442

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310379292.XA Active CN103439738B (en) 2013-08-27 2013-08-27 Seismic prospecting single shot record exception road recognition methods

Country Status (1)

Country Link
CN (1) CN103439738B (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104459779B (en) * 2014-12-09 2017-07-28 中国石油天然气集团公司 A kind of method and apparatus of the automatic identification in Unusual Seismic road
CN104865601B (en) * 2015-05-18 2017-08-25 中国石油集团川庆钻探工程有限公司地球物理勘探公司 The method for determining the type of the earthquake record gathered in seismic prospecting
CN109425894A (en) * 2017-08-23 2019-03-05 中国石油化工股份有限公司 A kind of seismic anomaly road detection method and device
CN107664778B (en) * 2017-08-24 2019-11-08 中国石油天然气股份有限公司 Method and device for determining seismic data quality
CN108181648B (en) * 2017-11-23 2019-07-09 中国石油天然气集团公司 The quality control method and device of continuous record aliased seismic acquisition data
CN107884826A (en) * 2017-12-13 2018-04-06 中国石油集团川庆钻探工程有限公司地球物理勘探公司 The method of dead trace recognition methods and optimization seismic data process in geological data
CN111123356B (en) * 2018-10-30 2022-07-12 中石化石油工程技术服务有限公司 Abnormal track intelligent identification method based on first arrival information
CN111596357B (en) * 2019-02-20 2023-11-28 中国石油天然气集团有限公司 Method and device for analyzing working state of submarine acquisition node
CN109932748A (en) * 2019-03-01 2019-06-25 中国石油天然气集团有限公司 A kind of earth surface consistency vibration amplitude compensation processing method, device and storage medium
CN111308551B (en) * 2020-03-12 2022-09-23 王美琪 Distance radius inner shot number mean value standard evaluation method
CN114114391B (en) * 2020-09-01 2023-08-22 中国石油天然气股份有限公司 Method and device for identifying negative samples of seismic data

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5265192A (en) * 1990-09-20 1993-11-23 Atlantic Richfield Company Method for the automated editing of seismic traces using an adaptive network
GB0015974D0 (en) * 2000-06-30 2000-08-23 Geco As Quality control of data
CN100349011C (en) * 2005-06-03 2007-11-14 中国石油集团东方地球物理勘探有限责任公司 Method for suppressing background noise not related to excitation source in seismic data enhancement
CN102193108B (en) * 2010-03-19 2013-07-03 中国石油天然气集团公司 Method for improving SNR (signal-to-noise ratio) of oil exploration data processing
CN102004264B (en) * 2010-10-18 2015-09-23 中国石油化工股份有限公司 A kind of earthquake-capturing data quality quantitative test and evaluation method

Also Published As

Publication number Publication date
CN103439738A (en) 2013-12-11

Similar Documents

Publication Publication Date Title
CN103439738B (en) Seismic prospecting single shot record exception road recognition methods
Poiata et al. Multiband array detection and location of seismic sources recorded by dense seismic networks
Fäh et al. ECOS-09 earthquake catalogue of Switzerland release 2011 report and database. Public catalogue, 17. 4. 2011. Swiss Seismological Service ETH Zurich
CN103365916B (en) Earthquake parameter estimates acquisition methods and system, seismic events search engine
CN103728659B (en) A kind of method improving detection of karst cave precision
CN100429527C (en) Three-dimensional seismic data processing quality monitoring technology
CN101900833B (en) Imaging velocity analysis method of seismic scattering P-P waves
CN101852863A (en) Method for processing seismic data by using high-precision single-channel spectrum analysis technology
US8639442B2 (en) Identifying invalid seismic data
CN104330823A (en) Method for determining observational parameters of vertical seismic profile
CN115327616B (en) Automatic positioning method for mine microseism focus driven by massive data
CN109581502B (en) Seismic exploration data compression method and device and compression storage structure
CN106405644B (en) Crack determination method and apparatus
CN101852864B (en) Method for processing mass seismic data by using surface consistent statistical spectrum analysis technology
CN107678057A (en) Determined in 3-D seismics or optimize wave detector or equipment investment amount and the method for arrangement
CN104133244B (en) A kind of near surface static corrections processing method of focus shear wave exploration
CN102692651A (en) Preliminary wave residual static correction method with space-variant velocity
CN105022091B (en) The far field focus method for rapidly positioning that a kind of nothing tests the speed in advance
CN113534259A (en) Vibroseis efficient acquisition real-time prestack time migration imaging method
CN102288995A (en) Seismic exploration collection quality quantitative analysis evaluation system
CN102053275B (en) Method for calculating relative statics correction value of combination in single-point earthquake chamber
CN103513279A (en) Lighting analysis calculation method and device based on wave equation of earthquake waves
CN105445788B (en) A kind of automatic means of interpretation of normal-moveout spectrum based on model and global optimizing
CN104865601B (en) The method for determining the type of the earthquake record gathered in seismic prospecting
CN103777244A (en) Method for quantitatively analyzing earthquake crack attribute volume

Legal Events

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

Effective date of registration: 20180208

Address after: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189

Patentee after: BGP INC., CHINA NATIONAL PETROLEUM Corp.

Address before: 610213 No. 1, No. 1, No. 1, Huayang Avenue, Huayang Town, Shuangliu County, Chengdu, Sichuan

Patentee before: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200918

Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No.

Patentee after: CHINA NATIONAL PETROLEUM Corp.

Patentee after: BGP Inc., China National Petroleum Corp.

Address before: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189

Patentee before: BGP Inc., China National Petroleum Corp.