VSP-BASED LEVEL CALIBRATION METHOD AND SEISMIC DEEP DOMAIN PROFILE DEVICE

Description

Field of the Invention

[0001] The present disclosure relates to the technical field of seismic exploration data interpretation, in particular to a VSP-based depth domain seismic profile horizon calibration method, a VSP-based depth domain seismic profile horizon calibration apparatus and a computer readable storage medium.Background of the Invention

[0002] With the continuous development of science and technology, more and more technologies have been applied to the technical field of seismic exploration, such as a borehole seismic exploration technology based on a vertical seismic profile (VSP). This technology accurately obtains seismic wavefields in a depth domain by placing a plurality of sensors in a well.

[0003] In the application process, since the sensors are placed in the well to receive data, which avoids the influence of near surface, and has significant advantages such as a high signal-to-noise ratio, a rich wavefield, and an accurate time-depth relationship, so VSPs are widely used to obtain formation velocity and other information, and also for accurate horizon calibration of surface seismic profiles.

[0004] However, in the actual application process, existing horizon calibration is generally carried out on a time domain profile, with the gradual application and promotion of a depth migration technology, people put forward higher requirements for the accuracy of depth calibration, but an existing technology often adopts a method of performing time-depth conversion based on a vertical seismic corridor stack profile, or utilizes a depth domain synthetic record from well-logging data for calibration, its calibration accuracy cannot meet the needs of current technicians.

[0005] CN-A-101071175 and US-A1-2021 / 103068 disclose a VSP-based depth domain seismic profile horizon calibration method as of the preamble of claim 1, and a VSP-based depth domain seismic profile horizon calibration apparatus as of the preamble of claim 4.Summary of the Invention

[0006] The invention is identified in the appended claims.

[0007] In order to overcome the above technical problems in the prior art, an embodiment of the present disclosure provides a VSP-based depth domain seismic profile horizon calibration method. Upgoing wave data in seismic data is extracted by analyzing the seismic data of the depth domain, accurate reflection information of the depth domain in a formation is determined according to the upgoing wave data, and a profile of the depth domain is further accurately calibrated, so that calibration accuracy is improved, and actual needs are met.

[0008] In order to achieve the above aims, the present invention provides a VSP-based depth domain seismic profile horizon calibration method. The method includes: obtaining VSP wavefield data of a depth domain; determining first arrival time information on the basis of the VSP wavefield data; processing the VSP wavefield data on the basis of a preset filtering method and the first arrival time information to obtain upgoing wave data; generating a zero offset profile of the depth domain on the basis of the upgoing wave data; and calibrating, on the basis of the zero offset profile, a depth-domain profile to be calibrated, so as to generate a corresponding calibration result. Generating the zero offset profile of the depth domain on the basis of the upgoing wave data includes: determining a plurality of preset output depth positions in the depth domain; determining a second time window for processing the upgoing wave on the basis of the first arrival time information; extracting the amplitude of the upgoing wave data on the plurality of preset output depth positions on the basis of the second time window to obtain the upgoing wave amplitude information; extracting the upgoing wave amplitude information to obtain a plurality of corresponding depth-domain data, wherein the upgoing wave amplitude information is in one-to-one correspondence to the plurality of preset output depth positions; and processing the plurality of depth-domain data on the basis of the preset multichannel filtering method to obtain the zero offset profile of the depth domain.

[0009] Optionally, the preset filtering method includes a median filtering method and a mean filtering method, and processing the VSP wavefield data on the basis of a preset filtering method and the first arrival time information to obtain upgoing wave data includes: processing the VSP wavefield data on the basis of the median filtering method or the mean filtering method to determine a corresponding upgoing wave and downgoing wave; determining a first time window on the basis of the first arrival time information; processing the downgoing wave on the basis of the first time window to obtain a downgoing wavelet; executing deconvolution processing on the upgoing wave in data after wavefield separation on the basis of the downgoing wavelet to obtain deconvoluted data; and executing amplitude compensation processing on the deconvoluted data to obtain the upgoing wave data.

[0010] Optionally, calibrating, on the basis of the zero offset profile, a depth-domain profile to be calibrated, so as to generate a corresponding calibration result includes: determining the to-be-calibrated depth-domain profile; determining a depth error between the zero offset profile and the to-be-calibrated depth-domain profile and a corresponding relationship of reflection characteristics; and analyzing the depth error and the corresponding relationship of the reflection characteristics to generate the corresponding calibration result, wherein the calibration result includes depth accuracy evaluation information of the profile of the depth domain and corresponding correction suggestion information.

[0011] Correspondingly, the present inventionfurther provides a VSP-based depth domain seismic profile horizon calibration apparatus, and the apparatus includes: a borehole seismic wavefield obtaining unit, configured to obtain VSP wavefield data of a depth domain; a first arrival time obtaining unit, configured to determine first arrival time information on the basis of the VSP wavefield data; a preprocessing unit, configured to process the VSP wavefield data on the basis of a preset filtering method and the first arrival time information to obtain upgoing wave data; a profile information generating unit, configured to generate a zero offset profile of the depth domain on the basis of the upgoing wave data; and a calibrating unit, configured to calibrate, on the basis of the zero offset profile, a depth-domain profile to be calibrated, so as to generate a corresponding calibration result. The profile information generating unit comprises: a depth domain amplitude value extracting module, configured to: determine a plurality of preset output depth positions in the depth domain; determine a second time window for processing the upgoing wave on the basis of the first arrival time information; and extract the amplitude of the upgoing wave data on the plurality of preset output depth positions on the basis of the second time window to obtain the upgoing wave amplitude information; and a filtering module, configured to: extract the upgoing wave amplitude information to obtain a plurality of corresponding depth-domain data, wherein the upgoing wave amplitude information is in one-to-one correspondence to the plurality of preset output depth positions; and process the plurality of depth-domain data on the basis of the preset multichannel filtering method to obtain the zero offset profile of the depth domain.

[0012] Optionally, the preset filtering method includes a median filtering method and a mean filtering method, and the preprocessing unit includes: a wavefield separating module, configured to process the VSP wavefield data on the basis of the median filtering method or the mean filtering method to determine a corresponding upgoing wave and downgoing wave; a first time window determining module, configured to determine a first time window on the basis of the first arrival time information; a wavelet obtaining module, configured to process the downgoing wave on the basis of the first time window to obtain a downgoing wavelet; a wavelet deconvolution module, configured to execute deconvolution processing on the upgoing wave in data after wavefield separation on the basis of the downgoing wavelet to obtain deconvoluted data; and an amplitude compensating module, configured to execute amplitude compensation processing on the deconvoluted data to obtain the upgoing wave data.

[0013] Optionally, the calibrating unit is specifically configured to: determine the to-be-calibrated depth-domain profile; determine a depth error between the zero offset profile and the to-be-calibrated depth-domain profile and a corresponding relationship of reflection characteristics; and analyzing the depth error and the corresponding relationship of the reflection characteristics to generate the corresponding calibration result, wherein the calibration result includes depth accuracy evaluation information of the depth-domain profile and corresponding correction suggestion information.

[0014] On the other hand, the present invention further provides a computer readable storage medium, storing a computer program, and the program, when executed by a processor, implements the VSP-based depth domain seismic profile horizon calibration method provided by the present disclosure.

[0015] Through the technical solution provided by the present disclosure, the present disclosure at least has the following technical effects: the accurate upgoing wave data is extracted from the seismic data by performing seismic data acquisition on the borehole depth domain, performing first arrival picking on the seismic data and performing wavefield separation and deconvolution, corresponding formation reflection information is determined according to the upgoing wave data so as to obtain a zero offset profile of the depth domain, therefore, the depth-domain profile is accurately calibrated, an error of a time-domain corridor in time-depth conversion in prior art is effectively overcome, and the calibration accuracy of the depth domain is improved.

[0016] Other features and advantages of the embodiments of the present disclosure will be illustrated in detail in the following part of specific implementations.Brief Description of Drawings

[0017] The accompanying drawings are used to provide further understanding of embodiments of the present disclosure, and form part of the specification, and are used to explain the embodiments of the present disclosure together with the following specific implementations, but they do not constitute a limitation to the embodiments of the present disclosure. In the accompanying drawings: Fig. 1 is a flow diagram of a specific implementation of a VSP-based depth domain seismic profile horizon calibration method provided by an embodiment of the present disclosure. Fig. 2 is a schematic diagram of VSP wavefield data in a VSP-based depth domain seismic profile horizon calibration method provided by an embodiment of the present disclosure. Fig. 3 is a flow diagram of a specific implementation of obtaining upgoing wave data in a VSP-based depth domain seismic profile horizon calibration method provided by an embodiment of the present disclosure. Fig. 4 is a schematic diagram of upgoing wave data in a VSP-based depth domain seismic profile horizon calibration method provided by an embodiment of the present disclosure. Fig. 5 is schematic contrast diagram of synthetic records (three in the middle) with different resolutions in a depth domain obtained by a VSP depth-domain zero offset profile (right) and well-logging data (left) in a VSP-based depth domain seismic profile horizon calibration method provided by an embodiment of the present disclosure. Fig. 6 is a flow diagram of a specific implementation of generating a calibration result in a VSP-based depth domain seismic profile horizon calibration method provided by an embodiment of the present disclosure. Fig. 7 is a schematic contrast diagram of embedding a well-logging depth-domain synthetic record and a VSP depth-domain zero offset profile into a surface seismic profile in a VSP-based depth domain seismic profile horizon calibration method provided by an embodiment of the present disclosure. Fig. 8 is a schematic structural diagram of a VSP-based depth domain seismic profile horizon calibration apparatus provided by an embodiment of the present disclosure. Detailed Description of the Embodiments

[0018] The specific implementations of embodiments of the present disclosure are illustrated in detail below in combination with the accompanying drawings. It should be understood that the specific implementations described herein are merely used to illustrate and explain the embodiments of the present disclosure and are not used to limit the embodiments of the present disclosure.

[0019] Terms "system" and "network" in the embodiments of the present disclosure may be used interchangeably. "A plurality of" means two or more, and in view of this, "a plurality of" may also be understood as "at least two" in the embodiments of the present disclosure. "And / or" describes an association relationship of associated objects, indicating that there may be three kinds of relationships, for example, A and / or B may indicate that: A exists alone, both A and B exist, and B exists alone. In addition, the character " / ", unless otherwise specified, generally indicates that the associated objects before and after the character is in an "or" relationship. In addition, it needs to be understood that in the description of the embodiments of the present disclosure, words "first", "second", etc. are merely for the purpose of distinguishing and are not to be understood as indicating or suggesting relative importance, nor as indicating or suggesting an order.

[0020] Please referring to Fig. 1, an embodiment of the present disclosure provides a VSP-based depth domain seismic profile horizon calibration method, including: S10) obtaining VSP wavefield data of a depth domain; S20) determining first arrival time information on the basis of the VSP wavefield data; S30) processing the VSP wavefield data on the basis of a preset filtering method and the first arrival time information to obtain upgoing wave data; S40) generating a zero offset profile of the depth domain on the basis of the upgoing wave data; and S50) calibrating, on the basis of the zero offset profile, a depth-domain profile to be calibrated, so as to generate a corresponding calibration result.

[0021] In a possible implementation, in order to accurately calibrate a formation of the depth domain, firstly, borehole seismic wavefield data in the depth domain is obtained, and data detected by a borehole acquisition instrument is transmitted and collected by arranging the borehole acquisition instrument in the formation of the depth domain and receiving and connecting the borehole acquisition instrument through a geophone or optical cable reception, and then, the corresponding seismic wavefield data is collected in the formation of the depth domain by using an excitation of an explosive seismic source or a vibroseis on the surface or in a well. Then, existing first arrival picking methods are used to analyze the above seismic wave data to obtain the corresponding first arrival time information, the existing first arrival picking method include but are not limited to a maximum energy method, a Short Term Average / Long Term Average (STA / LTA) method, a cross-correlation method, a phase tracking method and so on. Fig. 2 is a schematic diagram of VSP wavefield data provided by an embodiment of the present disclosure, and the corresponding first arrival time information may be obtained through the seismic wave data.

[0022] At this time, the seismic wave data is further analyzed and processed. Please referring to Fig. 3, in the embodiment of the present disclosure, the preset filtering method includes a median filtering method and a mean filtering method, and processing the VSP wavefield data on the basis of the preset filtering method and the first arrival time information to obtain the upgoing wave data includes: S31) processing the VSP wavefield data on the basis of the median filtering method or the mean filtering method to determine a corresponding upgoing wave and downgoing wave; S32) determining a first time window on the basis of the first arrival time information; S33) processing the downgoing wave on the basis of the first time window to obtain a downgoing wavelet; S34) executing deconvolution processing on the upgoing wave in data after wavefield separation on the basis of the downgoing wavelet to obtain deconvoluted data; and S35) executing amplitude compensation processing on the deconvoluted data to obtain the upgoing wave data.

[0023] In a possible implementation, firstly, the above seismic wave data is predicted according to the median filtering method or the mean filtering method, the corresponding upgoing wave and downgoing wave are determined, then, the first time window is determined according to the first arrival time information, corresponding data are intercepted in the downgoing wave through the first time window to obtain the corresponding downgoing wavelet, at this time, deconvolution processing is further performed on the upgoing wave in the data after wavefield separation according to the above downgoing wavelet to suppress the influence of multiple waves and obtain the corresponding deconvoluted data, and amplitude compensation processing is executed on the above deconvoluted data so as to remove the amplitude change influence caused by spherical spreading and obtain upgoing wavefield data. Please referring to Fig. 4 which is a schematic diagram of upgoing wave data provided by an embodiment of the present disclosure.

[0024] In the embodiment of the present disclosure, corresponding filtering and denoising processing are performed on both the obtained upgoing wave and downgoing wave of the seismic wave data, so that the more accurate upgoing wave data is obtained, accurate data support is provided for subsequent analysis of relevant information of the depth domain, and accuracy of subsequent depth-domain formation calibration is improved effectively.

[0025] In the embodiment of the present disclosure, generating the zero offset profile of the depth domain on the basis of the upgoing wave data includes: determining upgoing wave amplitude information of the depth domain on the basis of the upgoing wave data; and processing the upgoing wave amplitude information on the basis of a preset multichannel filtering method to obtain the zero offset profile of the depth domain.

[0026] In a possible implementation, after the above upgoing wave data is obtained, the upgoing wave amplitude information of the depth domain is further determined. And determining the upgoing wave amplitude information of the depth domain on the basis of the upgoing wave data includes: determining a plurality of preset output depth positions in the depth domain; determining a second time window on the basis of the first arrival time information; and processing the upgoing wave data on the plurality of preset output depth positions on the basis of the second time window to obtain the upgoing wave amplitude information.

[0027] After an upgoing wavefield is obtained, for each output depth position i in m depths of upgoing wave records, an upgoing wave amplitude A ij is further extracted by taking a short time window downward along the first arrival time of the upgoing wavefield, wherein j is a trace number relative to a depth position above the depth position i. According to a size of a depth interval of an input seismic trace and non-uniformity of a true vertical depth under an inclined well condition, an output depth may be deeply encrypted or interpolated to obtain a higher and more uniform depth sampling density.

[0028] Furthermore, in the embodiment of the present disclosure, processing the upgoing wave amplitude information on the basis of the preset multichannel filtering method to obtain the zero offset profile of the depth domain includes: extracting the upgoing wave amplitude information to obtain a plurality of corresponding depth-domain data, wherein the upgoing wave amplitude information is in one-to-one correspondence to the plurality of preset output depth positions; and processing the plurality of depth-domain data on the basis of the preset multichannel filtering method to obtain the zero offset profile of the depth domain.

[0029] In a possible implementation, after the upgoing wave amplitude information is obtained, n depth-domain record traces corresponding to m depth points are generated through the upgoing wave amplitude information, the depth-domain record traces are further processed through the preset multichannel filtering method so as to improve a signal to noise ratio, for example, the preset multichannel filtering method is a median filtering method or other multichannel filtering methods, and a corresponding zero offset profile of the depth domain is obtained. At this time, a calibration result for the current formation of the depth domain is generated according to the profile information, please referring to Fig. 5 which is a schematic contrast diagram of synthetic records (three in the middle) with different resolutions in a depth domain obtained by a VSP depth-domain zero offset profile(right) and well-logging data (left) provided by an embodiment of the present disclosure.

[0030] Please referring to Fig. 6, in the embodiment of the present disclosure, calibrating, on the basis of the zero offset profile, a depth-domain profile to be calibrated, so as to generate a corresponding calibration result includes: S51) determining the to-be-calibrated depth-domain profile; S52) determining a depth error between the zero offset profile and the to-be-calibrated depth-domain profile and a corresponding relationship of reflection characteristics; and S53) analyzing the depth error and the corresponding relationship of the reflection characteristics to generate the corresponding calibration result, wherein the calibration result includes depth accuracy evaluation information of the profile of the depth domain and corresponding correction suggestion information.

[0031] In a possible implementation, firstly, the to-be-calibrated depth-domain profile is determined, then, the zero offset profile of the depth domain is embedded in the to-be-calibrated depth-domain profile (see Fig. 7), so that a corresponding relationship of reflection characteristics of shallow, middle and deep horizons between the zero offset profile and the to-be-calibrated depth-domain profile and the depth error are analyzed, and depth accuracy evaluation information and depth correction suggestion information of the corresponding depth-domain profile are obtained.

[0032] In the embodiment of the present disclosure, formation reflection information is obtained by travel time along each depth upgoing wave of the depth domain, the zero offset depth profile of VSP is obtained, the to-be-calibrated depth profile is accurately calibrated, an error of a time-domain corridor in time-depth conversion in the prior art is effectively overcome, depth calibration accuracy is obviously improved, accurate data support is provided for subsequent interpretation of a reservoir depth domain, the depth domain interpretation accuracy is improved, and actual needs of technicists are met.

[0033] An apparatus for calibrating a formation of a depth domain provided by an embodiment of the present disclosure is illustrated below in combination with accompanying drawings.

[0034] Please referring to Fig. 8, based on the same inventive conception, an embodiment of the present disclosure provides a VSP-based depth domain seismic profile horizon calibration apparatus. The apparatus includes: a borehole seismic wavefield obtaining unit, configured to obtain VSP wavefield data of a depth domain; a first arrival time obtaining unit, configured to determine first arrival time information on the basis of the VSP wavefield data; a preprocessing unit, configured to process the VSP wavefield data on the basis of a preset filtering method and the first arrival time information to obtain upgoing wave data; a profile information generating unit, configured to generate a zero offset profile of the depth domain on the basis of the upgoing wave data; and a calibrating unit, configured to calibrate, on the basis of the zero offset profile, a depth-domain profile to be calibrated, so as to generate a corresponding calibration result.

[0035] In the embodiment of the present disclosure, the preset filtering method includes a median filtering method and a mean filtering method, and the preprocessing unit includes: a wavefield separating module, configured to process the VSP wavefield data on the basis of the median filtering method or the mean filtering method to determine a corresponding upgoing wave and downgoing wave; a first time window determining module, configured to determine a first time window on the basis of the first arrival time information; a wavelet obtaining module, configured to process the downgoing wave on the basis of the first time window to obtain a downgoing wavelet; a wavelet deconvolution module, configured to execute deconvolution processing on the upgoing wave in data after wavefield separation on the basis of the downgoing wavelet to obtain deconvoluted data; and an amplitude compensating module, configured to execute amplitude compensation processing on the deconvoluted data to obtain the upgoing wave data.

[0036] In the embodiment of the present disclosure, the profile information generating unit includes: a depth domain amplitude value extracting module, configured to determine upgoing wave amplitude information of the depth domain on the basis of the upgoing wave data; and a filtering module, configured to process the upgoing wave amplitude information on the basis of a preset multichannel filtering method to obtain the zero offset profile of the depth domain.

[0037] In the embodiment of the present disclosure, the depth domain amplitude value extracting module is specifically configured to: determine a plurality of preset output depth positions in the depth domain; determine a second time window on the basis of the first arrival time information; and process the upgoing wave data on the plurality of preset output depth positions on the basis of the second time window to obtain the upgoing wave amplitude information.

[0038] In the embodiment of the present disclosure, the filtering module is specifically configured to: analyze the upgoing wave amplitude information to obtain a plurality of corresponding depth-domain data, wherein the upgoing wave amplitude information is in one-to-one correspondence to the plurality of preset output depth positions; and process the plurality of depth-domain data on the basis of the preset multichannel filtering method to obtain the zero offset profile of the depth domain.

[0039] In the embodiment of the present disclosure, the calibrating unit is specifically configured to: determine the to-be-calibrated depth domain profile; determine a depth error between the zero offset profile and the to-be-calibrated depth-domain profile and a corresponding relationship of reflection characteristics; and analyze the depth error and the corresponding relationship of the reflection characteristics to generate the corresponding calibration result, wherein the calibration result includes depth accuracy evaluation information of the depth-domain profile and corresponding correction suggestion information.

[0040] Furthermore, an embodiment of the present disclosure further provides a computer readable storage medium, storing a computer program, and the program, when executed by a processor, implements the VSP-based depth domain seismic profile horizon calibration method of the present disclosure.

[0041] Those skilled in the art may understand that the implementation of all or part of steps in the methods of the above embodiments may be completed by instructing, through a program, relevant hardware, the program is stored in a storage medium, and includes a plurality of instructions to cause a single chip microcomputer, a chip, or a processor to execute all or part of the steps of the methods described in the various embodiments of the present application. The aforementioned storage mediums include: a USB flash disk, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a diskette or an optical disk and various other media that may store program codes.

Claims

1. A VSP-based depth domain seismic profile horizon calibration method, comprising: obtaining (S10) vertical seismic profile VSP wavefield data of a depth domain with a borehole seismic wavefield obtaining unit; determining (S20) first arrival time information on the basis of the VSP wavefield data; processing (S30) the VSP wavefield data on the basis of a preset filtering method and the first arrival time information to obtain upgoing wave data; generating (S40) a zero-offset profile of the depth domain on the basis of the upgoing wave data; and calibrating (S50), on the basis of the zero offset profile, a depth-domain profile to be calibrated, so as to generate a corresponding calibration result; characterized in that the generating (S40) the zero offset profile of the depth domain on the basis of the upgoing wave data comprises: determining a plurality of preset output depth positions in the depth domain; determining a second time window for processing the upgoing wave on the basis of the first arrival time information; extracting the amplitude of the upgoing wave data on the plurality of preset output depth positions on the basis of the second time window to obtain the upgoing wave amplitude information; extracting the upgoing wave amplitude information to obtain a plurality of corresponding depth-domain data, wherein the upgoing wave amplitude information is in one-to-one correspondence to the plurality of preset output depth positions; and processing the plurality of depth-domain data on the basis of the preset multichannel filtering method to obtain the zero offset profile of the depth domain.

2. The method according to claim 1, wherein the preset filtering method comprises a median filtering method and a mean filtering method; and processing the VSP wavefield data on the basis of a preset filtering method and the first arrival time information to obtain upgoing wave data comprises: processing (S31) the VSP wavefield data on the basis of the median filtering method or the mean filtering method to determine a corresponding upgoing wave and downgoing wave; determining (S32) a first time window for processing the downgoing wave on the basis of the first arrival time information; intercepting data from the downgoing wave on the basis of the first time window to obtain a downgoing wavelet; executing (S34) deconvolution processing on the upgoing wave in data after wavefield separation on the basis of the downgoing wavelet to obtain deconvoluted data; and executing (S35) amplitude compensation processing on the deconvoluted data to obtain the upgoing wave data.

3. The method according to claim 1, wherein calibrating, on the basis of the zero offset profile, a depth-domain profile to be calibrated, so as to generate a corresponding calibration result, comprises: determining (S51) the to-be-calibrated depth-domain profile; determining (S52) a depth error between the zero offset profile and the to-be-calibrated depth-domain profile and a corresponding relationship of reflection characteristics; and analyzing (S53) the depth error and the corresponding relationship of the reflection characteristics to generate the corresponding calibration result, wherein the calibration result comprises depth accuracy evaluation information of the profile of the depth domain and corresponding correction suggestion information.

4. A VSP-based depth domain seismic profile horizon calibration apparatus, comprising: a borehole seismic wavefield obtaining unit, configured to obtain VSP wavefield data of a depth domain; a first arrival time obtaining unit, configured to determine first arrival time information on the basis of the VSP wavefield data; a preprocessing unit, configured to process the VSP wavefield data on the basis of a preset filtering method and the first arrival time information to obtain upgoing wave data; a profile information generating unit, configured to generate a zero offset profile of the depth domain on the basis of the upgoing wave data; and a calibrating unit, configured to calibrate, on the basis of the zero offset profile, a depth-domain profile to be calibrated, so as to generate a corresponding calibration result; characterized in that the profile information generating unit comprises: a depth domain amplitude value extracting module, configured to: determine a plurality of preset output depth positions in the depth domain; determine a second time window for processing the upgoing wave on the basis of the first arrival time information; and extract the amplitude of the upgoing wave data on the plurality of preset output depth positions on the basis of the second time window to obtain the upgoing wave amplitude information; and a filtering module, configured to: extract the upgoing wave amplitude information to obtain a plurality of corresponding depth-domain data, wherein the upgoing wave amplitude information is in one-to-one correspondence to the plurality of preset output depth positions; and process the plurality of depth-domain data on the basis of the preset multichannel filtering method to obtain the zero offset profile of the depth domain.

5. The apparatus according to claim 4, wherein the preset filtering method comprises a median filtering method and a mean filtering method, and the preprocessing unit comprises: a wavefield separating module, configured to process the VSP wavefield data on the basis of the median filtering method or the mean filtering method to determine a corresponding upgoing wave and downgoing wave; a first time window determining module, configured to determine a first time window for processing the downgoing wave on the basis of the first arrival time information; a wavelet obtaining module, configured to intercept data from the downgoing wave on the basis of the first time window to obtain a downgoing wavelet; a wavelet deconvolution module, configured to executing deconvolution processing on the upgoing wave in data after wavefield separation on the basis of the downgoing wavelet to obtain deconvoluted data; and an amplitude compensating module, configured to execute amplitude compensation processing on the deconvoluted data to obtain the upgoing wave data.

6. The apparatus according to claim 4, wherein the calibrating unit is specifically configured to: determine the to-be-calibrated depth-domain profile; determine a depth error between the zero offset profile and the to-be-calibrated depth-domain profile and a corresponding relationship of reflection characteristics; and analyze the depth error and the corresponding relationship of the reflection characteristics to generate the corresponding calibration result, wherein the calibration result comprises depth accuracy evaluation information of the depth-domain profile and corresponding correction suggestion information.

7. A computer readable storage medium, storing a computer program, wherein the program, when executed by a processor, implements the VSP-based depth domain seismic profile horizon calibration method according to any of claims 1 to 3.