Velocity analysis method, device, electronic equipment and medium for pseudo-surface seismic data

Abstract

The present application discloses a velocity analysis method, device, electronic device and medium for pseudo-surface seismic data. The method may include: constructing pseudo-surface seismic data based on VSP downlink multiple wave data; constructing a common scattering point gather based on the initial seismic wave velocity and the pseudo-surface seismic data; performing velocity analysis on the common scattering point gather to obtain a final velocity analysis result. The present invention improves the velocity analysis accuracy of pseudo-surface seismic data, which is conducive to improving the accuracy of subsequent seismic imaging.

Description

Technical Field

[0001] The present invention relates to the field of seismic exploration, and more specifically, to a velocity analysis method, device, electronic equipment and medium for pseudo-surface seismic data. Background Art

[0002] Seismic interferometry is a seismic data reconstruction technology developed in recent years. Its essence is to construct a virtual source, that is, the cross-correlation of the records of two receiving points, just like the record received at another receiving point when one receiving point stimulates it. It is driven by measured data and does not require a velocity model, nor does it require information on the location of the source and the time of stimulation.

[0003] The current interferometric seismic method is widely used. The reconstruction of the acquisition datum between VSP, surface seismic, single-well seismic and inter-well seismic, as well as the surface wave recovery of seismic coda and the inversion of shear wave velocity have greatly promoted its progress. Among them, the VSP data acquisition datum reconstruction technology based on interferometric seismic is a technology that is widely studied and applied. The interferometric seismic method can reconstruct the multiple wave data of VSP into pseudo-surface seismic data through datum reconstruction, that is, extend the well detection points of VSP to the ground to form pseudo-surface seismic data, and then image it through the processing methods and processes of ground seismic data. This imaging technology effectively utilizes multiple wave data, expands the imaging range of VSP data, and enhances the application range of VSP methods. However, the reconstructed pseudo-surface seismic data generally has a low signal-to-noise ratio and a complex wave field, which poses a huge challenge to the existing ground processing velocity analysis technology.

[0004] Therefore, it is necessary to develop a velocity analysis method, device, electronic equipment and medium for VSP pseudo-surface seismic data.

[0005] The information disclosed in the background technology section of the present invention is only intended to deepen the understanding of the general background technology of the present invention, and should not be regarded as acknowledging or suggesting in any form that the information constitutes the prior art already known to those skilled in the art. Summary of the invention

[0006] The present invention proposes a velocity analysis method, device, electronic equipment and medium for pseudo-surface seismic data, which can improve the velocity analysis accuracy of pseudo-surface seismic data, thereby facilitating improving the accuracy of subsequent seismic imaging.

[0007] In a first aspect, an embodiment of the present disclosure provides a velocity analysis method for pseudo-surface seismic data, comprising:

[0008] Construct pseudo-surface seismic data based on VSP downlink multiple wave data;

[0009] constructing a common scattering point gather according to the initial seismic wave velocity and the pseudo-surface seismic data;

[0010] Velocity analysis is performed on the common scattering point gathers to obtain a final velocity analysis result.

[0011] Preferably, obtaining the final velocity analysis result includes:

[0012] The velocity obtained by velocity analysis on the common scattering point gather is used as the initial seismic wave velocity, and the common scattering point gather and velocity analysis are repeatedly constructed until the difference between the velocities obtained from the two velocity analyses is less than a set threshold, thereby obtaining a final velocity analysis result.

[0013] Preferably, the VSP downlink multiple wave data is downlink multiple wave data after wave field separation.

[0014] Preferably, constructing pseudo-surface seismic data includes:

[0015] Step 101: Determine the reconstructed pseudo-ground shot point position;

[0016] Step 102: Determine the virtual receiver position;

[0017] Step 103: Select a receiving point of the VSP, correlate the downlink direct wave of the receiving point at the virtual receiver position with the downlink wave multiple wave records of the receiving point at the pseudo-surface shot point position, and superimpose the correlation result on the pseudo-surface seismic trace;

[0018] Step 104: traverse all receiving points of the VSP and repeat step 103;

[0019] Step 105: repeating steps 102-104, traversing all virtual receiving points of the ground shot points to be reconstructed, and finally forming a reconstructed pseudo-ground shot collection;

[0020] Step 106: Repeat steps 101-105 to traverse all reconstructed pseudo-surface shot points to form all pseudo-surface shot collections, namely, pseudo-surface seismic data.

[0021] Preferably, it also includes:

[0022] Preprocessing is performed on the pseudo-surface seismic data, and then common scattering point gathers are constructed.

[0023] Preferably, the pretreatment comprises:

[0024] The linear noise of the pseudo-ground shot gather is attenuated by FK and median filtering;

[0025] Negatively correlated records in the quasi-ground gun set are removed.

[0026] As a specific implementation method of the embodiment of the present disclosure,

[0027] In a second aspect, the embodiment of the present disclosure further provides a velocity analysis device for pseudo-surface seismic data, comprising:

[0028] The pseudo-surface seismic data construction module constructs pseudo-surface seismic data based on VSP downlink multiple wave data;

[0029] A common scattering point gather construction module is used to construct a common scattering point gather according to the initial seismic wave velocity and the pseudo-surface seismic data;

[0030] The velocity analysis module performs velocity analysis on the common scattering point gathers to obtain a final velocity analysis result.

[0031] Preferably, obtaining the final velocity analysis result includes:

[0032] The velocity obtained by velocity analysis on the common scattering point gather is used as the initial seismic wave velocity, and the common scattering point gather and velocity analysis are repeatedly constructed until the difference between the velocities obtained from the two velocity analyses is less than a set threshold, thereby obtaining a final velocity analysis result.

[0033] Preferably, the VSP downlink multiple wave data is downlink multiple wave data after wave field separation.

[0034] Preferably, constructing pseudo-surface seismic data includes:

[0035] Step 101: Determine the reconstructed pseudo-ground shot point position;

[0036] Step 102: Determine the virtual receiver position;

[0037] Step 103: Select a receiving point of the VSP, correlate the downlink direct wave of the receiving point at the virtual receiver position with the downlink wave multiple wave records of the receiving point at the pseudo-surface shot point position, and superimpose the correlation result on the pseudo-surface seismic trace;

[0038] Step 104: traverse all receiving points of the VSP and repeat step 103;

[0039] Step 105: repeating steps 102-104, traversing all virtual receiving points of the ground shot points to be reconstructed, and finally forming a reconstructed pseudo-ground shot collection;

[0040] Step 106: Repeat steps 101-105 to traverse all reconstructed pseudo-surface shot points to form all pseudo-surface shot collections, namely, pseudo-surface seismic data.

[0041] Preferably, it also includes:

[0042] Preprocessing is performed on the pseudo-surface seismic data, and then common scattering point gathers are constructed.

[0043] Preferably, the pretreatment comprises:

[0044] The linear noise of the pseudo-ground shot gather is attenuated by FK and median filtering;

[0045] Negatively correlated records in the quasi-ground gun set are removed.

[0046] In a third aspect, an embodiment of the present disclosure further provides an electronic device, the electronic device comprising:

[0047] A memory storing executable instructions;

[0048] A processor runs the executable instructions in the memory to implement the velocity analysis method of pseudo-surface seismic data.

[0049] In a fourth aspect, an embodiment of the present disclosure further provides a computer-readable storage medium storing a computer program, which implements the velocity analysis method of pseudo-surface seismic data when executed by a processor.

[0050] The methods and apparatus of the present invention have other features and advantages that will be apparent from, or will be described in detail in, the accompanying drawings and subsequent detailed descriptions incorporated herein, which together serve to explain the specific principles of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] The above and other objects, features and advantages of the present invention will become more apparent through a more detailed description of exemplary embodiments of the present invention in conjunction with the accompanying drawings, wherein like reference numerals generally represent like components throughout the exemplary embodiments of the present invention.

[0052] Figure 1 A flow chart showing the steps of a velocity analysis method for pseudo-surface seismic data according to an embodiment of the present invention.

[0053] Figure 2 A schematic diagram showing downlink multiple wave data of a certain work area after wave field separation of input VSP data according to an embodiment of the present invention is shown.

[0054] Figure 3a , Figure 3b Schematic diagrams of a shot gather containing linear noise and a shot gather containing negatively correlated records according to an embodiment of the present invention are respectively shown.

[0055] Figure 4a , Figure 4b Schematic diagrams of a shot gather with linear noise removed and a shot gather with negative correlation channels removed are respectively shown according to an embodiment of the present invention.

[0056] Figure 5a , Figure 5b Schematic diagrams of a common-center point gather of pseudo-surface seismic data and its velocity spectrum according to an embodiment of the present invention are respectively shown.

[0057] Figure 6a , Figure 6b Schematic diagrams of common scattering point gathers and their velocity spectra reconstructed from pseudo-surface seismic data according to an embodiment of the present invention are respectively shown.

[0058] Figure 7 A schematic diagram showing a velocity profile of an output according to an embodiment of the present invention is shown.

[0059] Figure 8 A block diagram of a velocity analysis device for pseudo-surface seismic data according to an embodiment of the present invention is shown.

[0060] Description of reference numerals:

[0061] 201. Pseudo-surface seismic data construction module; 202. Common scattering point gather construction module; 203. Velocity analysis module. DETAILED DESCRIPTION

[0062] The preferred embodiments of the present invention will be described in more detail below. Although the preferred embodiments of the present invention are described below, it should be understood that the present invention can be implemented in various forms and should not be limited to the embodiments set forth herein.

[0063] The present invention provides a velocity analysis method for pseudo-surface seismic data, comprising:

[0064] Construct pseudo-surface seismic data based on VSP downlink multiple wave data;

[0065] Construct common scattering point gathers based on initial seismic wave velocity and pseudo-surface seismic data;

[0066] Velocity analysis is performed on the common scattering point gathers to obtain the final velocity analysis results.

[0067] In one example, obtaining the final velocity analysis results includes:

[0068] The velocity obtained by velocity analysis of the common scattering point gather is used as the initial seismic wave velocity. The common scattering point gather and velocity analysis are repeatedly constructed until the difference in the velocities obtained from the two velocity analyses is less than the set threshold, and the final velocity analysis result is obtained.

[0069] In one example, the VSP downlink multiple wave data is downlink multiple wave data after wavefield separation.

[0070] In one example, constructing pseudo-surface seismic data includes:

[0071] Step 101: Determine the reconstructed pseudo-ground shot point position;

[0072] Step 102: Determine the virtual receiver position;

[0073] Step 103: Select a receiving point of the VSP, correlate the downlink direct wave of the receiving point at the virtual receiver position with the downlink multiple wave records of the receiving point at the pseudo-ground shot point position, and superimpose the correlation results on the pseudo-ground seismic trace;

[0074] Step 104: traverse all receiving points of the VSP and repeat step 103;

[0075] Step 105: repeating steps 102-104, traversing all virtual receiving points of the ground shot points to be reconstructed, and finally forming a reconstructed pseudo-ground shot collection;

[0076] Step 106: Repeat steps 101-105 to traverse all reconstructed pseudo-surface shot points to form all pseudo-surface shot collections, namely, pseudo-surface seismic data.

[0077] In one example, it also includes:

[0078] The pseudo-surface seismic data are preprocessed and then the common scattering point gathers are constructed.

[0079] In one example, preprocessing includes:

[0080] The linear noise of the pseudo-ground shot gather is attenuated by FK and median filtering;

[0081] Negatively correlated records in the quasi-ground gun set are removed.

[0082] Specifically, the pseudo-surface seismic data are constructed based on the VSP downlink multiple wave data, wherein the VSP downlink multiple wave data are the downlink multiple wave data after wave field separation.

[0083] Constructing pseudo-surface seismic data includes: Step 101: Select a shot point S i (i=1,2,…N) will be used as the reconstructed pseudo-ground shot point position; Step 102: Select one shot point S among the remaining shot points j The position of (j=1,2,…N,j≠i) is used as a virtual receiver position R of the reconstructed gun j ; Step 103: Select a receiving point R of VSP k (i=1,2,…M), S j The receiving point R k The downlink direct wave and S i The receiving point R kThe downlink multiple wave records are correlated and the correlation results are superimposed on the output pseudo-surface seismic trace R i,j Step 104: traverse all receiving points R of VSP k Repeat step 103 to output the final pseudo ground channel R i,j Record; Step 105: Repeat steps 102-104 to traverse all virtual receiving points R of the ground shot point to be reconstructed j , and finally form a reconstructed pseudo-ground shot collection SG i (i=1,2,…N); Step 106: Repeat steps 101-105 to traverse all reconstructed pseudo-ground shot points S i , forming all pseudo-ground shot sets SG1, SG2, ...SG N , which is the pseudo-surface seismic data.

[0084] Preprocessing is performed on the pseudo-surface seismic data, including: attenuating linear noise of the pseudo-surface shot gathers through FK and median filtering; and removing negatively correlated records in the pseudo-surface shot gathers.

[0085] According to the initial seismic wave velocity and pseudo-surface seismic data, a common scattering point gather is constructed; velocity analysis is performed on the common scattering point gather, and a velocity spectrum is produced using the velocity scanning method. Velocity picking is performed through human-computer interaction, and the velocity obtained from the velocity analysis of the common scattering point gather is used as the initial seismic wave velocity. The common scattering point gather and velocity analysis are repeatedly constructed until the difference in the velocities obtained from the two velocity analyses is less than the set threshold, thereby obtaining the final velocity analysis result.

[0086] The present invention also provides a velocity analysis device for pseudo-surface seismic data, comprising:

[0087] The pseudo-surface seismic data construction module constructs pseudo-surface seismic data based on VSP downlink multiple wave data;

[0088] Common scattering point gather construction module, which constructs common scattering point gathers based on initial seismic wave velocity and pseudo-surface seismic data;

[0089] The velocity analysis module performs velocity analysis on the common scattering point gathers to obtain the final velocity analysis results.

[0090] In one example, obtaining the final velocity analysis results includes:

[0091] The velocity obtained by velocity analysis of the common scattering point gather is used as the initial seismic wave velocity. The common scattering point gather and velocity analysis are repeatedly constructed until the difference in the velocities obtained from the two velocity analyses is less than the set threshold, and the final velocity analysis result is obtained.

[0092] In one example, the VSP downlink multiple wave data is downlink multiple wave data after wavefield separation.

[0093] In one example, constructing pseudo-surface seismic data includes:

[0094] Step 101: Determine the reconstructed pseudo-ground shot point position;

[0095] Step 102: Determine the virtual receiver position;

[0096] Step 103: Select a receiving point of the VSP, correlate the downlink direct wave of the receiving point at the virtual receiver position with the downlink multiple wave records of the receiving point at the pseudo-ground shot point position, and superimpose the correlation results on the pseudo-ground seismic trace;

[0097] Step 104: traverse all receiving points of the VSP and repeat step 103;

[0098] Step 105: repeating steps 102-104, traversing all virtual receiving points of the ground shot points to be reconstructed, and finally forming a reconstructed pseudo-ground shot collection;

[0099] Step 106: Repeat steps 101-105 to traverse all reconstructed pseudo-surface shot points to form all pseudo-surface shot collections, namely, pseudo-surface seismic data.

[0100] In one example, it also includes:

[0101] The pseudo-surface seismic data are preprocessed and then the common scattering point gathers are constructed.

[0102] In one example, preprocessing includes:

[0103] The linear noise of the pseudo-ground shot gather is attenuated by FK and median filtering;

[0104] Negatively correlated records in the quasi-ground gun set are removed.

[0105] Specifically, the pseudo-surface seismic data are constructed based on the VSP downlink multiple wave data, wherein the VSP downlink multiple wave data are the downlink multiple wave data after wave field separation.

[0106] Constructing pseudo-surface seismic data includes: Step 101: Select a shot point S i (i=1,2,…N) will be used as the reconstructed pseudo-ground shot point position; Step 102: Select one shot point S among the remaining shot points j The position of (j=1,2,…N,j≠i) is used as a virtual receiver position R of the reconstructed gun j ; Step 103: Select a receiving point R of VSP k (i=1,2,…M), S jThe receiving point R k The downlink direct wave and S i The receiving point R k The downlink multiple wave records are correlated and the correlation results are superimposed on the output pseudo-surface seismic trace R i,j Step 104: traverse all receiving points R of VSP k Repeat step 103 to output the final pseudo ground channel R i,j Record; Step 105: Repeat steps 102-104 to traverse all virtual receiving points R of the ground shot point to be reconstructed j , and finally form a reconstructed pseudo-ground shot collection SG i (i=1,2,…N); Step 106: Repeat steps 101-105 to traverse all reconstructed pseudo-ground shot points S i , forming all pseudo-ground shot sets SG1, SG2, ...SG N , which is the pseudo-surface seismic data.

[0107] Preprocessing is performed on the pseudo-surface seismic data, including: attenuating linear noise of the pseudo-surface shot gathers through FK and median filtering; and removing negatively correlated records in the pseudo-surface shot gathers.

[0108] According to the initial seismic wave velocity and pseudo-surface seismic data, a common scattering point gather is constructed; velocity analysis is performed on the common scattering point gather, and a velocity spectrum is produced using the velocity scanning method. Velocity picking is performed through human-computer interaction, and the velocity obtained from the velocity analysis of the common scattering point gather is used as the initial seismic wave velocity. The common scattering point gather and velocity analysis are repeatedly constructed until the difference in the velocities obtained from the two velocity analyses is less than the set threshold, thereby obtaining the final velocity analysis result.

[0109] The present invention also provides an electronic device, which includes: a memory storing executable instructions; and a processor, which runs the executable instructions in the memory to implement the velocity analysis method of pseudo-ground seismic data.

[0110] The present invention also provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, the velocity analysis method of pseudo-surface seismic data is implemented.

[0111] To facilitate understanding of the solutions and effects of the embodiments of the present invention, four specific application examples are given below. Those skilled in the art should understand that the examples are only for facilitating understanding of the present invention, and any specific details thereof are not intended to limit the present invention in any way.

[0112] Example 1

[0113] Figure 1A flow chart showing the steps of a velocity analysis method for pseudo-surface seismic data according to an embodiment of the present invention.

[0114] like Figure 1 As shown, the velocity analysis method of the pseudo-surface seismic data includes: step 1, constructing the pseudo-surface seismic data according to the VSP downlink multiple wave data; step 2, constructing the common scattering point gather according to the initial seismic wave velocity and the pseudo-surface seismic data; step 3, performing velocity analysis on the common scattering point gather to obtain the final velocity analysis result.

[0115] Figure 2 A schematic diagram showing downlink multiple wave data of a certain work area after wave field separation of input VSP data according to an embodiment of the present invention is shown.

[0116] Figure 3a , Figure 3b Schematic diagrams of a shot gather containing linear noise and a shot gather containing negatively correlated records according to an embodiment of the present invention are respectively shown, both of which are reconstructed pseudo-surface seismic data shot gathers.

[0117] Figure 4a , Figure 4b Schematic diagrams of a shot gather with linear noise removed and a shot gather with negative correlation traces removed according to an embodiment of the present invention are respectively shown, both of which are pre-processed pseudo-surface seismic data.

[0118] Figure 5a , Figure 5b Schematic diagrams of a common-center point gather of pseudo-surface seismic data and its velocity spectrum according to an embodiment of the present invention are respectively shown. It can be seen that the energy cluster trend of the velocity spectrum of the common-center point gather is not obvious. With the increase of depth, the velocity of its energy cluster gradually decreases, which is inconsistent with general geological common sense.

[0119] Figure 6a , Figure 6b Schematic diagrams of common scattering point gathers and their velocity spectra reconstructed from pseudo-surface seismic data according to an embodiment of the present invention are respectively shown. It can be seen from the figure that the velocity spectrum of the common scattering point gathers is relatively concentrated in the shallow layer, and the velocity trend gradually increases from shallow to deep, which is consistent with geological knowledge.

[0120] The output velocity profile is as follows: Figure 7 shown.

[0121] Example 2

[0122] Figure 8 A block diagram of a velocity analysis device for pseudo-surface seismic data according to an embodiment of the present invention is shown.

[0123] like Figure 8As shown, the velocity analysis device for pseudo-surface seismic data comprises:

[0124] A pseudo-surface seismic data construction module 201 constructs pseudo-surface seismic data based on VSP downlink multiple wave data;

[0125] A common scattering point gather construction module 202 is used to construct a common scattering point gather according to the initial seismic wave velocity and the pseudo-surface seismic data;

[0126] The velocity analysis module 203 performs velocity analysis on the common scattering point gathers to obtain a final velocity analysis result.

[0127] As an option, obtain the final velocity analysis results including:

[0128] The velocity obtained by velocity analysis of the common scattering point gather is used as the initial seismic wave velocity. The common scattering point gather and velocity analysis are repeatedly constructed until the difference in the velocities obtained from the two velocity analyses is less than the set threshold, and the final velocity analysis result is obtained.

[0129] As an optional solution, the VSP downlink multiple wave data is the downlink multiple wave data after wave field separation.

[0130] As an optional solution, constructing pseudo-surface seismic data includes:

[0131] Step 101: Determine the reconstructed pseudo-ground shot point position;

[0132] Step 102: Determine the virtual receiver position;

[0133] Step 103: Select a receiving point of the VSP, correlate the downlink direct wave of the receiving point at the virtual receiver position with the downlink multiple wave records of the receiving point at the pseudo-ground shot point position, and superimpose the correlation results on the pseudo-ground seismic trace;

[0134] Step 104: traverse all receiving points of the VSP and repeat step 103;

[0135] Step 105: repeating steps 102-104, traversing all virtual receiving points of the ground shot points to be reconstructed, and finally forming a reconstructed pseudo-ground shot collection;

[0136] Step 106: Repeat steps 101-105 to traverse all reconstructed pseudo-surface shot points to form all pseudo-surface shot collections, namely, pseudo-surface seismic data.

[0137] As an option, it also includes:

[0138] The pseudo-surface seismic data are preprocessed and then the common scattering point gathers are constructed.

[0139] As an option, pre-processing includes:

[0140] The linear noise of the pseudo-ground shot gather is attenuated by FK and median filtering;

[0141] Negatively correlated records in the quasi-ground gun set are removed.

[0142] Example 3

[0143] The present disclosure provides an electronic device, which includes: a memory storing executable instructions; a processor running the executable instructions in the memory to implement the velocity analysis method of pseudo-ground seismic data.

[0144] An electronic device according to an embodiment of the present disclosure includes a memory and a processor.

[0145] The memory is used to store non-temporary computer-readable instructions. Specifically, the memory may include one or more computer program products, which may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include random access memory (RAM) and / or cache memory (cache), etc. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, etc.

[0146] The processor may be a central processing unit (CPU) or other forms of processing units having data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In one embodiment of the present disclosure, the processor is used to run the computer-readable instructions stored in the memory.

[0147] Those skilled in the art should be able to understand that in order to solve the technical problem of how to obtain a good user experience, the present embodiment may also include well-known structures such as a communication bus and an interface, and these well-known structures should also be included in the protection scope of the present disclosure.

[0148] For detailed description of this embodiment, reference may be made to the corresponding descriptions in the aforementioned embodiments, which will not be repeated here.

[0149] Example 4

[0150] An embodiment of the present disclosure provides a computer-readable storage medium storing a computer program, which implements the velocity analysis method of pseudo-surface seismic data when executed by a processor.

[0151] According to the computer-readable storage medium of the embodiment of the present disclosure, non-transitory computer-readable instructions are stored thereon. When the non-transitory computer-readable instructions are executed by a processor, all or part of the steps of the above-mentioned methods of each embodiment of the present disclosure are executed.

[0152] The above-mentioned computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROM and DVD), magneto-optical storage media (e.g., MO), magnetic storage media (e.g., magnetic tape or mobile hard disk), media with built-in rewritable non-volatile memory (e.g., memory card) and media with built-in ROM (e.g., ROM box).

[0153] Those skilled in the art should understand that the purpose of the above description of the embodiments of the present invention is only to exemplarily illustrate the beneficial effects of the embodiments of the present invention, and is not intended to limit the embodiments of the present invention to any given examples.

[0154] The embodiments of the present invention have been described above, and the above description is exemplary, not exhaustive, and is not limited to the disclosed embodiments. Many modifications and changes will be apparent to those skilled in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A velocity analysis method for pseudo-surface seismic data, characterized in that: include: Construct pseudo-surface seismic data based on VSP downlink multiple wave data; constructing a common scattering point gather according to the initial seismic wave velocity and the pseudo-surface seismic data; Performing velocity analysis on the common scattering point gathers to obtain a final velocity analysis result; The final speed analysis results include: The velocity obtained by performing velocity analysis on the common scattering point gather is used as the initial seismic wave velocity, and the common scattering point gather and velocity analysis are repeatedly constructed until the difference between the velocities obtained by the two velocity analyses is less than a set threshold value, thereby obtaining a final velocity analysis result; Wherein, the VSP downlink multiple wave data is the downlink multiple wave data after wave field separation; Among them, constructing pseudo-surface seismic data includes: Step 101: Determine the reconstructed pseudo-ground shot point position; Step 102: Determine the virtual receiver position; Step 103: Select a receiving point of the VSP, correlate the downlink direct wave of the receiving point at the virtual receiver position with the downlink wave multiple wave records of the receiving point at the pseudo-surface shot point position, and superimpose the correlation result on the pseudo-surface seismic trace; Step 104: traverse all receiving points of the VSP and repeat step 103; Step 105: repeating steps 102-104, traversing all virtual receiving points of the ground shot points to be reconstructed, and finally forming a reconstructed pseudo-ground shot collection; Step 106: repeating steps 101-105, traversing all reconstructed pseudo-surface shot points, forming all pseudo-surface shot collections, namely, pseudo-surface seismic data; Among them, it also includes: Preprocessing the pseudo-surface seismic data, and then constructing a common scattering point gather; Wherein, the preprocessing includes: The linear noise of the pseudo-ground shot gather is attenuated by FK and median filtering; Negatively correlated records in the quasi-ground gun set are removed.

2. A velocity analysis device for pseudo-surface seismic data, used to implement the velocity analysis method for pseudo-surface seismic data according to claim 1, characterized in that: include: The pseudo-surface seismic data construction module constructs pseudo-surface seismic data based on VSP downlink multiple wave data; A common scattering point gather construction module is used to construct a common scattering point gather according to the initial seismic wave velocity and the pseudo-surface seismic data; The velocity analysis module performs velocity analysis on the common scattering point gathers to obtain a final velocity analysis result.

3. The velocity analysis device for pseudo-surface seismic data according to claim 2, wherein: Obtaining the final speed analysis results includes: The velocity obtained by velocity analysis on the common scattering point gather is used as the initial seismic wave velocity, and the common scattering point gather and velocity analysis are repeatedly constructed until the difference between the velocities obtained from the two velocity analyses is less than a set threshold, thereby obtaining a final velocity analysis result.

4. An electronic device, characterized in that: The electronic device comprises: A memory storing executable instructions; A processor, wherein the processor runs the executable instructions in the memory to implement the velocity analysis method for pseudo-surface seismic data according to claim 1.

5. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the velocity analysis method for pseudo-surface seismic data according to claim 1 is implemented.

Images