CN108873067A - Diffraction coefficient method for solving and device - Google Patents
Diffraction coefficient method for solving and device Download PDFInfo
- Publication number
- CN108873067A CN108873067A CN201811128027.3A CN201811128027A CN108873067A CN 108873067 A CN108873067 A CN 108873067A CN 201811128027 A CN201811128027 A CN 201811128027A CN 108873067 A CN108873067 A CN 108873067A
- Authority
- CN
- China
- Prior art keywords
- gal
- diffraction coefficient
- diffracted wave
- domain
- primary domain
- 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.)
- Pending
Links
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. analysis, for interpretation, for correction
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)
- Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)
Abstract
The present invention provides a kind of diffraction coefficient method for solving and devices, belong to seismic prospecting and engineering geology technical field.Diffraction coefficient method for solving and device provided in an embodiment of the present invention obtain the diffracted wave under attenuation function effect;Time histories sample is combined with attenuation function, obtains the attenuation function of added-time window function;Gal uncle's transformation is carried out to the diffracted wave under attenuation function effect, obtains the diffracted wave in the primary domain of gal;According to the attenuation function of diffracted wave and added-time window function in the primary domain of gal, obtain the diffraction coefficient in the primary domain of gal, diffraction coefficient in the primary domain of gal is converted into the diffraction coefficient in time-domain, obtain the diffraction coefficient of subsurface structure geologic body in time-domain, diffraction coefficient is solved according to the primary shift theory of gal, and improved time histories sample is introduced into attenuation function, available high-precision diffraction coefficient, to identify that little structure geologic body provides more reliable foundation in seismic prospecting, effective ensure is provided to the safe working of mineral resources.
Description
Technical field
The present invention relates to geological prospecting and engineering geology technical fields, solve in particular to a kind of diffraction coefficient
Method and apparatus.
Background technique
In the exploitation of the mineral resources such as coal, petroleum, it is different usually to encounter the construction such as tomography, karst collapse col umn, structural trap
Conduit pipe in coal seam is often connected in Chang Ti, these anomalous bodys, causes Seam Roof And Floor water inrush accident, to the safety open of coal mine
It adopts and brings danger with personnel safety, if structural trap is generated near Faults, destructible trap causes the loss of oil gas, no
Conducive to the storage and exploitation of petroleum.Therefore, find out that little structure anomalous body has the exploitation of the mineral resources such as coal and petroleum
Highly important meaning.In seismic prospecting, when seismic wave is propagated to underground medium, these little structure geologic bodies are encountered often
It will form diffracted wave, little structure geologic body can be identified using diffracted wave.Can be obtained by diffracted wave numerical simulation principle, diffracted wave by
Diffraction coefficient is formed with corresponding attenuation function two parts.Attenuation function is determined by the attenuation by absorption effect on stratum, can only be influenced
Diffraction wave energy power situation of change, does not influence the variation characteristic of diffracted wave.Therefore, diffraction coefficient is to determine that diffracted wave passes
Broadcast the important symbol of feature and its Spatial Variation.
In the exploitation and seismic prospecting of the mineral resources such as coal, petroleum, effective mesh is usually identified with reflective wave method
Layer, and ignore carry little structure geologic body information diffracted wave.The method that routine seeks diffracted wave is directly in seismic wave
Separation back wave obtains, and the diffracted wave sought in this way is mingled with very more interference waves, is unfavorable for studying the generation of little structure geologic body
Diffracted wave variation characteristic, also, diffracted wave can be encountered in communication process underground medium attenuation by absorption effect, if not
By the attenuation by absorption process of diffracted wave in the earth formation in view of the diffraction coefficient in the calculating of diffraction coefficient, obtained can be not smart enough
Really, it cannot really reflect the Variation Features of diffracted wave, so that little structure geologic body cannot be identified effectively, and cannot be surveyed for earthquake
It finds and explains that little structure geologic body provides more reliable foundation in spy, and the safe working of mineral resources is provided and is effectively protected
Barrier.
Summary of the invention
For above-mentioned problems of the prior art, the present invention provides a kind of diffraction coefficient method for solving and device,
Available high-precision diffraction coefficient, so that the discrimination of little structure geologic body is improved, to identify little structure in seismic prospecting
Geologic body provides more reliable foundation, provides effective ensure to the safe working of mineral resources.
In a first aspect, the embodiment of the invention provides a kind of diffraction coefficient method for solving, the method includes:
Obtain the diffracted wave under attenuation function effect;
Time histories sample is combined with the attenuation function, obtains the attenuation function of added-time window function;
Gal uncle's transformation is carried out to the diffracted wave under attenuation function effect, obtains the diffracted wave in the primary domain of gal;
According to the attenuation function of diffracted wave and the added-time window function in the primary domain of the gal, the diffraction in the primary domain of gal is obtained
Coefficient;
Diffraction coefficient in the primary domain of the gal is converted into the diffraction coefficient in time-domain, obtains subsurface structure in time-domain
The diffraction coefficient of geologic body.
With reference to first aspect, the embodiment of the invention provides the first possible embodiments of first aspect, wherein institute
Stating time histories sample is the time histories sample for adding modifying factor.
With reference to first aspect, the embodiment of the invention provides second of possible embodiments of first aspect, wherein root
According to the attenuation function of diffracted wave and the added-time window function in the primary domain of the gal, the step of the diffraction coefficient in the primary domain of gal is obtained
Suddenly, including:
Using the diffracted wave in the primary domain of the gal as molecule, using the attenuation function of the added-time window function as denominator, meter
Calculation obtains the diffraction coefficient in the primary domain of the gal.
With reference to first aspect, the embodiment of the invention provides the third possible embodiments of first aspect, wherein obtains
Before the step of taking the diffracted wave under attenuation function effect, the method also includes:
Noise and random disturbances processing are carried out to the data of acquisition, the data that obtain that treated;
Treated that data are separated by described, obtains diffracted wave.
Second aspect, the embodiment of the invention also provides a kind of diffracted wave solving device, described device includes:
First acquisition unit, for obtaining the diffracted wave under attenuation function effect;
Second acquisition unit obtains the decaying of added-time window function for combining time histories sample with the attenuation function
Function;
Converter unit obtains the diffraction in the primary domain of gal for carrying out gal uncle's transformation to the diffracted wave under the attenuation
Wave;
First diffraction coefficient acquiring unit, after according to the diffracted wave and the added-time window function in the primary domain of the gal
Attenuation function obtains the diffraction coefficient in the primary domain of gal;
Second diffraction coefficient acquiring unit, for the diffraction coefficient in the primary domain of the gal to be converted to the diffraction in time-domain
Coefficient obtains the diffraction coefficient of subsurface structure geologic body in time-domain.
In conjunction with second aspect, the embodiment of the invention provides the first possible embodiments of second aspect, wherein institute
Stating time histories sample is the time histories sample for adding modifying factor.
In conjunction with second aspect, the embodiment of the invention provides second of possible embodiments of second aspect, wherein institute
The first diffraction coefficient acquiring unit is stated, is also used to:
Using the diffracted wave in the primary domain of the gal as molecule, using the attenuation function of the added-time window function as denominator, meter
Calculation obtains the diffraction coefficient in the primary domain of the gal.
In conjunction with second aspect, the embodiment of the invention provides the third possible embodiments of second aspect, wherein institute
Stating device further includes:
Diffracted wave acquiring unit, for carrying out noise to the data of acquisition and random disturbances are handled, the number that obtains that treated
According to;Treated that data are separated by described, obtains diffracted wave.
The third aspect, the embodiment of the invention also provides a kind of diffraction coefficients to solve equipment, including processor and memory,
The machine-executable instruction that can be executed by the processor is stored in the memory, the processor executes the machine
Executable instruction is to realize method provided by first aspect.
Fourth aspect, it is described computer-readable to deposit the embodiment of the invention also provides a kind of computer readable storage medium
Storage media is stored with machine-executable instruction, and the machine-executable instruction is when being called and being executed by processor, the machine
Executable instruction promotes the processor to realize method provided by first aspect.
The embodiment of the present invention brings following beneficial effect:
The present invention provides a kind of diffraction coefficient method for solving and devices, belong to seismic prospecting and engineering geology technology neck
Domain.Diffraction coefficient method for solving and device provided in an embodiment of the present invention, obtain attenuation function effect under diffracted wave, by when window
Function is combined with attenuation function, obtains the attenuation function of added-time window function, is carried out to the diffracted wave under attenuation function effect
Jia Bai transformation, obtains the diffracted wave in the primary domain of gal, according to the diffraction wave attenuation after the diffracted wave and added-time window function in the primary domain of gal
Function obtains the diffraction coefficient in the primary domain of gal, and the diffraction coefficient in the primary domain of gal is converted to the diffraction coefficient in time-domain, is obtained
The diffraction coefficient of subsurface structure geologic body in time-domain.Diffraction coefficient is solved according to the primary shift theory of gal, and attenuation function is drawn
Enter improved time histories sample, available high-precision diffraction coefficient, to identify that little structure geologic body provides in seismic prospecting
More reliable foundation provides effective ensure to the safe working of mineral resources.
Other features and advantages of the present invention will illustrate in the following description, also, partly become from specification
It obtains it is clear that understand through the implementation of the invention.The objectives and other advantages of the invention are in specification, claims
And specifically noted structure is achieved and obtained in attached drawing.
To enable the above objects, features and advantages of the present invention to be clearer and more comprehensible, preferred embodiment is cited below particularly, and cooperate
Appended attached drawing, is described in detail below.
Detailed description of the invention
It, below will be to specific in order to illustrate more clearly of the specific embodiment of the invention or technical solution in the prior art
Embodiment or attached drawing needed to be used in the description of the prior art be briefly described, it should be apparent that, it is described below
Attached drawing is some embodiments of the present invention, for those of ordinary skill in the art, before not making the creative labor
It puts, is also possible to obtain other drawings based on these drawings.
Fig. 1 is the flow chart of diffraction coefficient method for solving provided by the embodiment of the present invention;
Fig. 2 is the structural block diagram of diffraction coefficient solving device provided by the embodiment of the present invention;
Fig. 3 is the structural block diagram that diffraction coefficient provided by the embodiment of the present invention solves equipment.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with attached drawing to the present invention
Technical solution be clearly and completely described, it is clear that described embodiments are some of the embodiments of the present invention, rather than
Whole embodiments.The component of embodiments of the present invention, which are generally described and illustrated herein in the accompanying drawings can be matched with a variety of different
It sets to arrange and design.Therefore, the detailed description of the embodiment of the present invention provided in the accompanying drawings is not intended to limit below
The range of claimed invention, but it is merely representative of selected embodiment of the invention.Based on the embodiments of the present invention, originally
Field those of ordinary skill every other embodiment obtained without making creative work, belongs to the present invention
The range of protection.
The change of the diffracted wave of little structure geologic body generation can not be effectively studied for the existing method for seeking diffracted wave
Change Characteristic Problem, the embodiment of the invention provides a kind of diffraction coefficient method for solving and devices, below first to it is of the invention around
Coefficient method for solving is penetrated to describe in detail.
Embodiment one
A kind of diffraction coefficient method for solving is present embodiments provided, as shown in Figure 1, this method comprises the following steps:
Step S101 obtains the diffracted wave under attenuation function effect.
Specifically, according to the primary shift theory of gal, the diffracted wave under attenuation function effect is calculated.Under normal Q attenuation function
The diffracted wave expression of unstable state be:
Wherein, sQIt (t) is diffraction wave function, αQ(τ, f)=e-πfτ/Q+iH(πfτ/Q)For the normal Q attenuation function of underground medium, d (τ)
For diffraction coefficient.
Time histories sample is combined with attenuation function, obtains the attenuation function of added-time window function by step S102.
In the calculating of attenuation function introduce addition modifying factor time histories sample, can make diffracted wave calculating process with
The propagation time of diffracted wave does different decaying, obtains more accurate diffracted wave.Declining for the time histories sample of modifying factor is added
The expression of subtraction function is as follows:
αQ(τ-τ ', f)=e-πf(τ-τ')/Q+iH(πf(τ-τ')/Q)
Window when wherein τ is, τ ' are modifying factor,
Step S103 carries out gal uncle's transformation to the diffracted wave under attenuation function effect, obtains the diffracted wave in the primary domain of gal.
Lower diffracted wave is acted on to attenuation function and carries out gal uncle's transformation, is transformed into the primary domain of gal and is solved, obtain the primary domain of gal
In diffracted wave.The expression of diffracted wave in the primary domain of gal is as follows:
Vgs(τ-τ',f)≈W(f)αQ(τ-τ',f)Vgd(τ-τ',f)
In formula, VgS (τ-τ ', f) is the diffracted wave expression formula in the primary domain of gal, and W (f) is the Fourier of given seismic wavelet
Transformation, VgD (τ-τ ', f) is the diffraction coefficient in the gal uncle transformation i.e. primary domain of gal of diffraction coefficient.
Step S104, according to the attenuation function of diffracted wave and added-time window function in the primary domain of gal, obtain in the primary domain of gal around
Penetrate coefficient.
According to the attenuation function of diffracted wave, added-time window function in step S103 in the primary domain of gal, given seismic wavelet
The relationship of Fourier transformation and the diffraction coefficient in the primary domain of gal, obtaining diffraction coefficient expression in the primary domain of gal is:
Diffraction coefficient in the primary domain of gal is converted to the diffraction coefficient in time-domain by step S105, with obtaining in time-domain
The diffraction coefficient of lower tectonic geology body.
Diffraction coefficient in the primary domain of gal is converted into the diffraction coefficient in time-domain to get the subsurface structure into time-domain
The diffraction coefficient of plastid, and the diffraction coefficient of subsurface structure geologic body is in time-domain convenient for imaging.Diffraction coefficient in time-domain
Expression is as follows:
D (t)=G-1[Vgd(τ-τ'),f]
Wherein, G-1For the primary inverse transformation of gal.
In practical applications, before step S101, it should noise first is carried out to the data of acquisition and random disturbances are handled,
By treated, data are separated, and obtain diffracted wave, then execute step S101 again, calculate the diffraction under attenuation function effect
Wave.Wherein, random disturbances are generally the industrial noise of 50Hz, and noise refers generally to surface wave, sound wave or environmental factor and artificial
It is influenced caused by factor, the frequency of surface wave is lower, and generally between 10-30Hz, the frequency of sound wave is relatively high, generally in 100Hz
More than, these noises and interference different from seismic wave or diffraction wave frequency rate can be removed with a filter.
Diffraction coefficient method for solving provided in an embodiment of the present invention obtains the diffracted wave under attenuation function effect, will improve
Time histories sample afterwards is combined with attenuation function, obtains the attenuation function of added-time window function, to the diffraction under attenuation function effect
Wave, which is done, carries out gal uncle's transformation, obtains the diffracted wave in the primary domain of gal, according to after the diffracted wave and added-time window function in the primary domain of gal around
Ejected wave attenuation function obtains the diffraction coefficient in the primary domain of gal, and the diffraction coefficient in the primary domain of gal is converted to the diffraction in time-domain
Coefficient obtains the diffraction coefficient of subsurface structure geologic body in time-domain.Diffraction coefficient is solved according to the primary shift theory of gal, and will be changed
Time histories sample after is introduced into attenuation function, available high-precision diffraction coefficient, to identify little structure in seismic prospecting
Geologic body provides more reliable foundation, provides effective ensure to the safe working of mineral resources.
Embodiment two
A kind of diffraction coefficient solving device is present embodiments provided, as shown in Fig. 2, the device includes:
First acquisition unit 201, for obtaining the diffracted wave under attenuation function effect.
Second acquisition unit 202 obtains the decaying letter of added-time window function for combining time histories sample with attenuation function
Number.
Wherein, time histories sample be add modifying factor time histories sample, add the time histories sample of modifying factor, can make around
Ejected wave does different decaying with the propagation time of diffracted wave in calculating process, makes being accurately calculated for attenuation function, and make
Final diffraction coefficient is accurately calculated.
Converter unit 203 obtains the diffraction in the primary domain of gal for carrying out gal uncle's transformation to the diffracted wave under attenuation
Wave.
First diffraction coefficient acquiring unit 204, for according to the decaying after the diffracted wave and added-time window function in the primary domain of gal
Function obtains the diffraction coefficient in the primary domain of gal.
Specifically, it is calculated using the diffracted wave in the primary domain of gal as molecule using the attenuation function of added-time window function as denominator
Obtain the diffraction coefficient in the primary domain of gal.
Second diffraction coefficient acquiring unit 205 obtains subterranean in time-domain according to the diffraction coefficient in the primary domain of the gal
Make the diffraction coefficient of geologic body.
Optionally, which can also include:Diffracted wave acquiring unit, for the data progress noise to acquisition and at random
Interference processing, the data that obtain that treated;By treated, data are separated, and obtain diffracted wave.
Diffraction coefficient solving device provided in an embodiment of the present invention, obtain attenuation function effect under diffracted wave, by when window
Function is combined with attenuation function, obtains the attenuation function of added-time window function, is carried out to the diffracted wave under attenuation function effect
Jia Bai transformation, obtains the diffracted wave in the primary domain of gal, according to the diffraction wave attenuation after the diffracted wave and added-time window function in the primary domain of gal
Function obtains the diffraction coefficient in the primary domain of gal, and the diffraction coefficient in the primary domain of gal is converted to the diffraction coefficient in time-domain, is obtained
The diffraction coefficient of subsurface structure geologic body in time-domain.Diffraction coefficient is solved according to the primary shift theory of gal, and in attenuation function
Improved time histories sample is introduced, available high-precision diffraction coefficient mentions for identification little structure geologic body in seismic prospecting
For more reliable foundation, effective ensure is provided to the safe working of mineral resources.
Embodiment three
A kind of diffraction coefficient solution equipment is present embodiments provided, as shown in figure 3, the equipment includes processor 31 and storage
Device 32;Wherein, memory 32 is held for storing one or more computer instruction, one or more computer instruction by processor
Row, to realize above-mentioned distributed task scheduling processing method.Equipment shown in Fig. 3 further includes bus 33, processor 31 and memory 32
It is connected by bus 33.
Wherein, memory 32 may include high-speed random access memory (RAM, Random Access Memory),
It may further include non-labile memory (non-volatile memory), for example, at least a magnetic disk storage.Bus 33
It can be isa bus, pci bus or eisa bus etc..The bus can be divided into address bus, data/address bus, control bus
Deng.Only to be indicated with a four-headed arrow in Fig. 3 convenient for indicating, it is not intended that an only bus or a type of total
Line.
Processor 31 may be a kind of IC chip, the processing capacity with signal.During realization, above-mentioned side
Each step of method can be completed by the integrated logic circuit of the hardware in processor 31 or the instruction of software form.Above-mentioned
Processor 31 can be general processor, including central processing unit (Central Processing Unit, abbreviation CPU), network
Processor (Network Processor, abbreviation NP) etc.;It can also be digital signal processor (Digital Signal
Processing, abbreviation DSP), specific integrated circuit (Application Specific Integrated Circuit, referred to as
ASIC), ready-made programmable gate array (Field-Programmable Gate Array, abbreviation FPGA) or other are programmable
Logical device, discrete gate or transistor logic, discrete hardware components.It may be implemented or execute embodiment of the present invention
In disclosed each method, step and logic diagram.General processor can be microprocessor or the processor is also possible to
Any conventional processor etc..The step of method in conjunction with disclosed in embodiment of the present invention, can be embodied directly in hardware decoding
Processor executes completion, or in decoding processor hardware and software module combination execute completion.Software module can position
In random access memory, flash memory, read-only memory, programmable read only memory or electrically erasable programmable memory, register
In the storage medium of equal this fields maturation.The storage medium is located at memory 32, and processor 31 reads the information in memory 32,
The step of completing the method for aforementioned embodiments in conjunction with its hardware.
Further, the embodiment of the invention also provides a kind of machine readable storage medium, the machine readable storage mediums
It is stored with machine-executable instruction, when being called and being executed by processor, machine-executable instruction promotees the machine-executable instruction
Processor is set to realize above-mentioned diffraction coefficient method for solving.
Diffraction coefficient method for solving and device provided in an embodiment of the present invention technical characteristic having the same, so can also solve
Certainly identical technical problem reaches identical technical effect.
It should be noted that in embodiment provided by the present invention, it should be understood that disclosed system and method, it can
To realize by another way.The apparatus embodiments described above are merely exemplary, for example, the unit is drawn
Point, only a kind of logical function partition, there may be another division manner in actual implementation, in another example, multiple units or group
Part can be combined or can be integrated into another system, or some features can be ignored or not executed.It is described to be used as separation unit
The unit that part illustrates may or may not be physically separated, and component shown as a unit can be or can also
Not to be physical unit, it can it is in one place, or may be distributed over multiple network units.It can be according to reality
Needs some or all of the units may be selected to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in embodiment provided by the invention can integrate in one processing unit, it can also
To be that each unit physically exists alone, can also be integrated in one unit with two or more units.
It, can be with if the function is realized in the form of SFU software functional unit and when sold or used as an independent product
It is stored in a computer readable storage medium.Based on this understanding, technical solution of the present invention is substantially in other words
The part of the part that contributes to existing technology or the technical solution can be embodied in the form of software products, the meter
Calculation machine software product is stored in a storage medium, including some instructions are used so that a computer equipment (can be a
People's computer, server or network equipment etc.) it performs all or part of the steps of the method described in the various embodiments of the present invention.
And storage medium above-mentioned includes:USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited
The various media that can store program code such as reservoir (RAM, Random Access Memory), magnetic or disk.
In addition, term " first ", " second ", " third " are used for description purposes only, it is not understood to indicate or imply phase
To importance.
Finally it should be noted that:Embodiment described above, only a specific embodiment of the invention, to illustrate the present invention
Technical solution, rather than its limitations, scope of protection of the present invention is not limited thereto, although with reference to the foregoing embodiments to this hair
It is bright to be described in detail, those skilled in the art should understand that:Anyone skilled in the art
In the technical scope disclosed by the present invention, it can still modify to technical solution documented by previous embodiment or can be light
It is readily conceivable that variation or equivalent replacement of some of the technical features;And these modifications, variation or replacement, do not make
The essence of corresponding technical solution is detached from the spirit and scope of technical solution of the embodiment of the present invention, should all cover in protection of the invention
Within the scope of.Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (10)
1. a kind of diffraction coefficient method for solving, which is characterized in that the method includes:
Obtain the diffracted wave under attenuation function effect;
Added-time window function is combined with the attenuation function, obtains the attenuation function of added-time window function;
Gal uncle's transformation is carried out to the diffracted wave under attenuation function effect, obtains the diffracted wave in the primary domain of gal;
According to the diffracted wave and described plus modifying factor time histories sample attenuation function in the primary domain of the gal, obtain in the primary domain of gal
Diffraction coefficient;
Diffraction coefficient in the primary domain of the gal is converted into the diffraction coefficient in time-domain, obtains subsurface structure geology in time-domain
The diffraction coefficient of body.
2. the method according to claim 1, wherein the time histories sample is the when window letter for adding modifying factor
Number.
3. the method according to claim 1, wherein according to diffracted wave and the added-time window in the primary domain of the gal
Attenuation function after function, the step of obtaining the diffraction coefficient in the primary domain of gal, including:
It is calculated using the diffracted wave in the primary domain of the gal as molecule using the attenuation function after the added-time window function as denominator
Obtain the diffraction coefficient in the primary domain of the gal.
4. the method according to claim 1, wherein obtain attenuation function effect under diffracted wave the step of it
Before, the method also includes:
Noise and random disturbances processing are carried out to the data of acquisition, the data that obtain that treated;
Treated that data are separated by described, obtains diffracted wave.
5. a kind of diffracted wave solving device, which is characterized in that described device includes:
First acquisition unit, for obtaining the diffracted wave under attenuation function effect;
Second acquisition unit obtains the attenuation function of added-time window function for combining time histories sample with the attenuation function;
Converter unit obtains the diffracted wave in the primary domain of gal for carrying out gal uncle's transformation to the diffracted wave under the attenuation;
First diffraction coefficient acquiring unit, for according to the decaying after the diffracted wave and the added-time window function in the primary domain of the gal
Function obtains the diffraction coefficient in the primary domain of gal;
Diffraction coefficient in the primary domain of the gal is converted to the diffraction coefficient in time-domain, obtained by the second diffraction coefficient acquiring unit
The diffraction coefficient of subsurface structure geologic body into time-domain.
6. device according to claim 5, which is characterized in that the time histories sample is the when window letter for adding modifying factor
Number.
7. device according to claim 5, which is characterized in that the first diffraction coefficient acquiring unit is also used to:
It is calculated using the diffracted wave in the primary domain of the gal as molecule using the attenuation function of the added-time window function as denominator
To the diffraction coefficient in the primary domain of the gal.
8. device according to claim 5, which is characterized in that described device further includes:
Diffracted wave acquiring unit, for carrying out noise to the data of acquisition and random disturbances are handled, the data that obtain that treated;It will
Described treated that data are separated, and obtains diffracted wave.
9. a kind of diffraction coefficient solves equipment, which is characterized in that including processor and memory, be stored with energy in the memory
Enough machine-executable instructions executed by the processor, the processor execute the machine-executable instruction to realize right
It is required that method described in any one of 1 to 4.
10. a kind of computer readable storage medium, which is characterized in that the computer-readable recording medium storage has machine that can hold
Row instruction, for the machine-executable instruction when being called and being executed by processor, the machine-executable instruction promotes the place
It manages device and realizes method described in any one of Claims 1-4.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811128027.3A CN108873067A (en) | 2018-09-26 | 2018-09-26 | Diffraction coefficient method for solving and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811128027.3A CN108873067A (en) | 2018-09-26 | 2018-09-26 | Diffraction coefficient method for solving and device |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108873067A true CN108873067A (en) | 2018-11-23 |
Family
ID=64324872
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811128027.3A Pending CN108873067A (en) | 2018-09-26 | 2018-09-26 | Diffraction coefficient method for solving and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108873067A (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101073020A (en) * | 2004-09-13 | 2007-11-14 | 切夫里昂美国公司 | Methods for earth modeling and seismic imaging using interactive and selective updating |
US20120116703A1 (en) * | 2009-04-24 | 2012-05-10 | Nicolas Pavillon | Method and apparatus for enhanced spatial bandwidth wavefronts reconstructed from digital interferograms or holograms |
CN106772592A (en) * | 2016-11-10 | 2017-05-31 | 中国矿业大学(北京) | The analysis method and device of diffracted wave focus energy |
CN106772578A (en) * | 2016-12-07 | 2017-05-31 | 中国矿业大学(北京) | A kind of method and apparatus of synthetic seismogram |
-
2018
- 2018-09-26 CN CN201811128027.3A patent/CN108873067A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101073020A (en) * | 2004-09-13 | 2007-11-14 | 切夫里昂美国公司 | Methods for earth modeling and seismic imaging using interactive and selective updating |
US20120116703A1 (en) * | 2009-04-24 | 2012-05-10 | Nicolas Pavillon | Method and apparatus for enhanced spatial bandwidth wavefronts reconstructed from digital interferograms or holograms |
CN106772592A (en) * | 2016-11-10 | 2017-05-31 | 中国矿业大学(北京) | The analysis method and device of diffracted wave focus energy |
CN106772578A (en) * | 2016-12-07 | 2017-05-31 | 中国矿业大学(北京) | A kind of method and apparatus of synthetic seismogram |
Non-Patent Citations (2)
Title |
---|
彭涛: "动态反褶积方法的改进及其在储层预测中的应用", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
秦雪霏: "地震波多次散射波场属性分析", 《吉林大学硕士学位论文(2017年)》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lu et al. | Seismic spectral decomposition using deconvolutive short-time Fourier transform spectrogram | |
CN104345345B (en) | A kind of shale reservoir total organic carbon TOC content prediction methods | |
CN103257361B (en) | Based on oil gas forecasting method and the system of Zoeppritz equation approximate expression | |
Salem et al. | Interpretation of magnetic data using an enhanced local wavenumber (ELW) method | |
CN103713315B (en) | A kind of seismic anisotropy parameter full waveform inversion method and device | |
CN103758511B (en) | A kind of method and device of the down-hole reverse-time migration hidden reservoir of imaging identification | |
CN104101902B (en) | Seismic properties clustering method and device | |
CN103926623B (en) | Method for suppressing reverse time migration low frequency noise | |
CN105467460B (en) | electromagnetic exploration method and device | |
CN109188520A (en) | Thin reservoir thickness prediction method and device | |
CN106226841A (en) | A kind of fluvial facies three-dimensional sedimentary facies model Decided modelling method | |
CN105425299A (en) | Method and device for determining formation fracture distribution | |
Ni et al. | Lineament length and density analyses based on the segment tracing algorithm: a case study of the gaosong field in gejiu tin mine, China | |
CN106526678A (en) | Reflection acoustic logging wave field separation method and device | |
CN104345343A (en) | Prediction method for complex seabed related interbed multiples | |
CN103777241B (en) | Three dimensional seismic data rapid edge-detection method based on time domain generalized Hilbert transformation | |
CN110118993A (en) | Diffraction wave imaging method and device | |
Maity et al. | Reservoir characterization of an unconventional reservoir by integrating microseismic, seismic, and well log data | |
CN108873067A (en) | Diffraction coefficient method for solving and device | |
CN105319594B (en) | A kind of Fourier geological data reconstructing method based on least-squares parameter inverting | |
CN103364825A (en) | Three-dimensional Fourier transform amplitude-preserving pre-stack channel interpolation method | |
CN110058319A (en) | A kind of magnetotelluric data acquisition method, device and terminal device | |
CN109143344A (en) | Diffracted wave separation method, device and electronic equipment | |
CN112305591B (en) | Tunnel advanced geological prediction method and computer readable storage medium | |
Murty et al. | Mapping the sub-trappean Mesozoic sediments in the western part of Narmada–Tapti region of Deccan Volcanic Province, India |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181123 |