CN109100804B - Data reconstruction method, device and system for improving seismic data space sampling attribute - Google Patents
Data reconstruction method, device and system for improving seismic data space sampling attribute Download PDFInfo
- Publication number
- CN109100804B CN109100804B CN201810592971.8A CN201810592971A CN109100804B CN 109100804 B CN109100804 B CN 109100804B CN 201810592971 A CN201810592971 A CN 201810592971A CN 109100804 B CN109100804 B CN 109100804B
- Authority
- CN
- China
- Prior art keywords
- data
- space
- frequency
- seismic data
- task
- 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
Images
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
- 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
Abstract
The embodiment of the application provides a data reconstruction method, a device and a system for improving seismic data space sampling attributes, wherein the method comprises the following steps: acquiring a task from a slave node, wherein the task corresponds to one of a plurality of data blocks which are obtained by splitting a seismic data body of a time-space domain to be reconstructed by a master node in advance; carrying out fast Fourier transform on the data block along the time direction by the slave node to obtain a frequency-space domain data body; the slave node performs non-uniform Fourier transform on the frequency-space domain data volume according to the space position along the space direction to obtain a frequency-space wave number domain data volume; wave number spectrum components of all time frequency slices in the frequency space wave number domain data body are determined from the nodes and stored in an output matrix; the slave node inversely transforms the output matrix to a time space domain to obtain seismic data after data block reconstruction; and the slave node provides the reconstructed seismic data to the master node for result combination processing. According to the data reconstruction method and device, the universality and the reconstruction efficiency of data reconstruction can be improved.
Description
Technical Field
The application relates to the technical field of seismic data reconstruction, in particular to a data reconstruction method, a data reconstruction device and a data reconstruction system for improving seismic data space sampling attributes.
Background
In seismic exploration at the present stage, due to the influence of factors such as the limitation of exploration expenditure, field construction conditions and the like, the acquired seismic data are difficult to meet the requirements of subsequent imaging on the spatial regularity and spatial sampling density of the seismic data. Sparse or irregular spatial sampling properties can seriously affect the imaging effect of the seismic data, and the seismic data reconstruction technology can improve the spatial sampling properties of an observation system to a certain extent, so that the imaging effect of the seismic data is improved.
For example: the frequency-space domain channel interpolation technology, the frequency-wavenumber domain channel interpolation technology and the like can realize the encryption processing of the sparse space sampling interval. However, these methods first require that the input seismic data be regularly spatially sampled and must not have missing seismic traces, and then construct new seismic traces between each adjacent 2 traces. Moreover, due to field construction conditions, for example: due to the influence of land terrain conditions, ocean currents, etc., spatial sampling of raw seismic data often does not meet the requirements of regular sampling or no missing, and therefore such methods can generally only be used in the encryption process of post-stack data.
In addition, with the popularization and application of high-coverage and high-density acquisition technologies, data reconstruction for improving the spatial sampling attribute of seismic data needs to process mass data of dozens of TB levels, and if data reconstruction is realized in a traditional serial mode and a variable multi-operation mode, the computing efficiency will seriously affect the application of the data reconstruction in production.
Therefore, how to improve the universality and reconstruction efficiency of data reconstruction aiming at improving the spatial sampling attribute of seismic data is a technical problem to be solved urgently at present.
Disclosure of Invention
An object of the embodiments of the present application is to provide a data reconstruction method, apparatus, and system for improving spatial sampling attributes of seismic data, so as to improve universality and reconstruction efficiency of data reconstruction aiming at improving spatial sampling attributes of seismic data.
In order to achieve the above object, in one aspect, an embodiment of the present application provides a data reconstruction method for improving spatial sampling properties of seismic data, including:
acquiring a task from a slave node, wherein the task corresponds to one of a plurality of data blocks which are obtained by splitting a seismic data body of a time-space domain to be reconstructed by a master node in advance;
the slave node performs fast Fourier transform on the data block along the time direction to obtain a frequency-space domain data body;
the slave node performs non-uniform Fourier transform on the frequency-space domain data volume according to a space position along a space direction to obtain a frequency-space wave number domain data volume;
the slave node determines wave number spectrum components of all time frequency slices in the frequency space wave number domain data body and stores the wave number spectrum components into an output matrix;
the slave node inversely transforms the output matrix to a time space domain to obtain the seismic data after the data block is reconstructed;
and the slave node provides the reconstructed seismic data to the master node for result combination processing.
In the data reconstruction method for improving the seismic data spatial sampling attribute of the embodiment of the application, the number of the split data blocks is the product of the number of spatial windows in three specified directions in a data reconstruction range.
In the data reconstruction method for improving the spatial sampling property of the seismic data according to the embodiment of the application, the three designated directions include:
CMP Line direction, CMP direction, and Offset-x direction; alternatively, the first and second electrodes may be,
the CMP Line direction, the CMP direction, and the Offset bin direction.
According to the data reconstruction method for improving the seismic data spatial sampling attribute, a task request and data reading control thread, a calculation control thread and a result output control thread are established in the slave nodes.
In the data reconstruction method for improving the seismic data spatial sampling attribute according to the embodiment of the application, when the slave node is a multi-core node, each CPU core in the slave node establishes a task request and data reading control thread, a calculation control thread, and a result output control thread.
In the data reconstruction method for improving the spatial sampling property of the seismic data according to the embodiment of the application, the task acquisition from the node includes:
the slave node requests the task from the master node, transmits the last completed task ID and the processing sub-result to the master node, and the master node saves the last completed task ID and the processing sub-result as an information recording file;
the slave node receives the tasks distributed by the master node according to the task list; the task list includes data block IDs and the three specified directions.
In the data reconstruction method for improving the seismic data spatial sampling attribute of the embodiment of the application, the information recording file is also used as a proof for fault recovery; when the slave node fails, the master node distributes tasks which are not completed by the slave node to other slave nodes in the system for processing according to the information record file.
In the data reconstruction method for improving the spatial sampling property of the seismic data according to the embodiment of the application, the method further includes:
and after the operation is restarted, the main node reconstructs the task list according to the information recording file and checks the correctness of the processing sub-result.
In another aspect, an embodiment of the present application further provides a data reconstruction system for improving spatial sampling properties of seismic data, including a master node and a plurality of slave nodes, where each slave node includes:
the task acquisition module is used for acquiring a task, wherein the task corresponds to one of a plurality of data blocks which are obtained by splitting a seismic data body of a time space domain to be reconstructed by a main node in advance;
the first transformation module is used for performing fast Fourier transformation on the data block along the time direction to obtain a frequency-space domain data volume;
the second transformation module is used for carrying out non-uniform Fourier transformation on the frequency-space domain data volume according to the space position along the space direction to obtain a frequency-space wave number domain data volume;
the data storage module is used for determining wave number spectrum components of all time frequency slices in the frequency space wave number domain data body and storing the wave number spectrum components into an output matrix;
the third transformation module is used for inversely transforming the output matrix to a time space domain to obtain the seismic data after the data block is reconstructed;
and the result processing module is used for providing the reconstructed seismic data to the main node so as to carry out result combination processing.
In another aspect, an embodiment of the present application further provides a data reconstruction apparatus for improving spatial sampling properties of seismic data, including a memory, a processor, and a computer program stored on the memory, where the computer program is executed by the processor to perform the following steps:
acquiring a task, wherein the task corresponds to one of a plurality of data blocks which are obtained by splitting a seismic data body of a time space domain to be reconstructed by a main node in advance;
carrying out fast Fourier transform on the data block along the time direction to obtain a frequency-space domain data volume;
carrying out non-uniform Fourier transform on the frequency-space domain data volume according to the space position along the space direction to obtain a frequency-space wave number domain data volume;
determining wave number spectrum components of all time frequency slices in the frequency space wave number domain data body, and storing the wave number spectrum components in an output matrix;
inversely transforming the output matrix to a time space domain to obtain the seismic data after the data block is reconstructed;
and providing the reconstructed seismic data to the main node for result combination processing.
According to the technical scheme provided by the embodiment of the application, the seismic data are distributed to different nodes to be processed in parallel, the seismic data reconstruction time is greatly shortened on the basis of improving the imaging quality of the seismic data, and therefore the universality and reconstruction efficiency of data reconstruction aiming at improving the spatial sampling attribute of the seismic data are greatly improved.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, it is obvious that the drawings in the following description are only some embodiments described in the present application, and for those skilled in the art, other drawings can be obtained according to the drawings without any creative effort. In the drawings:
FIG. 1 is a schematic block diagram of a data reconstruction system for improving spatial sampling properties of seismic data according to some embodiments of the present application;
FIG. 2 is a schematic diagram of parallel processing by a data reconstruction system for improving spatial sampling properties of seismic data according to some embodiments of the present application;
FIG. 3 is a flow chart of a data reconstruction method for improving spatial sampling properties of seismic data according to some embodiments of the present application;
FIG. 4 is a block diagram of a data reconstruction device for improving spatial sampling properties of seismic data according to some embodiments of the present application;
FIG. 5 is a block diagram of a data reconstruction device for improving spatial sampling properties of seismic data according to further embodiments of the present application.
Detailed Description
In order to make those skilled in the art better understand the technical solutions in the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application. For example, in the following description, forming the second component over the first component may include embodiments in which the first and second components are formed in direct contact, embodiments in which the first and second components are formed in non-direct contact (i.e., additional components may be included between the first and second components), and so on.
Also, for ease of description, some embodiments of the present application may use spatially relative terms such as "above …," "below …," "top," "below," etc., to describe the relationship of one element or component to another (or other) element or component as illustrated in the various figures of the embodiments. It will be understood that the spatially relative terms are intended to encompass different orientations of the device in use or operation in addition to the orientation depicted in the figures. For example, if the device in the figures is turned over, elements or components described as "below" or "beneath" other elements or components would then be oriented "above" or "over" the other elements or components.
Referring to fig. 1, a data reconstruction system for improving spatial sampling properties of seismic data according to an embodiment of the present application may include a master node and a plurality of slave nodes. The plurality of slave nodes are configured to acquire tasks from the master node in parallel, process the tasks, and provide results obtained after the processing to the master node, for example, as shown in fig. 2. The master node can be used for splitting a seismic data body of a time-space domain to be reconstructed into a plurality of data blocks and generating a task list according to the data blocks; and after the tasks of the task list are processed and finished by the plurality of slave nodes, combining the results output by the plurality of slave nodes.
With reference to fig. 3, taking the slave node as an execution subject, the data reconstruction method for improving the spatial sampling property of seismic data according to the embodiment of the present application may include the following steps:
s301, a task is obtained from the slave node, and the task corresponds to one of a plurality of data blocks which are obtained by splitting a seismic data body of a time space domain to be reconstructed by the master node in advance.
In an embodiment of the application, before acquiring the task from the slave node, the master node splits a time-space domain seismic data volume to be reconstructed into a plurality of data blocks in advance, each data block is a small time-space domain seismic data volume, and establishes a task list according to the data blocks for task allocation.
It should be noted that, the splitting of the time-space domain seismic data volume to be reconstructed may be dividing data in a data reconstruction range into N data blocks according to the size of a space window, where the size of N is a product of the number of space windows in three specified directions in the data reconstruction range. In an exemplary embodiment, the three specified directions may be a CMPLine direction (a common midpoint line direction), a CMP direction (a common midpoint direction), and an Offset-x direction (an x-axis direction of an Offset distance); or may be a CMP Line direction, a CMP direction, and an Offset bin direction (Offset bin direction). Taking the CMP Line direction, the CMP direction, and the Offset-x direction as examples of the three designated directions, N is (the number of spatial windows in the CMP Line direction) × (the number of spatial windows in the CMP direction) × (the number of spatial windows in the Offset-x direction).
In an embodiment of the present application, the attributes of the generated task list may include the data block ID and the three specified directions.
In an embodiment of the present application, the acquiring a task from a node may include the following steps:
the slave node requests the task from the master node, transmits the last completed task ID and the processing sub-result to the master node, and the master node saves the last completed task ID and the processing sub-result as an information recording file; the slave node may then receive the tasks assigned by the master node according to the task list. Wherein, the information recording file can also be used as a proof for fault recovery. In addition, the information recording file can be named as the job name under the directory of the user sending the job, so that the job conflict among the users can be avoided.
Correspondingly, when the slave node fails, the master node can distribute tasks which are not completed by the slave node to other slave nodes in the system for processing according to the information record file. When a certain computing node fails, the master node checks the result according to the recorded file information and delivers the unfinished tasks to other slave nodes in the system to finish the tasks. Therefore, the method of the embodiment of the application can have the node fault tolerance capability.
S302, the slave node performs fast Fourier transform on the data block along the time direction to obtain a frequency-space domain data body.
In an embodiment of the present application, the data block may be transformed from the time space domain to the frequency space domain by fast fourier transform, and accordingly, the data block is transformed into a frequency space domain data volume for subsequent processing.
And S303, the slave node performs non-uniform Fourier transform on the frequency-space domain data volume according to the space position along the space direction to obtain the frequency-space wave number domain data volume.
In an embodiment of the present application, the frequency-space domain data volume is transformed from the frequency-space domain to the frequency-space wave number domain by non-uniform fourier transform, and accordingly, the frequency-space domain data volume is transformed into the frequency-space wave number domain data volume.
And S304, determining the wave number spectrum components of all time frequency slices in the frequency space wave number domain data body by the slave node, and storing the wave number spectrum components in an output matrix.
In an embodiment of the present application, in the frequency-space wavenumber domain, the wavenumber spectral component of each time frequency slice in the frequency-space wavenumber domain data volume can be obtained by calculation according to the frequency-space wavenumber domain data volume. For example, a constraint matrix may be obtained by calculation according to the amplitude spectrum of the frequency-space wavenumber domain data volume, and then the wavenumber spectrum component of each time frequency slice is obtained by calculation according to the constraint matrix.
S305, the slave node inversely transforms the output matrix to a time space domain to obtain the seismic data after the data block is reconstructed.
And S306, the slave node provides the reconstructed seismic data to the master node for result combination processing so as to obtain a final reconstruction result. In an embodiment of the present application, the result merging may be merging of information recording files, and it is necessary to count the header information of each node data again, and the physical storage file may not be changed and may be logically represented as one file. After generating the final reconstruction result, the master node may write the final reconstruction result into the database, for example, as shown in fig. 2.
In one embodiment of the present application, a task request and data read control thread, a computation control thread, and a result output control thread may be established within the slave node. Thus, the processing of the various nodes is parallel; in each node, the task request and data reading, calculation and output are parallel, so that the data reconstruction time is greatly shortened, and the data reconstruction efficiency is improved. In addition, when the slave node is a multi-core node, each CPU core in the slave node establishes a task request and data reading control thread, a computation control thread, and a result output control thread, so that resources of the node can be further fully utilized for parallel processing. In an exemplary embodiment, a thread may be created using a pthread method, for example, and thread control may be performed using a pthread _ mutex _ t function and a pthread _ cond _ t function.
In an embodiment of the present application, during job execution, when job execution fails due to various machine reasons and the like, the job needs to be restarted. And after the operation is restarted, the main node can rebuild the task list according to the information recording file and check the correctness of the processing sub-result. In this way, the data reconstruction of the present application can have a failure recovery function.
Therefore, the seismic data are distributed to different nodes to be processed in parallel, the seismic data reconstruction time is greatly shortened on the basis of improving the imaging quality of the seismic data, and the seismic data processing method has the node fault tolerance and fault recovery functions, and meets the requirement of seismic exploration data processing industrialization in a real sense.
Referring to fig. 4, in a logical structure, each slave node may include:
a task obtaining module 41, configured to obtain a task, where the task corresponds to one of a plurality of data blocks obtained by splitting a seismic data volume of a time-space domain to be reconstructed by a master node in advance;
a first transform module 42, configured to perform fast fourier transform on the data block along a time direction to obtain a frequency-space domain data volume;
a second transform module 43, configured to perform non-uniform fourier transform on the frequency-space domain data volume according to a spatial position along a spatial direction, so as to obtain a frequency-space wavenumber domain data volume;
a data storage module 44, configured to determine wave number spectrum components of all time frequency slices in the frequency-space wave number domain data volume, and store the wave number spectrum components in an output matrix;
a third transformation module 45, configured to inverse transform the output matrix to a time-space domain to obtain the seismic data after the data block is reconstructed;
a result processing module 46 may be configured to provide the reconstructed seismic data to the master node for result merging processing.
Referring to fig. 5, a data reconstruction apparatus for improving spatial sampling properties of seismic data according to an embodiment of the present application may include a memory, a processor, and a computer program stored on the memory, and when executed by the processor, the computer program performs the following steps:
acquiring a task, wherein the task corresponds to one of a plurality of data blocks which are obtained by splitting a seismic data body of a time space domain to be reconstructed by a main node in advance;
carrying out fast Fourier transform on the data block along the time direction to obtain a frequency-space domain data volume;
carrying out non-uniform Fourier transform on the frequency-space domain data volume according to the space position along the space direction to obtain a frequency-space wave number domain data volume;
determining wave number spectrum components of all time frequency slices in the frequency space wave number domain data body, and storing the wave number spectrum components in an output matrix;
inversely transforming the output matrix to a time space domain to obtain the seismic data after the data block is reconstructed;
and providing the reconstructed seismic data to the main node for result combination processing.
While the process flows described above include operations that occur in a particular order, it should be appreciated that the processes may include more or less operations that are performed sequentially or in parallel (e.g., using parallel processors or a multi-threaded environment).
For convenience of description, the above devices are described as being divided into various units by function, and are described separately. Of course, the functionality of the units may be implemented in one or more software and/or hardware when implementing the present application.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In a typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include forms of volatile memory in a computer readable medium, Random Access Memory (RAM) and/or non-volatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). Memory is an example of a computer-readable medium.
Computer-readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), Static Random Access Memory (SRAM), Dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), Read Only Memory (ROM), Electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), Digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, a computer readable medium does not include a transitory computer readable medium such as a modulated data signal and a carrier wave.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other like elements in a process, method or apparatus that comprises the element.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The application may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. The application may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the system embodiment, since it is substantially similar to the method embodiment, the description is simple, and for the relevant points, reference may be made to the partial description of the method embodiment.
The above description is only an example of the present application and is not intended to limit the present application. Various modifications and changes may occur to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the scope of the claims of the present application.
Claims (8)
1. A data reconstruction method for improving spatial sampling properties of seismic data, comprising:
acquiring a task from a slave node, wherein the task corresponds to one of a plurality of data blocks which are obtained by splitting a seismic data body of a time-space domain to be reconstructed by a master node in advance;
the slave node performs fast Fourier transform on the data block along the time direction to obtain a frequency-space domain data body;
the slave node performs non-uniform Fourier transform on the frequency-space domain data volume according to a space position along a space direction to obtain a frequency-space wave number domain data volume;
the slave node determines wave number spectrum components of all time frequency slices in the frequency space wave number domain data body and stores the wave number spectrum components into an output matrix;
the slave node inversely transforms the output matrix to a time space domain to obtain the seismic data after the data block is reconstructed;
the slave node provides the reconstructed seismic data to the master node for result combination processing; the number of the split data blocks is the product of the number of spatial windows in three specified directions in a data reconstruction range;
the three specified directions include:
CMP Line direction, CMP direction, and Offset-x direction; alternatively, the first and second electrodes may be,
the CMP Line direction, the CMP direction, and the Offset bin direction.
2. The method for reconstructing data to improve spatial sampling properties of seismic data according to claim 1, wherein a task request and data read control thread, a computation control thread and a result output control thread are established from within the nodes.
3. The method for reconstructing data to improve spatial sampling properties of seismic data according to claim 1, wherein when the slave node is a multi-core node, each CPU core in the slave node is established with a task request and data reading control thread, a computation control thread and a result output control thread.
4. The method of data reconstruction for improving spatial sampling properties of seismic data as claimed in claim 1 wherein said acquiring tasks from nodes comprises:
the slave node requests the task from the master node, transmits the last completed task ID and the processing sub-result to the master node, and the master node saves the last completed task ID and the processing sub-result as an information recording file;
the slave node receives the tasks distributed by the master node according to the task list; the task list includes data block IDs and the three specified directions.
5. The method of data reconstruction for improving spatial sampling properties of seismic data according to claim 4, wherein said information log file doubles as a proof of failure recovery; when the slave node fails, the master node distributes tasks which are not completed by the slave node to other slave nodes in the system for processing according to the information record file.
6. The method for data reconstruction that improves the spatial sampling properties of seismic data as set forth in claim 4, further including:
and after the operation is restarted, the main node reconstructs the task list according to the information recording file and checks the correctness of the processing sub-result.
7. A data reconstruction system for improving spatial sampling properties of seismic data, comprising a master node and a plurality of slave nodes, each of said slave nodes comprising:
the task acquisition module is used for acquiring a task, wherein the task corresponds to one of a plurality of data blocks which are obtained by splitting a seismic data body of a time space domain to be reconstructed by a main node in advance;
the first transformation module is used for performing fast Fourier transformation on the data block along the time direction to obtain a frequency-space domain data volume;
the second transformation module is used for carrying out non-uniform Fourier transformation on the frequency-space domain data volume according to the space position along the space direction to obtain a frequency-space wave number domain data volume;
the data storage module is used for determining wave number spectrum components of all time frequency slices in the frequency space wave number domain data body and storing the wave number spectrum components into an output matrix;
the third transformation module is used for inversely transforming the output matrix to a time space domain to obtain the seismic data after the data block is reconstructed;
the result processing module is used for providing the reconstructed seismic data to the main node so as to carry out result combination processing; the number of the split data blocks is the product of the number of spatial windows in three specified directions in a data reconstruction range;
the three specified directions include:
CMP Line direction, CMP direction, and Offset-x direction; alternatively, the first and second electrodes may be,
the CMP Line direction, the CMP direction, and the Offset bin direction.
8. A data reconstruction apparatus for improving spatial sampling properties of seismic data, comprising a memory, a processor, and a computer program stored on the memory, wherein the computer program when executed by the processor performs the steps of:
acquiring a task, wherein the task corresponds to one of a plurality of data blocks which are obtained by splitting a seismic data body of a time space domain to be reconstructed by a main node in advance;
carrying out fast Fourier transform on the data block along the time direction to obtain a frequency-space domain data volume;
carrying out non-uniform Fourier transform on the frequency-space domain data volume according to the space position along the space direction to obtain a frequency-space wave number domain data volume;
determining wave number spectrum components of all time frequency slices in the frequency space wave number domain data body, and storing the wave number spectrum components in an output matrix;
inversely transforming the output matrix to a time space domain to obtain the seismic data after the data block is reconstructed;
providing the reconstructed seismic data to the master node for result merging processing; wherein the content of the first and second substances,
the number of the split data blocks is the product of the number of spatial windows in three specified directions in the data reconstruction range;
the three specified directions include:
CMP Line direction, CMP direction, and Offset-x direction; alternatively, the first and second electrodes may be,
the CMP Line direction, the CMP direction, and the Offset bin direction.
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 CN109100804A (en) | 2018-12-28 |
| CN109100804B true 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) |
Families Citing this family (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 |
| CN115659579B (en) * | 2022-09-05 | 2023-07-04 | 中南大学 | Three-dimensional gravity field numerical simulation method and system based on 3D AS-FT |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9632193B2 (en) * | 2013-11-01 | 2017-04-25 | Conocophillips Company | Compressive sensing |
| CN104459781B (en) * | 2014-12-09 | 2017-09-01 | 中国石油天然气集团公司 | The random noise attenuation method of three-dimensional earthquake data before superposition |
| 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 |
| CN106371138B (en) * | 2016-08-17 | 2018-10-16 | 中国石油天然气集团公司 | Reconstruction of seismic data method and apparatus |
-
2018
- 2018-06-11 CN CN201810592971.8A patent/CN109100804B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN109100804A (en) | 2018-12-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ekanayake et al. | Mapreduce for data intensive scientific analyses | |
| Ji et al. | ispan: Parallel identification of strongly connected components with spanning trees | |
| CN109471171B (en) | Method, device and system for separating aliasing seismic data | |
| CN109100804B (en) | Data reconstruction method, device and system for improving seismic data space sampling attribute | |
| CN104123288A (en) | Method and device for inquiring data | |
| CN104423961B (en) | A kind of method and system for generating test script | |
| US9535948B1 (en) | Data access statement translation | |
| CN108268468B (en) | Big data analysis method and system | |
| CN106353816B (en) | A kind of earthquake-capturing footprint Noise Elimination method and system | |
| CN107561588B (en) | A kind of seismic data noise drawing method and device | |
| CN110825764B (en) | SQL script generation method, system, storage medium and processor | |
| CN112230284A (en) | Parallel random noise attenuation method, monitoring method and node for three-dimensional pre-stack data | |
| CN108304186B (en) | Method and equipment for executing multiple operations based on synthetic configuration file | |
| CN114968028A (en) | Method, apparatus, and medium for dynamically changing menu layout based on user behavior analysis | |
| CN108599773B (en) | Vibration signal data compression acquisition method based on deterministic measurement matrix | |
| CN112308317A (en) | Massive seismic observation data noise power spectrum calculation method and system based on distributed architecture | |
| CN112818367B (en) | File encryption method, device, storage medium and processor | |
| CN110990378A (en) | Block chain-based data consistency comparison method, device and medium | |
| CN108958732A (en) | A kind of data load method and equipment based on PHP | |
| CN117033526B (en) | Data storage method, data query method, device, equipment and storage medium | |
| CN109558221B (en) | Task processing method and device | |
| CN118035089A (en) | Code coverage rate calculating method and device based on multithreading | |
| CN114969200B (en) | Data synchronization method, device, electronic equipment and storage medium | |
| US20220171696A1 (en) | Techniques and Architectures for Providing an Isolated Runtime Context in a Shared Environment | |
| CN117092698A (en) | Seismic data reconstruction method, processor and storage medium based on binary tree structure |
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 |