CN109100804A - Improve data re-establishing method, the apparatus and system of seismic data spatial sampling attribute - Google Patents
Improve data re-establishing method, the apparatus and system of seismic data spatial sampling attribute Download PDFInfo
- Publication number
- CN109100804A CN109100804A CN201810592971.8A CN201810592971A CN109100804A CN 109100804 A CN109100804 A CN 109100804A CN 201810592971 A CN201810592971 A CN 201810592971A CN 109100804 A CN109100804 A CN 109100804A
- Authority
- CN
- China
- Prior art keywords
- data
- node
- task
- seismic data
- space
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000005070 sampling Methods 0.000 title claims abstract description 44
- 239000011159 matrix material Substances 0.000 claims abstract description 22
- 230000002123 temporal effect Effects 0.000 claims abstract description 11
- 230000003595 spectral effect Effects 0.000 claims abstract description 9
- 238000012545 processing Methods 0.000 claims description 24
- 238000004590 computer program Methods 0.000 claims description 13
- 238000003860 storage Methods 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 238000011084 recovery Methods 0.000 claims description 5
- 238000001228 spectrum Methods 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 238000007689 inspection Methods 0.000 claims description 2
- 230000006872 improvement Effects 0.000 description 14
- 230000006870 function Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005194 fractionation Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 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/32—Transforming one recording into another or one representation into another
- G01V1/325—Transforming one representation into another
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
- G01V2210/48—Other transforms
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 present application provides a kind of data re-establishing method, apparatus and system for improving seismic data spatial sampling attribute, this method comprises: obtaining task from node, task is corresponding with one of multiple data blocks in advance being split out time and space domain seismic data cube to be reconstructed by host node;Fast Fourier Transform (FFT) is carried out along time orientation to data block from node, obtains frequency space numeric field data body;From node along direction in space, Nonuniform fast Fourier transform is carried out to frequency space numeric field data body according to spatial position, obtains frequency space wave-number domain data volume;The wave number spectral component of all temporal frequency pieces in the wave-number domain data volume of frequency space is determined from node, and is stored into output matrix;Time and space domain is shifted to from node by output matrix contravariant, obtains the seismic data after data block is rebuild;It is supplied to host node from node by the seismic data after reconstruction, to carry out result merging treatment.The embodiment of the present application can be improved the versatility of data reconstruction and rebuild efficiency.
Description
Technical field
This application involves Reconstruction of seismic data technical fields, more particularly, to a kind of improvement seismic data spatial sampling attribute
Data re-establishing method, apparatus and system.
Background technique
In seismic prospecting at this stage, due to the influence of the factors such as limitation, the field construction condition of exploration funds, acquisition
To seismic data be often difficult to meet requirement of the subsequent imaging to earthquake data space systematicness and space sampling densities.It is dilute
Thin or irregular spatial sampling attribute can seriously affect the imaging effect of seismic data, and Reconstruction of seismic data technology can be with
Improve the spatial sampling attribute of observation system, to a certain extent so as to improve the imaging effect of seismic data.
Such as: Frequency-Space Domain trace interpolation technology, frequency wavenumber domain trace interpolation technology etc. can all be realized to evacuated space
The encryption in sampling interval.However, these methods require first input seismic data be rule space sampling and cannot
There is the seismic channel of missing, then, then between every adjacent 2 constructs new seismic channel.Moreover, because field construction item
Part, such as: the influence of subaerial relief condition, marine ocean current etc., the spatial sampling of original seismic data are often unsatisfactory for rule
Sampling or the requirement without missing, therefore such methods generally can be used only in the encryption of poststack data.
In addition, improving the number of seismic data spatial sampling attribute with high covering, the popularization and application of high-density acquisition technique
Need to handle tens TB grade of mass data according to rebuilding, if in a manner of conventional serial and the more operation modes realization data weights of variable
It builds, computational efficiency will seriously affect its application in production.
Therefore, how to improve the versatility for being intended to improve the data reconstruction of seismic data spatial sampling attribute and rebuild efficiency
It is a technical problem to be solved urgently.
Summary of the invention
A kind of data re-establishing method for being designed to provide improvement seismic data spatial sampling attribute of the embodiment of the present application,
Apparatus and system, to improve the versatility for the data reconstruction for being intended to improve seismic data spatial sampling attribute and rebuild efficiency.
In order to achieve the above objectives, on the one hand, the embodiment of the present application provides a kind of improvement seismic data spatial sampling attribute
Data re-establishing method, comprising:
From node obtain task, the task with time and space domain seismic data cube to be reconstructed is torn open by host node in advance
One of multiple data blocks separated are corresponding;
It is described that Fast Fourier Transform (FFT) is carried out along time orientation to the data block from node, obtain frequency space domain number
According to body;
It is described from node along direction in space, non-homogeneous Fu is carried out to the frequency space numeric field data body according to spatial position
In leaf transformation, obtain frequency space wave-number domain data volume;
The wave number spectral component that all temporal frequency pieces in the frequency space wave-number domain data volume are determined from node, and
It stores into output matrix;
It is described that the output matrix contravariant is shifted into time and space domain from node, obtain the earthquake after the data block is rebuild
Data;
It is described that the seismic data after the reconstruction is supplied to the host node from node, to carry out result merging treatment.
In the data re-establishing method of the improvement seismic data spatial sampling attribute of the embodiment of the present application, the data block that splits out
Number is the product of the spatial window number of three assigned directions within the scope of data reconstruction.
In the data re-establishing method of the improvement seismic data spatial sampling attribute of the embodiment of the present application, described three specified
Direction includes:
The direction CMP Line, the direction CMP and the direction Offset-x;Alternatively,
The direction CMP Line, the direction CMP and the direction Offset bin.
It is described built-in from node in the data re-establishing method of the improvement seismic data spatial sampling attribute of the embodiment of the present application
It is vertical to have task requests and reading data control thread, calculate control thread and result output control line journey.
In the data re-establishing method of the improvement seismic data spatial sampling attribute of the embodiment of the present application, it is from node when described
When multi-core node, being somebody's turn to do each CPU core out of node and establishing has respective task requests and reading data control thread, calculates control
Thread processed and result output control line journey.
It is described to be obtained from node in the data re-establishing method of the improvement seismic data spatial sampling attribute of the embodiment of the present application
Task, comprising:
It is described from node to the host node request task, and task ID that the last time is completed and handle sub- result and transmit
To the host node, the last task ID completed and the sub- result of processing are saved as into information record text by the host node
Part;
It is described that the task that the host node is distributed according to task list is received from node;The task list includes data block
ID and three assigned directions.
In the data re-establishing method of the improvement seismic data spatial sampling attribute of the embodiment of the present application, the information record text
The authority that part doubles as fault recovery;When it is described from node failure when, by the host node according to the information record file will
Not completing from node for the task distributes to other in system from node processing.
In the data re-establishing method of the improvement seismic data spatial sampling attribute of the embodiment of the present application, further includes:
After restarting operation, as host node task list according to information record file reconstruction, and check
The correctness of the sub- result of processing.
On the other hand, the embodiment of the present application also provides a kind of data reconstruction systems for improving seismic data spatial sampling attribute
System, including host node and multiple from node, it is each described to include: from node
Task acquisition module, for obtaining task, the task in advance by host node by time and space domain to be reconstructed
One of multiple data blocks that seismic data cube is split out are corresponding;
First conversion module obtains frequency for carrying out Fast Fourier Transform (FFT) along time orientation to the data block
Spatial domain data body;
Second conversion module, for being carried out to the frequency space numeric field data body according to spatial position along direction in space
Nonuniform fast Fourier transform obtains frequency space wave-number domain data volume;
Data memory module, for determining the wave-number spectrum of all temporal frequency pieces in the frequency space wave-number domain data volume
Component, and store into output matrix;
Third conversion module obtains the data block and rebuilds for the output matrix contravariant to be shifted to time and space domain
Seismic data afterwards;
Result treatment module, for the seismic data after the reconstruction to be supplied to the host node, to carry out result conjunction
And it handles.
On the other hand, the embodiment of the present application also provides a kind of data reconstruction dresses for improving seismic data spatial sampling attribute
It sets, including memory, processor and the computer program being stored on the memory, the computer program is described
Processor executes following steps when running:
Acquisition task, the task in advance time and space domain seismic data cube to be reconstructed is split out by host node
One of multiple data blocks are corresponding;
Fast Fourier Transform (FFT) is carried out along time orientation to the data block, obtains frequency space numeric field data body;
Along direction in space, Nonuniform fast Fourier transform is carried out to the frequency space numeric field data body according to spatial position,
Obtain frequency space wave-number domain data volume;
It determines the wave number spectral component of all temporal frequency pieces in the frequency space wave-number domain data volume, and stores to output
In matrix;
The output matrix contravariant is shifted into time and space domain, obtains the seismic data after the data block is rebuild;
Seismic data after the reconstruction is supplied to the host node, to carry out result merging treatment.
As can be seen from the technical scheme provided by the above embodiments of the present application, the embodiment of the present application is by the way that seismic data to be distributed to
Seismic data weight is greatly shortened on the basis of improving seismic data image quality in parallel processing on different nodes
The time built, to substantially increase the versatility for being intended to improve the data reconstruction of seismic data spatial sampling attribute and rebuild effect
Rate, and since spatial sampling of the embodiment of the present application for original seismic data is specific not as the prior art
It is required that therefore the data reconstruction of embodiments herein versatility it is stronger.
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 that the structure for the data reconstruction system for improving seismic data spatial sampling attribute in some embodiments of the application is shown
It is intended to;
Fig. 2 is the parallel place for improving the data reconstruction system of seismic data spatial sampling attribute in some embodiments of the application
Manage schematic diagram;
Fig. 3 is the process for improving the data re-establishing method of seismic data spatial sampling attribute in some embodiments of the application
Figure;
Fig. 4 is the structural frames for improving the data reconstruction device of seismic data spatial sampling attribute in some embodiments of the application
Figure;
Fig. 5 is the structure for improving the data reconstruction device of seismic data spatial sampling attribute in other embodiments of the application
Block diagram.
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.Such as in being described below, second component is formed above the first component, may include the first component and second component
It can also include the first component and second component in a manner of non-direct contact (i.e. the with the embodiment that the direct way of contact is formed
Can also include additional component between one component and second component) embodiment etc. that is formed.
Moreover, for ease of description, some embodiments of the application can be used such as " in ... top ", " ... under ",
The spatially relative terms such as " top ", " lower section ", with description such as each element shown in the drawings of embodiment or component and another
Relationship between (or other) element or component.It should be understood that space is opposite other than the orientation described in the attached drawing
Term also aims to the different direction including device in use or operation.If such as the device in attached drawing is reversed, and is described
For the element or component of " " other elements or component " below " or " under ", it will then be positioned as " " other elements or portion
Part " top " or " on ".
Refering to what is shown in Fig. 1, the data reconstruction system of the improvement seismic data spatial sampling attribute of the embodiment of the present application can be with
Including host node and multiple from node.Wherein, multiple to obtain task for parallel slave host node from node, handle the task
And the result obtained after processing is supplied to host node, such as shown in Fig. 2.Host node can be used for time and space to be reconstructed
Domain seismic data cube is split as multiple data blocks, and accordingly generates task list;Task list task by multiple from node
After the completion of processing, multiple results exported from node are merged into processing.
As shown in connection with fig. 3, using above-mentioned from node as executing subject, the improvement seismic data space of the embodiment of the present application is adopted
The data re-establishing method of sample attribute may comprise steps of:
S301, obtain task from node, the task in advance by host node by time and space domain earthquake number to be reconstructed
One of multiple data blocks split out according to body are corresponding.
In one embodiment of the application, before node acquisition task, the host node is in advance by the time to be reconstructed
Spatial domain seismic data cube splits into multiple data blocks, and each data block is a small time and space domain seismic data cube, and
Task list is established accordingly, to be used as task distribution.
It should be noted that the fractionation of above-mentioned time and space domain seismic data cube to be reconstructed, can be according to space
Data within the scope of data reconstruction are divided into N number of data block by the size of window, and the size of N is three within the scope of data reconstruction
The product of the spatial window number of assigned direction.In one exemplary embodiment, three assigned directions can be CMP
The direction Line (common point line direction), the direction CMP (common point direction) and the direction Offset-x (the x-axis side of offset distance
To);It or can be the direction CMP Line, the direction CMP and the direction Offset bin (offset distance face element direction).With CMP Line
For direction, the direction CMP and the direction Offset-x are as three assigned directions, N=(the spatial window number in the direction CMP Line) *
(the spatial window number in the direction CMP) * (the spatial window number in the direction Offset-x).
In one embodiment of the application, the attribute of the task list of generation may include that data block ID and described three specify
Direction.
It is described to may include steps of from node acquisition task in one embodiment of the application:
It is described from node to the host node request task, and task ID that the last time is completed and handle sub- result and transmit
To the host node, the last task ID completed and the sub- result of processing are saved as into information record text by the host node
Part;Then described to can receive the task that the host node is distributed according to task list from node.Wherein, the information records file
It can also be as the authority for fault recovery.Further, it is also possible to be to send under operation User Catalog by information record file designation
Job name, in order to avoid operation conflict between user.
It is corresponding, when it is described from node failure when, can by the host node according to the information record file will it is described from
The unfinished task of node distributes to other in system from node processing.Nodal test is carried out by host node, when some meter
Operator node breaks down, and host node carries out result inspection according to record the file information, and its unfinished task is transferred to system
In it is other from node complete.In this way, node fault-tolerance ability can be had by making the method for the embodiment of the present application.
S302, it is described from node to the data block along time orientation carry out Fast Fourier Transform (FFT), obtain frequency sky
Between numeric field data body.
In one embodiment of the application, data block can be converted by time and space domain by frequency by Fast Fourier Transform (FFT)
Spatial domain, correspondingly, data block is just transformed to frequency space numeric field data body, in order to subsequent processing.
S303, it is described from node along direction in space, the frequency space numeric field data body is carried out according to spatial position non-
Uniform Fourier transformation obtains frequency space wave-number domain data volume.
It, can be by frequency space numeric field data body by frequency space by Nonuniform fast Fourier transform in one embodiment of the application
Domain transforms to frequency space wave-number domain, correspondingly, frequency space numeric field data body is just transformed to frequency space wave-number domain data volume.
S304, the wave-number spectrum point that all temporal frequency pieces in the frequency space wave-number domain data volume are determined from node
Amount, and store into output matrix.
In one embodiment of the application, in frequency space wave-number domain, it can be calculated according to frequency space wave-number domain data volume
Obtain the wave number spectral component of each temporal frequency piece in the frequency space wave-number domain data volume.Such as it can be according to the frequency
Constraint matrix is calculated in the amplitude spectrum of space wave numeric field data body, then according to the constraint matrix, when being calculated each
Between frequency chip wave number spectral component.
S305, it is described the output matrix contravariant is shifted into time and space domain from node, obtain after the data block rebuilds
Seismic data.
S306, it is described the seismic data after the reconstruction is supplied to the host node from node, to carry out result merging
Processing, in order to obtain final reconstructed results.In one embodiment of the application, as a result merges and can be information record file
Merge, needs to count the trace header information of each node data again, physical stored file can be constant, can show as one in logic
A file.It generates after obtaining final reconstructed results, final reconstructed results, such as Fig. 2 can be written in host node into database
It is shown.
It is described that task requests and reading data control thread, meter can have been established out of node in one embodiment of the application
Calculate control thread and result output control line journey.In this way, the processing of each node is parallel;In each intra-node, task
Please summed data read, calculate, output be also it is parallel, to substantially reduce the time of data reconstruction, improve data weight
The efficiency built.In addition, when it is described from node be multi-core node when, each CPU core should establish out of node and have respective task
Request and reading data control thread calculate control thread and result output control line journey, so as to further make full use of
The resource of node carries out parallel processing.In one exemplary embodiment, such as using pthread method thread is created, and can
Thread control is carried out using pthread_mutex_t function and pthread_cond_t function.
In one embodiment of the application, during job execution, when causing job execution to lose because of various machine reasons etc.
Operation need to be restarted when losing.It, can be as host node task column according to the information record file reconstruction after restarting operation
Table, and check the correctness of the sub- result of processing.In this way, fault recovery can be had the function of by making the data reconstruction of the application.
It can be seen that the embodiment of the present application by the way that seismic data to be distributed to parallel processing on different nodes, is improving
On the basis of seismic data image quality, the time of seismic data reconstruction is greatly shortened, and the embodiment of the present application is also
Have the function of node fault-tolerance and fault recovery, truly meets seismic prospecting data and handle industrialized needs.
Refering to what is shown in Fig. 4, in logical construction, it is above-mentioned each to may include: from node
Task acquisition module 41, can be used for obtaining task, the task in advance by host node by the time to be reconstructed
One of multiple data blocks that spatial domain seismic data cube is split out are corresponding;
First conversion module 42 can be used for carrying out Fast Fourier Transform (FFT) along time orientation to the data block, obtain
To frequency space numeric field data body;
Second conversion module 43 can be used for along direction in space, according to spatial position to the frequency space numeric field data
Body carries out Nonuniform fast Fourier transform, obtains frequency space wave-number domain data volume;
Data memory module 44 is determined for all temporal frequency pieces in the frequency space wave-number domain data volume
Wave number spectral component, and store into output matrix;
Third conversion module 45 can be used for the output matrix contravariant shifting to time and space domain, obtain the data
Seismic data after block reconstruction;
Result treatment module 46 can be used for the seismic data after the reconstruction being supplied to the host node, to carry out
As a result merging treatment.
Refering to what is shown in Fig. 5, the data reconstruction device of the improvement seismic data spatial sampling attribute of the embodiment of the present application, it can be with
Including memory, processor and the computer program being stored on the memory, the computer program is by the processing
Device executes following steps when running:
Acquisition task, the task in advance time and space domain seismic data cube to be reconstructed is split out by host node
One of multiple data blocks are corresponding;
Fast Fourier Transform (FFT) is carried out along time orientation to the data block, obtains frequency space numeric field data body;
Along direction in space, Nonuniform fast Fourier transform is carried out to the frequency space numeric field data body according to spatial position,
Obtain frequency space wave-number domain data volume;
It determines the wave number spectral component of all temporal frequency pieces in the frequency space wave-number domain data volume, and stores to output
In matrix;
The output matrix contravariant is shifted into time and space domain, obtains the seismic data after the data block is rebuild;
Seismic data after the reconstruction is supplied to the host node, to carry out result merging treatment.
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).
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.
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 magnetic storage devices
Or any other non-transmission medium, can be used for storage can be accessed by a computing device information.As defined in this article, it calculates
Machine readable medium does not include temporary computer readable media (transitory media), such as the data-signal and carrier wave of modulation.
It should also be noted that, the terms "include", "comprise" or its any other variant are intended to nonexcludability
Include, so that process, method or equipment including a series of elements not only include those elements, but also including not having
There is the other element being expressly recited, or further includes for this process, method or the intrinsic element of equipment.Do not having more
In the case where more limitations, the element that is limited by sentence "including a ...", it is not excluded that including process, the side of the element
There is also other identical elements in method 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 (10)
1. a kind of data re-establishing method for improving seismic data spatial sampling attribute characterized by comprising
From node obtain task, the task with time and space domain seismic data cube to be reconstructed is split out by host node in advance
One of multiple data blocks it is corresponding;
It is described that Fast Fourier Transform (FFT) is carried out along time orientation to the data block from node, obtain frequency space numeric field data
Body;
It is described from node along direction in space, non-homogeneous Fourier is carried out to the frequency space numeric field data body according to spatial position
Transformation, obtains frequency space wave-number domain data volume;
The wave number spectral component that all temporal frequency pieces in the frequency space wave-number domain data volume are determined from node, and store
Into output matrix;
It is described that the output matrix contravariant is shifted into time and space domain from node, obtain the earthquake number after the data block is rebuild
According to;
It is described that the seismic data after the reconstruction is supplied to the host node from node, to carry out result merging treatment.
2. improving the data re-establishing method of seismic data spatial sampling attribute as described in claim 1, which is characterized in that split
Data block number out is the product of the spatial window number of three assigned directions within the scope of data reconstruction.
3. improving the data re-establishing method of seismic data spatial sampling attribute as claimed in claim 2, which is characterized in that described
Three assigned directions include:
The direction CMP Line, the direction CMP and the direction Offset-x;Alternatively,
The direction CMP Line, the direction CMP and the direction Offset bin.
4. improving the data re-establishing method of seismic data spatial sampling attribute as described in claim 1, which is characterized in that described
Establishing out of node has task requests and reading data control thread, calculates control thread and result output control line journey.
5. improving the data re-establishing method of seismic data spatial sampling attribute as claimed in claim 2, which is characterized in that work as institute
State from node be multi-core node when, each CPU core should establish out of node and there is respective task requests and reading data to control
Thread calculates control thread and result output control line journey.
6. improving the data re-establishing method of seismic data spatial sampling attribute as claimed in claim 2, which is characterized in that described
Task is obtained from node, comprising:
It is described from node to the host node request task, and by task ID that the last time completes and handle sub- result and send institute to
Host node is stated, the last task ID completed and the sub- result of processing are saved as into information record file by the host node;
It is described that the task that the host node is distributed according to task list is received from node;The task list include data block ID and
Three assigned directions.
7. improving the data re-establishing method of seismic data spatial sampling attribute as claimed in claim 6, which is characterized in that described
Information records file as the authority for fault recovery;When it is described from node failure when, by the host node according to the information
Not completing from node for the task is distributed to other in system from node processing by record file.
8. improving the data re-establishing method of seismic data spatial sampling attribute as claimed in claim 6, which is characterized in that also wrap
It includes:
After restarting operation, as host node task list according to information record file reconstruction, and described in inspection
Handle the correctness of sub- result.
9. a kind of data reconstruction system for improving seismic data spatial sampling attribute, which is characterized in that including host node and multiple
It is each described to include: from node from node
Task acquisition module, for obtaining task, the task in advance by host node by time and space domain earthquake to be reconstructed
One of multiple data blocks that data volume is split out are corresponding;
First conversion module obtains frequency space for carrying out Fast Fourier Transform (FFT) along time orientation to the data block
Numeric field data body;
Second conversion module, for being carried out according to spatial position to the frequency space numeric field data body non-equal along direction in space
Even Fourier transformation obtains frequency space wave-number domain data volume;
Data memory module, for determining the wave-number spectrum point of all temporal frequency pieces in the frequency space wave-number domain data volume
Amount, and store into output matrix;
Third conversion module, for the output matrix contravariant to be shifted to time and space domain, after obtaining the data block reconstruction
Seismic data;
Result treatment module, for the seismic data after the reconstruction to be supplied to the host node, to carry out at result merging
Reason.
10. a kind of data reconstruction device for improving seismic data spatial sampling attribute, including memory, processor and storage
Computer program on the memory, which is characterized in that executed such as when the computer program is run by the processor
Lower step:
Acquisition task, the task with time and space domain seismic data cube to be reconstructed split out by host node in advance it is multiple
One of data block is corresponding;
Fast Fourier Transform (FFT) is carried out along time orientation to the data block, obtains frequency space numeric field data body;
Along direction in space, Nonuniform fast Fourier transform is carried out to the frequency space numeric field data body according to spatial position, is obtained
Frequency space wave-number domain data volume;
It determines the wave number spectral component of all temporal frequency pieces in the frequency space wave-number domain data volume, and stores to output matrix
In;
The output matrix contravariant is shifted into time and space domain, obtains the seismic data after the data block is rebuild;
Seismic data after the reconstruction is supplied to the host node, to carry out result merging treatment.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810592971.8A CN109100804B (en) | 2018-06-11 | 2018-06-11 | Data reconstruction method, device and system for improving seismic data space sampling attribute |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810592971.8A CN109100804B (en) | 2018-06-11 | 2018-06-11 | Data reconstruction method, device and system for improving seismic data space sampling attribute |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109100804A true CN109100804A (en) | 2018-12-28 |
CN109100804B CN109100804B (en) | 2020-10-09 |
Family
ID=64796793
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810592971.8A Active CN109100804B (en) | 2018-06-11 | 2018-06-11 | Data reconstruction method, device and system for improving seismic data space sampling attribute |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109100804B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112230284A (en) * | 2019-07-15 | 2021-01-15 | 中国石油天然气集团有限公司 | Parallel random noise attenuation method, monitoring method and node for three-dimensional pre-stack data |
CN115659579A (en) * | 2022-09-05 | 2023-01-31 | 中南大学 | Three-dimensional gravity field numerical simulation method and system based on 3D AS-FT |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104459781A (en) * | 2014-12-09 | 2015-03-25 | 中国石油天然气集团公司 | Three-dimensional pre-stack seismic data random noise degeneration method |
CN105403913A (en) * | 2015-10-30 | 2016-03-16 | 中国石油天然气集团公司 | Pre-stack depth migration method and device |
CN106371138A (en) * | 2016-08-17 | 2017-02-01 | 中国石油天然气集团公司 | Seismic data reconstruction method and apparatus |
US9632193B2 (en) * | 2013-11-01 | 2017-04-25 | Conocophillips Company | Compressive sensing |
US10261207B2 (en) * | 2014-12-18 | 2019-04-16 | Pgs Geophysical As | Seismic noise mitigation system and method |
-
2018
- 2018-06-11 CN CN201810592971.8A patent/CN109100804B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9632193B2 (en) * | 2013-11-01 | 2017-04-25 | Conocophillips Company | Compressive sensing |
CN104459781A (en) * | 2014-12-09 | 2015-03-25 | 中国石油天然气集团公司 | Three-dimensional pre-stack seismic data random noise degeneration method |
US10261207B2 (en) * | 2014-12-18 | 2019-04-16 | Pgs Geophysical As | Seismic noise mitigation system and method |
CN105403913A (en) * | 2015-10-30 | 2016-03-16 | 中国石油天然气集团公司 | Pre-stack depth migration method and device |
CN106371138A (en) * | 2016-08-17 | 2017-02-01 | 中国石油天然气集团公司 | Seismic data reconstruction method and apparatus |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112230284A (en) * | 2019-07-15 | 2021-01-15 | 中国石油天然气集团有限公司 | Parallel random noise attenuation method, monitoring method and node for three-dimensional pre-stack data |
CN115659579A (en) * | 2022-09-05 | 2023-01-31 | 中南大学 | Three-dimensional gravity field numerical simulation method and system based on 3D AS-FT |
Also Published As
Publication number | Publication date |
---|---|
CN109100804B (en) | 2020-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chamberlain et al. | EQcorrscan: Repeating and near‐repeating earthquake detection and analysis in Python | |
Yang et al. | A nearly analytic discrete method for acoustic and elastic wave equations in anisotropic media | |
Krischer et al. | Large‐scale seismic inversion framework | |
CN110019835B (en) | Resource arranging method and device and electronic equipment | |
CN103455396A (en) | Test method and test device of hardware performance of electronic equipment | |
CN110082823B (en) | Seismic data interpolation method and device | |
Belonosov et al. | 3D numerical simulation of elastic waves with a frequency-domain iterative solver | |
CN107144877B (en) | Abnormal structure identification method and device | |
CN109100804A (en) | Improve data re-establishing method, the apparatus and system of seismic data spatial sampling attribute | |
Noble et al. | High-performance 3D first-arrival traveltime tomography | |
CN108445532A (en) | A kind of Depth Domain inverse migration method and device | |
Panetta et al. | Accelerating Kirchhoff migration by CPU and GPU cooperation | |
Perera et al. | Database scaling on Kubernetes | |
Okita et al. | High‐performance computing strategies for seismic‐imaging software on the cluster and cloud‐computing environments | |
CN107831541B (en) | Thin strate recognition methods and device based on high density VSP data | |
Li et al. | A 3D reflection ray-tracing method based on linear traveltime perturbation interpolation | |
CN112764105B (en) | HTI medium quasi-longitudinal wave forward simulation method and device, storage medium and processor | |
CN109597964A (en) | Theoretical flow through oil reservoir parameter value determines method, system, equipment, readable medium | |
Yang et al. | Unstructured mesh based elastic wave modelling on GPU: a double-mesh grid method | |
Wang et al. | Enhancement of computational efficiency for weighted total least squares | |
Xing et al. | Automated parallel data processing engine with application to large-scale feature extraction | |
US10082588B2 (en) | Adaptive structure-oriented operator | |
CN112230284A (en) | Parallel random noise attenuation method, monitoring method and node for three-dimensional pre-stack data | |
Saupi Teri et al. | GPU utilization in geoprocessing BIG Geodata: A review | |
CN113126162B (en) | Random noise attenuation calculation method and device |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |