Method and device for establishing initial model of continuous medium chromatography, and electronic equipment

Abstract

Embodiments of the present application provide a method and device for establishing a continuous medium tomography initial model and an electronic device. The method for establishing a continuous medium tomography initial model comprises: obtaining first arrival set control points in a plurality of preset regional units in a target work area; performing a series of transformations on offset-time information of each first arrival set control point to obtain a velocity-depth fitting function corresponding to the first arrival set control point; obtaining surface investigation control point parameters of the first arrival set control points according to the velocity-depth fitting function and the size of the grid unit of the continuous medium tomography initial model; and obtaining the continuous medium tomography initial model of the target work area by using a similarity coefficient modeling method according to the surface investigation control point parameters corresponding to the first arrival set control points in the plurality of preset regional units.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of geophysical exploration, and in particular to a method and device for establishing a continuous medium tomographic initial model and an electronic device. BACKGROUND

[0002] Static correction is a technique for eliminating the delay of seismic waves by low-velocity layers and velocity-drop layers near the surface, and correcting the reflected waves to the environment equivalent to no low-velocity band and velocity-drop band and to be excited and received in a horizontal plane. To achieve static correction, a geophysical model of low-velocity band and velocity-drop band needs to be established. The current common modeling method is to use the first arrival time on a single shot to achieve inversion.

[0003] In the related art, a geophysical model can be established by tomographic inversion modeling using first arrival information. Tomographic inversion modeling needs to be completed based on a ray tracing and first arrival time iterative process, and therefore has a greater dependence on an initial velocity model, so a more reasonable initial velocity model is necessary. A very effective method for establishing an initial velocity model has been developed for a layered structure near the surface, but for a continuous medium near-surface structure, there is a lack of a method for establishing a high-precision continuous medium tomographic initial model.

[0004] Therefore, how to establish a high-precision continuous medium tomographic initial model is a technical problem to be solved. SUMMARY

[0005] Embodiments of the present application provide a method and device for establishing a continuous medium tomographic initial model, an electronic device and a storage medium, for establishing a high-precision continuous medium tomographic initial model.

[0006] One of the embodiments of the present application provides a method for establishing a continuous medium tomography initial model, the method comprising: obtaining a set of first arrivals control points in a plurality of preset area units in a target work area; the set of first arrivals control points is a set of first arrival information of shot points in the preset area units, and the first arrival information at least includes a first arrival time and a moveout distance corresponding to the first arrival time, and each of the preset area units corresponds to a grid unit in the continuous medium tomography initial model; for each of the set of first arrivals control points, the following operations are performed: obtaining a time-moveout distance fitting function corresponding to the set of first arrivals control points according to the first arrival time and the moveout distance of each first arrival information in the set of first arrivals control points; performing a derivation operation on the time-moveout distance fitting function to obtain a velocity-moveout distance fitting function of the set of first arrivals control points; obtaining a velocity-depth fitting function corresponding to the set of first arrivals control points according to the velocity-moveout distance fitting function and a time-moveout distance function of refracted waves; obtaining a surface investigation control point parameter of the set of first arrivals control points according to the size of the grid unit of the continuous medium tomography initial model according to the velocity-depth fitting function; and obtaining the continuous medium tomography initial model of the target work area by using a similarity coefficient modeling method according to the surface investigation control point parameters corresponding to the set of first arrivals control points in the plurality of preset area units.

[0007] In some embodiments, the time-moveout distance fitting function of the set of first arrivals control points is obtained according to the first arrival time and the moveout distance of each first arrival information in the set of first arrivals control points, comprising: obtaining a plurality of coordinate points in a plane rectangular coordinate system according to the first arrival time and the moveout distance of each first arrival information; fitting the plurality of coordinate points to obtain a first fitting curve; determining a polynomial expression corresponding to the first fitting curve to obtain the time-moveout distance fitting function.

[0008] In some embodiments, the velocity-depth fitting function of the set of first arrivals control points is obtained according to the velocity-moveout distance fitting function and the time-moveout distance function of refracted waves, comprising: substituting the first arrival time and the moveout distance of each first arrival information in the set of first arrivals control points into the velocity-moveout distance fitting function and the time-moveout distance function of refracted waves to obtain a velocity-depth discrete data set; determining a coordinate point corresponding to each data in the velocity-depth discrete data set in a plane rectangular coordinate system; fitting a plurality of determined coordinate points to obtain a second fitting curve; determining a polynomial expression corresponding to the second fitting curve to obtain the velocity-depth fitting function of the set of first arrivals control points.

[0009] In some embodiments, the obtaining the surface investigation control point parameter of the set of first arrivals control points according to the size of the grid cell of the continuous medium tomography initial model and according to the velocity-depth fitting function comprises: obtaining a depth parameter of a surface investigation control point corresponding to the grid cell according to the size of the grid cell; and obtaining a velocity parameter of the surface investigation control point according to the size of the grid cell and by using the velocity-depth fitting function.

[0010] In some embodiments, the obtaining the depth parameter of the surface investigation control point corresponding to the grid cell according to the size of the grid cell comprises: taking a vertical distance between a position corresponding to a center point of the grid cell and the ground surface as the depth parameter of the surface investigation control point.

[0011] In some embodiments, the velocity parameter comprises an average velocity and a maximum velocity; and the obtaining the velocity parameter of the surface investigation control point according to the size of the grid cell and by using the velocity-depth fitting function comprises: dividing the grid cell into a plurality of sub-grid cells along a direction perpendicular to the ground surface, taking a distance between a center point of each of the sub-grid cells and the ground surface as a depth value; substituting a plurality of the depth values into the velocity-depth fitting function to obtain a plurality of velocities; averaging the plurality of velocities to obtain the average velocity in the velocity parameter; and taking a maximum value in the plurality of velocities as the maximum velocity in the velocity parameter.

[0012] One of the embodiments of the present application provides a device for establishing a continuous medium tomography initial model, the device comprising: a first obtaining module configured to obtain a set of first arrivals control points in a plurality of preset area units in a target work area; the set of first arrivals control points is a set of first arrival information of shot points in the preset area units, and the first arrival information at least includes a first arrival time and a moveout distance corresponding to the first arrival time; each of the preset area units corresponds to a grid unit in the continuous medium tomography initial model; a second obtaining module configured to, for each of the set of first arrivals control points, perform the following operations: obtaining a time-moveout distance fitting function corresponding to the set of first arrivals control points according to the first arrival time and the moveout distance of each of the first arrival information in the set of first arrivals control points; obtaining a velocity-moveout distance fitting function of the set of first arrivals control points by performing a derivation operation on the time-moveout distance fitting function; obtaining a velocity-depth fitting function corresponding to the set of first arrivals control points according to the velocity-moveout distance fitting function and a time-moveout distance function of refracted waves; and obtaining a surface investigation control point parameter of the set of first arrivals control points according to the velocity-depth fitting function and the size of the grid unit of the continuous medium tomography initial model; and a third obtaining module configured to obtain the continuous medium tomography initial model of the target work area by using a similarity coefficient modeling method according to the surface investigation control point parameters corresponding to the set of first arrivals control points in the plurality of preset area units.

[0013] The embodiments of the present application provide an electronic device, which comprises a memory and a processor, the memory stores a computer program, and the processor executes the program to perform the method described above.

[0014] The embodiments of the present application provide a storage medium for storing a computer readable program, which is executed to perform the method described above.

[0015] The above technical solutions provided by the embodiments of the present application have at least the following advantages compared with the prior art:

[0016] In the embodiments provided by the present application, the moveout distance-time information of the set of first arrivals control points in the plurality of preset area units in the target work area is transformed in a series of ways to obtain a velocity-depth fitting function corresponding to the set of first arrivals control points; the surface investigation control point parameter of the set of first arrivals control points is obtained according to the velocity-depth fitting function and the size of the grid unit of the continuous medium tomography initial model; and the continuous medium tomography initial model of the target work area is obtained by using a similarity coefficient modeling method according to the surface investigation control point parameters corresponding to the set of first arrivals control points in the plurality of preset area units. Therefore, a high-precision continuous medium tomography initial model can be established. BRIEF DESCRIPTION OF DRAWINGS

[0017] The application will be further described in the way of example embodiments, which will be described in detail with reference to the accompanying drawings. These embodiments are not restrictive, and in these embodiments, the same numbers represent the same structures, wherein:

[0018] Figure 1 is an example flow chart of a method for establishing a continuous medium tomography initial model according to some embodiments of the application;

[0019] Figure 2 is an example schematic diagram of an apparatus for establishing a continuous medium tomography initial model according to some embodiments of the application;

[0020] Figure 3 is an example structural schematic diagram of an electronic device according to some embodiments of the application. DETAILED DESCRIPTION

[0021] In order to more clearly illustrate the technical solutions of the embodiments of the application, the following will briefly introduce the drawings needed to be used in the embodiment description. Obviously, the drawings in the following description are only some examples or embodiments of the application, and for those skilled in the art, the application can also be applied to other similar scenarios without creative labor on the basis of these drawings. Unless it is clear from the language context or otherwise indicated, the same reference numbers in the drawings represent the same structures or operations.

[0022] It should be understood that the "system", "apparatus", "grid unit" and / or "module" used herein is a method for distinguishing different components, elements, parts, portions or assemblies at different levels. However, if other words can achieve the same purpose, the words can be replaced by other expressions.

[0023] As shown in the application and claims, unless the context clearly indicates otherwise, "one", "a", "an" and / or "the" do not refer to the singular, but also include the plural. Generally speaking, the terms "comprise" and "include" only indicate the inclusion of the steps and elements explicitly identified, and these steps and elements do not constitute an exclusive list, and the method or device can also include other steps or elements.

[0024] Flowcharts are used in the application to illustrate the operations performed by the system according to the embodiments of the application. It should be understood that the preceding or subsequent operations are not necessarily executed in sequence. On the contrary, each step can be processed in reverse order or simultaneously. At the same time, other operations can be added to these processes, or one or more steps of operation can be removed from these processes.

[0025] In order to facilitate understanding, the technical solutions of the application will be introduced below in combination with the drawings and embodiments.

[0026] Figure 1 is an exemplary flowchart of a method for establishing a continuous medium tomography initial model according to some embodiments of the present application. As shown in Figure 1 the method for establishing a continuous medium tomography initial model includes the following steps:

[0027] In step S110, a set of first arrival control points in a plurality of preset area units in a target work area are obtained.

[0028] Each preset area unit corresponds to a grid unit in the continuous medium tomography initial model. In a specific implementation process, the target work area can be divided into a plurality of preset area units according to the physical area corresponding to the grid unit in the continuous medium tomography initial model.

[0029] The set of first arrivals control points is a set of first arrival information of shot points in the preset area unit, and each preset area unit corresponds to a set of first arrival control points. The first arrival information at least includes a first arrival time and an offset distance corresponding to the first arrival time.

[0030] In seismic exploration, one-dimensional exploration is to observe the underground situation of a point; two-dimensional exploration is to observe the underground situation under a line; and three-dimensional exploration is to observe the underground situation under an area. Therefore, for three-dimensional exploration data, a set of first arrival control points is formed according to a plurality of shot points in the grid of the exploration area; and for two-dimensional data, a set of first arrival control points is formed according to a plurality of shot points in the exploration line segment. In a specific implementation process, the coincident shot points in the adjacent exploration area grid or the adjacent exploration line segment are selected, which can improve the rationality of the set of first arrival control points.

[0031] In a specific implementation process, the formation process of the set of first arrival control points can be realized by using professional software.

[0032] After obtaining the set of first arrival control points, the following steps can be performed for each set of first arrival control points:

[0033] In step S120, a time-offset distance fitting function corresponding to the set of first arrival control points is obtained according to the first arrival time and the offset distance of each first arrival information in the set of first arrival control points.

[0034] In some embodiments, the time-offset distance fitting function of the set of first arrival control points can be obtained by the following method: a plurality of coordinate points in a plane rectangular coordinate system are obtained according to the first arrival time and the offset distance of each first arrival information; the plurality of coordinate points are fitted to obtain a first fitting curve; a polynomial expression corresponding to the first fitting curve is determined to obtain the time-offset distance fitting function.

[0035] In the implementation process, the first arrival time and offset distance of each first arrival information can be placed in a plane rectangular coordinate system to obtain a plurality of discrete coordinate points. For example, the first arrival time can be taken as the horizontal coordinate and the offset distance can be taken as the vertical coordinate. For another example, the first arrival time can be taken as the vertical coordinate and the offset distance can be taken as the horizontal coordinate. The plurality of discrete points in the plane rectangular coordinate system are connected by a curve to obtain a first fitting curve. A polynomial fitting expression of the first fitting curve is obtained to obtain a time-offset distance function of the first arrival set control point.

[0036] In the implementation process, the drawing of the first fitting curve and the process of obtaining the polynomial fitting expression of the first fitting curve can be realized by using professional software.

[0037] In step S130, a derivative operation is performed on the time-offset distance fitting function to obtain a velocity-offset distance fitting function of the first arrival set control point.

[0038] In step S140, a velocity-depth fitting function corresponding to the first arrival set control point is obtained according to the velocity-offset distance fitting function and the time-offset distance function of the refracted wave.

[0039] In some embodiments, the velocity-depth fitting function of the first arrival set control point can be obtained in the following manner:

[0040] The first arrival time and offset distance of each first arrival information in the first arrival set control point are substituted into the velocity-offset distance fitting function obtained in step S130 and the time-offset distance function of the refracted wave to obtain a velocity-depth discrete data set. The coordinate points corresponding to each data in the velocity-depth discrete data set in the plane rectangular coordinate system are determined. The determined plurality of coordinate points are fitted to obtain a second fitting curve. A polynomial expression corresponding to the second fitting curve is determined to obtain the velocity-depth fitting function of the first arrival set control point.

[0041] The velocity-depth discrete data set includes a plurality of groups of discrete velocity-depth data pairs. In the implementation process, the plurality of groups of discrete velocity-depth data pairs can be placed in a plane rectangular coordinate system to form a plurality of coordinate points. The plurality of coordinate points in the plane rectangular coordinate system are connected by a curve to obtain a second fitting curve. A polynomial fitting expression of the second fitting curve is obtained to obtain the velocity-depth fitting function of the first arrival set control point.

[0042] In the implementation process, the drawing of the second fitting curve and the process of obtaining the polynomial fitting expression of the second fitting curve can be realized by using professional software.

[0043] In step S150, according to the size of the grid unit of the initial model of continuous medium tomography, the surface investigation control point parameters of the first arrival set control point are obtained according to the velocity-depth fitting function.

[0044] In some embodiments, the depth parameter of the near-surface survey control point corresponding to the grid cell can be obtained according to the size of the grid cell. For example only, the vertical distance between the position of the center point of the grid cell and the ground surface can be taken as the depth parameter of the near-surface survey control point.

[0045] In some embodiments, the velocity parameter of the near-surface survey control point can be obtained by using a velocity-depth fitting function according to the size of the grid cell.

[0046] For example only, the velocity parameter includes the average velocity and the maximum velocity, which can be obtained by the following method:

[0047] The grid cell is divided into a plurality of sub-grid cells along the direction perpendicular to the ground surface, and the distance between the center point of each sub-grid cell and the ground surface is taken as a depth value; the plurality of depth values are substituted into the velocity-depth fitting function to obtain a plurality of velocities; the average of the plurality of velocities is taken as the average velocity in the velocity parameter; and the maximum value in the plurality of velocities is taken as the maximum velocity in the velocity parameter.

[0048] In step S160, the initial model of the continuous medium tomography of the target work area is obtained by using the similarity coefficient modeling method according to the near-surface survey control point parameters corresponding to each of the first arrival set control points in the plurality of preset regional units.

[0049] In the specific implementation process, the initial model of the continuous medium tomography of the target work area can be established by using the similarity coefficient modeling method based on the near-surface survey control point parameters corresponding to each of the first arrival set control points in the plurality of preset regional units.

[0050] In the implementation provided in the present application, the first arrival set control points in the plurality of preset regional units in the target work area are obtained; a series of processing operations are performed on the offset-time information of each first arrival set control point to obtain the near-surface survey control point parameters of the first arrival set control points; and the initial model of the continuous medium tomography of the target work area is obtained according to the near-surface survey control point parameters corresponding to each of the first arrival set control points in the plurality of preset regional units. Therefore, the initial model of the continuous medium tomography obtained according to the embodiments of the present application can reflect the longitudinal and lateral changes of the continuous near-surface structure of the work area, thereby improving the modeling accuracy of the low-velocity zone in the longitudinal and lateral wave acquisition under the condition of the continuous medium structure.

[0051] The technical effects of the embodiments of the present application are described below through specific examples.

[0052] The eastern Qaidam Basin is a natural gas reservoir area, with a surface consisting of saline-alkali land filled by repeated evaporation and sandification from salt lakes. While shear wave exploration can provide a relatively complete description of the subsurface strata in gas-bearing areas, the near-surface structure in the saline-alkali region of eastern Qaidam exhibits significantly different shear wave and P-wave responses. It is not a layered structure like P-waves, but rather a continuous medium. Therefore, static correction results for near-surface structure models obtained using traditional initial modeling methods are unsatisfactory.

[0053] The method provided in this application was used for near-surface modeling in both the Sanhu 2D shear wave exploration area and the Taidong 3D shear wave exploration area in eastern Qaidam Basin. By establishing an initial continuous medium tomography model that conforms to the longitudinal and lateral variations in the Sanhu and Taidong areas, the tomography inversion results were effectively constrained, and high-precision near-surface structure models with shear wave continuity were established for these two areas, solving the serious static correction problem commonly found in shear wave exploration in the eastern Qaidam gas region.

[0054] Figure 2 This is an exemplary schematic diagram of an apparatus for establishing an initial model for continuous media chromatography according to some embodiments of this application.

[0055] like Figure 2 As shown, the apparatus for establishing the initial model of continuous medium chromatography includes: a first acquisition module 210, a second acquisition module 220, and a third acquisition module 230.

[0056] The first acquisition module 210 is used to acquire the initial arrival set control points in multiple preset area units within the target work area; the initial arrival set control points are a set of initial arrival information of shot points in the preset area units, and the initial arrival information includes at least the initial arrival time and the offset distance corresponding to the initial arrival time, and each preset area unit corresponds to a grid unit in the initial model of the continuous medium tomography.

[0057] The second acquisition module 220 is used to perform the following operations for each of the initial arrival set control points:

[0058] Based on the arrival time and offset of each arrival information in the initial arrival set control point, the time-offset fitting function corresponding to the initial arrival set control point is obtained; the derivative of the time-offset fitting function is performed to obtain the velocity-offset fitting function of the initial arrival set control point; based on the velocity-offset fitting function and the time-offset function of the refracted wave, the velocity-depth fitting function corresponding to the initial arrival set control point is obtained; based on the velocity-depth fitting function and according to the size of the grid cell of the initial model of the continuous medium tomography, the surface survey control point parameters of the initial arrival set control point are obtained.

[0059] The third obtaining module 230 is configured to obtain, according to the surface investigation control point parameters corresponding to each of the initial set control points in the plurality of preset area units, a continuous medium tomography initial model of the target work area by using a similarity coefficient modeling method.

[0060] In some embodiments, the step of obtaining the time-offset fitting function of each of the initial set control points according to the initial time and offset distance of the initial information of the initial set control point comprises: obtaining a plurality of coordinate points in a plane rectangular coordinate system according to the initial time and offset distance of the initial information of the initial set control point; fitting the plurality of coordinate points to obtain a first fitting curve; and determining a polynomial expression corresponding to the first fitting curve to obtain the time-offset fitting function.

[0061] In some embodiments, the step of obtaining the velocity-depth fitting function of the initial set control point according to the velocity-offset fitting function and the time-offset function of the refracted wave comprises: substituting the initial time and offset distance of the initial information of each of the initial set control points into the velocity-offset fitting function and the time-offset function of the refracted wave to obtain a velocity-depth discrete data set; determining a coordinate point corresponding to each data in the velocity-depth discrete data set in a plane rectangular coordinate system; fitting the plurality of determined coordinate points to obtain a second fitting curve; and determining a polynomial expression corresponding to the second fitting curve to obtain the velocity-depth fitting function of the initial set control point.

[0062] In some embodiments, the step of obtaining the surface investigation control point parameter of the initial set control point according to the size of the grid unit of the continuous medium tomography initial model comprises: obtaining a depth parameter of a surface investigation control point corresponding to the grid unit according to the size of the grid unit; and obtaining a velocity parameter of the surface investigation control point by using the velocity-depth fitting function according to the size of the grid unit.

[0063] In some embodiments, the step of obtaining the depth parameter of the surface investigation control point corresponding to the grid unit according to the size of the grid unit comprises: taking a vertical distance between a position corresponding to a center point of the grid unit and the ground surface as the depth parameter of the surface investigation control point.

[0064] In some embodiments, the velocity parameter comprises an average velocity and a maximum velocity; and the obtaining the velocity parameter of the surface survey control point according to the size of the grid cell and using the velocity-depth fitting function comprises: dividing the grid cell into a plurality of sub-grid cells along a direction perpendicular to the surface, taking a distance from the center point of each of the sub-grid cells to the surface as a depth value; substituting the plurality of depth values into the velocity-depth fitting function to obtain a plurality of velocities; and averaging the plurality of velocities to obtain the average velocity in the velocity parameter; and taking a maximum value in the plurality of velocities as the maximum velocity in the velocity parameter.

[0065] In the embodiment of the device for establishing the initial model of the continuous medium tomography, the specific processing of each module and the technical effects brought by the processing can be referred to the related description in the corresponding method embodiment, which will not be repeated here.

[0066] Figure 3 is an exemplary structural schematic diagram of an electronic device according to some embodiments of the present application.

[0067] As shown in Figure 3 The electronic device comprises at least one processor 301, at least one communication interface 302, at least one memory 303 and at least one communication bus 304. Optionally, the communication interface 302 can be an interface of a communication module, such as an interface of a GSM module. The processor 301 can be a processor CPU, or a specific integrated circuit ASIC (Application Specific Integrated Circuit), or one or more integrated circuits configured to implement embodiments of the present application. The memory 303 can include a high-speed RAM memory, and can also include a non-volatile memory, such as at least one disk memory. The memory 303 stores a program, and the processor 301 invokes the program stored in the memory 303 to execute part or all of the method embodiments described above.

[0068] The present application relates to a storage medium for storing a computer readable program, which, when executed, performs part or all of the method embodiments described above.

[0069] Optionally, the storage medium can be a non-transitory computer readable storage medium, for example, the non-transitory computer readable storage medium can be a ROM, a random access memory (RAM), a CD-ROM, a magnetic tape, a floppy disk and an optical data storage device, etc.

[0070] Based on the same inventive concept, the embodiment of the present application further provides a computer program product comprising a computer program, which, when executed by a processor, implements some or all of the method embodiments described above.

[0071] The foregoing merely illustrates the principles of the application. It will thus be appreciated that those skilled in the art will be able to devise various arrangements that, although not explicitly described or shown herein, embody the principles of the application and are thus within its spirit and scope. For example, although the above describes the system components as being implemented by hardware devices, the system components can also be implemented by software solutions, such as installing the described system on an existing server or mobile device.

[0072] In addition, the use of particular nomenclature in reference to certain features, aspects or characteristics of the application is not intended to limit the scope of the application, but rather is intended to highlight or emphasize the features, aspects or characteristics. As such, the foregoing description discloses only exemplary aspects of the application, and does not represent the only embodiments or implementations of the application. Accordingly, the scope of the application is indicated by the appended claims, rather than by the foregoing description, and all modifications that fall within the range of equivalents are intended to be embraced therein.

[0073] In addition, the order of processing elements or sequence of steps in the processes described is not limiting, except as can be expressly specified in the claims. Neither is the use of any of the terms in the description intended to limit the scope of the claims beyond the description accompanying the claims themselves. It will also be understood that any numerical range recited in the description is intended to include all derivatives within that range (e.g., endpoints are included and are divisible by both whole and fractional steps).

[0074] In addition, the use of particular nomenclature in reference to certain features, aspects or characteristics of the application is not intended to limit the scope of the application, but rather is intended to highlight or emphasize the features, aspects or characteristics. As such, the foregoing description discloses only exemplary aspects of the application, and does not represent the only embodiments or implementations of the application. Accordingly, the scope of the application is indicated by the appended claims, rather than by the foregoing description, and all modifications that fall within the range of equivalents are intended to be embraced therein.

[0075] In some embodiments, numbers that describe amounts, dimensions, and so forth, are used in the description of the embodiments. It should be understood that such numbers are used only to illustrate certain embodiments and that the application is not limited to the numbers. In some examples, such numbers are modified by the modifier "about" or "approximately." Unless otherwise specified, "about" or "approximately" means ±20% of the value of the measured quantity that the term describes. Accordingly, in some embodiments, the numerical parameters in the description and claims are approximations that can vary depending upon the desired properties sought to be obtained by the individual embodiments. In some embodiments, numerical parameters are determined by the use of standard techniques. Although the numerical ranges and parameters setting forth the broad scope of the application are approximations, the numerical values set forth in the specific examples are reported as precisely as practicable. The numerical values set forth in the specific examples are provided to be as precise as reasonably possible. However, some variations may occur depending on the choice of input used to develop or derive the numerical values in the examples.

[0076] Each patent, patent application, publication, and other material cited in this application is hereby incorporated by reference in its entirety. In the event of inconsistencies between the disclosure of this application and the materials incorporated by reference, the disclosure of this application is intended to prevail. Nothing in this application is intended to be forgo the common general knowledge of the art as of the priority date of this application. To the extent that any meaning or definition of a term in this application conflicts with the meaning or definition of the same term in a document incorporated by reference, the meaning or definition assigned to that term in this application shall control.

[0077] Finally, it should be understood that the embodiments described herein are merely exemplary of the application. Other variations of the embodiments can be devised by those skilled in the art without departing from the scope of the present application. Accordingly, the embodiments described herein are intended to be illustrative only and are not limiting of the scope of the present application.

Claims

1. A method for establishing a continuum chromatography initial model, characterized in that, The method comprises: acquiring first arrival set control points in a plurality of preset regional units in a target work area; the first arrival set control point is a set of first arrival information of a shot point in the preset regional unit, and the first arrival information at least comprises a first arrival time and a corresponding offset distance; each of the preset regional units corresponds to a grid unit in the continuous medium tomography initial model; for each of the first arrival set control points, the following operations are performed: obtaining a time-offset distance fitting function corresponding to the first arrival set control point according to the first arrival time and the offset distance of each first arrival information in the first arrival set control point; deriving the time-offset distance fitting function to obtain a velocity-offset distance fitting function of the first arrival set control point; obtaining a velocity-depth fitting function corresponding to the first arrival set control point according to the velocity-offset distance fitting function and a time-offset distance function of the refracted wave; obtaining a surface investigation control point parameter of the first arrival set control point according to the size of the grid unit of the continuous medium tomography initial model according to the velocity-depth fitting function; obtaining a continuous medium tomography initial model of the target work area by using a similarity coefficient modeling method according to the surface investigation control point parameters corresponding to the first arrival set control points in the plurality of preset regional units; obtaining a surface investigation control point parameter of the first arrival set control point according to the size of the grid unit of the continuous medium tomography initial model according to the velocity-depth fitting function, comprising: obtaining a depth parameter of the surface investigation control point corresponding to the grid unit according to the size of the grid unit; obtaining a velocity parameter of the surface investigation control point by using the velocity-depth fitting function according to the size of the grid unit; obtaining a depth parameter of the surface investigation control point corresponding to the grid unit according to the size of the grid unit, comprising: taking the vertical distance between the position corresponding to the center point of the grid unit and the ground surface as the depth parameter of the surface investigation control point; the velocity parameter comprises an average velocity and a maximum velocity; obtaining a velocity parameter of the surface investigation control point by using the velocity-depth fitting function according to the size of the grid unit, comprising: dividing the grid unit into a plurality of sub-grid units along the direction perpendicular to the ground surface, and taking the distance between the center point of each sub-grid unit and the ground surface as a depth value; substituting a plurality of the depth values into the velocity-depth fitting function to obtain a plurality of velocities; averaging the plurality of velocities to obtain the average velocity in the velocity parameter; taking the maximum value in the plurality of velocities as the maximum velocity in the velocity parameter.

2. The method of claim 1, wherein, obtaining a time-offset distance fitting function of the first arrival set control point according to the first arrival time and the offset distance of each first arrival information in the first arrival set control point, comprising: obtaining a plurality of coordinate points in a plane rectangular coordinate system according to the first arrival time and the offset distance of each first arrival information; fitting the plurality of coordinate points to obtain a first fitting curve; determine a polynomial expression corresponding to the first fitting curve to obtain the time-offset fitting function.

3. The method of claim 1, wherein, The velocity-depth fitting function of the initial arrival set control point is obtained according to the velocity-offset fitting function and a time-offset function of the refracted wave. The initial arrival time and offset of each initial arrival information in the initial arrival set control point are substituted into the velocity-offset fitting function and the time-offset function of the refracted wave to obtain a velocity-depth discrete data set. A coordinate point corresponding to each data in the velocity-depth discrete data set in a plane rectangular coordinate system is determined. A second fitting curve is obtained by fitting the determined multiple coordinate points. A polynomial expression corresponding to the second fitting curve is determined to obtain the velocity-depth fitting function of the initial arrival set control point.

4. An apparatus for establishing a continuous medium chromatography initial model, characterized by The device comprises: A first obtaining module is configured to obtain initial arrival set control points in a plurality of preset regional units in a target work area. The initial arrival set control point is a set of initial arrival information of a shot point in the preset regional unit. The initial arrival information at least includes an initial arrival time and an offset corresponding to the initial arrival time. Each preset regional unit corresponds to a grid unit in the continuous medium tomography initial model. A second obtaining module is configured to perform the following operations for each initial arrival set control point: A time-offset fitting function corresponding to the initial arrival set control point is obtained according to the initial arrival time and offset of each initial arrival information in the initial arrival set control point. A velocity-offset fitting function of the initial arrival set control point is obtained by performing a derivative operation on the time-offset fitting function. A velocity-depth fitting function corresponding to the initial arrival set control point is obtained according to the velocity-offset fitting function and a time-offset function of the refracted wave. A surface investigation control point parameter of the initial arrival set control point is obtained according to the size of the grid unit of the continuous medium tomography initial model based on the velocity-depth fitting function. A third obtaining module is configured to obtain a continuous medium tomography initial model of the target work area by using a similarity coefficient modeling method based on the surface investigation control point parameters corresponding to the initial arrival set control points in the plurality of preset regional units. The surface investigation control point parameter of the initial arrival set control point is obtained according to the size of the grid unit of the continuous medium tomography initial model based on the velocity-depth fitting function, which comprises: A depth parameter of a surface investigation control point corresponding to the grid unit is obtained according to the size of the grid unit. A velocity parameter of the surface investigation control point is obtained by using the velocity-depth fitting function according to the size of the grid unit. The depth parameter of the surface investigation control point corresponding to the grid unit is obtained according to the size of the grid unit, which comprises: A vertical distance between a position corresponding to a center point of the grid unit and the ground surface is taken as the depth parameter of the surface investigation control point. The velocity parameter comprises an average velocity and a maximum velocity. The velocity parameter of the surface investigation control point is obtained by using the velocity-depth fitting function according to the size of the grid unit, which comprises: dividing the grid cell into a plurality of sub-grid cells in a direction perpendicular to the earth surface, and taking a distance from the center point of each of the sub-grid cells to the earth surface as a depth value; substituting the plurality of depth values into the velocity-depth fitting function to obtain a plurality of velocities; averaging the plurality of velocities to obtain an average velocity in the velocity parameter; taking a maximum value in the plurality of velocities as a highest velocity in the velocity parameter. 5.An electronic device, comprising a memory and a processor, the memory storing a computer program, and the processor executing the program to perform the method of any one of claims 1 to 3. 6.A storage medium for storing a computer-readable program, the computer-readable program being executed to perform the method of any one of claims 1 to 3.

Images