Three-dimensional variable speed drawing method, device, equipment, storage medium and computer program

Abstract

The application relates to the technical field of seismic data acquisition and processing, and discloses a three-dimensional variable-speed mapping method, device, equipment, storage medium and computer program, which comprises the following steps: obtaining underground geological data, performing model construction on the underground geological data to obtain a velocity model; obtaining seismic data, adjusting model parameters of the velocity model according to the seismic data to obtain a mixed velocity model; obtaining a migration time structure map, converting the time structure map into a depth domain structure map by using the mixed velocity model; and correcting the depth domain structure map according to the seismic data to obtain a mixed velocity model structure map. The application can improve the three-dimensional velocity modeling precision and the structure mapping precision, and especially can improve the micro-amplitude structure mapping precision.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of seismic data interpretation, and in particular to a three-dimensional variable velocity mapping method, device, equipment, storage medium and computer program. BACKGROUND

[0002] Seismic data interpretation is a process of converting seismic data after acquisition and computer processing into geological information, that is, a process of inverting underground geological features from known actual observed seismic data. If reliable seismic data interpretation is to be obtained, the primary step is to obtain accurate structure maps. Specifically, ensuring the accuracy of structure maps mainly depends on obtaining accurate seismic velocity field data.

[0003] Currently, obtaining accurate seismic velocity field is usually achieved by using a variable velocity mapping method. However, the current variable velocity mapping technology mainly focuses on two-dimensional fields such as velocity field modeling and pre-stack migration imaging, and in the application of three-dimensional fields, the model accuracy of the variable velocity mapping method is still insufficient to meet the needs of three-dimensional fine interpretation structure mapping, thereby causing the problem of low precision of three-dimensional structure fine interpretation. SUMMARY

[0004] To solve the above problems, embodiments of the present application provide a three-dimensional variable velocity mapping method, device, equipment, storage medium and computer program.

[0005] In a first aspect, the embodiments of the present application provide a three-dimensional variable velocity mapping method, comprising:

[0006] obtaining underground geological data, constructing a model based on the underground geological data to obtain a velocity model;

[0007] obtaining seismic data, adjusting model parameters of the velocity model based on the seismic data to obtain a hybrid velocity model;

[0008] obtaining a migration time structure map, converting a time-depth relationship of the time structure map based on the hybrid velocity model to obtain a depth domain structure map;

[0009] correcting results of the depth domain structure map based on the seismic data to obtain a hybrid velocity model structure map.

[0010] According to the embodiments of the present application, the model is constructed based on the underground geological data to obtain a velocity model, comprising:

[0011] generating a target layer seismic record fine calibration map based on the underground geological data;

[0012] performing time-depth correlation picking on the target layer seismic record fine calibration map to obtain a time-depth relationship curve;

[0013] The time-depth curve is interpolated by using the Kriging interpolation method to obtain a velocity model of the time-depth curve.

[0014] According to an embodiment of the present application, the model parameter adjustment of the velocity model according to the seismic data is to obtain a mixed velocity model, which comprises:

[0015] The velocity model is parameter-constrained by using the vertical seismic profile logging data in the seismic data to obtain a vertical seismic profile constrained velocity model.

[0016] The vertical seismic profile constrained velocity model is updated and iterated by using the layered data in the seismic data to obtain a layered iterative velocity model.

[0017] The layered iterative velocity model is model-constrained by using the small least square tomography according to the pre-stack depth migration layer velocity body data in the seismic data to obtain a mixed velocity model. According to an embodiment of the present application,

[0018] According to an embodiment of the present application, the velocity model is parameter-constrained by using the vertical seismic profile logging data in the seismic data to obtain a vertical seismic profile constrained velocity model, which comprises:

[0019] A fitting function is established according to the vertical seismic profile logging data and the velocity model.

[0020] The residual sum of squares of the fitting function and the vertical seismic profile logging data is calculated.

[0021] The constrained velocity model is calculated by using the least square method according to the residual sum of squares.

[0022] The above steps are performed on the constrained velocity model for a preset number of times to obtain a vertical seismic profile constrained velocity model.

[0023] According to an embodiment of the present application, the vertical seismic profile constrained velocity model is updated and iterated by using the layered data in the seismic data to obtain a layered iterative velocity model, which comprises:

[0024] Well list data is obtained, and the seismic wave velocity of the well list data is calculated by using the vertical seismic profile constrained model.

[0025] The seismic wave data is time-depth converted to obtain a simulated depth.

[0026] The error between the depth of the geological layer of the well list data and the simulated depth is calculated.

[0027] Screen out the well with an error value less than a preset constant as a to-be-trained well, update the vertical seismic profile constrained velocity model by using the error corresponding to the to-be-trained well, and obtain a layered iterative velocity model.

[0028] According to the embodiment of the application, the time-depth relationship conversion of the time structure map by using the mixed velocity model to obtain a depth domain structure map comprises:

[0029] According to the time structure map and the velocity model, the corresponding velocity on the time structure map is calculated;

[0030] The time-depth relationship conversion of the corresponding velocity on the time structure map and the time structure map is performed to obtain a depth value;

[0031] According to the corresponding velocity on the time structure map and the time structure map, a depth structure map is constructed;

[0032] Geological horizon data is acquired, and the depth domain structure map is corrected by using the geological horizon data to obtain a depth domain structure map.

[0033] In a second aspect, the embodiment of the application provides a three-dimensional variable velocity mapping device, characterized in that a model construction module is used to acquire underground geological data, construct a model by using the underground geological data, and obtain a velocity model;

[0034] A model optimization module is used to acquire seismic data, adjust model parameters of the velocity model according to the seismic data, and obtain a mixed velocity model;

[0035] A structure map construction module is used to acquire a migration time structure map, perform time-depth relationship conversion of the time structure map by using the mixed velocity model, and obtain a depth domain structure map;

[0036] A structure map correction module is used to correct the depth domain structure map according to the seismic data, and obtain a mixed velocity model structure map.

[0037] In a third aspect, the embodiment of the application provides an electronic device, which comprises:

[0038] A processor;

[0039] A memory for storing instructions executable by the processor;

[0040] The processor is configured to execute the instructions to implement the three-dimensional variable velocity mapping method according to the first aspect.

[0041] In a fourth aspect, the embodiments of the present application provide a computer readable storage medium, which stores a computer program. When the computer program is executed by a processor, the three-dimensional variable velocity mapping method according to the first aspect is implemented.

[0042] In a fifth aspect, the embodiments of the present application provide a computer program. When the computer program is executed by a processor, the three-dimensional variable velocity mapping method according to the first aspect is implemented.

[0043] Compared with the prior art, the above technical solution of the present application has the following beneficial effects:

[0044] The embodiments of the present application constrain the velocity model by using the vertical seismic profile logging data in the seismic data, obtain the vertical seismic profile constrained velocity model, update and iterate the vertical seismic profile constrained velocity model by using the layered data in the seismic data, obtain the layered iterative velocity model, constrain the layered iterative velocity model by using the pre-stack depth migration layer velocity body data in the seismic data, obtain the mixed velocity model, update and constrain the model from the vertical seismic profile time-depth curve, the geological layering and the seismic horizon and the pre-stack depth migration layer velocity body, so that the accuracy of the finally obtained mixed velocity model is further improved; calculate the corresponding velocity on the time structure map according to the time structure map and the velocity model, convert the corresponding velocity and the time structure map to obtain the depth value; construct the depth structure map according to the depth value, the corresponding velocity and the time structure map; obtain the geological horizon data, correct the depth domain structure map by using the geological horizon data to obtain the depth domain structure map, and the time-depth-velocity three-dimensional depth domain structure map can be constructed, which is convenient for obtaining accurate seismic velocity field data and further meeting the demand of three-dimensional fine interpretation structure mapping. Therefore, the three-dimensional variable velocity mapping method, device, electronic equipment, computer readable storage medium and computer program provided by the present application can improve the three-dimensional velocity modeling accuracy and the structure mapping accuracy, especially the micro-amplitude structure mapping accuracy. BRIEF DESCRIPTION OF DRAWINGS

[0045] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the following will briefly introduce the drawings needed to be used in the embodiments or prior art description. Obviously, the drawings in the following description are only some embodiments of the present application, and for those skilled in the art, other drawings can also be obtained without creative labor.

[0046] Figure 1 The working flowchart of the three-dimensional variable velocity mapping method of the embodiment one of the present application is shown;

[0047] Figure 2 The W1 well synthetic seismic record calibration graph of the embodiments of the present application is shown.

[0048] Figure 3a Figure 1 shows a three-dimensional velocity model of an embodiment of the present application;

[0049] Figure 3b Figure 2 shows a three-dimensional velocity model of an embodiment of the present application;

[0050] Figure 4a Figure 3 shows a well-to-well velocity profile of W1 and W2 wells obtained by the velocity model of an embodiment of the present application;

[0051] Figure 4b Figure 4 shows a well-to-well velocity profile of W1 and W2 wells obtained by the velocity model of an embodiment of the present application;

[0052] Figure 5a Figure 5 shows a vertical seismic profile constrained velocity model of an embodiment of the present application;

[0053] Figure 5b Figure 6 shows a layered iterative velocity model of an embodiment of the present application;

[0054] Figure 5c Figure 7 shows a depth domain layer velocity constrained velocity model of an embodiment of the present application;

[0055] Figure 6a Figure 8 shows a time-depth conversion structure map of a target layer of a certain work area obtained by the conventional velocity model of an embodiment of the present application;

[0056] Figure 6b Figure 9 shows a time-depth conversion structure map of a target layer of a certain work area obtained by the mixed velocity model of an embodiment of the present application;

[0057] Figure 7 Figure 10 shows a functional module diagram of a three-dimensional variable velocity mapping device of an embodiment of the present application;

[0058] Figure 8 Figure 11 shows a composition structure diagram of an electronic device for implementing the three-dimensional variable velocity mapping method of an embodiment of the present application. DETAILED DESCRIPTION

[0059] The present disclosure will be further described below with reference to the embodiments shown in the drawings.

[0060] It should be noted that the embodiments in the present application and the features in the embodiments can be combined with each other without conflict. The present application will be described in detail below with reference to the embodiments and in combination with the drawings.

[0061] The application proposes a three-dimensional variable velocity mapping method based on model construction, which is based on the synthesis of the fine calibration map of the seismic record, combined with the least square method, builds a three-dimensional velocity model, and constructs a three-dimensional fine interpretation map with higher accuracy.

[0062] Embodiment one

[0063] As shown in Figure 1 The application proposes a three-dimensional variable velocity mapping method, which comprises the following steps:

[0064] S1, obtain underground geological data, and perform model construction on the underground geological data to obtain a velocity model.

[0065] In the embodiment of the application, the underground geological data comprises acoustic logging data, density logging data, pre-stack depth migration seismic data volume, and structure interpretation migration time horizon data.

[0066] In the embodiment of the application, the three-dimensional velocity model is established by performing model construction on the underground geological data, and the velocity model obtained by performing velocity interpolation on the model can reduce data redundancy and improve the analysis and utilization efficiency of geological data.

[0067] In the embodiment of the application, the Kriging interpolation method is used to perform interpolation processing on the underground geological data to obtain a velocity model; since the Kriging interpolation method can adaptively adjust the interpolation result according to the density and distribution of sample points, the obtained velocity model is more accurate and has higher accuracy.

[0068] In the embodiment of the application, the model construction on the underground geological data to obtain a velocity model comprises:

[0069] generating a fine calibration map of the seismic record of the target layer according to the underground geological data;

[0070] performing time-depth correlation picking on the fine calibration map of the seismic record of the target layer to obtain a time-depth relationship curve;

[0071] performing interpolation processing on the time-depth relationship curve by using the Kriging interpolation method to obtain a velocity model of the time-depth curve.

[0072] In the embodiment of the application, the grid method can be used to perform time-depth correlation picking to obtain a fine calibration map of the seismic record, and since the grid method can accurately obtain full-wavefield seismic records, the fine calibration map of the seismic record of the target layer obtained has higher accuracy.

[0073] In another optional embodiment of the present application, ray tracing method can also be used to pick up the time-depth correlation, to obtain the fine calibration map of the seismic record, since the ray tracing method only extracts part of the information of the wave field, thereby reducing the calculation amount of the synthesis, and therefore the efficiency of the fine calibration map of the seismic record of the target layer is higher.

[0074] In the embodiment of the present application, the time-depth relationship curve can be obtained by using the software calculation method, since the software and algorithm directly obtain the time-depth relationship curve according to the acoustic integral, and therefore the accuracy of the obtained time-depth curve is higher.

[0075] In another optional embodiment of the present application, the time-depth relationship curve can also be obtained by using the contrast method, since the contrast method compares the fine calibration map of the seismic record with the underground geological data to obtain the correspondence between the time axis and the depth range, and then obtains the time-depth relationship, and therefore the efficiency of the time-depth relationship curve is higher.

[0076] In the embodiment of the present application, in the step of obtaining the fine calibration map of the seismic record of the target layer according to the underground geological data, the fine calibration map of the seismic record includes the time axis, the seismic waveform, the epicenter distance, the acoustic wave, and the density data, as shown in Figure 2 , which is the fine calibration map of the seismic record generated by the underground geological data of well W1 in a certain work area.

[0077] In the embodiment of the present application, the time-depth relationship curve can be obtained by generating the fine calibration map of the seismic record of the target layer, and the time-depth curve with high accuracy can be quickly obtained, the model can be obtained by interpolation according to the time-depth curve, and the efficiency of the construction of the velocity model can be improved.

[0078] S2, obtaining seismic data, adjusting the model parameters of the velocity model according to the seismic data, and obtaining a mixed velocity model.

[0079] In the embodiment of the present application, the seismic data refers to the vertical seismic profile logging data, the logging layering data, the pre-stack depth migration layer velocity body data, the seismic horizon data, the depth domain horizon data, and the depth domain imaging structure.

[0080] In the embodiment of the present application, the model updating of the velocity model from the aspect of the seismic data can make the accuracy of the finally obtained mixed velocity model higher.

[0081] In the embodiment of the present application, the adjusting of the model parameters of the velocity model according to the seismic data to obtain the mixed velocity model comprises:

[0082] The vertical seismic profile logging data in the seismic data is used to constrain the parameters of the velocity model, to obtain a vertical seismic profile constrained velocity model.

[0083] updating and iterating the vertical seismic profile constraint velocity model by using the layered data in the seismic data, to obtain a layered iterative velocity model;

[0084] model-constraining the layered iterative velocity model by using a small least square tomography according to the pre-stack depth migration layer velocity body data in the seismic data, to obtain a mixed velocity model.

[0085] In the embodiment of the present application, the least square method can be used to constrain the velocity model parameters, to obtain a vertical seismic profile constraint velocity model, since the least square method does not increase new features and can also solve the error between the vertical seismic profile logging time-depth and the time-depth relationship curve, so that the precision of the vertical seismic profile constraint velocity model is higher.

[0086] In the embodiment of the present application, the final obtained mixed velocity model includes time, velocity and layered three aspects, referring to FIG. 3, it is a 3D mixed velocity model example of a certain work area established by the method, wherein, Figure 3a is a three-dimensional velocity whole model, Figure 3b is a three-dimensional velocity body internal and target layer velocity model.

[0087] In the embodiment of the present application, the velocity model is parameter-constrained by using the vertical seismic profile logging data in the seismic data, to obtain a vertical seismic profile constraint velocity model, comprising:

[0088] establishing a fitting function according to the vertical seismic profile logging data and the velocity model;

[0089] calculating the residual sum of squares of the fitting function and the vertical seismic profile logging data;

[0090] calculating a constraint velocity model by using the least square method according to the residual sum of squares;

[0091] performing the above steps on the constraint velocity model for a preset number of times to obtain a vertical seismic profile constraint velocity model.

[0092] In the embodiment of the present application, the vertical seismic profile constraint velocity model obtained by constraining the velocity model by the profile logging data can eliminate part of the abnormal data in the velocity model and improve the precision of the model, referring to FIG. 4, it is a comparison diagram of the velocity model and the vertical seismic profile time-depth curve constraint velocity model of W1 and W2 wells in a certain work area based on the time-depth curve, wherein, Figure 4a is a W1 well and W2 well connecting well velocity profile diagram obtained by the velocity model, Figure 4bThe well connecting velocity profile of W1 well and W2 well based on the vertical seismic profile constrained velocity model is obtained, and it can be seen from the figure that after the time-depth curve of the vertical seismic profile is constrained, the velocity model anomaly is eliminated and the accuracy is improved.

[0093] In the embodiment of the present application, the layered iterative velocity model is obtained by updating and iterating the vertical seismic profile constrained velocity model using the layered data in the seismic data.

[0094] The well list data is obtained, and the seismic wave velocity of the well list data is calculated using the vertical seismic profile constrained model;

[0095] The time-depth conversion is performed on the seismic wave data to obtain the simulated depth;

[0096] The error between the depth of the geological layer of the well list data and the simulated depth is calculated;

[0097] The well with an error value less than a preset constant is selected as a training well, and the vertical seismic profile constrained velocity model is updated using the error corresponding to the training well to obtain a layered iterative velocity model.

[0098] In the embodiment of the present application, the well list refers to a list of drilled wells used to obtain seismic data and underground geological data, and the well in the well list refers to a geological exploration well drilled on the earth's surface.

[0099] In another optional embodiment of the present application, the layered iterative velocity model is constrained by the prestack depth migration layer velocity body data in the seismic data using the least square tomography method to obtain a hybrid velocity model, which comprises:

[0100] Forward simulation is performed using the layered iterative velocity model to obtain forward data;

[0101] The prestack depth migration layer velocity body data is subtracted from the forward data to obtain a residual error;

[0102] The layered iterative model is updated using the residual error according to the inversion method to obtain a hybrid velocity model.

[0103] In the embodiment of the present application, the layered iterative model can be updated using regularization inversion, and since regularization inversion can effectively control noise interference in the inversion process, when the layered iterative model is updated using regularization inversion, the stability and reliability of updating the layered iterative model can be improved.

[0104] In another optional embodiment of the present application, the layered iterative model can also be updated by using a quasi-Newton method; since the quasi-Newton method is used to update parameters by using a second-order derivative, the calculation process is simple, and therefore when the layered iterative model is updated by using the quasi-Newton method, the efficiency of updating the layered iterative model can be improved.

[0105] In the embodiment of the present application, the layered iterative velocity model is model-constrained by using the vertical seismic profile logging data, the layered data and the pre-stack depth migration layer velocity body data, and the mixed velocity model can gradually improve the accuracy of the model; as shown in FIG. 5, it is a well tie velocity profile comparison chart of a certain work area, wherein, Figure 5a the vertical seismic profile constraint velocity model, Figure 5b the layered iterative velocity model, Figure 5c the depth domain layer velocity constraint velocity model, and

[0106] In the embodiment of the present application, the model parameter adjustment of the velocity model according to the seismic data can update and constrain the velocity model from the aspect of seismic data, and improve the accuracy of the final mixed model.

[0107] S3, a time structure map is obtained, the time structure map is converted by using the mixed velocity model to obtain a depth domain structure map.

[0108] In the embodiment of the present application, the time structure map is converted by using the mixed velocity model to obtain the corresponding velocity corresponding to the time on the time structure map, and then more accurate depth values are obtained through time-depth conversion, and the accuracy of the depth domain structure map is improved.

[0109] In the embodiment of the present application, the time structure map is converted by using the mixed velocity model to obtain the depth domain structure map, including:

[0110] the corresponding velocity on the time structure map is calculated according to the time structure map and the velocity model;

[0111] the corresponding velocity on the time structure map and the time structure map are converted to obtain a depth value;

[0112] the depth structure map is constructed according to the corresponding velocity on the time structure map and the time structure map;

[0113] geological horizon data is obtained, the depth domain structure map is corrected by using the geological horizon data to obtain a depth domain structure map.

[0114] In the embodiment of the present application, the depth domain structure map obtained by using the mixed velocity model can improve the accuracy of three-dimensional structure interpretation, as shown in FIG. 6, Figure 6a is a conventional velocity model time-depth conversion structure map of a target layer in a certain work area, Figure 6b is a mixed velocity model time-depth conversion structure map, as shown in the figure, the depth domain structure map obtained by using the mixed velocity model to convert the time-depth relationship of the time structure map can improve the accuracy of the structure map.

[0115] In the embodiment of the present application, the depth domain structure map obtained by using the mixed velocity model to convert the time-depth relationship of the time structure map can improve the accuracy of the depth domain structure map, and thus can meet the actual exploration and development needs.

[0116] S4, according to the seismic data data, the result of the depth domain structure map is corrected, and a mixed velocity model structure map is obtained.

[0117] In the embodiment of the present application, the result of the depth domain structure map is corrected by using the seismic data data, which can improve the reliability of the finally obtained mixed velocity model structure map.

[0118] In the embodiment of the present application, the result of the depth domain structure map is corrected according to the seismic data data, and a mixed velocity model structure map is obtained, comprising:

[0119] Using the well tie profile analysis method to construct the seismic data data and the time structure map, a time-depth structure map is obtained;

[0120] Projecting the time-depth structure map onto the depth domain structure map obtains a projection result;

[0121] Matching and analyzing the projection result and the depth domain structure map obtains a mixed velocity model structure map.

[0122] In the embodiment of the present application, by correcting the result of the depth domain structure map, the comparison and matching of the seismic data data and the depth domain structure map can be realized, and then the data difference is obtained, the reliability of the depth domain structure map is verified, and finally the accurate and reliable mixed velocity model structure map is obtained.

[0123] Embodiment two

[0124] As shown in FIG. Figure 7 The embodiment also provides a functional module diagram of a three-dimensional variable speed mapping device.

[0125] The three-dimensional variable velocity mapping device 100 described in the embodiment can be installed in a device. According to the implemented functions, the three-dimensional variable velocity mapping device 100 can include a model construction module 101, a model optimization module 102, a model optimization module 103, and a structure map construction module 104. The modules described in the present application can also be referred to as units, which refer to a series of computer program segments that can be executed by a device processor and can complete a fixed function, which are stored in the memory of the device.

[0126] In the embodiment, the functions of each module / unit are as follows:

[0127] The model construction module 101 is configured to obtain underground geological data, construct a model based on the underground geological data, and obtain a velocity model.

[0128] The model optimization module 102 is configured to obtain seismic data, adjust model parameters of the velocity model based on the seismic data, and obtain a hybrid velocity model.

[0129] The structure map construction module 103 is configured to obtain a migration time structure map, convert the time structure map to a depth domain structure map using the hybrid velocity model.

[0130] The structure map correction module 104 is configured to correct the depth domain structure map based on the seismic data, and obtain a hybrid velocity model structure map.

[0131] In detail, each module in the three-dimensional variable velocity mapping device 100 in the embodiment uses the same technical means as the three-dimensional variable velocity mapping method described in Embodiment One and Embodiment Two, and can produce the same technical effects, which will not be described here.

[0132] Embodiment Three

[0133] As shown in Figure 8 The embodiment also provides a computer electronic device, which can include a processor 10, a memory 11, a communication bus 12, and a communication interface 13, and can also include a computer program stored in the memory 11 and executable on the processor 10, such as a product recommendation program based on user portrait.

[0134] The processor 10 may, in some embodiments, be composed of integrated circuits, for example, a single packaged integrated circuit, or a plurality of packaged integrated circuits of the same or different functions, including one or more combinations of central processing units (CPUs), microprocessors, digital processing chips, graphics processors, and various control chips. The processor 10 is the control unit of the electronic device, connects various components of the electronic device through various interfaces and lines, executes programs or modules stored in the memory 11 (for example, a product recommendation program based on user portrait), and calls data stored in the memory 11 to perform various functions of the electronic device and process data.

[0135] The memory 11 includes at least one type of readable storage medium, including flash memories, mobile hard disks, multimedia cards, card-type memories (for example, SD or DX memories, etc.), magnetic memories, magnetic disks, optical disks, etc. The memory 11 may, in some embodiments, be an internal storage unit of the electronic device, for example, a mobile hard disk of the electronic device. The memory 11 may, in other embodiments, also be an external storage device of the electronic device, for example, a plug-in mobile hard disk, a smart media card (SMC), a secure digital (SD) card, a flash card, etc. Further, the memory 11 may include both an internal storage unit and an external storage device of the electronic device. The memory 11 can be used not only to store application software and various data installed in the electronic device, for example, codes of a product recommendation program based on user portrait, but also to temporarily store data that has been or will be output.

[0136] The communication bus 12 can be a peripheral component interconnect (PCI) bus or an extended industry standard architecture (EISA) bus, etc. The bus can be divided into an address bus, a data bus, a control bus, etc. The bus is configured to realize the connection and communication between the memory 11 and at least one processor 10, etc.

[0137] The communication interface 13 is used for communication between the electronic device and other devices, including a network interface and a user interface. Optionally, the network interface can include a wired interface and / or a wireless interface (such as a WI-FI interface, a Bluetooth interface, etc.), which is usually used to establish a communication connection between the electronic device and other electronic devices. The user interface can be a display, an input unit (such as a keyboard), and optionally, the user interface can also be a standard wired interface, a wireless interface. Optionally, in some embodiments, the display can be an LED display, a liquid crystal display, a touch liquid crystal display, an OLED (Organic Light-Emitting Diode) touch, etc. Among them, the display can also be appropriately called a display screen or a display unit, which is used to display information processed in the electronic device and to display a visualized user interface.

[0138] Only the electronic device with components is shown in the figure, and those skilled in the art can understand that the structure shown in the figure does not constitute a limitation on the electronic device, and can include fewer or more components than the figure, or combine certain components, or different component arrangements.

[0139] For example, although not shown, the electronic device can also include a power supply (such as a battery) for powering each component. Preferably, the power supply can be logically connected to the at least one processor 10 through a power management device, so as to realize functions such as charge management, discharge management, and power consumption management through the power management device. The power supply can also include one or more direct current or alternating current power supplies, recharging devices, power supply fault detection circuits, power supply converters or inverters, power supply status indicators, and any other components. The electronic device can also include various sensors, Bluetooth modules, Wi-Fi modules, etc., which are not described here.

[0140] It should be understood that the embodiments are only for illustration and are not limited in the scope of the patent application by this structure.

[0141] The product recommendation program based on the user portrait stored in the memory 11 in the electronic device is a combination of multiple instructions, which, when running in the processor 10, can realize:

[0142] Obtaining underground geological data, constructing a model based on the underground geological data, and obtaining a velocity model;

[0143] Obtaining seismic data, adjusting model parameters of the velocity model based on the seismic data, and obtaining a hybrid velocity model;

[0144] An offset time structure map is acquired, and the time structure map is converted into a depth domain structure map by using the mixed velocity model;

[0145] The depth domain structure map is corrected according to the seismic data to obtain a mixed velocity model structure map.

[0146] Specifically, the specific implementation method of the processor 10 to the above instructions can refer to the description of the related steps in the corresponding embodiment of the drawings, which will not be described here.

[0147] Further, the modules / units of the electronic device, if realized in the form of software function units and sold or used as independent products, can be stored in a computer readable storage medium. The computer readable storage medium can be volatile or non-volatile. For example, the computer readable medium can include any entity or device capable of carrying the computer program code, recording medium, U disk, mobile hard disk, magnetic disk, optical disk, computer memory, read-only memory (ROM).

[0148] Embodiment four

[0149] The embodiment provides a storage medium storing a computer program, and the computer program is executed by a processor to realize the steps of the three-dimensional variable speed mapping method described above.

[0150] These program codes can also be loaded into a computer or other programmable data processing device to make a series of operation steps executed on the computer or other programmable device to produce a computer implemented process, so that the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in the flow Figure One In one flow or multiple flows.

[0151] The storage medium includes permanent and non-permanent, removable and non-removable media, and can be realized by any method or technology. Information can be computer readable instructions, data structures, program modules or other data. Examples of storage media can include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassette, magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information accessible by a computing device.

[0152] Embodiment five

[0153] An embodiment of the present invention provides a computer program, which, when executed by a processor, implements the steps of the three-dimensional variable speed mapping method as described in the first aspect.

[0154] These computer-readable program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, thereby producing a machine such that when these instructions are executed by the processor of the computer or other programmable data processing device, a device is generated that implements the functions / actions specified in one or more blocks in the flowchart and / or block diagram. These computer programs can also be stored in a computer-readable storage medium. These computer programs cause the computer, programmable data processing device, and / or other equipment to operate in a specific manner. Thus, the computer-readable medium storing the computer program comprises an article of manufacture that includes instructions for implementing various aspects of the functions / actions specified in one or more blocks in the flowchart and / or block diagram.

[0155] It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit the exemplary embodiments according to the present application. When the terms "include" and / or "include" are used in this specification, they indicate the presence of features, steps, operations, devices, components and / or their combinations.

[0156] It is to be understood that the terms so used are interchangeable under appropriate circumstances such that the embodiments of the application described herein are, for example, capable of operation in sequences other than those illustrated or described herein.

[0157] In the several embodiments provided by the present invention, it should be understood that the disclosed devices, apparatuses, and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative. For example, the module division is merely a logical function division, and other division methods may be used in actual implementation.

[0158] The modules described as separate components may or may not be physically separate, and the components shown as modules may or may not be physical units, that is, they may be located in one place or distributed across multiple network elements. Some or all of the modules may be selected to achieve the purpose of the solution of this embodiment according to actual needs.

[0159] In addition, the functional modules in various embodiments of the present invention may be integrated into a single processing unit, each unit may exist physically separately, or two or more units may be integrated into a single unit. The aforementioned integrated units may be implemented in the form of hardware or hardware plus software functional modules.

[0160] It will be obvious to a person skilled in the art that the application is not limited to the details of the above-described exemplary embodiments but can be implemented in other embodiments without departing from the scope of the application.

[0161] The embodiments should therefore be considered in all respects as illustrative and not restrictive, the scope of the application being indicated by the appended claims rather than by the description given above, so that all changes coming within the meaning and equivalency range of the claims are intended to be embraced therein. Any reference signs in the claims should not be construed as limiting the scope of the claims.

[0162] Embodiments of the present application can acquire and process related data based on artificial intelligence technology. Wherein, artificial intelligence (AI) is to use digital computers or digital computer controlled machine to simulate, extend and expand human intelligence, to perceive the environment, obtain knowledge and use knowledge to obtain the best results.

[0163] Furthermore, it is clear that the word "comprising" does not exclude other elements or steps, and the singular does not exclude the plural. A plurality of units or devices can be implemented by one unit or device by means of software or hardware. The terms first, second and the like do not denote any order, quantity, number or importance, but are used to identify different components.

[0164] Finally, it should be noted that the above embodiments are merely used to illustrate the technical solutions of the present application, rather than limit the present application. Even though the present application has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that the technical solutions of the present application can be modified or replaced equivalently without departing from the spirit and scope of the present application.

Claims

1. A three-dimensional variable speed mapping method, characterized in that: The method comprises: Acquiring underground geological data, constructing a model for the underground geological data, and obtaining a velocity model; Acquiring seismic data, and adjusting model parameters of the velocity model according to the seismic data to obtain a hybrid velocity model; Acquire a migration time structure map, and use the mixed velocity model to perform time-depth relationship conversion on the time structure map to obtain a depth domain structure map; The depth domain structural map is corrected according to the seismic data to obtain a mixed velocity model structural map.

2. The three-dimensional variable speed mapping method according to claim 1, wherein: The step of constructing a model for the underground geological data to obtain a velocity model includes: generating a fine calibration map of seismic records of a target layer according to the underground geological data; Performing time-depth correlation picking on the fine calibration map of the target layer seismic record to obtain a time-depth relationship curve; The time-depth relationship curve is interpolated using the Kriging interpolation method to obtain a velocity model of the time-depth curve.

3. The three-dimensional variable speed mapping method according to claim 1, wherein: The adjusting the model parameters of the velocity model according to the seismic data to obtain a hybrid velocity model includes: Using vertical seismic profile logging data in the seismic data to constrain parameters of the velocity model, to obtain a vertical seismic profile constrained velocity model; Using the layered data in the seismic data to iterate the vertical seismic profile constrained velocity model, to obtain a layered iterative velocity model; The layered iterative velocity model is constrained by using a least squares tomography method according to the pre-stack depth migration layer velocity volume data in the seismic data to obtain a mixed velocity model.

4. The three-dimensional speed-varying mapping method according to claim 3, wherein: The method of using the vertical seismic profile logging data in the seismic data to constrain the parameters of the velocity model to obtain a vertical seismic profile constrained velocity model includes: Establishing a fitting function based on the vertical seismic profile logging data and the velocity model; Calculating the residual sum of squares of the fitting function and the vertical seismic profile logging data; A constrained velocity model is obtained by calculating the residual square sum using the least square method; The above steps are performed for a preset number of times on the constrained velocity model to obtain a vertical seismic profile constrained velocity model.

5. The three-dimensional speed-varying mapping method according to claim 3, wherein: The method of updating and iterating the vertical seismic profile constrained velocity model using the layered data in the seismic data to obtain a layered iterative velocity model includes: Acquire well list data, and calculate seismic wave velocity of the well list data using the vertical seismic profile constraint model; Performing time-depth conversion on the seismic wave data to obtain a simulated depth; Calculating the error between the depth of the geological layer of the well list data and the simulated depth; Wells with error values ​​less than a preset constant are selected as training wells, and the vertical seismic profile constrained velocity model is updated using the errors corresponding to the training wells to obtain a layered iterative velocity model.

6. The three-dimensional speed-varying mapping method according to claim 1, wherein: The step of converting the time-depth relationship of the time structure map using the mixed velocity model to obtain a depth domain structure map includes: Calculating a corresponding speed on the time structure diagram according to the time structure diagram and the speed model; Performing time-depth relationship conversion on the corresponding speed on the time structure graph and the time structure graph to obtain a depth value; constructing a depth structure map according to the corresponding speed on the time structure map and the time structure map; Acquire geological stratum data, and use the geological stratum data to perform geological stratum correction on the depth domain structural map to obtain the depth domain structural map.

7. A three-dimensional variable speed mapping device, characterized in that: The device comprises: A model building module is used to acquire underground geological data, build a model for the underground geological data, and obtain a velocity model; A model optimization module is used to obtain seismic data, adjust model parameters of the velocity model according to the seismic data, and obtain a hybrid velocity model; a structure map construction module, configured to obtain a migration time structure map, and perform time-depth relationship conversion on the time structure map using the mixed velocity model to obtain a depth domain structure map; The structural map correction module is used to correct the depth domain structural map according to the seismic data to obtain a mixed velocity model structural map.

8. An electronic device, characterized in that: include: processor; a memory for storing instructions executable by the processor; The processor is configured to execute the instructions to implement the steps of the method for three-dimensional speed-varying mapping according to any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that A computer program is stored thereon, and when the program is executed by a processor, the steps of the three-dimensional variable speed mapping method according to any one of claims 1 to 6 are implemented.

10. A computer program, characterized in that When the program is executed by a processor, the steps of the three-dimensional variable speed mapping method according to any one of claims 1 to 6 are implemented.

Images