Method and device for establishing initial velocity model of well-seismic joint full waveform inversion

Abstract

The present invention provides a method and device for establishing an initial velocity model of well-seismic combined full waveform inversion. The method comprises: acquiring seismic data and acoustic logging data of a study area; acquiring the seismic velocity between wells, the seismic velocity of well points and the corresponding logging velocity according to the seismic data and the acoustic logging data; correcting the seismic velocity between wells by using the logging velocity ratio corresponding to the seismic velocity of the well points; and constructing an initial velocity model according to a transformation matrix and the corrected seismic velocity between wells, wherein the lateral variation trend of the seismic velocity can be maintained when the inter-well velocity is interpolated by the transformation matrix, the logging and seismic velocity information are integrated, and a high-precision initial velocity model of waveform inversion is established.

Description

Technical Field

[0001] The present invention relates to the field of seismic exploration technology, and in particular to a method and device for establishing an initial velocity model of well-seismic combined full waveform inversion. Background Art

[0002] Full waveform inversion is one of the research hotspots in the field of seismic exploration technology. This technology combines the kinematic and dynamic information in the seismic wave field to invert high-precision velocity, density, and anisotropy models of underground media, thereby predicting structural, reservoir, and fluid information in complex geological settings.

[0003] Full waveform inversion uses an optimization algorithm to iteratively solve the residual of the observed data and the wave equation forward data to solve the underground medium parameters. Since there is a nonlinear relationship between the underground medium parameters and the seismic wave field, this method is a mathematically ill-posed problem with multiple solutions and convergence problems. At present, the common method to solve nonlinear inversion is to use the Born approximation to transform the nonlinear problem into a linear problem, and then use the local optimization method to iteratively invert the velocity model. In the theory of full waveform inversion, to meet the Born approximation, the error between the travel time of the wave in the forward data of the initial model and the travel time of the actual observed data must be less than half a period in the time domain and less than half a wavelength in the depth domain. If the initial model is not accurate enough, the inversion objective function will converge to the local minimum and the "cycle jump" phenomenon will occur. Therefore, initial velocity modeling is a core issue of full waveform inversion. Usually, the initial velocity model of full waveform inversion is obtained using ray tomography, offset velocity analysis or sonic logging interpolation. Ray tomography and offset velocity analysis use the travel time information and coaxial curvature of seismic waves, and the velocity model obtained is of low accuracy and can only reflect relative velocity changes. Velocity modeling using acoustic logging velocity interpolation usually uses layer control, and linear or bilinear interpolation is used within the layer, which cannot reflect the lateral change of velocity. Moreover, after the velocity model is established using the above method, it is necessary to perform overall spatial smoothing according to the full waveform inversion multi-scale framework, which will also introduce velocity errors. Summary of the invention

[0004] In view of the problems in the prior art, the present invention provides a method and device for establishing an initial velocity model of well-seismic combined full waveform inversion, an electronic device and a computer-readable storage medium, which can at least partially solve the problems in the prior art.

[0005] In order to achieve the above object, the present invention adopts the following technical solution:

[0006] In a first aspect, a method for establishing an initial velocity model by well-seismic combined full waveform inversion is provided, comprising:

[0007] Obtain seismic data and sonic logging data for the study area;

[0008] Acquire the seismic velocity between wells, the seismic velocity of well points and the corresponding logging velocity according to the seismic data and the acoustic logging data;

[0009] Correcting the seismic velocity between the wells using the logging velocity ratio corresponding to the seismic velocity of the well point;

[0010] The initial velocity model is constructed based on the transformation matrix and the corrected inter-well seismic velocity.

[0011] Furthermore, the logging data includes: the distance between wells, and the transformation matrix is:

[0012]

[0013] Among them, d is the distance between the wells, g1 is the logging velocity difference between the two wells, g2 is the differential velocity between the wells, v1 is the seismic velocity at the corresponding position of the first well, and V1 is the logging velocity at the corresponding position of the first well.

[0014] Furthermore, the seismic data includes: seismic velocity field and stratigraphic data;

[0015] The obtaining of the inter-well seismic velocity, the seismic velocity of the well point and the corresponding logging velocity according to the seismic data and the acoustic logging data comprises:

[0016] performing median filtering on the seismic velocity field;

[0017] Matching the seismic data and the acoustic logging data using well-seismic calibration technology;

[0018] Extracting the seismic velocity between wells and the seismic velocity of the well point in the seismic velocity field after median filtering by using the layer data and the well positions in the acoustic logging data;

[0019] The logging velocity corresponding to each seismic velocity is extracted from the acoustic logging data after well-seismic matching.

[0020] Furthermore, the extracting of the seismic velocity between wells and the seismic velocity of the well point in the seismic velocity field after median filtering by using the well positions in the horizon data and the acoustic logging data comprises:

[0021] Establishing a stratigraphic framework model using the stratigraphic data;

[0022] The stratigraphic framework model and the well positions in the acoustic logging data are used to extract the seismic velocities between wells and the seismic velocities of the well points in the seismic velocity field after median filtering.

[0023] Furthermore, the initial velocity model is constructed according to the transformation matrix and the corrected inter-well seismic velocity, including:

[0024] Multiplying the transformation matrix with the corrected inter-well seismic velocity to obtain inter-well differential velocity;

[0025] The initial velocity model is constructed according to the differential velocities between the wells.

[0026] In a second aspect, a device for establishing an initial velocity model of well-seismic combined full waveform inversion is provided, comprising:

[0027] Data acquisition module, which acquires seismic data and acoustic logging data of the study area;

[0028] A velocity extraction module, for acquiring the seismic velocity between wells, the seismic velocity of well points and the corresponding logging velocity according to the seismic data and the acoustic logging data;

[0029] A velocity correction module, which corrects the seismic velocity between the wells using the logging velocity ratio corresponding to the seismic velocity of the well point;

[0030] The modeling module constructs an initial velocity model based on the transformation matrix and the corrected inter-well seismic velocity.

[0031] Furthermore, the logging data includes: the distance between wells, and the transformation matrix is:

[0032]

[0033] Among them, d is the distance between the wells, g1 is the logging velocity difference between the two wells, g2 is the differential velocity between the wells, v1 is the seismic velocity at the corresponding position of the first well, and V1 is the logging velocity at the corresponding position of the first well.

[0034] Furthermore, the seismic data includes: seismic velocity field and stratigraphic data;

[0035] The speed extraction module comprises:

[0036] A median filtering unit, performing median filtering on the seismic velocity field;

[0037] A well seismic calibration unit, which uses well seismic calibration technology to match the seismic data and the acoustic logging data;

[0038] A seismic velocity extraction unit, which extracts the seismic velocity between wells and the seismic velocity of the well point in the seismic velocity field after median filtering by using the layer data and the well positions in the acoustic logging data;

[0039] The logging velocity extraction unit extracts the logging velocity corresponding to each seismic velocity from the acoustic logging data after well-seismic matching.

[0040] In a third aspect, an electronic device is provided, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the steps of the method for establishing an initial velocity model of the above-mentioned well-seismic combined full waveform inversion are implemented.

[0041] In a fourth aspect, a computer-readable storage medium is provided, on which a computer program is stored, and when the computer program is executed by a processor, the steps of the method for establishing an initial velocity model of the above-mentioned well-seismic joint full waveform inversion are implemented.

[0042] The present invention provides a method and device for establishing an initial velocity model of well-seismic combined full waveform inversion. The method comprises: acquiring seismic data and acoustic logging data of a study area; acquiring the seismic velocity between wells, the seismic velocity of well points and the corresponding logging velocity according to the seismic data and the acoustic logging data; correcting the seismic velocity between wells using the logging velocity ratio corresponding to the seismic velocity of the well points; and constructing an initial velocity model according to a transformation matrix and the corrected seismic velocity between wells, wherein the lateral variation trend of the seismic velocity can be maintained when the inter-well velocity is interpolated by the transformation matrix, the logging and seismic velocity information are integrated, and a high-precision waveform inversion initial velocity model is established.

[0043] In order to make the above and other purposes, features and advantages of the present invention more obvious and easy to understand, preferred embodiments are specifically cited below and described in detail with reference to the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0044] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the prior art descriptions. Obviously, the drawings described below are some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work. In the drawings:

[0045] Figure 1 is a schematic diagram of the architecture between the server S1 and the client device B1 in an embodiment of the present invention;

[0046] Figure 2 Schematic diagram of the architecture between the server S1, the client device B1 and the database server S2 in an embodiment of the present invention;

[0047] Figure 3 It is a flow chart of a method for establishing an initial velocity model by well-seismic combined full waveform inversion in an embodiment of the present invention;

[0048] Figure 4 A schematic diagram of the extreme values ​​of the interwell velocity before and after correction in an embodiment of the present invention is shown;

[0049] Figure 5 A schematic diagram of matrix transformation interpolation in an embodiment of the present invention is shown;

[0050] Figure 6 The specific steps of step S200 in the embodiment of the present invention are shown;

[0051] Figure 7 The specific steps of step S230 in the embodiment of the present invention are shown;

[0052] Figure 8 A schematic diagram of an inter-well formation framework in an embodiment of the present invention is shown;

[0053] Fig. 9 The specific steps of step S400 in the embodiment of the present invention are shown;

[0054] Fig.10 It is a structural block diagram of a device for establishing an initial velocity model by well-seismic combined full waveform inversion in an embodiment of the present invention;

[0055] Fig.11 The specific structure of the speed extraction module in the embodiment of the present invention is shown;

[0056] Fig.12 4 is a structural diagram of an electronic device according to an embodiment of the present invention. DETAILED DESCRIPTION

[0057] In order to enable those skilled in the art to better understand the solution of the present application, the technical solution in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of the present application.

[0058] Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0059] It should be noted that the terms "including" and "having" and any variations thereof in the specification and claims of the present application and the above-mentioned drawings are intended to cover non-exclusive inclusions. For example, a process, method, system, product or apparatus comprising a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or apparatuses.

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

[0061] Velocity modeling using acoustic logging velocity interpolation usually uses layer control, and linear or bilinear interpolation is used within the layer, which cannot reflect the lateral change of velocity. Moreover, after the velocity model is established using the above method, it is necessary to perform overall spatial smoothing according to the full waveform inversion multi-scale framework, which will also introduce velocity errors.

[0062] In view of this, the present application provides a device for establishing an initial velocity model for well-seismic joint full waveform inversion, which can be a server S1, see Figure 1 , the server S1 can be connected to at least one client device B1 in communication, the client device B1 can send the seismic data and acoustic logging data of the study area to the server S1, and the server S1 can receive the seismic data and acoustic logging data of the study area online. The server S1 can pre-process the acquired seismic data and acoustic logging data of the study area online or offline, obtain the seismic velocity between wells, the seismic velocity of the well points and the corresponding logging velocity according to the seismic data and the acoustic logging data; correct the seismic velocity between wells using the logging velocity ratio corresponding to the seismic velocity of the well points; and construct an initial velocity model according to the transformation matrix and the corrected seismic velocity between wells. Then, the server S1 can send the initial velocity model to the client device B1 online. The client device B1 can receive the initial velocity model online.

[0063] Also, see Figure 2 The server S1 may also be connected to at least one database server S2 for communication, and the database server S2 is used to store the transformation matrix. The database server S2 sends the historical transformation matrix to the server S1 online, and the server S1 may receive the transformation matrix online.

[0064] It is understandable that the client device B1 may include a smart phone, a tablet electronic device, a network set-top box, a portable computer, a desktop computer, a personal digital assistant (PDA), a vehicle-mounted device, a smart wearable device, etc. Among them, the smart wearable device may include smart glasses, a smart watch, a smart bracelet, etc.

[0065] In practical applications, the part of establishing the initial velocity model of well-seismic joint full waveform inversion can be executed on the server S1 side as described above, that is, Figure 1 The architecture shown in the figure can also be that all operations are completed in the client device B1, and the client device B1 can be directly connected to the database server S2 for communication. The specific selection can be made according to the processing capability of the client device B1 and the limitations of the user's usage scenario. This application does not limit this. If all operations are completed in the client device B1, the client device B1 can also include a processor for performing specific processing of establishing an initial velocity model for well-seismic joint full waveform inversion.

[0066] The server and the client device may communicate using any suitable network protocol, including network protocols that have not yet been developed on the date of filing this application. The network protocol may include, for example, TCP / IP protocol, UDP / IP protocol, HTTP protocol, HTTPS protocol, etc. Of course, the network protocol may also include, for example, RPC protocol (Remote Procedure Call Protocol) and REST protocol (Representational State Transfer) used on top of the above protocols.

[0067] Figure 3 is a flow chart of a method for establishing an initial velocity model by well-seismic combined full waveform inversion in an embodiment of the present invention; Figure 3 As shown, the method for establishing the initial velocity model of the well-seismic combined full waveform inversion may include the following contents:

[0068] Step S100: Acquire seismic data and acoustic logging data of the study area;

[0069] Among them, well logging data includes: distance between wells, etc. Seismic data includes seismic velocity field and layer data, etc.

[0070] Step S200: Acquire the seismic velocity between wells, the seismic velocity of well points and the corresponding logging velocity according to the seismic data and the acoustic logging data;

[0071] Step S300: Correcting the seismic velocity between the wells using the logging velocity ratio corresponding to the seismic velocity of the well point; wherein: Figure 4 A schematic diagram of the extreme values ​​of the interwell velocity before and after correction in an embodiment of the present invention is shown;

[0072] Step S400: constructing an initial velocity model according to the transformation matrix and the corrected inter-well seismic velocity.

[0073] The transformation matrix is:

[0074]

[0075] Among them, d is the distance between the wells, g1 is the logging velocity difference between the two wells, g2 is the differential velocity between the wells, v1 is the seismic velocity at the corresponding position of the first well, and V1 is the logging velocity at the corresponding position of the first well.

[0076] The embodiment of the present invention is applicable to well-seismic joint velocity modeling of full waveform inversion, and provides a transformation matrix to add the lateral variation trend of seismic velocity to the inter-well velocity interpolation. Figure 5 The schematic diagram of matrix transformation interpolation in the embodiment of the present invention is shown to realize well-seismic joint velocity modeling. The model test results show the effectiveness of this method.

[0077] In an alternative embodiment, see Figure 6 , the step S200 may include the following contents:

[0078] Step S210: performing median filtering on the seismic velocity field;

[0079] Among them, the singular values ​​in the seismic velocity field are removed by median filtering.

[0080] Step S220: matching the seismic data and the acoustic logging data using well seismic calibration technology;

[0081] Step S230: extracting the seismic velocity between wells and the seismic velocity of the well point in the seismic velocity field after median filtering by using the layer data and the well positions in the acoustic logging data;

[0082] Step S240: extracting the logging velocity corresponding to each seismic velocity from the acoustic logging data after well-seismic matching.

[0083] In an alternative embodiment, participating Figure 7 , the step S230 may include the following contents:

[0084] Step S231: Establish a stratigraphic framework model using the stratigraphic data.

[0085] in, Figure 8 A schematic diagram of an inter-well formation framework in an embodiment of the present invention is shown.

[0086] Step S232: extracting the seismic velocity between wells and the seismic velocity of well points in the seismic velocity field after median filtering by using the stratigraphic framework model and the well positions in the acoustic logging data.

[0087] Among them, each seismic velocity has a corresponding logging velocity.

[0088] In an alternative embodiment, see Fig. 9 , the step S400 may include the following contents:

[0089] Step S410: multiplying the transformation matrix by the corrected inter-well seismic velocity to obtain inter-well differential velocity;

[0090] Step S420: constructing the initial velocity model according to the differential velocities between the wells.

[0091] Based on the same inventive concept, the embodiments of the present application also provide a device for establishing an initial velocity model of well-seismic combined full waveform inversion, which can be used to implement the method described in the above embodiments, as described in the following embodiments. Since the principle of solving the problem by the device for establishing an initial velocity model of well-seismic combined full waveform inversion is similar to that of the above methods, the implementation of the device for establishing an initial velocity model of well-seismic combined full waveform inversion can refer to the implementation of the above methods, and the repeated parts will not be repeated. As used below, the term "unit" or "module" can be a combination of software and / or hardware that implements a predetermined function. Although the device described in the following embodiments is preferably implemented in software, the implementation of hardware, or a combination of software and hardware, is also possible and conceived.

[0092] Fig.10 1 is a structural block diagram of a device for establishing an initial velocity model by well-seismic combined full waveform inversion in an embodiment of the present invention. Fig.10 As shown, the apparatus for establishing an initial velocity model by well-seismic combined full waveform inversion specifically comprises: a data acquisition module 10 , a velocity extraction module 20 , a velocity correction module 30 and a modeling module 40 .

[0093] Data acquisition module 10, to acquire seismic data and acoustic logging data of the study area;

[0094] A velocity extraction module 20, for acquiring the seismic velocity between wells, the seismic velocity of well points and the corresponding logging velocity according to the seismic data and the acoustic logging data;

[0095] A velocity correction module 30, which corrects the seismic velocity between the wells using the logging velocity ratio corresponding to the seismic velocity of the well point;

[0096] The modeling module 40 constructs an initial velocity model according to the transformation matrix and the corrected inter-well seismic velocity.

[0097] In an alternative embodiment, see Fig.11 , the seismic data includes: seismic velocity field and layer data;

[0098] The velocity extraction module includes: a median filtering unit 21 , a well seismic calibration unit 22 , a seismic velocity extraction unit 23 and a well logging velocity extraction unit 24 .

[0099] A median filtering unit 21 performs median filtering on the seismic velocity field;

[0100] A well seismic calibration unit 22, which uses a well seismic calibration technology to match the seismic data and the acoustic logging data;

[0101] A seismic velocity extraction unit 23 extracts the seismic velocity between wells and the seismic velocity of the well point in the seismic velocity field after median filtering by using the layer data and the well positions in the acoustic logging data;

[0102] The logging velocity extraction unit 24 extracts the logging velocity corresponding to each seismic velocity from the acoustic logging data after well-seismic matching.

[0103] The devices, modules or units described in the above embodiments may be implemented by computer chips or entities, or by products with certain functions. A typical implementation device is an electronic device, and the electronic device may be, for example, a personal computer, a laptop computer, a cellular phone, a camera phone, a smart phone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or a combination of any of these devices.

[0104] In a typical example, the electronic device specifically includes a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, the steps of the method for establishing the initial velocity model of the well-seismic combined full waveform inversion are implemented.

[0105] Reference below Fig.12 , which shows a schematic structural diagram of an electronic device 600 suitable for implementing an embodiment of the present application.

[0106] like Fig.12As shown, electronic device 600 includes a central processing unit (CPU) 601, which can perform various appropriate operations and processes according to a program stored in a read-only memory (ROM) 602 or a program loaded from a storage portion 608 into a random access memory (RAM) 603. In RAM 603, various programs and data required for the operation of system 600 are also stored. CPU 601, ROM 602, and RAM 603 are connected to each other via a bus 604. An input / output (I / O) interface 605 is also connected to bus 604.

[0107] The following components are connected to the I / O interface 605: an input section 606 including a keyboard, a mouse, etc.; an output section 607 including a cathode ray tube (CRT), a liquid crystal display (LCD), etc., and a speaker, etc.; a storage section 608 including a hard disk, etc.; and a communication section 609 including a network interface card such as a LAN card, a modem, etc. The communication section 609 performs communication processing via a network such as the Internet. A drive 610 is also connected to the I / O interface 605 as needed. A removable medium 611, such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, etc., is installed on the drive 610 as needed, so that a computer program read therefrom is installed as needed as the storage section 608.

[0108] In particular, according to an embodiment of the present invention, the process described above with reference to the flowchart can be implemented as a computer software program. For example, an embodiment of the present invention includes a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the above-mentioned method for establishing an initial velocity model of well-seismic combined full waveform inversion.

[0109] In such an embodiment, the computer program may be downloaded and installed from a network via the communication section 609 , and / or installed from the removable medium 611 .

[0110] Computer readable media include permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. Information can be computer readable instructions, data structures, program modules or other data. Examples of computer storage media 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 technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer readable media does not include temporary computer readable media (transitory media), such as modulated data signals and carrier waves.

[0111] For the convenience of description, the above device is described in various units according to their functions. Of course, when implementing the present application, the functions of each unit can be implemented in the same or multiple software and / or hardware.

[0112] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0113] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0114] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0115] It should also be noted that the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, commodity or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, commodity or device. In the absence of more restrictions, the elements defined by the sentence "comprises a ..." do not exclude the existence of other identical elements in the process, method, commodity or device including the elements.

[0116] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that contain computer-usable program code.

[0117] The present application may be described in the general context of computer-executable instructions executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform specific tasks or implement specific abstract data types. The present application may also be practiced in distributed computing environments where tasks are performed by remote processing devices connected through a communication network. In a distributed computing environment, program modules may be located in local and remote computer storage media, including storage devices.

[0118] Each embodiment in this specification is described in a progressive manner, and the same or similar parts between the embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the system embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment.

[0119] The above is only an embodiment of the present application and is not intended to limit the present application. For those skilled in the art, the present application may have various changes and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the scope of the claims of the present application.

Claims

1. A method for establishing an initial velocity model by well-seismic combined full waveform inversion, characterized in that: include: Obtain seismic data and sonic logging data for the study area; Acquire the seismic velocity between wells, the seismic velocity of well points and the corresponding logging velocity according to the seismic data and the acoustic logging data; Correcting the seismic velocity between the wells using the logging velocity ratio corresponding to the seismic velocity of the well point; An initial velocity model is constructed according to the transformation matrix and the corrected inter-well seismic velocity; Wherein, the logging data includes: the distance between wells, and the transformation matrix is: Wherein, d is the distance between the wells, g1 is the logging velocity difference between the two wells, g2 is the differential velocity between the wells, v1 is the seismic velocity at the corresponding position of the first well, V1 is the logging velocity at the corresponding position of the first well, and x1 is the distance between the corresponding position of the extreme velocity and the corresponding position of the first well; The step of constructing an initial velocity model based on the transformation matrix and the corrected inter-well seismic velocity includes: Multiplying the transformation matrix with the corrected inter-well seismic velocity to obtain inter-well differential velocity; The initial velocity model is constructed according to the differential velocities between the wells.

2. The method for establishing an initial velocity model by well-seismic combined full waveform inversion according to claim 1 is characterized in that: The seismic data include: seismic velocity field and layer data; The obtaining of the inter-well seismic velocity, the seismic velocity of the well point and the corresponding logging velocity according to the seismic data and the acoustic logging data comprises: performing median filtering on the seismic velocity field; Matching the seismic data and the acoustic logging data using well-seismic calibration technology; Extracting the seismic velocity between wells and the seismic velocity of the well point in the seismic velocity field after median filtering by using the layer data and the well positions in the acoustic logging data; The logging velocity corresponding to each seismic velocity is extracted from the acoustic logging data after well-seismic matching.

3. The method for establishing an initial velocity model by well-seismic combined full waveform inversion according to claim 2 is characterized in that: The method of extracting the seismic velocity between wells and the seismic velocity of well points in the seismic velocity field after median filtering by using the layer data and the well positions in the acoustic logging data comprises: Establishing a stratigraphic framework model using the stratigraphic data; The stratigraphic framework model and the well positions in the acoustic logging data are used to extract the seismic velocities between wells and the seismic velocities of the well points in the seismic velocity field after median filtering.

4. A device for establishing an initial velocity model by well-seismic combined full waveform inversion, characterized in that: include: Data acquisition module, which acquires seismic data and acoustic logging data of the study area; A velocity extraction module, for acquiring the seismic velocity between wells, the seismic velocity of well points and the corresponding logging velocity according to the seismic data and the acoustic logging data; A velocity correction module, which corrects the seismic velocity between the wells using the logging velocity ratio corresponding to the seismic velocity of the well point; A modeling module, which constructs an initial velocity model based on the transformation matrix and the corrected inter-well seismic velocity; Wherein, the logging data includes: the distance between wells, and the transformation matrix is: Wherein, d is the distance between the wells, g1 is the logging velocity difference between the two wells, g2 is the differential velocity between the wells, v1 is the seismic velocity at the corresponding position of the first well, V1 is the logging velocity at the corresponding position of the first well, and x1 is the distance between the corresponding position of the extreme velocity and the corresponding position of the first well; Wherein, the modeling module is further configured as follows: Multiplying the transformation matrix with the corrected inter-well seismic velocity to obtain inter-well differential velocity; The initial velocity model is constructed according to the differential velocities between the wells.

5. The device for establishing initial velocity model of well-seismic combined full waveform inversion according to claim 4 is characterized in that: The seismic data include: seismic velocity field and layer data; The speed extraction module comprises: A median filtering unit, performing median filtering on the seismic velocity field; A well seismic calibration unit, which uses well seismic calibration technology to match the seismic data and the acoustic logging data; A seismic velocity extraction unit, which extracts the seismic velocity between wells and the seismic velocity of the well point in the seismic velocity field after median filtering by using the layer data and the well positions in the acoustic logging data; The logging velocity extraction unit extracts the logging velocity corresponding to each seismic velocity from the acoustic logging data after well-seismic matching.

6. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the program, the steps of the method for establishing an initial velocity model by well-seismic combined full waveform inversion according to any one of claims 1 to 3 are implemented.

7. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method for establishing an initial velocity model by well-seismic combined full waveform inversion as described in any one of claims 1 to 3 are implemented.

Images