Method, apparatus and storage medium for determining continuous pulse seismic record data

Abstract

Embodiments of the present application provide a continuous pulse seismic record data determination method, device, equipment and storage medium, the method comprises: obtaining the seismic wave of the target work area continuously excited by the pulse source in the preset excitation interval duration, wherein the preset excitation interval duration is the time difference between two continuous excitation events;According to the preset excitation interval duration, the continuously excited seismic wave is separated to obtain the independent seismic response generated by each excitation;The independent seismic response generated by each excitation is time corrected;The independent seismic response generated by each excitation after time correction is weighted and superimposed to generate the synthetic seismic data of the target work area continuously excited, more high-quality seismic data can be obtained in the same time, the acquisition efficiency is significantly improved, which is suitable for complex terrain and high-density exploration, provides high-precision and efficient data support for seismic exploration engineering, and has important economic benefits and broad application prospect.

Description

Technical Field

[0001] This application relates to the field of petroleum seismic exploration technology, specifically to a method, apparatus, equipment, and storage medium for determining continuous pulse seismic record data. Background Technology

[0002] In seismic exploration, traditional data acquisition often employs intermittent excitation and recording, i.e., excitation followed by waiting for wave decay before re-excitation. While widely used, this method is inefficient, especially in complex terrain requiring high-density data acquisition. Intermittent excitation may not fully utilize the source energy, and the time interval between each excitation can affect data quality and continuity.

[0003] In related technologies, the industry has explored more efficient pulse seismic recording techniques. However, this brings the challenge of effectively separating and synthesizing single-shot data, which has become a key aspect of technological development. Although traditional methods have achieved significant success in seismic exploration, their limitations in efficiency and signal continuity have spurred the development of pulse technology. It is hoped that the research and application of new technologies will further promote the advancement of seismic exploration technology. Summary of the Invention

[0004] This application provides a method, apparatus, device, and storage medium for determining continuous pulse seismic record data.

[0005] The first aspect of this application provides a method for determining continuous pulse seismic record data, including:

[0006] The seismic waves continuously excited by the pulse source on the target work area within a preset excitation interval are obtained, wherein the preset excitation interval is the time difference between two consecutive excitation events;

[0007] Based on the preset excitation interval, the continuously excited seismic waves are separated to obtain the independent seismic response generated by each excitation;

[0008] Time correction is performed on the independent seismic response generated by each excitation;

[0009] The independent seismic responses generated by each excitation after time correction are weighted and superimposed to generate synthetic seismic data for continuous excitation of the target work area.

[0010] In an optional embodiment of this application, the continuously excited seismic waves are separated according to a preset excitation interval to obtain the independent seismic response generated by each excitation, including:

[0011] Based on the preset excitation interval, the independent seismic response generated by each excitation is obtained according to the difference between the seismic waves of two adjacent excitations.

[0012] In an optional embodiment of this application, the independent seismic response generated by each excitation is obtained based on the difference between the seismic waves of two adjacent excitations, using the following expression:

[0013] D(t)=R(t)-R(t-Δt)

[0014] Where Δt is the preset excitation interval, D(t) is the difference between two adjacent excitations of seismic waves, R(t) is the amplitude of the excitation of seismic waves at time t, and R(t-Δt) is the amplitude of the excitation of seismic waves at time t-Δt.

[0015] In an optional embodiment of this application, time correction is performed on the independent seismic response generated by each excitation, including:

[0016] The time offset of the i-th excitation relative to the first shot is determined based on the start time recorded by the first shot and the actual start time of the independent seismic response generated by the i-th excitation.

[0017] Based on the time offset, the independent seismic response generated by each excitation is time-corrected to obtain the adjusted data start time;

[0018] Compare whether the independent seismic response waveforms before and after the data start time adjustment are consistent under each excitation.

[0019] If the independent seismic response waveforms before and after the data start time adjustment are inconsistent under each excitation, the time offset is adjusted until the independent seismic response waveforms before and after the data start time adjustment are consistent under each excitation.

[0020] In an optional embodiment of this application, the following expression is used to compare whether the independent seismic response waveforms before and after the data start time adjustment are consistent under each excitation:

[0021] C(τ)=∫R(t)R(t+τ)dt

[0022] Where C(τ) is the cross-correlation function, R(t) is the independent seismic response waveform before the data start time adjustment under each excitation, R(t+τ) is the independent seismic response waveform after the data start time adjustment under each excitation, and τ is the time delay. The time delay corresponding to the peak value of the cross-correlation function is close to 0, which determines that the independent seismic response waveforms before and after the data start time adjustment under each excitation are consistent.

[0023] In an optional embodiment of this application, the independent seismic responses generated by each time-corrected excitation are weighted and superimposed to generate synthetic seismic data continuously excited for the target work area, including:

[0024] Squaring the independent seismic responses generated by each time-corrected excitation;

[0025] According to pre-set weights, the squared data points are weighted and averaged to generate synthetic seismic data continuously excited for the target work area.

[0026] The pre-set weights are determined based on the signal-to-noise ratio (SNR) of the seismic waves generated in each excitation. The SNR of the seismic waves generated in each excitation is determined using the following expression:

[0027]

[0028]

[0029] Where d is depth or time, R i For the seismic wave generated by the i-th excitation, SNR(d, R) i () represents the adjusted signal-to-noise ratio corresponding to the seismic wave generated by the i-th excitation. Let be the unadjusted signal-to-noise ratio of the seismic wave generated by the i-th excitation. It is a depth- or time-dependent adjustment coefficient. It is a coefficient adjustment function, SNR local This represents the local signal-to-noise ratio.

[0030] In an optional embodiment of this application, the independent seismic response generated by each time-corrected excitation is squared using the following expression:

[0031]

[0032] Among them, W tj Let Yj be the weighting coefficient for the j-th channel at time t, L be the length of the time window considered, Ysj be the amplitude value of the j-th channel at time s, and t be the current time point under consideration.

[0033] Using the following expression, the squared data points are weighted and averaged according to pre-defined weights to generate synthetic seismic data continuously excited for the target work area:

[0034]

[0035] Among them, X t For the synthetic seismic data continuously excited at time t for the target work area, x tj Let N be the amplitude value of the j-th channel at time t, and N be the total number of channels.

[0036] A second aspect of this application provides a device for determining continuous pulse seismic record data, comprising:

[0037] The acquisition module is used to acquire the seismic waves continuously excited by the pulse source on the target work area within a preset excitation interval, wherein the preset excitation interval is the time difference between two consecutive excitation events;

[0038] The separation module is used to separate continuously excited seismic waves according to a preset excitation interval to obtain the independent seismic response generated by each excitation.

[0039] The correction module is used to perform time correction on the independent seismic response generated by each excitation;

[0040] The overlay module is used to weight and overlay the independent seismic responses generated by each excitation after time correction, to generate synthetic seismic data for continuous excitation of the target work area.

[0041] A third aspect of this application provides a computer device, including a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the steps of any of the above methods for determining continuous pulse seismic record data.

[0042] A fourth aspect of the present application provides a computer-readable storage medium having a computer program stored thereon, characterized in that the computer program, when executed by a processor, implements the steps of the continuous pulse seismic record data determination method as described in any of the preceding claims.

[0043] Compared with the prior art, the technical solutions provided in this application have at least some or all of the following advantages:

[0044] The continuous pulse seismic record data determination method described in this application acquires seismic waves continuously excited by a pulse source on a target work area within a preset excitation interval, wherein the preset excitation interval is the time difference between two consecutive excitation events; based on the preset excitation interval, the continuously excited seismic waves are separated to obtain the independent seismic response generated by each excitation; time correction is performed on the independent seismic response generated by each excitation; weighted and superimposed processing is performed on the time-corrected independent seismic responses generated by each excitation to generate synthetic seismic data continuously excited on the target work area. This method can acquire more high-quality seismic data in the same time, significantly improving acquisition efficiency. It is suitable for complex terrain and high-density exploration, providing high-precision and efficient data support for seismic exploration engineering, and has significant economic benefits and broad application prospects. Attached Figure Description

[0045] The accompanying drawings, which are included to provide a further understanding of this application and form part of this application, illustrate exemplary embodiments and are used to explain this application, but do not constitute an undue limitation of this application. In the drawings:

[0046] Figure 1 A flowchart illustrating a method for determining continuous pulse seismic record data according to an embodiment of this application;

[0047] Figure 2 A flowchart illustrating a method for determining continuous pulse seismic record data according to another embodiment of this application;

[0048] Figure 3 One embodiment of this application provides a method for determining continuous pulse seismic record data by forward modeling shot records;

[0049] Figure 4 The continuous pulse seismic record data determination method provided in one embodiment of this application is a forward modeling simulation of the seismic record of the second shot separated from the blasting record;

[0050] Figure 5 The continuous pulse seismic record data determination method provided in one embodiment of this application is a forward modeling simulation of the seismic record of the third shot separated from the blasting record;

[0051] Figure 6 The continuous pulse seismic record data determination method provided in one embodiment of this application is a forward modeling simulation of shot records, which combines the three separated shots into a composite seismic record.

[0052] Figure 7 A schematic diagram of a continuous pulse seismic record data determination device provided in one embodiment of this application;

[0053] Figure 8 This is a schematic diagram of a computer device structure provided in one embodiment of this application. Detailed Implementation

[0054] To make the technical solutions and advantages of the embodiments of this application clearer, the exemplary embodiments of this application will be described in further detail below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of this application, and not an exhaustive list of all embodiments. It should be noted that, unless otherwise specified, the embodiments and features in the embodiments of this application can be combined with each other.

[0055] Please see Figure 1 The method for determining continuous pulse seismic record data provided in this application includes the following steps 100 to 400:

[0056] Step 100: Obtain the seismic waves continuously excited by the pulse source on the target work area within a preset excitation interval, wherein the preset excitation interval is the time difference between two consecutive excitation events;

[0057] Step 200: According to the preset excitation interval, the continuously excited seismic waves are separated to obtain the independent seismic response generated by each excitation;

[0058] Step 300: Perform time correction on the independent seismic response generated by each excitation;

[0059] Step 400: Weight and superimpose the independent seismic responses generated by each excitation after time correction to generate synthetic seismic data continuously excited for the target work area.

[0060] In an optional embodiment of this application, step 100, acquiring the seismic waves continuously excited by the pulse source on the target work area within a preset excitation interval, includes:

[0061] A pulse source is used to continuously generate seismic waves within a predetermined excitation interval, and a multichannel seismograph is used to continuously record the seismic waves generated by the pulse source on the target work area within a predetermined excitation interval. The pulse source includes, but is not limited to, an electromagnetic pulse source, a water jet source, and an air jet source. The excitation parameters include excitation energy and repetition frequency. The multichannel seismograph recording parameters involve sampling rate and filter settings. By exciting with a pulse source and recording with a seismograph, an aliased pulse seismic record is obtained, which serves as the seismic wave continuously generated by the pulse source on the target work area within a predetermined excitation interval.

[0062] In an optional embodiment of this application, step 100 involves the seismic instrument continuously recording seismic waves excited by the pulse source within a preset interval, including the following steps:

[0063] First, a pulse source continuously excites seismic waves at preset intervals. The pulse source can be an electromagnetic pulse source, a water jet source, or an air jet source, each generating seismic waves according to its specific physical mechanism. Excitation parameters, such as excitation energy and repetition frequency, are optimized based on the purpose and specific conditions of the geological exploration. Excitation energy determines the penetration depth of the seismic waves, while the repetition frequency affects the speed and efficiency of data acquisition. Continuous excitation by the pulse source generates a series of seismic waves that penetrate underground structures and are received by seismic instruments at the surface or underground. The arrangement and configuration of the seismic instruments are designed to maximize the capture of seismic wave information. To accommodate different exploration targets, seismic instruments are typically deployed in multichannel configurations to cover a wider geographical area and improve the spatial resolution of the data. The seismic instruments continuously record the reflected, refracted, or scattered signals of the seismic waves. Recording parameters, including sampling rate and filter settings, are adjusted according to the exploration target and the expected signal-to-noise ratio. The sampling rate must be high enough to ensure that all frequency components of the seismic signal are captured according to the Nyquist criterion, avoiding aliasing.

[0064] In an optional embodiment of this application, step 200 involves separating the continuously excited seismic waves according to a preset excitation interval to obtain the independent seismic response generated by each excitation, including:

[0065] Based on the preset excitation interval, the independent seismic response generated by each excitation is obtained according to the difference between the seismic waves of two adjacent excitations.

[0066] In an optional embodiment of this application, step 200 involves identifying the specific time points of each excitation and segmenting the continuously recorded data into separate segments based on these time points. Each segment represents a record between one excitation event and the next. By applying a difference operation to these segments, the influence of the previous excitation on the record of the next excitation is eliminated, thereby achieving the separation of aliased data.

[0067] In an optional embodiment of this application, the independent seismic response generated by each excitation is obtained based on the difference between the seismic waves of two adjacent excitations, using the following expression:

[0068] D(t)=R(t)-R(t-Δt)

[0069] Where Δt is the preset excitation interval, D(t) is the difference between two adjacent excitations of seismic waves, R(t) is the amplitude of the excitation of seismic waves at time t, and R(t-Δt) is the amplitude of the excitation of seismic waves at time t-Δt.

[0070] In an optional embodiment of this application, a differential method is used to process the shot gather data, thereby eliminating aliasing interference and separating individual shot gather data. The differential method is based on the time difference between two adjacent consecutive excitation events, and cancels or reduces mutual interference by calculating the difference between the records of adjacent excitation events.

[0071] In an optional embodiment of this application, step 300, performing time correction on the independent seismic response generated by each excitation, includes:

[0072] The time offset of the i-th excitation relative to the first shot is determined based on the start time recorded by the first shot and the actual start time of the independent seismic response generated by the i-th excitation.

[0073] Based on the time offset, the independent seismic response generated by each excitation is time-corrected to obtain the adjusted data start time;

[0074] Compare whether the independent seismic response waveforms before and after the data start time adjustment are consistent under each excitation.

[0075] If the independent seismic response waveforms before and after the data start time adjustment are inconsistent under each excitation, the time offset is adjusted until the independent seismic response waveforms before and after the data start time adjustment are consistent under each excitation.

[0076] In an optional embodiment of this application, in step 300, the exact start time of the first shot in the continuous pulse recording is determined and used as a reference point for aligning all other individual shot records. Subsequently, the start time of each separated individual shot record is adjusted accordingly by calculating the time offset of each individual shot record relative to the first shot record, so that the start time of all individual shot records is consistent with that of the first shot.

[0077] In an optional embodiment of this application, the following expression is used to compare whether the independent seismic response waveforms before and after the data start time adjustment are consistent under each excitation:

[0078] C(τ)=∫R(t)R(t+τ)dt

[0079] Where C(τ) is the cross-correlation function, R(t) is the independent seismic response waveform before the data start time adjustment under each excitation, R(t+τ) is the independent seismic response waveform after the data start time adjustment under each excitation, and τ is the time delay. The time delay corresponding to the peak value of the cross-correlation function is close to 0, which determines that the independent seismic response waveforms before and after the data start time adjustment under each excitation are consistent.

[0080] In this application, the time offset of the i-th firing relative to the first firing is determined based on the start time of the first firing record and the actual start time of the single firing record generated by the i-th firing. Based on the time offset, time correction is performed on each single firing record, and the adjusted data start time Ti′ can be obtained using the following formula: T i ′=T0+ΔT i In practice, time correction may be achieved through interpolation or resampling to maintain data continuity and integrity. Assuming the sampling interval of the original data is δt, the interpolated time point can be expressed by the following formula: t new =t old +δt, where t old The time point before correction, t newThe corrected time point is used for data verification. The accuracy of the time correction is verified by comparing the waveform consistency before and after correction. This is done by calculating the cross-correlation function of the data before and after correction: C(τ)=∫R(t)R(t+τ)dt, where R(t) is the record before correction, R(t+τ) is the record after correction, and τ is the time delay. Ideally, the time delay corresponding to the peak of the cross-correlation function should be close to 0, indicating a high degree of consistency between the data before and after correction. Through the above steps and methods, it is ensured that the starting time point of all single-shot records is consistent with the first shot in the continuous pulse record.

[0081] In an optional embodiment of this application, step 400 involves weighting and superimposing the independent seismic responses generated by each time-corrected excitation to generate synthetic seismic data continuously excited for the target work area, including:

[0082] Squaring the independent seismic responses generated by each time-corrected excitation;

[0083] According to pre-set weights, the squared data points are weighted and averaged to generate synthetic seismic data continuously excited for the target work area.

[0084] The pre-set weights are determined based on the signal-to-noise ratio (SNR) of the seismic waves generated in each excitation. The SNR of the seismic waves generated in each excitation is determined using the following expression:

[0085]

[0086] Where d is depth or time, R i For the seismic wave generated by the i-th excitation, SNR(d, R) i () represents the adjusted signal-to-noise ratio corresponding to the seismic wave generated by the i-th excitation. Let be the unadjusted signal-to-noise ratio of the seismic wave generated by the i-th excitation. It is a depth- or time-dependent adjustment coefficient. It is a coefficient adjustment function, SNR local This represents the local signal-to-noise ratio.

[0087] In the method for determining continuous pulse seismic record data in this application, in order to adapt to possible continuous excitation situations and optimize the specific depth interval of each excitation, a general signal-to-noise ratio function SNR(d, R) is introduced. i This method dynamically adjusts the weight of each excitation and optimizes the processing effect, with higher signal-to-noise ratio excitations receiving greater weights. Through the design and implementation of this universal signal-to-noise ratio function, it is possible to handle any number of continuous excitation events more flexibly and provide optimal weight adjustments for each depth interval of each excitation, significantly improving the reliability of continuous pulse seismic waveform separation.

[0088] In an optional embodiment of this application, the independent seismic response generated by each time-corrected excitation is squared using the following expression:

[0089]

[0090] Among them, W tj Let Yj be the weighting coefficient for the j-th channel at time t, L be the length of the time window considered, Ysj be the amplitude value of the j-th channel at time s, and t be the current time point under consideration.

[0091] Using the following expression, the squared data points are weighted and averaged according to pre-defined weights to generate synthetic seismic data continuously excited for the target work area:

[0092]

[0093] Among them, X t For the synthetic seismic data continuously excited at time t for the target work area, X tj Let N be the amplitude value of the j-th channel at time t, and N be the total number of channels.

[0094] This application employs square diversity weighting to stack the separated single-shot records, resulting in synthesized data. This step utilizes square diversity weighting techniques from signal processing to improve the signal-to-noise ratio and data quality of the stacked seismic data.

[0095] In an optional embodiment of this application, in step 400, the weighting factor in this process can be dynamically adjusted according to the signal intensity and noise level in the single-shot record to ensure that the weighting process can effectively highlight the useful signal while suppressing noise. After completing the above sum-of-squares weighting processing, all processed single-shot records are superimposed to form the final synthetic seismic data.

[0096] like Figure 2 As shown below, the process of applying the continuous pulse seismic record data determination method of this application to complex terrain includes the following steps:

[0097] The first step, data acquisition in complex terrain, involves using pulse seismic sources (including electromagnetic pulses, water jets, or air jets) to continuously generate seismic waves at preset time intervals. Excitation parameters such as excitation energy and repetition frequency need to be optimized according to the exploration target. Multiple seismic instruments are then deployed to continuously record seismic waves. Recording parameters such as sampling rate and filter settings need to be adjusted based on the signal-to-noise ratio and frequency components to ensure the capture of all frequency signals. Through continuous excitation and recording, aliased pulse seismic record data is obtained, such as... Figure 3 The forward modeling record shown is as follows. Multiple seismic waves are severely superimposed on the record, making it difficult to distinguish individual responses.

[0098] The second step, data separation for complex terrain, involves applying a differential method to the aliased data based on a preset excitation interval Δt. This process eliminates mutual interference by calculating the difference between adjacent excitation records, thus separating the independent seismic responses for each event. Figure 4 and Figure 5 As shown, the separated records of the second and third shots are displayed respectively. Furthermore, a universal signal-to-noise ratio function SNR(d,R) is introduced. i The function dynamically adjusts the weights of each excitation record to optimize the processing effect in different depth ranges. This design allows for flexible handling of any number of excitations, providing optimal weights for each depth in each excitation.

[0099] The third step involves data correction and overlay for complex terrain. This includes time correction of each isolated single-shot record to ensure that all records are consistent with the start time of the first shot. This involves steps such as time offset calculation, interpolation correction, and verification of correction results.

[0100] The fourth step involves using square diversity weighting to weight and superimpose the time-corrected single-shot records for synthesis.

[0101] a) Perform a square operation on each single-shot recorded data point;

[0102] b) The squared values ​​are weighted averaged according to the preset weight allocation scheme, wherein the weights can be dynamically adjusted according to the signal-to-noise ratio;

[0103] c) Overlay all weighted individual shot records to obtain the final synthetic seismic data, such as Figure 6 As shown, the signal-to-noise ratio and signal continuity of the final synthesized seismic data were significantly improved.

[0104] The continuous pulse seismic record data determination method of this application significantly improves the efficiency and quality of seismic data acquisition. In particular, in the application of the continuous pulse seismic record data separation and synthesis method, by continuously excitation within a preset interval, not only is the acquisition speed of seismic waves accelerated, but the differential method is also used to effectively separate accurate single-shot data from the aliased seismic records. This process significantly improves the efficiency of data processing. At the same time, precise time correction and square diversity weighted superposition processing further optimize the data quality, improve the signal-to-noise ratio of the seismic data, and make the final synthesized seismic data more accurate in interpretation and reduce the impact of noise.

[0105] The method for determining continuous pulse seismic record data in this application is highly adaptable and can be effectively applied to electromagnetic pulse sources, water gun sources, and air gun sources. It provides reliable data support for the exploration of complex geological structures. Compared with traditional intermittent excitation and recording methods, this application has significant improvements in data acquisition efficiency, data quality, and adaptability, providing a new solution for the development of seismic exploration technology.

[0106] It should be understood that although the steps in the flowchart are shown sequentially according to the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless explicitly stated herein, there is no strict order constraint on the execution of these steps, and they can be executed in other orders. Moreover, at least some steps in the diagram may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily completed at the same time, but can be executed at different times. The execution order of these sub-steps or stages is not necessarily sequential, but can be performed alternately or in turn with other steps or at least some of the sub-steps or stages of other steps.

[0107] Please see Figure 7 One embodiment of this application provides a continuous pulse seismic record data determination device 700, comprising:

[0108] The acquisition module 710 is used to acquire the seismic waves continuously excited by the pulse source on the target work area within a preset excitation interval, wherein the preset excitation interval is the time difference between two consecutive excitation events.

[0109] The separation module 720 is used to separate continuously excited seismic waves according to a preset excitation interval to obtain the independent seismic response generated by each excitation.

[0110] The correction module 730 is used to perform time correction on the independent seismic response generated by each excitation;

[0111] The overlay module 740 is used to weight and overlay the independent seismic responses generated by each excitation after time correction, and generate synthetic seismic data for continuous excitation of the target work area.

[0112] The continuous pulse seismic record data determination device of this application uses a pulse source to continuously excite seismic waves within a preset interval and uses a multi-channel seismograph for continuous recording. According to the preset excitation interval, the aliased data is separated using the differential method to obtain the independent response of each excitation. Time correction is performed on the separated single-shot records to make the time baseline of all records consistent. The corrected records are then subjected to square-weighted superposition processing to synthesize the final seismic data. This device can solve the problems of low efficiency and insufficient signal continuity in traditional intermittent seismic data acquisition methods, reduce the time loss caused by waiting for seismic wave attenuation, and achieve the goal of collecting more seismic data in the same amount of time.

[0113] The continuous pulse seismic record data determination device of this application can effectively separate single-shot data from continuously aliased seismic records, improve the continuity and consistency of data, and ensure high-quality seismic signals for subsequent interpretation. Especially for complex terrain or situations requiring high-density acquisition, it overcomes the limitations of traditional methods in these cases and enhances the signal-to-noise ratio and interpretation accuracy of the final data by improving the acquisition and processing methods of seismic data.

[0114] For specific limitations regarding the aforementioned device 700, please refer to the limitations on the pseudo-random sequence generation method described above, which will not be repeated here. Each module in the aforementioned device 700 can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in hardware or independently of the processor in a computer device, or stored in software in the memory of a computer device, so that the processor can call and execute the operations corresponding to each module.

[0115] In one embodiment, a computer device is provided, the internal structure of which can be as follows: Figure 8 As shown. The computer device includes a processor, memory, network interface, and database connected via a system bus. The processor provides computing and control capabilities. The memory includes a non-volatile storage medium and internal memory. The non-volatile storage medium stores the operating system, computer programs, and the database. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage medium. The database stores data. The network interface communicates with external terminals via a network connection. When the computer program is executed by the processor, it implements the pseudo-random sequence generation method described above. It includes: memory and a processor; the memory stores the computer program; and the processor executes the computer program to implement any step of the pseudo-random sequence generation method described above.

[0116] In one embodiment, a computer-readable storage medium is provided having a computer program stored thereon, which, when executed by a processor, can perform any of the steps in the pseudo-random sequence generation method described above.

[0117] Those skilled in the art will understand that embodiments of this application can be provided as methods, systems, or computer program products. Therefore, this application can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, this application can 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, etc.) containing computer-usable program code.

[0118] This application is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of this application. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can 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, generate instructions for implementing the flowchart... Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0119] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0120] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0121] Although preferred embodiments of this application have been described, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments as well as all changes and modifications falling within the scope of this application.

[0122] Obviously, those skilled in the art can make various modifications and variations to this application without departing from the spirit and scope of this application. Therefore, if such modifications and variations fall within the scope of the claims of this application and their equivalents, this application also intends to include such modifications and variations.

Claims

1. A method for determining continuous pulse seismic record data, characterized in that, include: The seismic waves continuously excited by the pulse source on the target work area within a preset excitation interval are obtained, wherein the preset excitation interval is the time difference between two consecutive excitation events; Based on the preset excitation interval, the continuously excited seismic waves are separated to obtain the independent seismic response generated by each excitation; Time correction is performed on the independent seismic response generated by each excitation; The independent seismic responses generated by each excitation after time correction are weighted and superimposed to generate synthetic seismic data for continuous excitation of the target work area.

2. The method according to claim 1, characterized in that, Based on a preset excitation interval, the continuously excited seismic waves are separated to obtain the independent seismic response generated by each excitation, including: Based on the preset excitation interval, the independent seismic response generated by each excitation is obtained according to the difference between the seismic waves of two adjacent excitations.

3. The method according to claim 2, characterized in that, The independent seismic response generated by each excitation is obtained using the following expression, based on a preset excitation interval and the difference between the seismic waves from two adjacent excitations: D(t)=R(t)-R(t-Δt) Where Δt is the preset excitation interval, D(t) is the difference between two adjacent excitations of seismic waves, R(t) is the amplitude of the excitation of seismic waves at time t, and R(t-Δt) is the amplitude of the excitation of seismic waves at time t-Δt.

4. The method according to claim 1, characterized in that, Time correction is performed on the independent seismic response generated by each excitation, including: The time offset of the i-th excitation relative to the first shot is determined based on the start time recorded by the first shot and the actual start time of the independent seismic response generated by the i-th excitation. Based on the time offset, the independent seismic response generated by each excitation is time-corrected to obtain the adjusted data start time; Compare whether the independent seismic response waveforms before and after the data start time adjustment are consistent under each excitation. If the independent seismic response waveforms before and after the data start time adjustment are inconsistent under each excitation, the time offset is adjusted until the independent seismic response waveforms before and after the data start time adjustment are consistent under each excitation.

5. The method according to claim 4, characterized in that, The following expression is used to compare whether the independent seismic response waveforms before and after data start time adjustment are consistent under each excitation: C(τ)=∫R(t)R(t+τ)dt Where C(τ) is the cross-correlation function, R(t) is the independent seismic response waveform before the data start time adjustment under each excitation, R(t+τ) is the independent seismic response waveform after the data start time adjustment under each excitation, and τ is the time delay. Where the time delay corresponding to the peak value of the cross-correlation function is close to 0, it is determined that the independent seismic response waveforms before and after the data start time adjustment under each excitation are consistent.

6. The method according to claim 1, characterized in that, The independent seismic responses generated by each time-corrected excitation are weighted and superimposed to generate synthetic seismic data for the target work area, including: Squaring the independent seismic responses generated by each time-corrected excitation; According to pre-set weights, the squared data points are weighted and averaged to generate synthetic seismic data continuously excited for the target work area. The pre-set weights are determined based on the signal-to-noise ratio (SNR) of the seismic waves generated in each excitation. The SNR of the seismic waves generated in each excitation is determined using the following expression: Where d is depth or time, R i For the seismic wave generated by the i-th excitation, SNR(d, R) i () represents the adjusted signal-to-noise ratio corresponding to the seismic wave generated by the i-th excitation. Let be the unadjusted signal-to-noise ratio of the seismic wave generated by the i-th excitation. It is a depth- or time-dependent adjustment coefficient. It is a coefficient adjustment function, SNR local This represents the local signal-to-noise ratio.

7. The method according to claim 6, characterized in that, The independent seismic responses generated by each time-corrected excitation are squared using the following expression: Among them, W tj Let Yj be the weighting coefficient for the j-th channel at time t, L be the length of the time window considered, Ysj be the amplitude value of the j-th channel at time s, and t be the current time point under consideration. Using the following expression, the squared data points are weighted and averaged according to pre-defined weights to generate synthetic seismic data continuously excited for the target work area: Among them, X t For the synthetic seismic data continuously excited at time t for the target work area, x tj Let N be the amplitude value of the j-th channel at time t, and N be the total number of channels.

8. A device for determining continuous pulse seismic record data, characterized in that, include: The acquisition module is used to acquire the seismic waves continuously excited by the pulse source on the target work area within a preset excitation interval, wherein the preset excitation interval is the time difference between two consecutive excitation events; The separation module is used to separate continuously excited seismic waves according to a preset excitation interval to obtain the independent seismic response generated by each excitation. The correction module is used to perform time correction on the independent seismic response generated by each excitation; The overlay module is used to weight and overlay the independent seismic responses generated by each excitation after time correction, to generate synthetic seismic data for continuous excitation of the target work area.

9. A computer device, comprising: A memory and a processor, the memory storing a computer program, characterized in that the processor, when executing the computer program, implements the steps of the continuous pulse seismic record data determination method according to any one of claims 1 to 7.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the steps of the continuous pulse seismic record data determination method according to any one of claims 1 to 7.

Images