CN110058306B

CN110058306B - Method, apparatus and computer storage medium for determining three-dimensional velocity volume

Info

Publication number: CN110058306B
Application number: CN201910290017.8A
Authority: CN
Inventors: 周连敏; 王晶晶; 董晓伟; 刘天鹤; 张卫江; 赵玉琳; 黄强强; 娄小娟; 段英豪
Original assignee: Petrochina Co Ltd
Current assignee: Petrochina Co Ltd
Priority date: 2019-04-11
Filing date: 2019-04-11
Publication date: 2020-10-09
Anticipated expiration: 2039-04-11
Also published as: CN110058306A

Abstract

The application discloses a method and a device for determining a three-dimensional velocity body and a computer storage medium, and belongs to the technical field of oil and gas field development. In the application, when the depth domain structure map of the research area is determined, the hierarchical data of each of the M wells is used, so that the determined depth domain structure map of the research area can more accurately reflect the actual depth condition of each of a plurality of horizons in the stratum of the research area. And because the horizon data are accurate, the time domain structure diagram determined according to the horizon data is also accurate, so that the velocity data of each of the M wells determined according to the depth domain structure diagram and the time domain structure diagram can accurately reflect the propagation velocity of the seismic waves in the well, and finally, the three-dimensional velocity volume of the research area determined according to the velocity data of each of the M wells can accurately reflect the actual propagation velocity of the seismic waves at each stratum position in the research area.

Description

Method, apparatus and computer storage medium for determining three-dimensional velocity volume

Technical Field

The application relates to the technical field of oil and gas field development, in particular to a method and a device for determining a three-dimensional velocity body and a computer storage medium.

Background

In the process of oil and gas field exploration or development of a research area, a three-dimensional velocity body of the research area is generally required to be determined, so that characteristics such as depth at any stratum position in the research area or lithology characteristics at the stratum position are determined according to the three-dimensional velocity body of the research area. Wherein the three-dimensional velocity volume of the study region is used to indicate the velocity of seismic waves propagating at any of the formation locations of the study region.

In a related technique, a seismic synthetic record is acquired for each of a plurality of wells located in a study area while determining a three-dimensional velocity volume for the study area. The stratum corresponding to any well comprises a plurality of layers, the seismic synthetic record comprises a plurality of propagation velocities, the propagation velocities are in one-to-one correspondence with the layers, and each propagation velocity is used for indicating the propagation velocity of seismic waves at the junction of the corresponding layer and the adjacent layer. And for any two adjacent wells in the plurality of wells, determining the propagation velocity of the seismic waves in each horizon in the stratum between the two wells in a numerical interpolation mode according to the propagation velocities in the two seismic synthetic records corresponding to the two wells respectively. For example, if well A and well B are two adjacent wells, and the propagation velocity of the seismic wave is a at horizon 1 of well A and B at horizon 1 of well B, then the propagation velocity of the seismic wave at horizon 1 in the formation between well A and well B can be determined according to a and B. After the plurality of wells are processed according to the method, the propagation speed of seismic waves at any stratum position of the research area can be obtained, and a three-dimensional speed body of the research area is determined.

The above process determines a three-dimensional velocity volume from the seismic synthetic recordings for each of a plurality of wells in the area of interest, resulting in the determined three-dimensional velocity volume not being able to accurately represent the actual propagation velocity of seismic waves at various formation locations in the area of interest.

Content of application

The embodiment of the application provides a method and a device for determining a three-dimensional velocity body and a computer storage medium, wherein the determined three-dimensional velocity body of a research area can more accurately represent the actual propagation velocity of seismic waves at each stratum position in the research area. The technical scheme is as follows:

in a first aspect, a method of determining a three-dimensional velocity volume is provided, the method comprising:

determining a depth domain structure map of a research area according to horizon data of the research area, a seismic synthetic record of each well of N wells of the research area, and hierarchical data of each well of M wells of the research area, wherein the stratum of the research area comprises a plurality of horizons, the horizon data is used for indicating the time used by seismic waves to propagate from a surface position to any two adjacent horizon junctions of the research area, the seismic synthetic record comprises a plurality of propagation velocities, the propagation velocities correspond to the horizons in a one-to-one manner, each propagation velocity is used for indicating the propagation velocity of the seismic waves at the junctions of the corresponding horizon and the adjacent horizons, the hierarchical data is used for indicating the depth of any two adjacent horizon junctions of the stratum corresponding to each well, and N and M are positive integers greater than or equal to 1, and said N is less than said M;

determining a time domain structure map of the research region according to the horizon data of the research region;

determining velocity data for each of the M wells from the depth domain configuration map of the area of interest and the time domain configuration map of the area of interest, the velocity data being indicative of the propagation velocity of the seismic waves at any formation location in the respective well;

determining a three-dimensional velocity volume for the region of interest from the velocity data for each of the M wells.

Optionally, the determining a depth domain tectonic map of the investigation region from the horizon data of the investigation region, the seismic synthetic logs for each of the N wells of the investigation region, and the stratigraphic data for each of the M wells of the investigation region comprises:

determining an initial depth domain formation map for the area of interest from the horizon data and the seismic synthetic logs for each of the N wells;

and correcting the initial depth domain structure map according to the hierarchical data of each well in the M wells to obtain the depth domain structure map of the research area.

Optionally, after the correcting the initial depth domain configuration map of the research region according to the stratification data of each of the M wells, the method further includes:

and if the corrected initial depth domain structure map has interface discontinuity of two adjacent layers, performing data compensation on depth values of discontinuous interfaces in the corrected initial depth domain structure map, and determining the depth domain structure map after data compensation as the depth domain structure map of the research area.

and if the interface discontinuity of two adjacent horizons does not exist in the corrected initial depth domain structure map, determining the corrected initial depth domain structure map as the depth domain structure map of the research area.

Optionally, the determining a time domain configuration map of the research region according to the horizon data of the research region includes:

determining an initial time domain configuration map of the region of interest from the horizon data;

and if the interface discontinuity of two adjacent layers exists in the initial time domain structure diagram, performing data compensation on the time value of the discontinuous interface in the initial time domain structure diagram, and determining the time domain structure diagram after the data compensation as the time domain structure diagram of the research area.

Optionally, after determining the initial time domain configuration map of the research region according to the horizon data, the method further includes:

and if the interface discontinuity of the two adjacent horizons does not exist in the initial time domain structure map, determining the initial time domain structure map as the time domain structure map of the research area.

Optionally, the determining velocity data for each of the M wells from the depth-domain formation map of the area of interest and the time-domain formation map of the area of interest comprises:

for a first well of the M wells, determining the propagation velocity of the seismic wave of the first well at the junction of every two adjacent horizons in the plurality of horizons of the research area according to the depth domain structure diagram of the research area and the time domain structure diagram of the research area, wherein the first well is any one well of the M wells;

and determining the velocity data of the first well according to the propagation velocity of the seismic wave of the first well at the junction of every two adjacent horizons in the plurality of horizons of the research area, the depth domain structure diagram of the research area and the time domain structure diagram of the research area.

Optionally, said determining a three-dimensional velocity volume for the region of interest from the velocity data for each of the M wells comprises:

dividing the depth domain configuration map of the study region or the time domain configuration map of the study region into a plurality of meshes;

determining a velocity of each grid of a first type in the plurality of grids from the velocity data for each of the M wells, the velocity of each grid being indicative of a propagation velocity of seismic waves at a location of a geological formation in the corresponding grid, the grid of the first type being a grid of the plurality of grids in the M wells;

and determining the speed of each grid in a second type of grid according to the speed of each grid in the first type of grid, wherein the second type of grid is a grid except the first type of grid in the plurality of grids.

In a second aspect, there is provided an apparatus for determining a three-dimensional velocity volume, the apparatus comprising:

a first determining module, configured to determine a depth domain structure map of a research area according to horizon data of the research area, a seismic synthetic record of each of N wells of the research area, and hierarchical data of each of M wells of the research area, where a stratum of the research area includes a plurality of horizons, the horizon data is used to indicate time used for seismic waves to propagate from a surface position to a boundary between any two adjacent horizons of the research area, the seismic synthetic record includes a plurality of propagation velocities, the plurality of propagation velocities are in one-to-one correspondence with the plurality of horizons, each propagation velocity is used to indicate a propagation velocity of the seismic waves at a boundary between a corresponding horizon and an adjacent horizon, the hierarchical data is used to indicate a depth of any two adjacent layer boundaries in the stratum corresponding to each well, and N and M are positive integers greater than or equal to 1, and said N is less than said M;

a second determination module, configured to determine a time domain configuration map of the research region according to the horizon data of the research region;

a third determination module for determining velocity data for each of the M wells from the depth domain configuration map of the area of interest and the time domain configuration map of the area of interest, the velocity data being indicative of the propagation velocity of the seismic waves at any formation location in the respective well;

a fourth determination module to determine a three-dimensional velocity volume for the region of interest based on the velocity data for each of the M wells.

Optionally, the first determining module includes:

a first determination unit for determining an initial depth domain configuration map for the area of interest from the horizon data and the seismic synthetic logs for each of the N wells;

and the correcting unit is used for correcting the initial depth domain structure map according to the hierarchical data of each well in the M wells to obtain the depth domain structure map of the research area.

Optionally, the first determining module further includes:

and the first compensation unit is used for performing data compensation on depth values of discontinuous interfaces in the corrected initial depth domain structural map if the discontinuous interfaces of two adjacent layers exist in the corrected initial depth domain structural map, and determining the depth domain structural map after the data compensation as the depth domain structural map of the research area.

Optionally, the first determining module further includes:

a second determining unit, configured to determine the corrected initial depth domain structure map as the depth domain structure map of the research region if there is no interface discontinuity between two adjacent horizons in the corrected initial depth domain structure map.

Optionally, the second determining module includes:

a third determination unit for determining an initial time domain configuration map of the investigation region from the horizon data;

and a third compensation unit, configured to, if there is interface discontinuity between two adjacent horizons in the initial time domain configuration diagram, perform data compensation on a time value of a discontinuous interface in the initial time domain configuration diagram, and determine the time domain configuration diagram after data compensation as the time domain configuration diagram of the study region.

Optionally, the second determining module further includes:

a fourth determining unit, configured to determine the initial time domain configuration map as the time domain configuration map of the research region if there is no interface discontinuity between two adjacent horizons in the initial time domain configuration map.

Optionally, the third determining module includes:

a fifth determining unit, configured to determine, for a first well of the M wells, a propagation velocity of a seismic wave of the first well at a boundary of every two adjacent horizons of a plurality of horizons of the study area according to the depth domain structure diagram of the study area and the time domain structure diagram of the study area, where the first well is any one of the M wells;

a sixth determining unit, configured to determine velocity data of the first well according to a propagation velocity of the seismic wave of the first well at a boundary of every two adjacent horizons in the plurality of horizons of the study area, the depth domain structural map of the study area, and the time domain structural map of the study area.

Optionally, the fourth determining module includes:

a dividing unit, configured to divide the depth domain structure map of the research region or the time domain structure map of the research region into a plurality of grids;

a seventh determining unit, configured to determine, according to the velocity data of each of the M wells, a velocity of each of a first type of meshes in the multiple meshes, where the velocity of each mesh is used to indicate a propagation velocity of a seismic wave at a location layer position of the corresponding mesh, and the first type of meshes are meshes in the M wells in the multiple meshes;

an eighth determining unit, configured to determine a speed of each mesh in a second type of mesh according to the speed of each mesh in the first type of mesh, where the second type of mesh is a mesh other than the first type of mesh in the multiple meshes.

In a third aspect, there is provided an apparatus for determining a three-dimensional velocity volume, the apparatus comprising:

a processor;

a memory for storing processor-executable instructions;

wherein the processor is configured to perform the steps of any of the methods of the first aspect described above.

In a fourth aspect, there is provided a computer readable storage medium having stored thereon instructions which, when executed by a processor, carry out the steps of any of the methods of the first aspect described above.

In a fifth aspect, there is provided a computer program product comprising instructions which, when run on a computer, cause the computer to perform the steps of any of the methods of the first aspect described above.

The technical scheme provided by the embodiment of the application has the following beneficial effects:

in the application, when the depth domain structure map of the research area is determined, the determined depth domain structure map of the research area can more accurately reflect the actual depth situation of each horizon in a plurality of horizons in a stratum of the research area according to the seismic synthetic record of each well in N wells of the research area, the horizon data of the research area and the hierarchical data of each well in M wells of the research area, wherein M is larger than N. And because the horizon data are accurate, the time domain structure diagram determined according to the horizon data is also accurate, so that the velocity data of each of the M wells determined according to the depth domain structure diagram and the time domain structure diagram can accurately reflect the propagation velocity of the seismic waves in the well, and finally, the three-dimensional velocity volume of the research area determined according to the velocity data of each of the M wells can accurately reflect the actual propagation velocity of the seismic waves at each stratum position in the research area.

Drawings

In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present application, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.

FIG. 1 is a flow chart of a method for determining a three-dimensional velocity volume according to an embodiment of the present application;

FIG. 2 is a schematic diagram of a layer A having a fault phenomenon according to an embodiment of the present disclosure;

FIG. 3 is a schematic diagram of horizon data for horizon A in a formation of a study area provided by an embodiment of the present application;

FIG. 4 is a schematic illustration of a first well including a plurality of horizons provided by an embodiment of the present application;

FIG. 5 is a schematic illustration of the propagation velocities of seismic waves included in a grid provided by an embodiment of the present application;

FIG. 6 is a schematic plan view of a grid provided by an embodiment of the present application;

FIG. 7 is a schematic structural diagram of an apparatus for determining a three-dimensional velocity volume according to an embodiment of the present application;

fig. 8 is a schematic structural diagram of a terminal according to an embodiment of the present application.

Detailed Description

To make the objects, technical solutions and advantages of the present application more clear, embodiments of the present application will be described in further detail below with reference to the accompanying drawings.

Fig. 1 is a flowchart of a method for determining a three-dimensional velocity volume according to an embodiment of the present application, and as shown in fig. 1, the method includes the following steps:

step 101: a depth domain tectonic map for the study region is determined from the horizon data for the study region, the seismic synthetic logs for each of the N wells of the study region, and the stratigraphic data for each of the M wells of the study region.

The stratum of the research area comprises a plurality of layers, layer data is used for indicating the time of seismic waves transmitted to any two adjacent layer junctions of the research area from an earth surface position, the seismic synthetic record comprises a plurality of propagation speeds, the propagation speeds are in one-to-one correspondence with the layers, each propagation speed is used for indicating the propagation speed of the seismic waves at the corresponding layer junction and the adjacent layer junction, layered data is used for indicating the depth of any two adjacent layer junctions of the stratum corresponding to each well, N and M are positive integers greater than or equal to 1, and N is smaller than M.

In a possible implementation manner, step 101 may specifically be: an initial depth domain configuration map for the area of interest is determined from the horizon data and the seismic synthetic logs for each of the N wells. And correcting the initial depth domain structure diagram according to the hierarchical data of each well in the M wells to obtain the depth domain structure diagram of the research area.

The seismic synthetic record comprises a plurality of propagation speeds, wherein the propagation speeds correspond to the plurality of layers one to one, and each propagation speed is used for indicating the propagation speed of seismic waves at the junction of the corresponding layer and the adjacent layer. Thus, from the horizon data and the seismic synthetic logs for each of the N wells, an implementation in determining an initial depth domain configuration map for the region of interest may be: the stratum corresponding to the first well in the N wells comprises a plurality of layers, and the depth of the junction of the first layer and the second layer in the stratum corresponding to the first well can be obtained by multiplying the time for the seismic waves to propagate to the junction of the first layer and the second layer of the research area by the propagation speed of the seismic waves at the junction of the first layer and the second layer in the stratum corresponding to the first well. And then the depth of the junction of the first layer and the second layer corresponding to the first well is used as the depth of the junction of the first layer and the second layer in the stratum close to the first well. For other layers except the first layer and the second layer in the research area, the depths of the junctions of all adjacent layers in the other layers can be determined according to the method, the depths of the junctions of the first layer and the second layer and the depths of the junctions of all adjacent layers in the other layers can be determined, the depths of the junctions of all adjacent layers in the research area can be determined, and the depths of the junctions of all adjacent layers in the research area are determined, namely the initial depth domain structure diagram of the research area is determined. The first well is any one of N wells, the first layer is any one of a plurality of layers, and the second layer is a layer adjacent to the first layer in the plurality of layers.

For example, N is 1, and the stratum corresponding to this well includes five layers of a layer a, a layer B, a layer C, a layer D, and a layer E which are adjacent in sequence, the time used for acquiring seismic waves in the layer data of the research area to propagate to the boundary of the layer a and the layer B is 200 milliseconds, the time used for seismic waves to propagate to the boundary of the layer B and the layer C is 300 milliseconds, the time used for seismic waves to propagate to the boundary of the layer C and the layer D is 350 milliseconds, the time used for seismic waves to propagate to the boundary of the layer D and the layer E is 400 milliseconds, the propagation speed of seismic waves at the boundary of the layer a and the layer B is 15 m/millisecond, the propagation speed of seismic waves at the boundary of the layer B and the layer C is 17 m/millisecond, and the propagation speed of seismic waves at the boundary of the layer C and the layer D is 18 m/millisecond, the seismic wave has a propagation velocity at the interface of horizon D and horizon E of 20 m/ms. Multiplying 200 milliseconds and 15 m/millisecond, 300 milliseconds and 17 m/millisecond, 350 milliseconds and 18 m/millisecond, 400 milliseconds and 20 m/millisecond respectively to obtain 3000 m, 5100 m, 6300 m and 8000 m which are the depth of the well at the junction of the horizon A and the horizon B, the depth of the junction of the horizon B and the horizon C, the depth of the junction of the horizon C and the horizon D and the depth of the junction of the horizon D and the horizon E respectively. And then 3000 meters are taken as the depth of the junction of the level A and the level B of the research area, 5100 meters are taken as the depth of the junction of the level B and the level C of the research area, 6300 meters are taken as the depth of the junction of the level C and the level D of the research area, and 8000 meters are taken as the depth of the junction of the level D and the level E of the research area. The depth of every two adjacent horizons in the five horizons is determined, and the depth of every two adjacent horizons is taken as the depth of every two adjacent horizons in the five horizons in the research area, namely, an initial depth domain structure map of the research area is determined.

In addition, because the depth of each of the plurality of horizons of the stratigraphic layer of the study area may vary, the initial depth domain configuration map of the study area may not correctly reflect the actual depth of each of the horizons of the stratigraphic layer of the study area. Therefore, the initial depth domain structure map of the research region needs to be corrected, so that the corrected initial depth domain structure map can accurately reflect the actual depth of each horizon in the stratum of the research region. For the above reasons, the initial depth domain structure map of the research region needs to be corrected to obtain the depth domain structure map of the research region. In a possible implementation manner, the initial depth domain structure map may be corrected according to the hierarchical data of each of the M wells, so as to obtain a depth domain structure map of the research area.

The implementation manner of correcting the initial depth domain structure diagram according to the hierarchical data of each well of the M wells to obtain the depth domain structure diagram of the research area may specifically be: and finding the position of each well in the M wells in the initial depth domain structure diagram, and replacing the depth of the junction of the first layer and the second layer at the position of the well in the initial depth domain structure diagram with the depth of the junction of the first layer and the second layer at the position of the well in the formation corresponding to the well Z for the well Z in the M wells so as to realize the correction of the depth at the position. For other layers except the first layer and the second layer at the well Z position in the depth domain structure diagram, the depth of the junction of all adjacent layers in other layers can be determined according to the depth of the junction of the first layer and the second layer, and after the depth of the junction of all adjacent layers at the well Z position in the depth domain structure diagram is determined, the correction of the well Z position in the depth domain structure diagram is completed. The position of each of the wells other than well Z in the M wells in the depth domain map may be corrected in the manner described above. After correction for the location of each of the M wells in the depth-domain map, a depth-domain map of the area of interest is obtained. Wherein, well Z is any one well in M wells.

Because the initial depth domain structure map is corrected by using the hierarchical data of each well in the M wells, and the hierarchical data of each well in the M wells is accurate, the corrected depth domain structure map can accurately reflect the actual depth of each horizon in a plurality of horizons of the research area.

It should be noted that, because there are a plurality of horizons in the stratum of the research area, the depth of each horizon is different, so that the horizon data corresponding to each horizon is also different, and because each horizon of the research area is in a three-dimensional form in space, that is, a plurality of horizons of the research area are in a three-dimensional form in space. Thus, the initial depth domain map determined from the horizon data for the area of interest and the seismic synthetic logs for each of the N wells is also a three-dimensional stereogram, and the corrected initial depth domain map is also a three-dimensional stereogram.

In addition, there may be a case where any one of the plurality of horizons in the stratum of the study region is discontinuous, that is, a fault phenomenon occurs in the horizon, and if a fault phenomenon occurs in any one of the horizons in the stratum of the study region, there is no time taken for the seismic waves to propagate to the fault at the fault, that is, the horizon data at the fault is missing at the fault. For example, fig. 2 is a schematic diagram of a fault phenomenon occurring at a layer a, a fault dip is θ, and the fault dip is an included angle between a fault plane and a horizontal plane of the layer a. Generally, locations in the horizon where faults occur are also referred to as stratigraphic fault zones. FIG. 3 is a schematic diagram of horizon data of horizon A in a stratum of a study area, where the horizon data of horizon A is missing due to a fault.

If a fault phenomenon occurs in a horizon in a stratum of a study area, the horizon data of the study area is lost at a fault of a corresponding horizon in the study area, so that a depth value of the corresponding horizon in an initial depth domain structure map determined according to the horizon data of the corresponding horizon in the study area and a seismic synthetic record of each of N wells is also lost, and the depth value of the corresponding horizon in the corrected initial depth domain structure map is also lost, namely, an interface between the horizon and an adjacent horizon is discontinuous. Therefore, if the interface discontinuity of two adjacent layers exists in the corrected initial depth domain structure map, performing data compensation on the depth values of the discontinuous interface in the corrected initial depth domain structure map, and determining the depth domain structure map after the data compensation as the depth domain structure map of the study area.

The implementation manner of performing data compensation on the depth values of the discontinuous interfaces in the corrected initial depth domain configuration diagram may be: determining the starting position and the ending position of the discontinuous interface discontinuity, and performing data compensation between the depth of the starting position of the discontinuous interface and the depth of the ending position of the discontinuous interface. The method specifically comprises the following steps: determining the depth of the starting position of the discontinuity of the interface, determining the depth of the ending position of the discontinuity of the interface, and performing data compensation on the discontinuity of the interface according to a preset depth compensation rule by taking the depth of the starting position as a starting point until the last data in the compensated data is equal to the depth of the ending position.

For example, the depth of the start position is 3500 m, the depth of the end position is 3700 m, and when data compensation is performed, one data is compensated between 3500 m and 3700 m at an interval of 2 m, that is, the compensated data is 3502 and 3504 …. Data compensation ends until the compensated data is 3700.

And if the corrected initial depth domain structure map does not have the condition that the interface of the two adjacent layers is discontinuous, determining the corrected initial depth domain structure map as the depth domain structure map of the research area. That is, if the corrected initial depth is not continuous with the interface of two adjacent horizons of the structural map, the fault phenomenon does not exist in any one of the horizons of the stratum of the research area.

It should be noted that, the horizon data of the research area is obtained after seismic interpretation is performed on the seismic data of the research area, and when seismic interpretation is performed, there is a certain accuracy requirement in the seismic interpretation, for example, seismic interpretation is performed on the research area according to 20 meters by 20 meters, so the horizon data of the research area also has a certain accuracy. The accuracy of the horizon data of the research area refers to the size of each grid when the interface between the horizon indicated by the horizon data corresponding to each horizon in the horizon data and the adjacent horizon is divided into planar grids.

Step 102: and determining a time domain structure map of the research region according to the horizon data of the research region.

The implementation manner of determining the time domain structure diagram of the research region according to the horizon data of the research region may be: and determining an initial time domain structural diagram of the research area according to the horizon data, if the interface of two adjacent horizons is discontinuous in the initial time domain structural diagram, performing data compensation on the time value of the discontinuous interface in the initial time domain structural diagram, and determining the time domain structural diagram after the data compensation as the time domain structural diagram of the research area.

The horizon data is used for indicating the time used by seismic waves propagating from the surface position to the junction of any two adjacent horizons of the research area, so that the implementation mode of determining the initial time domain structural diagram of the research area according to the horizon data of the research area can be as follows: the time used by the junction of any two adjacent horizons of the horizon data is determined, the time is formed into a time set, the time set is determined into a time surface, and the time surface can represent the space form of the seismic wave in the research area which is propagated to the horizon at the corresponding horizon. After the time planes at the boundaries of all adjacent horizons in the horizon data have been determined, an initial time domain configuration map for the region of interest is determined.

In addition, if a fault occurs in a horizon in the stratum of the study region, the horizon data of the study region is lost in the study region at the fault of the corresponding horizon, and therefore, the initial time domain structure map determined based on the horizon data is also lost in the corresponding horizon. When the initial time domain structure map has the absence at the corresponding horizon, the interface discontinuity of two adjacent horizons in the initial time domain structure map is indicated. And if the interface discontinuity of two adjacent layers exists in the initial time domain structure diagram, performing data compensation on the time value of the discontinuous interface in the initial time domain structure diagram, and determining the time domain structure diagram after the data compensation as the time domain structure diagram of the research area.

The implementation manner of performing data compensation on the time values of the discontinuous interfaces in the initial time domain configuration diagram may be as follows: determining the starting position and the ending position of the discontinuous interface discontinuity, and performing data compensation between the time of the starting position of the discontinuous interface and the time of the ending position of the discontinuous interface. The method specifically comprises the following steps: determining the time of the starting position of the discontinuity of the interface, determining the time of the ending position of the discontinuity of the interface, and performing data compensation on the discontinuity of the interface according to a preset time compensation rule by taking the time of the starting position as a starting point until the time of the last data in the compensated data is equal to the time of the ending position.

For example, the time of the start position is 25 ms, the depth of the end position is 40 ms, and when data compensation is performed, one data is compensated between 25 ms and 40 ms at 1 ms intervals, that is, the compensated data is 26, 27 …. Data compensation is ended until the compensated data is 40.

In addition, if the initial time domain structure map does not have the condition that the interface of the two adjacent layers is discontinuous, the initial time domain structure map is determined as the time domain structure map of the research area. If the initial time is not continuous with the interface of two adjacent horizons in the structural map, the fault phenomenon does not exist in any one of the horizons of the research area.

It should be noted that, because there are a plurality of horizons in the stratum of the research area, and the depth of each horizon is different, the horizon data corresponding to each horizon is also different, and because each horizon of the research area is three-dimensional in space, that is, a plurality of horizons of the research area are three-dimensional in space, the initial time domain structural map determined according to the horizon data of the research area is also a three-dimensional stereo map.

Step 103: from the depth domain configuration map of the area of interest and the time domain configuration map of the area of interest, velocity data is determined for each of the M wells, the velocity data being indicative of the propagation velocity of seismic waves at any of the formation locations in the respective well.

In a possible implementation manner, determining the velocity data of each of the M wells according to the depth domain structure diagram of the study area and the time domain structure diagram of the study area may specifically be: for a first well of the M wells, determining a propagation velocity of the seismic wave of the first well at a boundary of every two adjacent horizons in the plurality of horizons of the study area according to the depth domain structure map of the study area and the time domain structure map of the study area. And determining velocity data of the first well according to the propagation velocity of the seismic wave of the first well at the junction of every two adjacent horizons in the plurality of horizons of the research area, the depth domain structure diagram of the research area and the time domain structure diagram of the research area. Wherein the first well is any one of the M wells.

For a first well of the M wells, according to the depth domain structure diagram of the study area and the time domain structure diagram of the study area, an implementation manner of determining the propagation velocity of the seismic wave of the first well at the boundary of every two adjacent horizons in the multiple horizons of the study area may be: the method comprises the steps of finding the depth of a junction of a first layer and a second layer of a first well in a plurality of layers of a research area in a depth domain structure diagram, finding the time of the junction of the first well in the first layer and the second layer of the plurality of layers of the research area in a time domain structure diagram, and dividing the depth of the junction of the first well in the first layer and the second layer of the plurality of layers of the research area by the time of the junction of the first well in the first layer and the second layer of the plurality of layers of the research area to obtain the propagation velocity of seismic waves of the first well at the junction of the first layer and the second layer of the plurality of layers of the research area. The first horizon is any one of a plurality of horizons. The second horizon is a horizon of the plurality of horizons that is adjacent to the first horizon.

In addition, according to the propagation velocity of the seismic wave of the first well at the boundary of every two adjacent horizons in the plurality of horizons of the research area, the depth domain structural diagram of the research area and the time domain structural diagram of the research area, the implementation manner of determining the velocity data of the first well may be: determining a first intersection point of a first well at a junction of a first layer and a second layer in a plurality of layers of a research area, determining a second intersection point of the first well at a junction of the second layer and a third layer in the plurality of layers of the research area, determining a third intersection point of the first well at a junction of a third layer and a fourth layer in the plurality of layers of the research area, and determining a midpoint of the first intersection point and the second intersection point and a midpoint of the second intersection point and the third intersection point. And determining the propagation velocity of the seismic wave at the midpoint between the first intersection point and the second intersection point according to the depth of the first intersection point, the time of the seismic wave propagating to the first intersection point, the depth of the second intersection point and the time of the seismic wave propagating to the second intersection point. And determining the propagation velocity of the seismic wave at the midpoint between the second intersection point and the third intersection point according to the depth of the second intersection point, the time of the seismic wave propagating to the second intersection point, the depth of the third intersection point and the time of the seismic wave propagating to the third intersection point. And determining the coefficient of the propagation acceleration of the seismic waves when the seismic waves propagate in the first well according to the velocity of the seismic waves at the midpoint between the first intersection point and the second intersection point, the velocity of the seismic waves at the midpoint between the second intersection point and the third intersection point, the depth of the seismic waves at the midpoint between the first intersection point and the second intersection point and the depth of the seismic waves at the midpoint between the second intersection point and the third intersection point. And determining the propagation velocity of the seismic wave at any position between the first intersection point and the midpoint between the first intersection point and the second intersection point according to the propagation acceleration coefficient when the seismic wave propagates in the first well and the propagation velocity of the seismic wave at the first intersection point. And determining the propagation velocity of the seismic wave at any position between the midpoint between the first intersection point and the second intersection point according to the propagation acceleration coefficient when the seismic wave propagates in the first well and the propagation velocity of the seismic wave at the midpoint between the first intersection point and the second intersection point. The propagation velocity of the seismic wave at any position between the second intersection point and the third intersection point can be determined in the above manner. The propagation velocity of the seismic wave at any position in each layer of the other layers except the first layer, the second layer and the third layer in the first well can be determined in the above mode, and after the propagation velocity of the seismic wave at any position in all the layers in the first well is determined, the velocity data of the first well is determined.

The determination of the propagation velocity of the seismic wave at the midpoint between the first intersection point and the second intersection point according to the depth of the first intersection point, the time of the seismic wave propagating to the first intersection point, the depth of the second intersection point and the time of the seismic wave propagating to the second intersection point can be realized by the following formula:

V₃＝(Z₂-Z₁)/(T₂-T₁)

in the above formula Z₁Is the depth of the first intersection point, Z₂Is the depth of the second intersection, T₁Time of propagation of seismic waves to the first intersection, T₂Is the time, V, of the seismic wave propagation to the first intersection₃Is the velocity of propagation of the seismic wave at the midpoint between the first intersection point and the second intersection point.

In addition, the coefficient for determining the propagation acceleration of the seismic wave when propagating in the first well according to the velocity of the seismic wave at the midpoint between the first intersection point and the second intersection point, the velocity of the seismic wave at the midpoint between the second intersection point and the third intersection point, the depth of the seismic wave at the midpoint between the first intersection point and the second intersection point, and the depth of the seismic wave at the midpoint between the second intersection point and the third intersection point can be realized by the following formula:

K＝(V₄-V₃)/(Z₄-Z₃)

in the above formula V₃Is the velocity, V, of the seismic wave propagating to the midpoint between the first and second intersection points₄Is the velocity, Z, of the seismic wave propagating to a midpoint between the second and third intersection points₃Is the depth at the midpoint between the first and second intersection points, Z₄And K is the propagation acceleration coefficient of seismic waves when the seismic waves propagate in the first well.

In addition, determining the propagation velocity of the seismic wave at any position between the first intersection point and the midpoint between the first intersection point and the second intersection point according to the propagation acceleration coefficient when the seismic wave propagates in the first well and the propagation velocity of the seismic wave at the first intersection point can be realized by the following formula:

V＝V₁+K×(Z-Z₁)

in the above formula V₁Is the propagation velocity of the seismic wave at the first intersection point, Z is the depth of any position between the first intersection point and the midpoint between the first intersection point and the second intersection point, Z₁And K is the coefficient of the propagation acceleration of the seismic wave when the seismic wave propagates in the first well, and V is the propagation velocity of the seismic wave at any position between the first intersection point and the midpoint between the first intersection point and the second intersection point.

For example, as shown in fig. 4, the intersection point of the first well and the junction of level 1 and level 2 is a, i.e., the first intersection point is a, the intersection point of the first well and the junction of level 2 and level 3 is B, i.e., the second intersection point is B, the intersection point of the first well and the junction of level 3 and level 4 is C, i.e., the third intersection point is C, the midpoint of the first intersection point a and the second intersection point B is D, and the midpoint of the second intersection point B and the third intersection point C is E. The propagation velocity of the seismic wave at point D and the propagation velocity of the seismic wave at point E can be determined according to the following equations:

V_D＝(Z_B-Z_A)/(T_B-T_A)

V_E＝(Z_C-Z_B)/(T_C-T_B)

in the above formula Z_A、Z_BAnd Z_CThe depth of the first intersection point A, the depth of the second intersection point B and the depth of the third intersection point C, T_A、T_B、T_CRespectively the time of seismic wave propagation to the first intersection point AThe time of the seismic wave to propagate to the second intersection point B and the time of the seismic wave to propagate to the third intersection point C, V_DAnd V_EThe propagation velocity of the seismic wave at point D and the propagation velocity of the seismic wave at point E, respectively.

The propagation acceleration coefficient of the seismic wave as it propagates in the first well may be determined according to the following equation:

K＝(V_E-V_D)/(Z_E-Z_D)

in the above formula Z_EAnd Z_DDepth of point E and depth of point D, V, respectively_DAnd V_EThe propagation velocity of the seismic wave at the point D and the propagation velocity of the seismic wave at the point E are respectively, and K is the propagation acceleration coefficient when the seismic wave propagates in the first well.

The propagation velocity of the seismic wave at any position between the first intersection point a and the point D can be determined according to the following equation:

V＝V_A+K×(Z-Z_A)

in the above formula V_AThe propagation velocity of the seismic wave at a first intersection point A, K is the propagation acceleration coefficient of the seismic wave when propagating in the first well, Z is the depth of any position of the first intersection point A and a point D, and Z is the propagation velocity of the seismic wave at the first intersection point A_AThe depth of the first intersection point A, and V is the propagation speed of the seismic waves at any position of the first intersection point A and the point D.

Step 104: from the velocity data for each of the M wells, a three-dimensional velocity volume for the area of interest is determined.

In one possible implementation, determining the three-dimensional velocity volume of the study area according to the velocity data of each of the M wells may specifically be: the method comprises the steps of dividing a depth domain structure diagram of a research area or a time domain structure diagram of the research area into a plurality of grids, determining the speed of each grid in a first type of grid in the plurality of grids according to the speed data of each well in M wells, and determining the speed of each grid in a second type of grid according to the speed of each grid in the first type of grid. The velocity of each grid is used for indicating the propagation velocity of seismic waves at the position of a stratum where the corresponding grid is located, the first type of grid is a grid in M wells in the multiple grids, and the second type of grid is a grid except the first type of grid in the multiple grids.

The depth domain structure diagram of the research area or the time domain structure diagram of the research area is divided into a plurality of grids, and the implementation manner may be as follows: and dividing the depth domain structural diagram or the time domain structural diagram into a plurality of grids according to a preset rule. The preset rule may be to divide a planar grid in a plane parallel to the ground surface and to divide a vertical grid in a plane perpendicular to the ground surface, that is, each grid of the multiple grids after the depth domain structural diagram or the time domain structural diagram is divided according to the preset rule is a three-dimensional cubic grid.

It should be noted that, when dividing the plane parallel to the ground surface into planar meshes, the division accuracy may be consistent with the accuracy of the horizon data.

Additionally, determining the velocity of each of the first type of mesh in the plurality of meshes from the velocity data for each of the M wells may be accomplished by: determining the number of the propagation speeds of the seismic waves contained in each grid in the first type of grid in the first well according to the speed data of the first well in the M wells, determining the average value of the propagation speeds of the seismic waves contained in each grid according to the number of the propagation speeds of the seismic waves contained in each grid in the first grid, and taking the average value of the propagation speeds of the seismic waves contained in each grid as the speed of the grid.

For example, as shown in fig. 5, the grid includes three seismic wave propagation velocities V1, V2, and V3, and the velocity in the grid is (V1+ V2+ V3)/3.

Optionally, according to the velocity data of each of the M wells, the implementation of determining the velocity of each of the first type of mesh in the plurality of meshes may further be: and determining the maximum value of the propagation velocity of the plurality of seismic waves in the first well, wherein the plurality of seismic waves are contained in each grid in the first type of grid, and taking the maximum value of the propagation velocity as the velocity of each grid in the first type of grid.

In addition, according to the speed of each mesh in the first type of mesh, the implementation manner of determining the speed of each mesh in the second type of mesh may be: and determining the speed of each grid in the second type of grids surrounding the first grid in the first type of grids according to the weight of the speed of each grid and the speed of the first grid. The first grid is any grid in the first type of grid. Each mesh of the second type of mesh surrounding the first mesh refers to each mesh of the second type of mesh that is at the same level as the first mesh.

Wherein the weighting for determining the velocity of each mesh in the second type of mesh around the first mesh in the first type of mesh may be implemented according to the following formula:

in the above formula_iAs a weight of the velocity of each mesh in the second type of mesh surrounding the first mesh, d_i0P is the distance between each mesh in the second type of mesh around the first mesh and the first mesh, and is an index value.

In addition, determining the velocity of each mesh in the second type of mesh around the first mesh based on the weight of the velocity of each mesh and the velocity of the first mesh may be implemented according to the following formula:

in the above formula V (S)₀) Is the velocity of the first grid, λ_iWeight for the velocity of each mesh in the second type of mesh surrounding the first mesh, V (S)_i) The velocity of each mesh in the second type of mesh around the first mesh.

For example, as shown in fig. 6, there are 15 second type meshes at the same level as the first mesh, of which 3 are located at the left side of the first mesh and 12 are located at the right side of the first mesh.

Fig. 7 is a schematic structural diagram of an apparatus for determining a three-dimensional velocity body according to an embodiment of the present application, and as shown in fig. 7, the apparatus 700 includes:

a first determining module 701, configured to determine a depth domain structure diagram of a research region according to horizon data of the research region, a seismic synthetic record of each of N wells of the research region, and hierarchical data of each of M wells of the research region, where a stratum of the research region includes a plurality of horizons, where the horizon data is used to indicate time used for seismic waves to propagate from an earth surface position to a junction of any two adjacent horizons of the research region, the seismic synthetic record includes a plurality of propagation velocities, the propagation velocities are in one-to-one correspondence with the horizons, each propagation velocity is used to indicate a propagation velocity of seismic waves at a junction of a corresponding horizon and an adjacent horizon, the hierarchical data is used to indicate a depth of a junction of any two adjacent horizons in the stratum corresponding to each well, N and M are positive integers greater than or equal to 1, and N is less than M;

a second determining module 702, configured to determine a time domain configuration map of the research region according to the horizon data of the research region;

a third determining module 703, configured to determine, according to the depth domain structural diagram of the research area and the time domain structural diagram of the research area, velocity data of each of the M wells, where the velocity data is used to indicate a propagation velocity of the seismic wave at any formation position in the corresponding well;

a fourth determination module 704 for determining a three-dimensional velocity volume for the area of interest based on the velocity data for each of the M wells.

Optionally, the first determining module 701 includes:

a first determination unit for determining an initial depth domain configuration map for the area of interest based on the horizon data and the seismic synthetic logs for each of the N wells;

and the correcting unit is used for correcting the initial depth domain structure diagram according to the hierarchical data of each well in the M wells to obtain the depth domain structure diagram of the research area.

Optionally, the first determining module 701 further includes:

and the first compensation unit is used for performing data compensation on depth values of discontinuous interfaces in the corrected initial depth domain structural map if the discontinuous interfaces of two adjacent layers exist in the corrected initial depth domain structural map, and determining the depth domain structural map after the data compensation as the depth domain structural map of the study area.

Optionally, the first determining module 701 further includes:

and a second determining unit, configured to determine the corrected initial depth domain structure map as the depth domain structure map of the research region if the corrected initial depth domain structure map does not have interface discontinuity of two adjacent horizons.

Optionally, the second determining module 702 includes:

and the third compensation unit is used for performing data compensation on time values of discontinuous interfaces in the initial time domain structure diagram if the discontinuous interfaces of two adjacent layers exist in the initial time domain structure diagram, and determining the time domain structure diagram after the data compensation as the time domain structure diagram of the research area.

Optionally, the second determining module 702 further includes:

and a fourth determining unit, configured to determine the initial time domain structure map as the time domain structure map of the study region if the interface discontinuity of the two adjacent horizons does not exist in the initial time domain structure map.

Optionally, the third determining module 703 includes:

a fifth determining unit, configured to determine, for a first well of the M wells, a propagation velocity of the seismic wave of the first well at a boundary between every two adjacent horizons in the multiple horizons of the study area according to the depth domain structural diagram of the study area and the time domain structural diagram of the study area, where the first well is any one of the M wells;

and the sixth determining unit is used for determining the velocity data of the first well according to the propagation velocity of the seismic wave of the first well at the junction of every two adjacent horizons in the plurality of horizons of the research area, the depth domain structure diagram of the research area and the time domain structure diagram of the research area.

Optionally, the fourth determining module 704 includes:

a dividing unit for dividing the depth domain structure map of the research area or the time domain structure map of the research area into a plurality of grids;

a seventh determining unit, configured to determine, according to the velocity data of each of the M wells, a velocity of each of a first type of mesh in the multiple meshes, where the velocity of each mesh is used to indicate a propagation velocity of the seismic wave at a location layer position of the corresponding mesh, and the first type of mesh is a mesh in the M wells in the multiple meshes;

and an eighth determining unit, configured to determine a speed of each mesh in the second type of mesh according to the speed of each mesh in the first type of mesh, where the second type of mesh is a mesh other than the first type of mesh in the plurality of meshes.

It should be noted that: the apparatus for determining a three-dimensional velocity volume provided in the above embodiment is only illustrated by the division of the above functional modules when determining the three-dimensional velocity volume, and in practical applications, the above function distribution may be completed by different functional modules according to needs, that is, the internal structure of the device is divided into different functional modules, so as to complete all or part of the above described functions. In addition, the apparatus for determining a three-dimensional velocity body and the method for determining a three-dimensional velocity body provided in the above embodiments belong to the same concept, and specific implementation processes thereof are described in the method embodiments, and are not described herein again.

Fig. 8 shows a block diagram of a terminal 800 according to an exemplary embodiment of the present application. The terminal 800 may be: a smart phone, a tablet computer, an MP3 player (Moving Picture Experts Group Audio Layer III, motion video Experts compression standard Audio Layer 3), an MP4 player (Moving Picture Experts Group Audio Layer iv, motion video Experts compression standard Audio Layer 4), a notebook computer, or a desktop computer. The terminal 800 may also be referred to by other names such as user equipment, portable terminal, laptop terminal, desktop terminal, etc.

In general, the terminal 800 includes: a processor 801 and a memory 802.

The processor 801 may include one or more processing cores, such as a 4-core processor, an 8-core processor, and so forth. The processor 801 may be implemented in at least one hardware form of a DSP (Digital Signal Processing), an FPGA (Field-Programmable Gate Array), and a PLA (Programmable Logic Array). The processor 801 may also include a main processor and a coprocessor, where the main processor is a processor for processing data in an awake state, and is also called a Central Processing Unit (CPU); a coprocessor is a low power processor for processing data in a standby state. In some embodiments, the processor 801 may be integrated with a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content required to be displayed on the display screen. In some embodiments, the processor 801 may further include an AI (Artificial Intelligence) processor for processing computing operations related to machine learning.

Memory 802 may include one or more computer-readable storage media, which may be non-transitory. Memory 802 may also include high speed random access memory, as well as non-volatile memory, such as one or more magnetic disk storage devices, flash memory storage devices. In some embodiments, a non-transitory computer readable storage medium in memory 802 is used to store at least one instruction for execution by processor 801 to implement the method of determining a three-dimensional velocity volume provided by the method embodiments herein.

In some embodiments, the terminal 800 may further include: a peripheral interface 803 and at least one peripheral. The processor 801, memory 802 and peripheral interface 803 may be connected by bus or signal lines. Various peripheral devices may be connected to peripheral interface 803 by a bus, signal line, or circuit board. Specifically, the peripheral device includes: at least one of a radio frequency circuit 804, a touch screen display 805, a camera assembly 806, an audio circuit 807, a positioning assembly 808, and a power supply 809.

The peripheral interface 803 may be used to connect at least one peripheral related to I/O (Input/Output) to the processor 801 and the memory 802. In some embodiments, the processor 801, memory 802, and peripheral interface 803 are integrated on the same chip or circuit board; in some other embodiments, any one or two of the processor 801, the memory 802, and the peripheral interface 803 may be implemented on separate chips or circuit boards, which are not limited by this embodiment.

The Radio Frequency circuit 804 is used for receiving and transmitting RF (Radio Frequency) signals, also called electromagnetic signals. The radio frequency circuitry 804 communicates with communication networks and other communication devices via electromagnetic signals. The rf circuit 804 converts an electrical signal into an electromagnetic signal to be transmitted, or converts a received electromagnetic signal into an electrical signal. Optionally, the radio frequency circuit 804 includes: an antenna system, an RF transceiver, one or more amplifiers, a tuner, an oscillator, a digital signal processor, a codec chipset, a subscriber identity module card, and so forth. The radio frequency circuit 804 may communicate with other terminals via at least one wireless communication protocol. The wireless communication protocols include, but are not limited to: metropolitan area networks, various generation mobile communication networks (2G, 3G, 4G, and 5G), Wireless local area networks, and/or WiFi (Wireless Fidelity) networks. In some embodiments, the radio frequency circuit 804 may further include NFC (Near Field Communication) related circuits, which are not limited in this application.

The display screen 805 is used to display a UI (User Interface). The UI may include graphics, text, icons, video, and any combination thereof. When the display 805 is a touch display, the display 805 also has the ability to capture touch signals on or above the surface of the display 805. The touch signal may be input to the processor 801 as a control signal for processing. At this point, the display 805 may also be used to provide virtual buttons and/or a virtual keyboard, also referred to as soft buttons and/or a soft keyboard. In some embodiments, the display 805 may be one, providing the front panel of the terminal 800; in other embodiments, the display 805 may be at least two, respectively disposed on different surfaces of the terminal 800 or in a folded design; in still other embodiments, the display 805 may be a flexible display disposed on a curved surface or a folded surface of the terminal 800. Even further, the display 805 may be arranged in a non-rectangular irregular pattern, i.e., a shaped screen. The Display 805 can be made of LCD (liquid crystal Display), OLED (Organic Light-Emitting Diode), and the like.

The camera assembly 806 is used to capture images or video. Optionally, camera assembly 806 includes a front camera and a rear camera. Generally, a front camera is disposed at a front panel of the terminal, and a rear camera is disposed at a rear surface of the terminal. In some embodiments, the number of the rear cameras is at least two, and each rear camera is any one of a main camera, a depth-of-field camera, a wide-angle camera and a telephoto camera, so that the main camera and the depth-of-field camera are fused to realize a background blurring function, and the main camera and the wide-angle camera are fused to realize panoramic shooting and VR (Virtual Reality) shooting functions or other fusion shooting functions. In some embodiments, camera assembly 806 may also include a flash. The flash lamp can be a monochrome temperature flash lamp or a bicolor temperature flash lamp. The double-color-temperature flash lamp is a combination of a warm-light flash lamp and a cold-light flash lamp, and can be used for light compensation at different color temperatures.

The audio circuit 807 may include a microphone and a speaker. The microphone is used for collecting sound waves of a user and the environment, converting the sound waves into electric signals, and inputting the electric signals to the processor 801 for processing or inputting the electric signals to the radio frequency circuit 804 to realize voice communication. For the purpose of stereo sound collection or noise reduction, a plurality of microphones may be provided at different portions of the terminal 800. The microphone may also be an array microphone or an omni-directional pick-up microphone. The speaker is used to convert electrical signals from the processor 801 or the radio frequency circuit 804 into sound waves. The loudspeaker can be a traditional film loudspeaker or a piezoelectric ceramic loudspeaker. When the speaker is a piezoelectric ceramic speaker, the speaker can be used for purposes such as converting an electric signal into a sound wave audible to a human being, or converting an electric signal into a sound wave inaudible to a human being to measure a distance. In some embodiments, the audio circuitry 807 may also include a headphone jack.

The positioning component 808 is used to locate the current geographic position of the terminal 800 for navigation or LBS (location based Service). The positioning component 808 may be a positioning component based on the GPS (global positioning System) in the united states, the beidou System in china, the graves System in russia, or the galileo System in the european union.

Power supply 809 is used to provide power to various components in terminal 800. The power supply 809 can be ac, dc, disposable or rechargeable. When the power source 809 comprises a rechargeable battery, the rechargeable battery may support wired or wireless charging. The rechargeable battery may also be used to support fast charge technology.

In some embodiments, terminal 800 also includes one or more sensors 810. The one or more sensors 810 include, but are not limited to: acceleration sensor 811, gyro sensor 812, pressure sensor 813, fingerprint sensor 814, optical sensor 815 and proximity sensor 816.

The acceleration sensor 811 may detect the magnitude of acceleration in three coordinate axes of the coordinate system established with the terminal 800. For example, the acceleration sensor 811 may be used to detect the components of the gravitational acceleration in three coordinate axes. The processor 801 may control the touch screen 805 to display the user interface in a landscape view or a portrait view according to the gravitational acceleration signal collected by the acceleration sensor 811. The acceleration sensor 811 may also be used for acquisition of motion data of a game or a user.

The gyro sensor 812 may detect a body direction and a rotation angle of the terminal 800, and the gyro sensor 812 may cooperate with the acceleration sensor 811 to acquire a 3D motion of the user with respect to the terminal 800. From the data collected by the gyro sensor 812, the processor 801 may implement the following functions: motion sensing (such as changing the UI according to a user's tilting operation), image stabilization at the time of photographing, game control, and inertial navigation.

Pressure sensors 813 may be disposed on the side bezel of terminal 800 and/or underneath touch display 805. When the pressure sensor 813 is disposed on the side frame of the terminal 800, the holding signal of the user to the terminal 800 can be detected, and the processor 801 performs left-right hand recognition or shortcut operation according to the holding signal collected by the pressure sensor 813. When the pressure sensor 813 is disposed at a lower layer of the touch display screen 805, the processor 801 controls the operability control on the UI interface according to the pressure operation of the user on the touch display screen 805. The operability control comprises at least one of a button control, a scroll bar control, an icon control and a menu control.

The fingerprint sensor 814 is used for collecting a fingerprint of the user, and the processor 801 identifies the identity of the user according to the fingerprint collected by the fingerprint sensor 814, or the fingerprint sensor 814 identifies the identity of the user according to the collected fingerprint. Upon identifying that the user's identity is a trusted identity, the processor 801 authorizes the user to perform relevant sensitive operations including unlocking a screen, viewing encrypted information, downloading software, paying for and changing settings, etc. Fingerprint sensor 814 may be disposed on the front, back, or side of terminal 800. When a physical button or a vendor Logo is provided on the terminal 800, the fingerprint sensor 814 may be integrated with the physical button or the vendor Logo.

The optical sensor 815 is used to collect the ambient light intensity. In one embodiment, the processor 801 may control the display brightness of the touch screen 805 based on the ambient light intensity collected by the optical sensor 815. Specifically, when the ambient light intensity is high, the display brightness of the touch display screen 805 is increased; when the ambient light intensity is low, the display brightness of the touch display 805 is turned down. In another embodiment, the processor 801 may also dynamically adjust the shooting parameters of the camera assembly 806 based on the ambient light intensity collected by the optical sensor 815.

A proximity sensor 816, also known as a distance sensor, is typically provided on the front panel of the terminal 800. The proximity sensor 816 is used to collect the distance between the user and the front surface of the terminal 800. In one embodiment, when the proximity sensor 816 detects that the distance between the user and the front surface of the terminal 800 gradually decreases, the processor 801 controls the touch display 805 to switch from the bright screen state to the dark screen state; when the proximity sensor 816 detects that the distance between the user and the front surface of the terminal 800 becomes gradually larger, the processor 801 controls the touch display 805 to switch from the screen-on state to the screen-on state.

Those skilled in the art will appreciate that the configuration shown in fig. 8 is not intended to be limiting of terminal 800 and may include more or fewer components than those shown, or some components may be combined, or a different arrangement of components may be used.

Embodiments of the present application further provide a non-transitory computer-readable storage medium, and when instructions in the storage medium are executed by a processor of a terminal, the terminal is enabled to execute the method for determining a three-dimensional velocity volume provided in the embodiment shown in fig. 1.

Embodiments of the present application also provide a computer program product containing instructions, which when run on a computer, cause the computer to perform the method for determining a three-dimensional velocity volume provided in the embodiment shown in fig. 1.

It will be understood by those skilled in the art that all or part of the steps for implementing the above embodiments may be implemented by hardware, or may be implemented by a program instructing relevant hardware, where the program may be stored in a computer-readable storage medium, and the above-mentioned storage medium may be a read-only memory, a magnetic disk or an optical disk, etc.

In summary, the present application is only a preferred embodiment and is not intended to be limited by the scope of the present application, and any modifications, equivalents, improvements, etc. made within the spirit and principle of the present application should be included in the protection scope of the present application.

Claims

1. A method of determining a three-dimensional velocity volume, the method comprising:

determining an initial depth domain structure map of a research area according to horizon data of the research area and seismic synthetic records of each well of N wells of the research area, correcting the initial depth domain structure map according to hierarchical data of each well of M wells of the research area to obtain the depth domain structure map of the research area, wherein the strata of the research area comprise a plurality of horizons, the horizon data are used for indicating the time of seismic waves propagating from a surface position to the junction of any two adjacent horizons of the research area, the seismic synthetic records comprise a plurality of propagation velocities, the plurality of propagation velocities are in one-to-one correspondence with the plurality of horizons, each propagation velocity is used for indicating the propagation velocity of the seismic waves at the junction of the corresponding horizon and adjacent horizon, and the hierarchical data are used for indicating the depth of any two adjacent horizons in the strata corresponding to each well, the N and the M are positive integers which are more than or equal to 1, and the N is less than the M;

2. The method of claim 1, wherein after correcting the initial depth domain configuration map for the area of interest based on the stratigraphic data for each of the M wells, further comprising:

3. The method of claim 1, wherein after correcting the initial depth domain configuration map for the area of interest based on the stratigraphic data for each of the M wells, further comprising:

4. The method of claim 1, wherein determining a time domain configuration map for the region of interest from the horizon data for the region of interest comprises:

5. The method of claim 4, wherein after determining an initial time domain configuration map for the region of interest from the horizon data, further comprising:

6. The method of any one of claims 1 to 5, wherein determining the velocity data for each of the M wells from the depth domain map of the area of interest and the time domain map of the area of interest comprises:

7. The method of any one of claims 1 to 5, wherein said determining a three-dimensional velocity volume for the area of interest from the velocity data for each of the M wells comprises:

8. An apparatus for determining a three-dimensional velocity volume, the apparatus comprising:

a fourth determination module for determining a three-dimensional velocity volume for the region of interest based on the velocity data for each of the M wells;

the first determining module includes:

9. The apparatus of claim 8, wherein the first determining module further comprises:

10. The apparatus of claim 8, wherein the first determining module further comprises:

11. The apparatus of claim 8, wherein the second determining module comprises:

12. The apparatus of claim 11, wherein the second determining module further comprises:

13. The apparatus of any of claims 8 to 12, wherein the third determining module comprises:

14. The apparatus of any of claims 8 to 12, wherein the fourth determining module comprises:

15. An apparatus for determining a three-dimensional velocity volume, the apparatus comprising:

a processor;

a memory for storing processor-executable instructions;

wherein the processor is configured to perform the steps of any one of the method of claim 1 to claim 7.

16. A computer readable storage medium having stored thereon instructions which, when executed by a processor, carry out the steps of the method of any one of claims 1 to 7.