A full-space variable-speed mapping method and system based on morphological inversion

Abstract

The present invention relates to a full-space variable-velocity mapping method and system based on morphological inversion, the method comprising: obtaining a three-dimensional time-domain seismic data volume of a target structure, establishing a velocity mathematical model for morphological inversion mapping, and constructing a velocity field according to the velocity mathematical model; performing time-depth conversion on the three-dimensional time-domain seismic data volume according to the velocity field, and calculating a morphological correction factor; substituting the morphological correction factor into the velocity mathematical model to correct the velocity field; judging whether the morphological correction factor is lower than a threshold value: if so, substituting the morphological correction factor into the velocity mathematical model to obtain a final velocity field; and determining a depth structural map of the target structure according to the final velocity field. The present invention applies the full-space variable-velocity mapping technology of the morphological inversion method, which can realize high-precision mapping of structures of different complexity, and provide guidance for accurate implementation of structures, well site deployment, exploration and development, and oilfield production and reserve increase.

Description

Technical Field

[0001] The invention belongs to seismic exploration and development technology in the field of petroleum exploration, and in particular relates to a full-space variable-speed mapping method and system based on morphological inversion. Background Art

[0002] Two crucial parameters in geophysics are reflection time and velocity. As oil and gas exploration continues to deepen, higher requirements are placed on the development of seismic technology. From imaging accuracy to lithology analysis, from acquisition to data processing and interpretation, velocity information analysis and utilization are particularly important.

[0003] In seismic exploration, the concept of velocity is more often used in seismic data processing and interpretation. There are many concepts of velocity due to different application links. For example, in the process of seismic data processing, the dynamic correction velocity is the stacking velocity, and the migration processing uses the migration velocity. In the process of data interpretation, the time-depth conversion involves the average velocity, and the interpretation of stratum lithology involves the concept of layer velocity.

[0004] In the process of establishing the mapping velocity, the layer velocity of each layer or the average velocity of each seismic interface is usually obtained. When the lateral changes of the formation attitude and physical properties are not large, the imaging velocity (stack velocity) obtained by the common center point gather is close to the root mean square velocity. At the location where the formation dip is large, the stacking velocity must be corrected for the dip, and the corrected velocity is also close to the root mean square velocity of each layer at that point. The stacking velocity dip correction formula is:

[0005] V r =V s cosθ

[0006] Where: V r is the root mean square velocity, V s is the stacking velocity, and θ is the formation dip. The formation dip is obtained from the stacking velocity and the structural map.

[0007] The layer velocity of each layer can be obtained by converting the root mean square velocity through the Dix formula. The formula for calculating the layer velocity from the root mean square velocity is:

[0008]

[0009] Where: t 0,n is the travel time of the earthquake from the 1st floor to the nth floor, is the RMS velocity of the nth layer.

[0010] The average velocity is narrowly defined as the average velocity of the medium above a certain interface in a group of horizontal layered media, which is the ratio of the total thickness of the layers above the interface vertically through the seismic wave to the total propagation time. The average velocity of N layers of horizontal layered media is:

[0011]

[0012] Where: V a is the average speed, h i , V i are the thickness and speed of each layer respectively.

[0013] The propagation of seismic waves in the medium follows the principle of "following the shortest path in the required time", that is, the Fermat principle. According to the characteristics of the work area and the difference in the velocity type used for time-depth conversion, the variable velocity mapping method is very different. At present, the most widely used methods are the average velocity method, the curved line method and the refraction line method. In addition, in data processing, the offset velocity field can be directly used to perform prestack depth migration to obtain depth domain data, and through depth domain interpretation, the depth domain structural map can be directly obtained.

[0014] The average velocity method is the most commonly used time-to-depth conversion method. This method directly extracts the average velocity of a certain reflection interface from the superposition velocity field and uses the time-to-depth conversion method to obtain the depth map of the reflection interface.

[0015] Using the velocity field formed by stacking velocity to perform time-depth conversion is based on the assumption that the horizontal layered medium and the CMP gather reflection wave time-distance curve are hyperbolic to achieve velocity analysis. When the velocity changes greatly laterally, the above method will cause large errors.

[0016] In summary, there is an urgent need for a more effective, accurate and universal method to construct maps. Summary of the invention

[0017] In order to solve the problem that in the existing process of establishing mapping velocity, when the target area structure is complex and the velocity changes laterally, there are large errors, in the first aspect of the present invention, a full-space variable velocity mapping method based on morphological inversion is provided, including: obtaining a three-dimensional time domain seismic data volume of the target structure, establishing a velocity mathematical model for morphological inversion mapping, and constructing a velocity field according to the velocity mathematical model; performing time-to-depth conversion on the three-dimensional time domain seismic data volume according to the velocity field, and calculating a morphological correction factor; substituting the morphological correction factor into the velocity mathematical model to correct the velocity field; judging whether the morphological correction factor is lower than a threshold value: if so, substituting the morphological correction factor into the velocity mathematical model to obtain a final velocity field; and determining a depth structure map of the target structure according to the final velocity field.

[0018] In some embodiments of the present invention, performing time-to-depth conversion on the three-dimensional time-domain seismic data volume according to the velocity field and calculating the morphological correction factor includes: performing time-to-depth conversion on the three-dimensional time-domain seismic data volume according to the velocity field to obtain the time-domain seismic profile morphology, depth-domain seismic profile morphology, time structural morphology and depth structural morphology corresponding to multiple survey lines; comparing the time-domain and depth-domain seismic profile morphologies of each survey line to obtain the profile morphology difference of each survey line; comparing the morphology of the time structural map and the depth structural map of each survey line to obtain the structural morphology difference of each survey line; and calculating the morphological correction factor based on the profile morphology difference and structural morphology difference of each survey line.

[0019] Furthermore, the morphology correction factor is calculated by the following method:

[0020]

[0021] Among them, δ n is the morphological correction factor of the nth sample point, a is the relative error accuracy of the structural mapping, and its value range is 0-0.003; H is the structural depth; h0 is the base surface altitude; t0 is the reflection time; V is the mapping speed; n is the sequence number of the sample point, and numsamples(V) is the number of sample points of the mapping speed.

[0022] In some embodiments of the present invention, the velocity mathematical model includes: extracting the average velocity of the reflection interface of the target structure from the stacking velocity field; and performing pre-stack depth migration using the migration velocity field to obtain depth domain data.

[0023] In some embodiments of the present invention, determining the depth structure map of the target structure according to the final velocity field includes: calculating the apparent depth structure map of the target structure according to the final velocity field; and correcting the apparent depth structure map to obtain the depth structure map of the target structure.

[0024] In the above embodiment, constructing the velocity field according to the velocity mathematical model includes: inputting the stacking velocity field of the target structure; obtaining the iso-T0 structure map, and obtaining the inclination correction factor according to the inclination correction factor; performing inclination correction on the stacking velocity field according to the inclination correction factor; and calculating the root mean square velocity field, the layer velocity field and the average velocity field in sequence according to the corrected stacking velocity field to obtain the average velocity field of the target structure.

[0025] The second aspect of the present invention provides a full-space variable-velocity mapping system based on morphological inversion, comprising: a construction module for acquiring a three-dimensional time-domain seismic data volume of a target structure, establishing a velocity mathematical model for morphological inversion mapping, and constructing a velocity field based on the velocity mathematical model; a calculation module for performing time-to-depth conversion on the three-dimensional time-domain seismic data volume based on the velocity field, and calculating a morphological correction factor; a correction module for substituting the morphological correction factor into the velocity mathematical model to correct the velocity field; a determination module for judging whether the morphological correction factor is lower than a threshold value: if so, substituting the morphological correction factor into the velocity mathematical model to obtain a final velocity field; and determining a depth structure map of the target structure based on the final velocity field.

[0026] Furthermore, the calculation module includes: a conversion unit, which is used to perform time-depth conversion on the three-dimensional time-domain seismic data volume according to the velocity field to obtain the corresponding time-domain seismic profile morphology, depth-domain seismic profile morphology, time structure morphology and depth structure morphology of multiple survey lines; a comparison unit, which is used to compare the time-domain and depth-domain seismic profile morphologies of each survey line to obtain the profile morphology difference of each survey line; compare the morphology of the time structure map and the depth structure map of each survey line to obtain the structural morphology difference of each survey line; and a calculation unit, which is used to calculate the morphology correction factor according to the profile morphology difference and structural morphology difference of each survey line.

[0027] The third aspect of the present invention provides an electronic device, comprising: one or more processors; a storage device for storing one or more programs, when the one or more programs are executed by the one or more processors, the one or more processors implement the full-space variable-speed mapping method based on morphological inversion provided by the first aspect of the present invention.

[0028] A fourth aspect of the present invention provides a computer-readable medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the full-space variable-speed mapping method based on morphological inversion provided by the first aspect of the present invention.

[0029] The beneficial effects of the present invention are:

[0030] 1. The full-space variable-speed mapping technology using the morphological inversion method can achieve high-precision mapping of structures of different complexity, providing guidance for the accurate implementation of structures, well site deployment, exploration and development, and oilfield production and reserve increase;

[0031] 2. Through the iteration of the morphological correction factor, the morphological difference between the deep structural map and the t0 structural map is eliminated; the deep structural map obtained by the velocity field of the morphological inversion method has a gentle structural change, more delicate fault characterization, obvious structural characteristics, and a relatively reasonable result with high accuracy. BRIEF DESCRIPTION OF THE DRAWINGS

[0032] Figure 1 It is a basic flow chart of a full-space variable-speed mapping method based on morphological inversion in some embodiments of the present invention;

[0033] Figure 2 It is a schematic diagram of a specific process of a full-space variable-speed mapping method based on morphological inversion in some embodiments of the present invention;

[0034] Figure 3a isot0 configuration diagram of the target configuration in some embodiments of the present invention;

[0035] Figure 3b A target structure diagram predicted by using prior art in some embodiments of the present invention;

[0036] Figure 3c A target structure diagram predicted by the present invention in some embodiments of the present invention;

[0037] Figure 3d A velocity cloud diagram obtained by using the prior art in some embodiments of the present invention;

[0038] Figure 3e The velocity cloud diagram obtained by using the present invention in some embodiments of the present invention;

[0039] Figure 4a A deep structure map of a target structure predicted by the present invention in some embodiments of the present invention;

[0040] Figure 4b A depth structure map of a target structure predicted by using the prior art in some embodiments of the present invention;

[0041] Figure 5 Schematic diagram of the specific structure of the full-space variable speed mapping system based on morphological inversion in some embodiments of the present invention;

[0042] Figure 6 It is a schematic diagram of the structure of an electronic device in some embodiments of the present invention. DETAILED DESCRIPTION

[0043] The principles and features of the present invention are described below in conjunction with the accompanying drawings. The examples given are only used to explain the present invention and are not used to limit the scope of the present invention.

[0044] refer to Figure 1 and Figure 2In a first aspect of the present invention, a full-space variable-velocity mapping method based on morphological inversion is provided, comprising: S100. Acquire a three-dimensional time-domain seismic data volume of a target structure, establish a velocity mathematical model for morphological inversion mapping, and construct a velocity field according to the velocity mathematical model; S200. Perform time-to-depth conversion on the three-dimensional time-domain seismic data volume according to the velocity field, and calculate a morphological correction factor; S300. Substitute the morphological correction factor into the velocity mathematical model to correct the velocity field; S400. Determine whether the morphological correction factor is lower than a threshold value: if so, substitute the morphological correction factor into the velocity mathematical model to obtain a final velocity field; determine a depth structure map of the target structure according to the final velocity field.

[0045] In step S200 of some embodiments of the present invention, the time-depth conversion of the three-dimensional time-domain seismic data volume according to the velocity field and the calculation of the morphological correction factor include: S201. According to the velocity field, the time-depth conversion of the three-dimensional time-domain seismic data volume is performed to obtain the corresponding time-domain seismic profile morphology, depth-domain seismic profile morphology, time-tectonic morphology and depth-tectonic morphology of multiple survey lines; S202. Comparing the time-domain and depth-domain seismic profile morphology of each survey line to obtain the profile morphology difference of each survey line; S203. Comparing the morphology of the time-tectonic map and the depth-tectonic map of each survey line to obtain the structural morphology difference of each survey line; S204. Calculating the morphological correction factor according to the profile morphology difference and structural morphology difference of each survey line. It should be noted that the length and shape of the survey line can be set as needed, for example, the survey line can be composed of one or more straight lines, curves or multi-segment lines, etc.

[0046] It can be understood that the use of seismic processing stacking velocity to construct a velocity field for time-to-depth conversion is a velocity analysis based on the assumption that the horizontal layered medium and the CMP (Common Middle Point) gather reflection wave time-distance curve are hyperbolic. However, when the target area is complex and the velocity changes laterally, the existing technical methods will cause large errors and their limitations will become apparent. In view of this, the morphology correction factor δ is introduced, which covers the structural undulation and velocity change information, thereby overcoming this defect. Further, the morphology correction factor is calculated by the following method:

[0047]

[0048] Among them, δ nis the morphological correction factor of the nth sample point, a is the relative error accuracy of structural mapping in the seismic exploration industry, and its value range is 0-0.003; H is the structural depth in meters; h0 is the reference surface altitude in meters; t0 is the reflection time; V is the mapping speed in meters per second; n represents the sequence number of the sample point, and numsamples(V) is the number of sample points of the mapping speed.

[0049] Indicatively, Figure 3a - Figure 3e The relevant structural diagrams of the target structure predicted by the prior art and the present invention are shown respectively, wherein the target structure predicted by the prior art is Figure 3b The circled area in the lower left corner is Figure 3a The area indicated by the circle in the lower left corner of the figure is obviously deformed, and the area after applying the method of the embodiment of the present application is Figure 3c and Figure 3a The structural morphology is basically the same. Figure 3d is the speed of the existing method. The local speed in the circle is reversed, and the distortion is obvious. Figure 3e This is the mapping speed of the method in the embodiment of the present application, the longitudinal layering trend is obvious, and the lateral change is uniform and reasonable.

[0050] Figure 4 (a) and Figure 4 (b) are schematic diagrams comparing the results of prediction using the method of the embodiment of the present application and the results of prediction using the prior art. It can be clearly seen from the comparison that the existing method predicts large changes in the plane and deformation of the structural morphology, while the results of the full-space variable speed mapping method of the morphological inversion method show that the structural changes are gentle, the fractures are more delicate, and the "flower-shaped" structural features in the southeast corner of the target area are obvious. The results are more reasonable and more accurate. Please see the table below for specific measurement errors:

[0051]

[0052] In the above embodiment, the velocity mathematical model includes: extracting the average velocity of the reflection interface of the target structure from the stacking velocity field; and performing prestack depth migration using the migration velocity field to obtain depth domain data.

[0053] In the above embodiment, determining the depth structure map of the target structure according to the final velocity field includes: calculating the apparent depth structure map of the target structure according to the final velocity field; and correcting the apparent depth structure map to obtain the depth structure map of the target structure.

[0054] refer to Figure 2In the above embodiment, constructing the velocity field according to the velocity mathematical model includes: inputting the stacking velocity field of the target structure; obtaining the isot0 structure map, and obtaining the inclination correction factor according to the inclination correction factor; performing inclination correction on the stacking velocity field according to the inclination correction factor; and calculating the root mean square velocity field, the layer velocity field and the average velocity field in sequence according to the corrected stacking velocity field to obtain the average velocity field of the target structure.

[0055] refer to Figure 2 In a specific embodiment, the first step is to establish a mathematical model of the velocity of the morphological inversion method; the second step is to calculate a full-space velocity field, and use the velocity field to perform time-depth conversion on the three-dimensional interpretation seismic data volume. Compare the morphology of the seismic profile in the time domain and the depth domain of any survey line, and the morphology of the time structure map and the depth structure map; the third step is to calculate the morphological difference (morphological correction factor); the fourth step is to bring the morphological correction factor into the mathematical model of the velocity of the morphological inversion method to correct the velocity field; the fifth step is to repeat the second, third, and fourth processes until the morphological correction factor is less than the specified value. This iterative method eliminates the morphological difference between the depth structure map and the t0 structure map, and the depth structure map obtained by the velocity field of the morphological inversion method can meet very high accuracy requirements.

[0056] Example 2

[0057] refer to Figure 5 In the second aspect of the present invention, a full-space variable velocity mapping system 1 based on morphological inversion is provided, comprising: a construction module 11, used to obtain a three-dimensional time-domain seismic data volume of a target structure, establish a velocity mathematical model for morphological inversion mapping, and construct a velocity field according to the velocity mathematical model; a calculation module 12, used to perform time-depth conversion on the three-dimensional time-domain seismic data volume according to the velocity field, and calculate a morphological correction factor; a correction module 13, used to substitute the morphological correction factor into the velocity mathematical model to correct the velocity field; a determination module 14, used to determine whether the morphological correction factor is lower than a threshold value: if so, substitute the morphological correction factor into the velocity mathematical model to obtain a final velocity field; and determine the depth structure map of the target structure according to the final velocity field. It can be understood that the morphological correction factor reflects the relative error of the apparent structural depth of the target structure.

[0058] Furthermore, the calculation module 12 includes: a conversion unit, which is used to perform time-depth conversion on the three-dimensional time-domain seismic data volume according to the velocity field to obtain the corresponding time-domain seismic profile morphology, depth-domain seismic profile morphology, time structure morphology and depth structure morphology of multiple survey lines; a comparison unit, which is used to compare the time-domain and depth-domain seismic profile morphologies of each survey line to obtain the profile morphology difference of each survey line; compare the morphology of the time structure map and the depth structure map of each survey line to obtain the structural morphology difference of each survey line; and a calculation unit, which is used to calculate the morphology correction factor according to the profile morphology difference and structural morphology difference of each survey line.

[0059] Example 3

[0060] refer to Figure 6 According to a third aspect of the present invention, an electronic device is provided, comprising: one or more processors; a storage device for storing one or more programs, wherein when the one or more programs are executed by the one or more processors, the one or more processors implement the method of the first aspect of the present invention.

[0061] The electronic device 500 may include a processing device (e.g., a central processing unit, a graphics processing unit, etc.) 501, which can perform various appropriate actions and processes according to a program stored in a read-only memory (ROM) 502 or a program loaded from a storage device 508 into a random access memory (RAM) 503. In the RAM 503, various programs and data required for the operation of the electronic device 500 are also stored. The processing device 501, the ROM 502, and the RAM 503 are connected to each other via a bus 504. An input / output (I / O) interface 505 is also connected to the bus 504.

[0062] Typically, the following devices may be connected to the I / O interface 505: an input device 506 including, for example, a touch screen, a touch pad, a keyboard, a mouse, a camera, a microphone, an accelerometer, a gyroscope, etc.; an output device 507 including, for example, a liquid crystal display (LCD), a speaker, a vibrator, etc.; a storage device 508 including, for example, a hard disk, etc.; and a communication device 509. The communication device 509 may allow the electronic device 500 to communicate with other devices wirelessly or by wire to exchange data. Although Figure 6 An electronic device 500 having various devices is shown, but it should be understood that it is not required to implement or have all the devices shown. More or fewer devices may be implemented or have alternatively. Each box shown in the figure may represent one device, or may represent multiple devices as required.

[0063] In particular, according to an embodiment of the present disclosure, the process described above with reference to the flowchart can be implemented as a computer software program. For example, an embodiment of the present disclosure includes a computer program product, which includes a computer program carried on a computer-readable medium, and the computer program includes a program code for executing the method shown in the flowchart. In such an embodiment, the computer program can be downloaded and installed from the network through a communication device 509, or installed from a storage device 508, or installed from a ROM 502. When the computer program is executed by the processing device 501, the above functions defined in the method of the embodiment of the present disclosure are executed. It should be noted that the computer-readable medium described in the embodiment of the present disclosure can be a computer-readable signal medium or a computer-readable storage medium or any combination of the above two. The computer-readable storage medium can be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device or device, or any combination of the above. More specific examples of computer-readable storage media may include, but are not limited to, an electrical connection with one or more conductors, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above. In an embodiment of the present disclosure, a computer-readable storage medium may be any tangible medium containing or storing a program that can be used by or in conjunction with an instruction execution system, apparatus, or device. In an embodiment of the present disclosure, a computer-readable signal medium may include a data signal propagated in a baseband or as part of a carrier wave, which carries a computer-readable program code. Such propagated data signals may take a variety of forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination of the above. A computer-readable signal medium may also be any computer-readable medium other than a computer-readable storage medium, which may send, propagate, or transmit a program for use by or in conjunction with an instruction execution system, apparatus, or device. The program code embodied on the computer readable medium may be transmitted using any appropriate medium, including but not limited to: wire, optical cable, RF (radio frequency), etc., or any suitable combination of the foregoing.

[0064] The computer-readable medium may be included in the electronic device, or may exist independently without being installed in the electronic device. The computer-readable medium carries one or more computer programs. When the one or more programs are executed by the electronic device, the electronic device:

[0065] Computer program code for performing the operations of embodiments of the present disclosure may be written in one or more programming languages ​​or a combination thereof, including object-oriented programming languages, such as Java, Smalltalk, C++, Python, and conventional procedural programming languages, such as "C" or similar programming languages. The program code may be executed entirely on a user's computer, partially on a user's computer, as a separate software package, partially on a user's computer and partially on a remote computer, or entirely on a remote computer or server. In cases involving a remote computer, the remote computer may be connected to the user's computer via any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computer (e.g., via the Internet using an Internet service provider).

[0066] The flow chart and block diagram in the accompanying drawings illustrate the possible architecture, function and operation of the system, method and computer program product according to various embodiments of the present disclosure. In this regard, each square box in the flow chart or block diagram can represent a module, a program segment or a part of a code, and the module, the program segment or a part of the code contains one or more executable instructions for realizing the specified logical function. It should also be noted that in some implementations as replacements, the functions marked in the square box can also occur in a sequence different from that marked in the accompanying drawings. For example, two square boxes represented in succession can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each square box in the block diagram and / or flow chart, and the combination of the square boxes in the block diagram and / or flow chart can be implemented with a dedicated hardware-based system that performs a specified function or operation, or can be implemented with a combination of dedicated hardware and computer instructions.

[0067] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims

1. A full-space variable-speed mapping method based on morphological inversion, characterized in that: include: Acquire a three-dimensional time-domain seismic data volume of the target structure, establish a velocity mathematical model for morphological inversion mapping, and construct a velocity field based on the velocity mathematical model; The three-dimensional time-domain seismic data volume is subjected to time-depth conversion according to the velocity field, and a morphological correction factor is calculated: the three-dimensional time-domain seismic data volume is subjected to time-depth conversion according to the velocity field to obtain the time-domain seismic profile morphology, depth-domain seismic profile morphology, time-tectonic morphology and depth-tectonic morphology corresponding to a plurality of survey lines; the time-domain and depth-domain seismic profile morphologies of each survey line are compared to obtain the profile morphological difference of each survey line; the morphology of the time-tectonic map and the depth-tectonic map of each survey line are compared to obtain the structural morphological difference of each survey line; the morphological correction factor is calculated according to the profile morphological difference and the structural morphological difference of each survey line; the morphological correction factor is calculated by the following method: ； in, For the n The morphology correction factor for each sample point is a To construct the relative error accuracy of the drawing, its value range is 0-0.003; H is the structural depth; h 0 is the base elevation; t 0 is the reflection time; V is the drawing speed, V n Indicates n The speed of mapping samples, H n To indicate the n The structural depth of each sample point; n Indicates the sequence number of the sample point, is the number of sample points for the mapping speed; Substituting the morphology correction factor into the velocity mathematical model to correct the velocity field; Determine whether the morphology correction factor is lower than a threshold value: if so, substitute the morphology correction factor into the velocity mathematical model to obtain a final velocity field; and determine a depth structure map of the target structure according to the final velocity field.

2. The full-space variable speed mapping method based on morphological inversion according to claim 1 is characterized in that: The speed mathematical model includes: Extract the average velocity of the reflection interface of the target structure from the stacked velocity field; The prestack depth migration is performed using the migration velocity field to obtain depth domain data.

3. The full-space variable speed mapping method based on morphological inversion according to claim 1 is characterized in that: Determining the depth structure map of the target structure according to the final velocity field includes: Calculating an apparent depth structure map of the target structure according to the final velocity field; The apparent depth structure map is corrected to obtain a depth structure map of the target structure.

4. The full-space variable speed mapping method based on morphological inversion according to any one of claims 1 to 3, characterized in that: The constructing of the velocity field according to the velocity mathematical model comprises: Input the stacking velocity field of the target structure; Obtain the iso-T0 structural map and calculate the inclination correction factor based on it; Performing inclination correction on the stacking velocity field according to the inclination correction factor; According to the corrected stacking velocity field, the root mean square velocity field, layer velocity field and average velocity field are calculated in sequence to obtain the average velocity field of the target structure.

5. A full-space variable-speed mapping system based on morphological inversion, characterized in that: include: A construction module is used to obtain a three-dimensional time-domain seismic data volume of a target structure, establish a velocity mathematical model for morphological inversion and mapping, and construct a velocity field based on the velocity mathematical model; A calculation module is used to perform time-depth conversion on the three-dimensional time-domain seismic data volume according to the velocity field, and calculate a morphological correction factor: perform time-depth conversion on the three-dimensional time-domain seismic data volume according to the velocity field to obtain the time-domain seismic profile morphology, depth-domain seismic profile morphology, time-tectonic morphology and depth-tectonic morphology corresponding to multiple survey lines; compare the time-domain and depth-domain seismic profile morphologies of each survey line to obtain the profile morphological difference of each survey line; compare the morphology of the time-tectonic map and the depth-tectonic map of each survey line to obtain the structural morphological difference of each survey line; calculate the morphological correction factor according to the profile morphological difference and the structural morphological difference of each survey line; the morphological correction factor is calculated by the following method: ； in, For the n The morphology correction factor for each sample point is a To construct the relative error accuracy of the drawing, its value range is 0-0.003; H is the structural depth; h 0 is the base elevation; t 0 is the reflection time; V is the drawing speed, V n Indicates n The speed of mapping samples, H n To indicate the n The structural depth of each sample point; n Indicates the sequence number of the sample point, is the number of sample points for the mapping speed; A correction module, used for substituting the morphology correction factor into the velocity mathematical model to correct the velocity field; The determination module is used to determine whether the morphology correction factor is lower than a threshold value. If so, the morphology correction factor is substituted into the velocity mathematical model to obtain a final velocity field. The depth structure map of the target structure is determined according to the final velocity field.

6. An electronic device comprising: one or more processors; A storage device for storing one or more programs, when the one or more programs are executed by the one or more processors, the one or more processors implement the full-space variable-speed mapping method based on morphological inversion as described in any one of claims 1 to 4.

7. A computer readable medium having a computer program stored thereon, wherein: When the computer program is executed by a processor, the full-space variable-speed mapping method based on morphological inversion as claimed in any one of claims 1 to 4 is implemented.

Images