CN106970418B - Method and device for acquiring earthquake synthetic record - Google Patents
Method and device for acquiring earthquake synthetic record Download PDFInfo
- Publication number
- CN106970418B CN106970418B CN201710213246.0A CN201710213246A CN106970418B CN 106970418 B CN106970418 B CN 106970418B CN 201710213246 A CN201710213246 A CN 201710213246A CN 106970418 B CN106970418 B CN 106970418B
- Authority
- CN
- China
- Prior art keywords
- reflection
- reflecting interface
- time
- domain
- wave
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 230000005540 biological transmission Effects 0.000 claims abstract description 23
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 230000000750 progressive effect Effects 0.000 claims description 6
- 238000003786 synthesis reaction Methods 0.000 claims description 6
- 238000010304 firing Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 22
- 238000003860 storage Methods 0.000 description 10
- 238000004590 computer program Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 8
- 239000002131 composite material Substances 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000005389 magnetism Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001615 p wave Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The embodiment of the application provides a method and a device for acquiring a seismic synthetic record, wherein the method comprises the following steps: determining the time domain primary wave reflection coefficient of each reflection interface in a target work area; determining downlink reflection control coefficients of all the reflection interfaces; the downlink reflection control coefficient is used for controlling whether each reflection interface generates multiple reflection or not; acquiring a time domain second reflection coefficient containing transmission loss and multiple reflection according to the downlink reflection control coefficient and the time domain primary wave reflection coefficient; and determining a seismic wavelet sequence, and performing convolution on the seismic wavelet sequence and the time domain second reflection coefficient to obtain a seismic synthetic record. The method and the device can improve the precision of the earthquake synthetic record and help to identify the multiple waves in the actual earthquake record.
Description
Technical field
This application involves the seismic forward simulation technical fields more particularly to a kind of acquisition earthquake conjunction in geophysical exploration
At the method and device of record.
Background technique
Synthetic seismic record is the intermediary for geological model being converted into earthquake information, i.e., with acoustic logging or vertically
It shakes sectional data synthetic seismic record (seismic channel).Synthetic seismic record is very widely used one in seismic modeling technique
The basis of the work such as kind and horizon calibration, reservoir description.Therefore, the precision of synthetic seismic record is with directly influencing earthquake
The accurate calibration of matter layer position.In general, synthetic seismic record is seismic forward simulation as a result, so-called seismic forward simulation is
Simplification appropriate is done to specific geologic body, forms simplified a data model or physical model, the side calculated using data
Method or physical simulating method obtain the process of seismic response.
In seismic prospecting, multiple wave often generates many negative effects.Such as multiple wave will affect final earthquake number
According to offset mass；Multiple wave is also possible to perplex to subsequent explanation work belt, leads to the layer position explanation results and mistake of mistake
Seismic inversion.However current synthetic seismic record method does not account for the influence of multiple wave, so as to cause synthesis
Earthquake record precision it is not high, be unfavorable for subsequent horizon calibration and Forecast Oil Reservoir Distribution.
Summary of the invention
The embodiment of the present application is designed to provide a kind of method and device for obtaining synthetic seismic record, to improve earthquake
The precision of composite traces.
In order to achieve the above objectives, on the one hand, the embodiment of the present application provides a kind of method for obtaining synthetic seismic record, packet
Include following steps：
Determine the time-domain primary wave reflection coefficient of each reflecting interface in target work area；And determine each reflection circle
The downgoing reflection control coefrficient in face；Whether to generate multiple wave anti-for controlling each reflecting interface for the downgoing reflection control coefrficient
It penetrates；
According to the downgoing reflection control coefrficient and the time-domain primary wave reflection coefficient, obtain comprising transmission loss and
The second reflection coefficient of time-domain of multiple wave reflection；
It determines seismic wavelet sequence, and the seismic wavelet sequence and second reflection coefficient of time-domain is subjected to pleat
Product obtains synthetic seismic record.
The embodiment of the present application acquisition synthetic seismic record method, the determining each reflecting interface in target work area when
Between domain primary wave reflection coefficient, including：
According to the velocity logging data and density log data in target work area, each reflection circle in the target work area is determined
The wave impedance in face；
The time-domain primary wave reflection coefficient of each reflecting interface is determined according to the wave impedance of each reflecting interface.
The method of the acquisition synthetic seismic record of the embodiment of the present application, the downlink of determination each reflecting interface are anti-
Control coefrficient is penetrated, including：
According to preset rules, the downgoing reflection control coefrficient { c of each reflecting interface is initializedj| j=0,1 ... Nr- 1 },
Wherein, cjIt is the downgoing reflection control coefrficient of j-th of reflecting interface, works as cjWhen=0, indicate synthetic seismic record in do not include with
The related multiple wave reflection of j-th of reflecting interface, works as cjWhen=1, indicate synthetic seismic record in include and j-th of reflecting interface
Related multiple wave reflection.
The method of the acquisition synthetic seismic record of the embodiment of the present application, it is described according to the downgoing reflection control coefrficient and institute
Time-domain primary wave reflection coefficient is stated, the second reflection coefficient of time-domain comprising transmission loss and multiple wave reflection is obtained, including：
According to formula dj+1=-cjrju+(1-rj)djAnd u=rjdj+(1+rj) u each reflection circle of progressive updating from the bottom to top
The downlink wave energy d in facej+1With uplink wave energy u, and using the uplink wave energy u on reflection circle of top layer as including transmission loss
With the second time-domain reflection coefficient of multiple wave reflection
Wherein,For the second time-domain reflection coefficient of k-th of reflecting interface, rjFor the time-domain one of j-th of reflecting interface
Subwave reflection coefficient, cjFor the downgoing reflection control coefrficient of j-th of reflecting interface.
The method of the acquisition synthetic seismic record of the embodiment of the present application, the time-domain primary wave of each reflecting interface are anti-
It penetrates whens coefficient waits and samples.
On the other hand, the embodiment of the present application also provides a kind of devices for obtaining synthetic seismic record, including：
First coefficient determination module, for determining wave reflection system of time-domain of each reflecting interface in target work area
Number；And determine the downgoing reflection control coefrficient of each reflecting interface；The downgoing reflection control coefrficient is each for controlling
Whether reflecting interface generates multiple wave reflection；
Second coefficient determination module, for according to wave reflection system of the downgoing reflection control coefrficient and the time-domain
Number obtains the second reflection coefficient of time-domain comprising transmission loss and multiple wave reflection；
Synthesis module is recorded, for determining seismic wavelet sequence, and by the seismic wavelet sequence and the time-domain the
Two reflection coefficients carry out convolution, obtain synthetic seismic record.
The embodiment of the present application acquisition synthetic seismic record device, the determining each reflecting interface in target work area when
Between domain primary wave reflection coefficient, including：
According to the velocity logging data and density log data in target work area, each reflection circle in the target work area is determined
The wave impedance in face；
The time-domain primary wave reflection coefficient of each reflecting interface is determined according to the wave impedance of each reflecting interface.
The device of the acquisition synthetic seismic record of the embodiment of the present application, the downlink of determination each reflecting interface are anti-
Control coefrficient is penetrated, including：
According to preset rules, the downgoing reflection control coefrficient { c of each reflecting interface is initializedj| j=0,1 ... Nr- 1 },
Wherein, cjIt is the downgoing reflection control coefrficient of j-th of reflecting interface, works as cjWhen=0, indicate synthetic seismic record in do not include with
The related multiple wave reflection of j-th of reflecting interface, works as cjWhen=1, indicate synthetic seismic record in include and j-th of reflecting interface
Related multiple wave reflection.
The device of the acquisition synthetic seismic record of the embodiment of the present application, it is described according to the downgoing reflection control coefrficient and institute
Time-domain primary wave reflection coefficient is stated, the second reflection coefficient of time-domain comprising transmission loss and multiple wave reflection is obtained, including：
According to formula dj+1=-cjrju+(1-rj)djAnd u=rjdj+(1+rj) u each reflection circle of progressive updating from the bottom to top
The downlink wave energy d in facej+1With uplink wave energy u, and using the uplink wave energy u on reflection circle of top layer as including transmission loss
With the second time-domain reflection coefficient of multiple wave reflection
Wherein,For the second time-domain reflection coefficient of k-th of reflecting interface, rjFor the time-domain one of j-th of reflecting interface
Subwave reflection coefficient, cjFor the downgoing reflection control coefrficient of j-th of reflecting interface.
The device of the acquisition synthetic seismic record of the embodiment of the present application, the time-domain primary wave of each reflecting interface are anti-
It penetrates whens coefficient waits and samples.
As can be seen from the technical scheme provided by the above embodiments of the present application, the embodiment of the present application can be according to actual needs by changing
Become downgoing reflection control coefrficient, flexibly can analyze and position the source that multiple wave generates in synthetic seismic record, for complexity
The earthquake multiple reflective analysis on stratum provides a kind of simple, flexible new tool, so that being based on the embodiment of the present application
Synthetic seismic record the obtained synthetic seismic record of method, the earthquake synthesis under the influence of any reflecting interface can be obtained
Record, thus the precision of synthetic seismic record is improved, and can help to identify the multiple wave in real seismic record.
Detailed description of the invention
In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
The some embodiments recorded in application, for those of ordinary skill in the art, in the premise of not making the creative labor property
Under, it is also possible to obtain other drawings based on these drawings.In the accompanying drawings：
Fig. 1 is the method flow diagram that one embodiment of the application obtains synthetic seismic record；
Fig. 2 a~2f be respectively p-wave impedance schematic diagram in one embodiment of the application, primary wave reflection coefficient schematic diagram,
Downgoing reflection control coefrficient schematic diagram, the time-domain primary wave reflection coefficient schematic diagram comprising transmission loss and multiple wave reflection,
Synthetic seismic record based on primary wave reflection coefficient 10 repetition display schematic diagrams transversely, comprising transmission loss and repeatedly
The 10 repetition display schematic diagrams of the synthetic seismic record of the time-domain primary wave reflection coefficient of wave reflection transversely；
Fig. 3 is the seismic wavelet sequence diagram of one embodiment of the application；
Fig. 4 a~4f be respectively in Fig. 2 a~2f corresponding part 1800 milliseconds of downgoing reflection control coefrficients below be 1,
And other downgoing reflection control coefrficients be 0 when composite result enlarged diagram；
Fig. 4 g is the 10 repetition display schematic diagrams of the real seismic record of one embodiment of the application transversely；
Fig. 5 a~5f be respectively in Fig. 2 a~2f corresponding part 0~1230 millisecond downgoing reflection control coefrficient be 1, and
Composite result enlarged diagram when other downgoing reflection control coefrficients are 0；
Fig. 6 a~6f is respectively when 1800 milliseconds of corresponding part downgoing reflection control coefrficients below are 1 in Fig. 5 a~5f
Composite result enlarged diagram；
Fig. 7 a~7f be respectively in Fig. 2 a~2f corresponding part 0~1600 millisecond of downgoing reflection control coefrficient be 1, and its
Composite result enlarged diagram when its downgoing reflection control coefrficient is 0；
Fig. 8 a~8f be respectively in Fig. 7 a~7f corresponding part 1800 milliseconds of downgoing reflection control coefrficients below be 1 when
Composite result enlarged diagram；
Fig. 9 is the apparatus structure block diagram that one embodiment of the application obtains synthetic seismic record.
Specific embodiment
In order to make those skilled in the art better understand the technical solutions in the application, below in conjunction with the application reality
The attached drawing in example is applied, the technical scheme in the embodiment of the application is clearly and completely described, it is clear that described implementation
Example is merely a part but not all of the embodiments of the present application.Based on the embodiment in the application, this field is common
The application protection all should belong in technical staff's every other embodiment obtained without creative efforts
Range.
Refering to what is shown in Fig. 1, the method for the acquisition synthetic seismic record of the embodiment of the present application may comprise steps of：
S101, the time-domain primary wave reflection coefficient for determining each reflecting interface in target work area；And determination is described each
The downgoing reflection control coefrficient of reflecting interface；The downgoing reflection control coefrficient for control each reflecting interface whether generate it is more
Secondary wave reflection.
In the embodiment of the present application, the time-domain primary wave reflection coefficient of each reflecting interface can in the determining target work area
To include：First according to the velocity logging data and density log data in target work area, determine each anti-in the target work area
The wave impedance in firing area face, such as shown in Fig. 2 a；Then each reflecting interface is determined according to the wave impedance of each reflecting interface
Time-domain primary wave reflection coefficient, such as shown in Fig. 2 b.In addition, in the embodiment of the present application, each reflecting interface when
Between domain primary wave reflection coefficient sample whens can be equal.
In the exemplary embodiment of the application, the wave impedance according to each reflecting interface determines each
The time-domain primary wave reflection coefficient of reflecting interface, such as following formula the time-domain of each reflecting interface can be determined according to
Primary wave reflection coefficient：
Wherein, rjFor the time-domain primary wave reflection coefficient of j-th of reflecting interface, ρiviIt is i-th
The wave impedance of reflecting interface, ρiFor the velocity of longitudinal wave of i-th of reflecting interface, viFor the density of i-th of reflecting interface, j=0,
1 ... Nr- 1, i=0,1 ... Nr- 1, NrFor the length of seismic wavelet sequence.
In the application one embodiment, the downgoing reflection control coefrficient for control each reflecting interface whether generate it is more
Secondary wave reflection may include following：
According to preset rules, the downgoing reflection control coefrficient { c of each reflecting interface is initializedj| j=0,1 ... Nr- 1 }, such as
Shown in Fig. 2 c, downgoing reflection control coefrficient is all assigned a value of 0 by initializing.Wherein, cjBe j-th of reflecting interface downlink it is anti-
Control coefrficient is penetrated, the value of the downgoing reflection control coefrficient is 0 or 1.Work as cjWhen=0, indicate not wrap in synthetic seismic record
Containing multiple wave reflection related with j-th of reflecting interface, work as cjWhen=1, indicate in synthetic seismic record comprising being reflected with j-th
The related multiple wave reflection in interface.It can be seen that passing through the downgoing reflection control for changing each reflecting interface (under preset rules)
Coefficient processed flexibly can analyze and position the source that multiple wave generates in synthetic seismic record, to be the ground of bad ground
Shake multiple wave reflective analysis provides a kind of simple, flexible new tool.
S102, according to the downgoing reflection control coefrficient and the time-domain primary wave reflection coefficient, obtain comprising transmission
The second reflection coefficient of time-domain of loss and multiple wave reflection.
It is described according to the wave reflection of the downgoing reflection control coefrficient and the time-domain in the application one embodiment
Coefficient obtains the second reflection coefficient of time-domain comprising transmission loss and multiple wave reflection, including：
Intermediate variable u and sequence { d needed for pre-defined calculating processj| j=0,1 ... Nr, variable u is for saving
With update uplink wave energy, sequence { dj| j=0,1 ... NrFor saving and updating downlink wave energy.
Reflecting interface k each for underground, successively sets k=0,1 ..., N from shallow to deepr- 1, as k=0, setting ground
Table initial downlink wave energy d0=1, as k ≠ 0, set earth's surface initial downlink wave energy d0=0.
Second comprising transmission loss and multiple wave reflection gradually can be sought according to the following steps after completing above-mentioned definition
Time-domain reflection coefficient
Set initial uplink wave energy u=0；For subsurface reflective boundary j, j=k, k-1 ... are successively shallowly set by being deep to,
0, according to formula dj+1=-cjrju+(1-rj)djAnd u=rjdj+(1+rj) downlink wave energy on each reflecting interface of u progressive updating
Measure dj+1With uplink wave energy u, and finally using j=0 when uplink wave energy u as containing transmission loss and multiple wave reflection
The second time-domain reflection coefficient
The second reflection coefficient of time-domain comprising transmission loss and multiple wave reflection is for example as shown in Figure 2 d.
S103, determine seismic wavelet sequence, and by the seismic wavelet sequence and second reflection coefficient of time-domain into
Row convolution obtains synthetic seismic record.
In the application one embodiment, seismic wavelet sequence { w for example shown in Fig. 3 can be chosen according to actual needsk|k
=0,1 ... Nw- 1 }, seismic wavelet sequence and the second reflection coefficient of time-domain are then subjected to convolution, to generate earthquake synthesis
Record { sj| j=0,1 ... Nr, such as Fig. 2 f -1 }.Obviously, with shown in Fig. 2 e directly by seismic wavelet sequence and time-domain one
The synthetic seismic record that subwave reflection coefficient obtains is compared, and the precision of Fig. 2 f is higher.
In one exemplary embodiment of the application, the convolution operation can for example pass through following formulaIt realizes, wherein sjFor synthetic seismic record,For the time-domain comprising transmission loss and multiple wave reflection
Second reflection coefficient,For seismic wavelet sequence, wherein MwThe center of seismic wavelet sequence, j=0,1 ... Nr- 1, and
And work as j-k+Mw< 0 or j-k+Mw> NwWhen -1, take
In the application one embodiment, refer to Fig. 4 a~4f, wherein comparison diagram 4e, Fig. 4 f and Fig. 4 g it can be found that
In target zone 1 and target zone 2, Fig. 4 f and Fig. 4 g is increasingly similar, shows to contain multiple wave reflection in real seismic record.
In order to further analyze multiple wave Producing reason, as shown in Fig. 5 a~5f, only changing above-described embodiment step
Downgoing reflection control coefrficient (i.e. { c in S101j=0 | j=0~1230 }, { cj=1 | j=1231~Nr- 1 } after), Fig. 5 f phase
Than being changed in Fig. 2 f；Especially in conjunction with shown in Fig. 6 a~6f, it can be seen that in target zone 1, Fig. 6 e and Fig. 6 f are
It is more similar, show that downgoing reflection caused by 1230 milliseconds or more of reflecting interface causes larger impact to target zone 1.
In order to further analyze multiple wave Producing reason, as shown in Fig. 7 a~7f, only changing above-described embodiment step
Downgoing reflection control coefrficient (i.e. { c in S101j=0 | j=0~1600 }, { cj=1 | j=1601~Nr- 1 } after), Fig. 7 f phase
Than being changed in Fig. 2 f；Especially in conjunction with shown in Fig. 8 a~8f, it can be seen that in target zone 2, the phase of Fig. 8 e and Fig. 8 f
Similarity degree like degree compared with Fig. 6 e and Fig. 6 f increased, and but without similar enough, show 1600 milliseconds or more of reflection
Downgoing reflection caused by interface causes larger impact to target zone 2, moreover, also showing 1600 milliseconds of reflecting interfaces below
Generated downgoing reflection still can impact target zone 2.
In conclusion the embodiment of the present application can be according to actual needs by changing downgoing reflection control coefrficient, it can be flexible
Analysis and positioning synthetic seismic record in multiple wave generate source, provided for the earthquake multiple reflective analysis of bad ground
A kind of simple, flexible new tool, so that the method for the synthetic seismic record based on the embodiment of the present application is obtained
Synthetic seismic record can obtain the synthetic seismic record under the influence of any reflecting interface, thus improve synthetic seismic record
Precision.
Although procedures described above process includes the multiple operations occurred with particular order, it should however be appreciated that understand,
These processes may include more or fewer operations, these operations can be executed sequentially or be executed parallel (such as using parallel
Processor or multi-thread environment).
Refering to what is shown in Fig. 9, the device of the acquisition synthetic seismic record of the embodiment of the present application, may include：
First coefficient determination module 91, the time-domain primary wave for being determined for each reflecting interface in target work area are anti-
Penetrate coefficient；And determine the downgoing reflection control coefrficient of each reflecting interface；The downgoing reflection control coefrficient is for controlling
Whether each reflecting interface generates multiple wave reflection；
Second coefficient determination module 92 can be used for according to the downgoing reflection control coefrficient and the time-domain primary wave
Reflection coefficient obtains the second reflection coefficient of time-domain comprising transmission loss and multiple wave reflection；
Record synthesis module 93, be determined for seismic wavelet sequence, and by the seismic wavelet sequence and it is described when
Between the second reflection coefficient of domain carry out convolution, obtain synthetic seismic record.
The device of the acquisition synthetic seismic record of the embodiment of the present application and above-mentioned acquisition synthetic seismic record shown in FIG. 1
Method is corresponding, therefore, about the details of the device of the acquisition synthetic seismic record of the embodiment of the present application, refers to above-mentioned Fig. 1
Shown in obtain synthetic seismic record method, details are not described herein.
For convenience of description, it is divided into various units when description apparatus above with function to describe respectively.Certainly, implementing this
The function of each unit can be realized in the same or multiple software and or hardware when application.
It should be understood by those skilled in the art that, the embodiment of the present invention can provide as method, system or computer program
Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the present invention
Apply the form of example.Moreover, it wherein includes the computer of computer usable program code that the present invention, which can be used in one or more,
The computer program implemented in usable storage medium (including but not limited to magnetic disk storage, CD~ROM, optical memory etc.) produces
The form of product.
The present invention be referring to according to the method for the embodiment of the present invention, the process of equipment (system) and computer program product
Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions
The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs
Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real
The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
In a typical configuration, calculating equipment includes one or more processors (CPU), input/output interface, net
Network interface and memory.
Memory may include the non-volatile memory in computer-readable medium, random access memory (RAM) and/or
The forms such as Nonvolatile memory, such as read-only memory (ROM) or flash memory (flash RAM).Memory is computer-readable medium
Example.
Computer-readable medium includes permanent and non-permanent, removable and non-removable media can be by any method
Or technology come realize information store.Information can be computer readable instructions, data structure, the module of program or other data.
The example of the storage medium of computer includes, but are not limited to phase change memory (PRAM), static random access memory (SRAM), moves
State random access memory (DRAM), other kinds of random access memory (RAM), read-only memory (ROM), electric erasable
Programmable read only memory (EEPROM), flash memory or other memory techniques, read-only disc read only memory (CD-ROM) (CD~
ROM), digital versatile disc (DVD) or other optical storage, magnetic cassettes, tape magnetic disk storage or other magnetism are deposited
Equipment or any other non-transmission medium are stored up, can be used for storage can be accessed by a computing device information.According to boundary herein
Fixed, computer-readable medium does not include temporary computer readable media (transitory media), such as the data-signal of modulation and
Carrier wave.
It should also be noted that, the terms "include", "comprise" or its any other variant are intended to nonexcludability
It include so that the process, method, commodity or the equipment that include a series of elements not only include those elements, but also to wrap
Include other elements that are not explicitly listed, or further include for this process, method, commodity or equipment intrinsic want
Element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that including described want
There is also other identical elements in the process, method of element, commodity or equipment.
It will be understood by those skilled in the art that embodiments herein can provide as method, system or computer program product.
Therefore, complete hardware embodiment, complete software embodiment or embodiment combining software and hardware aspects can be used in the application
Form.It is deposited moreover, the application can be used to can be used in the computer that one or more wherein includes computer usable program code
The shape for the computer program product implemented on storage media (including but not limited to magnetic disk storage, CD~ROM, optical memory etc.)
Formula.
The application can describe in the general context of computer-executable instructions executed by a computer, such as program
Module.Generally, program module includes routines performing specific tasks or implementing specific abstract data types, programs, objects, group
Part, data structure etc..The application can also be practiced in a distributed computing environment, in these distributed computing environments, by
Task is executed by the connected remote processing devices of communication network.In a distributed computing environment, program module can be with
In the local and remote computer storage media including storage equipment.
All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment
Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for system reality
For applying example, since it is substantially similar to the method embodiment, so being described relatively simple, related place is referring to embodiment of the method
Part explanation.
The above description is only an example of the present application, is not intended to limit this application.For those skilled in the art
For, various changes and changes are possible in this application.All any modifications made within the spirit and principles of the present application are equal
Replacement, improvement etc., should be included within the scope of the claims of this application.
Claims (8)
1. a kind of method for obtaining synthetic seismic record, which is characterized in that include the following steps：
Determine the time-domain primary wave reflection coefficient of each reflecting interface in target work area；And determine each reflecting interface
Downgoing reflection control coefrficient；The downgoing reflection control coefrficient is for controlling whether each reflecting interface generates multiple wave reflection；
According to the downgoing reflection control coefrficient and the time-domain primary wave reflection coefficient, acquisition is comprising transmission loss and repeatedly
The second reflection coefficient of time-domain of wave reflection；
It determines seismic wavelet sequence, and the seismic wavelet sequence and second reflection coefficient of time-domain is subjected to convolution, obtain
Obtain synthetic seismic record；Wherein,
It is described according to the downgoing reflection control coefrficient and the time-domain primary wave reflection coefficient, obtain comprising transmission loss and
The second reflection coefficient of time-domain of multiple wave reflection, including：
According to formula dj+1=-cjrju+(1-rj)djAnd u=rjdj+(1+rj) u progressive updating each reflecting interface from the bottom to top
Downlink wave energy dj+1With uplink wave energy u, and using top layer reflection circle uplink wave energy u be used as comprising transmission loss with it is more
Second time-domain reflection coefficient of secondary wave reflection
Wherein,For the second time-domain reflection coefficient of k-th of reflecting interface, rjFor the time-domain primary wave of j-th of reflecting interface
Reflection coefficient, cjFor the downgoing reflection control coefrficient of j-th of reflecting interface.
2. the method according to claim 1 for obtaining synthetic seismic record, which is characterized in that the determining target work area is each
The time-domain primary wave reflection coefficient of a reflecting interface, including：
According to the velocity logging data and density log data in target work area, each reflecting interface in the target work area is determined
Wave impedance；
The time-domain primary wave reflection coefficient of each reflecting interface is determined according to the wave impedance of each reflecting interface.
3. the method according to claim 1 for obtaining synthetic seismic record, which is characterized in that the determination is described each anti-
The downgoing reflection control coefrficient in firing area face, including：
According to preset rules, the downgoing reflection control coefrficient { c of each reflecting interface is initializedj| j=0,1 ... Nr- 1 }, wherein
cjFor the downgoing reflection control coefrficient of j-th of reflecting interface, work as cjWhen=0, indicate not including and j-th in synthetic seismic record
The related multiple wave reflection of reflecting interface, works as cjWhen=1, indicate in synthetic seismic record comprising related with j-th of reflecting interface
Multiple wave reflection.
4. the method according to claim 1 for obtaining synthetic seismic record, which is characterized in that each reflecting interface
Time-domain primary wave reflection coefficient samples whens equal.
5. a kind of device for obtaining synthetic seismic record, which is characterized in that including：
First coefficient determination module, for determining the time-domain primary wave reflection coefficient of each reflecting interface in target work area；And
Determine the downgoing reflection control coefrficient of each reflecting interface；The downgoing reflection control coefrficient is for controlling each reflection circle
Whether face generates multiple wave reflection；
Second coefficient determination module is used for according to the downgoing reflection control coefrficient and the time-domain primary wave reflection coefficient,
Obtain the second reflection coefficient of time-domain comprising transmission loss and multiple wave reflection；
Synthesis module is recorded, for determining seismic wavelet sequence, and the seismic wavelet sequence and the time-domain second is anti-
It penetrates coefficient and carries out convolution, obtain synthetic seismic record；Wherein,
It is described according to the downgoing reflection control coefrficient and the time-domain primary wave reflection coefficient, obtain comprising transmission loss and
The second reflection coefficient of time-domain of multiple wave reflection, including：
According to formula dj+1=-cjrju+(1-rj)djAnd u=rjdj+(1+rj) u progressive updating each reflecting interface from the bottom to top
Downlink wave energy dj+1With uplink wave energy u, and using top layer reflection circle uplink wave energy u be used as comprising transmission loss with it is more
Second time-domain reflection coefficient of secondary wave reflection
Wherein,For the second time-domain reflection coefficient of k-th of reflecting interface, rjFor the time-domain primary wave of j-th of reflecting interface
Reflection coefficient, cjFor the downgoing reflection control coefrficient of j-th of reflecting interface.
6. the device according to claim 5 for obtaining synthetic seismic record, which is characterized in that the determining target work area is each
The time-domain primary wave reflection coefficient of a reflecting interface, including：
According to the velocity logging data and density log data in target work area, each reflecting interface in the target work area is determined
Wave impedance；
The time-domain primary wave reflection coefficient of each reflecting interface is determined according to the wave impedance of each reflecting interface.
7. the device according to claim 5 for obtaining synthetic seismic record, which is characterized in that the determination is described each anti-
The downgoing reflection control coefrficient in firing area face, including：
According to preset rules, the downgoing reflection control coefrficient { c of each reflecting interface is initializedj| j=0,1 ... Nr- 1 }, wherein
cjIt is the downgoing reflection control coefrficient of j-th of reflecting interface, works as cjWhen=0, indicate not including and j-th in synthetic seismic record
The related multiple wave reflection of reflecting interface, works as cjWhen=1, indicate in synthetic seismic record comprising related with j-th of reflecting interface
Multiple wave reflection.
8. the device according to claim 5 for obtaining synthetic seismic record, which is characterized in that each reflecting interface
Time-domain primary wave reflection coefficient samples whens equal.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710213246.0A CN106970418B (en) | 2017-04-01 | 2017-04-01 | Method and device for acquiring earthquake synthetic record |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710213246.0A CN106970418B (en) | 2017-04-01 | 2017-04-01 | Method and device for acquiring earthquake synthetic record |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106970418A CN106970418A (en) | 2017-07-21 |
| CN106970418B true CN106970418B (en) | 2018-11-16 |
Family
ID=59336693
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710213246.0A Active CN106970418B (en) | 2017-04-01 | 2017-04-01 | Method and device for acquiring earthquake synthetic record |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106970418B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109143335B (en) * | 2018-07-27 | 2020-04-14 | 中国地质调查局沈阳地质调查中心 | Method, system, medium and apparatus for producing synthetic seismic record |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7336561B2 (en) * | 2004-09-07 | 2008-02-26 | Pgs Americas, Inc. | System for attenuation of water bottom multiples in seismic data recorded by pressure sensors and particle motion sensors |
| CN102169189B (en) * | 2011-01-19 | 2015-01-14 | 中国海洋石油总公司 | Method for eliminating deep water interlamination multiple |
| US9513390B2 (en) * | 2014-11-03 | 2016-12-06 | Westerngeco L.L.C. | Compensating for space and slowness/angle blurring of reflectivity |
| CN105510976B (en) * | 2016-01-11 | 2016-11-23 | 中国海洋大学 | A kind of many subwaves combination adaptive attenuation method |
-
2017
- 2017-04-01 CN CN201710213246.0A patent/CN106970418B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106970418A (en) | 2017-07-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA3122986C (en) | Automated seismic interpretation-guided inversion | |
| CN106896406B (en) | Method and device for predicting high-quality reservoir stratum based on impedance domain attributes | |
| AU722750B2 (en) | Presentation and interpretation of seismic data | |
| CN108181654B (en) | AVAF analogy method and device based on multi-scale rock physical model | |
| EP2539844B1 (en) | Method and system for using multiple-point statistics simulation to model reservoir property trends | |
| CN108535775B (en) | Non-stationary seismic data sound impedance inversion method and device | |
| CN104375188A (en) | Seismic wave transmission attenuation compensation method and device | |
| CN113919219A (en) | Stratum evaluation method and system based on logging big data | |
| WO2021194524A1 (en) | Method and system for automated velocity model updating using machine learning | |
| Grechka | Shear-wave group-velocity surfaces in low-symmetry anisotropic media | |
| CN105425289A (en) | Method and device of determining low frequency wave impedance | |
| CN109164491A (en) | A kind of seismic facies recognition methods and system based on category support vector machines | |
| CN109633742A (en) | A kind of full waveform inversion method and device | |
| CN108445532B (en) | A kind of Depth Domain inverse migration method and device | |
| US20210311223A1 (en) | System and method for seismic inversion | |
| CN106970418B (en) | Method and device for acquiring earthquake synthetic record | |
| CN104459774B (en) | Geological lithology difference identification method and system | |
| US10454713B2 (en) | Domain decomposition using a multi-dimensional spacepartitioning tree | |
| CN106772599A (en) | A kind of method and device for calculating formation shear speed | |
| CN107831541B (en) | Thin strate recognition methods and device based on high density VSP data | |
| CN106125133A (en) | A kind of based on the fine velocity modeling method under the constraint of gas cloud district | |
| US20090043549A1 (en) | Methods, apparatus, and products for seismic ray tracing | |
| CN114114421A (en) | Deep learning-based guided self-learning seismic data denoising method and device | |
| Kazei et al. | Acquisition and near-surface impacts on VSP mini-batch FWI and RTM imaging in desert environment | |
| US10620341B2 (en) | System and method for modifying an earth model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |