Seismic wave travel time calculation method, device, electronic device and storage medium

Abstract

The present invention provides a method, apparatus, electronic device and storage medium for calculating seismic wave travel time. The method for calculating seismic wave travel time includes: obtaining an analytical travel time field of a homogeneous medium in a target work area, using a first-order solution of the factored eikonal equation to generate a first-order accurate smooth scaling factor according to the analytical travel time field; and using a high-order accurate fast marching solution of the factored eikonal equation to determine a high-order accurate smooth scaling factor and a travel time field of the target work area according to the first-order accurate smooth scaling factor. The method of the present invention can accelerate the calculation speed while ensuring high-order accuracy in calculating seismic wave travel time, and can be used to improve imaging accuracy and near-surface velocity modeling.

Description

Technical Field

[0001] The present invention relates to the field of oil exploration, especially to the field of seismic data processing technology, and particularly relates to a method, device, electronic device and storage medium for calculating seismic wave travel time. Background Art

[0002] Through exploration in theory and actual work, the importance of near-surface velocity has become a consensus: Although the near-surface velocity is only the shallow layer of the velocity model, its impact on imaging is huge. Modeling of the near-surface velocity is inseparable from the calculation of first arrival travel time. However, while improving the calculation accuracy of travel time, the calculation efficiency is often reduced.

[0003] Specifically, the existing technology for calculating seismic wave travel time has the following problems: The ray theory method has low accuracy and causes problems in the shadow area; The conventional finite difference solution of the eikonal equation has problems with source singularities and limited accuracy; When using high-order finite differences, the calculation efficiency will be reduced. Summary of the Invention

[0004] According to the method, device, electronic device and storage medium for calculating seismic wave travel time provided by the present invention, the eikonal equation can be directly solved to solve the problem of the shadow area; and by using the factored eikonal equation, the problem of source singularities is solved; while improving the accuracy, the calculation speed is accelerated.

[0005] To achieve the above object, a method for calculating seismic wave travel time is provided, including:

[0006] Obtain the analytical travel time field of the homogeneous medium in the target work area,

[0007] Using the first-order solution of the factored eikonal equation, generate a first-order accuracy smooth scaling factor according to the analytical travel time field;

[0008] Using the high-order accuracy fast sweeping solution of the factored eikonal equation, determine the high-order accuracy smooth scaling factor and travel time field of the target work area according to the first-order accuracy smooth scaling factor.

[0009] In one embodiment, the obtaining the analytical travel time field of the homogeneous medium in the target work area includes:

[0010] Obtain the velocity model of the target work area;

[0011] Generate the analytical travel time field corresponding to the velocity model according to the pre-generated homogeneous velocity field and the velocity model.

[0012] In one embodiment, the step of generating the homogeneous velocity field includes:

[0013] Determine the homogeneous velocity field according to the background velocity of the target work area;

[0014] The error between the analytical travel time calculated according to the uniform velocity field and the true travel time of the velocity model is within a preset range.

[0015] In one embodiment, before generating the analytical travel time field corresponding to the velocity model according to the preset uniform velocity field and the velocity model, it further includes:

[0016] Determine the undulating surface in the velocity model;

[0017] In the velocity model, perform velocity filling on the strata above the undulating surface;

[0018] In the velocity model, perform velocity smoothing on the strata below the undulating surface.

[0019] In one embodiment, the method for generating the first-order accuracy smooth scaling factor according to the analytical travel time field by using the first-order solution of the factored eikonal equation includes:

[0020] Using the first-order accuracy fast marching solution of the factored eikonal equation, calculate the first arrival travel time according to the initial smooth scaling factor, the source position of the target block, and the positions of the geophones, so as to generate the first-order accuracy smooth scaling factor.

[0021] In one embodiment, the method for determining the high-order accuracy smooth scaling factor and travel time field of the target work area according to the first-order accuracy smooth scaling factor by using the high-order accuracy fast marching solution of the factored eikonal equation includes:

[0022] Set the first-order accuracy smooth scaling factor as the initial smooth scaling factor of the high-order accuracy factored eikonal equation;

[0023] Input the initial smooth scaling factor into the high-order accuracy factored eikonal equation to determine the high-order accuracy smooth scaling factor and travel time field of the target work area.

[0024] In one embodiment, the initial smooth scaling factor, the first-order accuracy smooth scaling factor, and the high-order accuracy smooth scaling factor are scale coefficient fields.

[0025] In a second aspect, the present invention provides a seismic wave travel time calculation device, and the device includes:

[0026] A travel time field acquisition unit, configured to acquire the analytical travel time field of a homogeneous medium in a target work area,

[0027] A first-order factor generation unit, configured to generate a first-order accuracy smooth scaling factor according to the analytical travel time field by using the first-order solution of the factored eikonal equation;

[0028] A travel-time field determination unit, which is used to utilize the high-order accurate fast sweeping method for the factorized eikonal equation to determine the high-order accurate smooth scaling factor and the travel-time field of the target work area according to the first-order accurate smooth scaling factor.

[0029] In one embodiment, the travel-time field acquisition unit includes:

[0030] A velocity model acquisition module, which is used to acquire the velocity model of the target work area;

[0031] A travel-time field generation module, which is used to generate the analytical travel-time field corresponding to the velocity model according to the pre-generated uniform velocity field and the velocity model.

[0032] In one embodiment, the seismic wave travel-time calculation device further includes: a velocity field generation unit, which is used to generate the uniform velocity field, and the velocity field generation unit includes:

[0033] A uniform velocity field determination module, which is used to determine the uniform velocity field according to the background velocity of the target work area;

[0034] The error between the analytical travel-time calculated according to the uniform velocity field and the true travel-time of the velocity model is within a preset range.

[0035] In one embodiment, the seismic wave travel-time calculation device further includes:

[0036] An undulating surface determination unit, which is used to determine the undulating surface in the velocity model;

[0037] A velocity filling unit, which is used to fill the velocity of the strata above the undulating surface in the velocity model;

[0038] A velocity smoothing unit, which is used to smooth the velocity of the strata below the undulating surface in the velocity model.

[0039] In one embodiment, the first-order factor generation unit is specifically used to utilize the first-order accurate fast sweeping method for the factorized eikonal equation to calculate the first arrival travel-time according to the initial smooth scaling factor, the source position of the target block, and the position of the geophone, so as to generate the first-order accurate smooth scaling factor.

[0040] In one embodiment, the travel-time field determination unit includes:

[0041] A scaling factor setting module, which is used to set the first-order accurate smooth scaling factor as the initial smooth scaling factor of the high-order accurate factorized eikonal equation;

[0042] A travel-time field determination single module is used to input the initial smooth scaling factor into the high-order accurate factorization eikonal equation to determine the high-order accurate smooth scaling factor and travel-time field of the target work area.

[0043] In one embodiment, the initial smooth scaling factor, the first-order accurate smooth scaling factor, and the high-order accurate smooth scaling factor are scale coefficient fields.

[0044] In a third aspect, the present invention provides an electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the program, the steps of the seismic wave travel-time calculation method are implemented.

[0045] In a fourth aspect, the present invention provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the steps of the seismic wave travel-time calculation method are implemented.

[0046] As can be seen from the above description, for the seismic wave travel-time calculation method and device provided by the embodiments of the present invention, first, an analytical travel-time field of a homogeneous medium in the target work area is obtained. Then, according to the analytical travel-time field, a first-order accurate smooth scaling factor is generated by using the first-order solution of the factorization eikonal equation. Finally, by using the high-order accurate fast sweeping solution of the factorization eikonal equation, the high-order accurate smooth scaling factor and travel-time field of the target work area are determined based on the first-order accurate smooth scaling factor. This application directly solves the eikonal equation to solve the problem of the shadow area; and by adopting the factorization eikonal equation, the problem of source singularity is solved; and a method combining high-order and low-order solutions is adopted to improve the accuracy while accelerating the calculation speed. BRIEF DESCRIPTION OF THE DRAWINGS

[0047] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required for the description of the embodiments or the prior art. Obviously, the following drawings are some embodiments of the present invention. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0048] Figure 1 It is a schematic flowchart of the seismic wave travel-time calculation method provided in the embodiments of the present invention Figure 1 ;

[0049] Figure 2 It is a schematic flowchart of step 100 of the seismic wave travel-time calculation method in the embodiments of the present invention;

[0050] Figure 3 It is a schematic flowchart of generating the homogeneous velocity field in the embodiments of the present invention;

[0051] Figure 4 Schematic flowchart of the seismic wave travel time calculation method provided in the embodiment of the present invention Figure 2 ;

[0052] Figure 5 Schematic flowchart of step 200 of the seismic wave travel time calculation method in the embodiment of the present invention;

[0053] Figure 6 Schematic flowchart of step 300 of the seismic wave travel time calculation method in the embodiment of the present invention;

[0054] Figure 7 Schematic flowchart of the seismic wave travel time calculation method in the specific application example of the present invention;

[0055] Figure 8 Schematic diagram of the Marmousi velocity model in the specific application example of the present invention;

[0056] Figure 9 Schematic diagram of the travel time field of the uniform velocity model in the specific application example of the present invention;

[0057] Figure 10 Schematic diagram of the smooth scaling factor in the specific application example of the present invention;

[0058] Figure 11 Schematic diagram of the high-order scaling factor in the specific application example of the present invention;

[0059] Figure 12 Schematic diagram of the scaling factor update amount in the specific application example of the present invention;

[0060] Figure 13 Schematic diagram of the high-precision travel time of the Marmousi model in the specific application example of the present invention;

[0061] Figure 14 Schematic structure of the seismic wave travel time calculation device in the embodiment of the present invention Figure 1 ;

[0062] Figure 15 Schematic diagram of the structure of the time field acquisition unit in the embodiment of the present invention;

[0063] Figure 16 Schematic structure of the seismic wave travel time calculation device in the embodiment of the present invention Figure 2 ;

[0064] Figure 17 Schematic diagram of the structure of the velocity field generation unit in the embodiment of the present invention;

[0065] Figure 18 Schematic structure of the seismic wave travel time calculation device in the embodiment of the present invention Figure 3 ;

[0066] Figure 19 This is a schematic structural diagram of the travel-time field determination unit during travel in the embodiments of the present invention;

[0067] Figure 20 This is a schematic structural diagram of the electronic device in the embodiments of the present invention. Specific embodiments

[0068] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are some, but not all, of the embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0069] The embodiments of the present invention also provide a specific implementation manner of a seismic wave travel-time calculation method. Refer to Figure 1 , and the method specifically includes the following content:

[0070] Step 100: Obtain the analytical travel-time field of the homogeneous medium in the target work area.

[0071] Step 200: Use the first-order solution of the factored eikonal equation to generate a first-order accurate smooth scaling factor according to the analytical travel-time field.

[0072] Specifically, apply a first-order accurate operator to the factored eikonal equation, perform first-arrival travel-time calculation according to the positions of the source and the geophone, and output a first-order accurate smooth scaling factor.

[0073] Step 300: Use the high-order accurate fast marching method of the factored eikonal equation to determine the high-order accurate smooth scaling factor and travel-time field of the target work area according to the first-order accurate smooth scaling factor.

[0074] In step 300, regard the iterative process of solving the factored eikonal equation as an inversion process, reduce the number of iterations by improving the accuracy of the initial model; use the first-order accurate smooth scaling factor as the input of the iterative process of the high-order accurate fast marching method, so as to accelerate the iterative convergence and generate the high-order accurate smooth scaling factor and travel-time field of the target work area.

[0075] As can be seen from the above description, the seismic wave travel time calculation method provided by the embodiments of the present invention first obtains the analytical travel time field of the homogeneous medium in the target work area. Then, according to the analytical travel time field, a first-order accurate smooth scaling factor is generated by using the first-order solution of the factored eikonal equation. Finally, by using the high-order accurate fast marching method of the factored eikonal equation, the high-order accurate smooth scaling factor and the travel time field of the target work area are determined according to the first-order accurate smooth scaling factor. This application directly solves the eikonal equation to solve the problem of the shadow area; and by adopting the factored eikonal equation, the problem of source singularity is solved; a method combining high-order and low-order solutions is adopted to improve the accuracy while accelerating the calculation speed.

[0076] In one embodiment, referring to Figure 2 , step 100 further includes:

[0077] Step 101: Obtain the velocity model of the target work area;

[0078] Step 102: Generate the analytical travel time field corresponding to the velocity model according to the pre-generated homogeneous velocity field and the velocity model.

[0079] First, obtain a suitable homogeneous velocity field in the target work area; based on this homogeneous velocity field, obtain the analytical travel time field corresponding to this velocity field.

[0080] In one embodiment, referring to Figure 3 , the steps of generating the homogeneous velocity field include:

[0081] Step S100: Determine the homogeneous velocity field according to the background velocity of the target work area;

[0082] It can be understood that the error between the analytical travel time calculated according to the homogeneous velocity field in step S100 and the true travel time of the velocity model needs to be within a preset range.

[0083] In one embodiment, referring to Figure 4 , before step 300, it further includes:

[0084] Step 400: Determine the undulating ground surface in the velocity model;

[0085] Step 500: Fill the velocity of the strata above the undulating ground surface in the velocity model;

[0086] Step 600: Smooth the velocity of the strata below the undulating ground surface in the velocity model.

[0087] It can be understood that between step 200 and step 300, the velocity model also needs to be processed, including the establishment of the undulating ground surface, the velocity filling above the undulating ground surface, and the velocity smoothing below the undulating ground surface.

[0088] In one embodiment, referring to Figure 5 , step 200 specifically includes:

[0089] Step 201: Using the first-order accuracy fast scanning solution of the factorized eikonal equation, calculate the first arrival travel time according to the initial smooth scaling factor, the source position of the target block, and the position of the geophone, so as to generate the first-order accuracy smooth scaling factor.

[0090] In the prior art, when obtaining the travel time using the eikonal equation, due to the general use of point sources, the problem of source singularities will occur, that is, the travel time is non-differentiable at the source, and there is a large truncation error in the directly obtained result, which may even lead to non-convergence. The factorized eikonal equation then came into being. The so-called factorized eikonal equation is to consider the travel time as the product of the constant-velocity analytical travel time and the smooth scaling factor. The constant-velocity travel time can be analytically obtained, and the scaling factor is obtained by iteratively solving the factorized eikonal equation. Multiplying the constant-velocity travel time by the scaling factor can obtain the travel time field. The factorized eikonal equation can obtain a more accurate travel time field and effectively solve the problem of source singularity.

[0091] In one embodiment, referring to Figure 6 , step 300 specifically includes:

[0092] Step 301: Set the first-order accuracy smooth scaling factor as the initial smooth scaling factor of the high-order accuracy factorized eikonal equation;

[0093] It can be understood that the factorized eikonal equation decomposes the travel time field into the product of the homogeneous model analytical travel time field and the smooth scaling factor, and iteratively improves the smooth scaling factor through the fast scanning method to obtain an accurate travel time field;

[0094] Step 302: Input the initial smooth scaling factor into the high-order accuracy factorized eikonal equation to determine the high-order accuracy smooth scaling factor and travel time field of the target work area.

[0095] It should be noted that the output smooth scaling factor and travel time field are of high order, and both the smooth scaling factor field and the travel time field have the application conditions for undulating ground;

[0096] In one embodiment, the initial smooth scaling factor, the first-order accuracy smooth scaling factor, and the high-order accuracy smooth scaling factor are scale coefficient fields.

[0097] It can be understood that multiplying the smooth scaling factor of the complex medium by the travel time field of the homogeneous medium can obtain the corresponding travel time field of the complex medium.

[0098] In view of the following problems existing in calculating the travel time of seismic waves in the prior art: the ray theory method has low accuracy and causes the problem of shadow zones; the conventional finite-difference solution method of the eikonal equation has the problem of source singularity and limited accuracy; when using high-order finite differences, the calculation efficiency will be reduced. The seismic wave travel time calculation method provided by the embodiments of the present application can directly solve the eikonal equation and solve the problem of shadow zones; adopt the factored eikonal equation to solve the problem of source singularity; adopt a method combining high-order and low-order solutions to accelerate the calculation speed while improving the accuracy.

[0099] To further illustrate the present solution, the present application provides a specific application example of the seismic wave travel time calculation method, which specifically includes the following content. See Figure 7 .

[0100] It can be understood that the factored eikonal equation (Formula (1)), that is, considering the travel time as the product of the constant-velocity analytical travel time and the smooth scaling factor. For the constant-velocity travel time, it can be analytically obtained. By iteratively solving the factored eikonal equation, the scaling factor is obtained, and multiplying it by the constant-velocity travel time can obtain the travel time field. The factored eikonal equation can obtain a more accurate travel time field and effectively solve the problem of source singularity.

[0101]

[0102] In the above formula, s is the velocity of the underground medium, s0 is the uniform velocity field, α is the velocity scaling factor; τ is the travel time, τ0 is the travel time corresponding to s0, and w is the smooth scaling factor of the travel time.

[0103] S1: Generate a first-order accurate smooth scaling factor.

[0104] Specifically, based on the input velocity model ( Figure 8 ), a suitable uniform velocity field is obtained; based on the uniform velocity field, the corresponding analytical travel time field ( Figure 9 ) of this velocity field is obtained; based on the initial smooth scaling factor, using the first-order accurate fast sweeping solution method of the factored eikonal equation, a first-order accurate smooth scaling factor ( Figure 10 ) is obtained;

[0105] S2: Obtain the initial input of the high-order accurate iterative solution method.

[0106] Based on the first-order accurate smooth scaling factor, after processing, the initial input of the high-order accurate iterative solution method is obtained (see Formula (2));

[0107]

[0108] Solving the eikonal equation of the factorization type involves a large amount of computation and requires multiple iterations to converge stably. To improve the computational speed, the iterative process is regarded as an inversion process, and the initial smooth scaling factor is improved to accelerate the iterative convergence. The first-order scaling factor ( Figure 10 ) is used as the initial input for the high-order calculation unit in the next stage. Figure 10 ).

[0109] S3: Obtain the high-order accuracy smooth scaling factor and the travel time field.

[0110] Specifically, based on the first-order accuracy smooth scaling factor ( Figure 10 ), the high-order accuracy fast sweeping method for the eikonal equation of the factorization type is used to obtain the high-order accuracy smooth scaling factor ( Figure 11 ) and the travel time field ( Figure 13 ). Figure 10 ) Figure 11 ) and the travel time field ( Figure 13 ).

[0111] From Figures 11 to 13 it can be seen that by substituting the smooth scaling factor obtained by the first-order accuracy method into the high-order accuracy method, high-precision travel times can be obtained. At the same time, Table 1 shows the comparison of the number of iterations and computational time for the first-order method, high-order method, and the hybrid fast solution of this application.

[0112] Table 1 Comparison of computational efficiency

[0113] Number of iterations Error Time First-order solution 26 1e-4 2s Higher-order solution 1702 1e-6 18s Mixed-strategy solution 1013 1e-6 10s

[0114] From the above process, it can be concluded (Table 1) that the seismic wave travel time calculation method provided by this specific application example can effectively improve the computational efficiency while ensuring the accuracy.

[0115] From the above description, it can be seen that the seismic wave travel time calculation method provided by the specific application example of the present invention accelerates the calculation process while ensuring the calculation accuracy. Specifically, first, a velocity model is obtained; then, the first-order accuracy fast sweeping method for the eikonal equation of the factorization type is used to output the first-order accuracy travel time and the first-order accuracy smooth scaling factor; then, based on the first-order accuracy smooth scaling factor, it is processed and used as the input for the high-order accuracy fast sweeping method; finally, the eikonal equation of the factorization type is used, and the high-order fast sweeping method is used for calculation to output the high-order accuracy travel time and the high-order accuracy smooth scaling factor. The seismic wave travel time calculation method provided by the specific application example of the present invention can accelerate the calculation speed while ensuring high-order accuracy, and can be used to improve the imaging accuracy and near-surface velocity modeling.

[0116] Based on the same inventive concept, an embodiment of the present application further provides a seismic wave travel time calculation device, which can be used to implement the method described in the above embodiment, as in the following embodiment. Since the principle of solving problems by the seismic wave travel time calculation device is similar to that of the seismic wave travel time calculation method, the implementation of the seismic wave travel time calculation device can refer to the implementation of the seismic wave travel time calculation method, and the repeated parts will not be elaborated. As used hereinafter, the term "unit" or "module" can be a combination of software and / or hardware that can achieve a predetermined function. Although the system described in the following embodiments is preferably implemented in software, implementation in hardware, or a combination of software and hardware is also possible and contemplated.

[0117] An embodiment of the present invention provides a specific implementation manner of a seismic wave travel time calculation device capable of implementing the seismic wave travel time calculation method, see Figure 14 , the seismic wave travel time calculation device specifically includes the following content:

[0118] A time field acquisition unit 10, configured to acquire an analytical travel time field of a homogeneous medium in a target work area,

[0119] A first-order factor generation unit 20, configured to generate a first-order accuracy smooth scaling factor according to the analytical travel time field by using a first-order solution of the factored eikonal equation;

[0120] A travel time field determination unit 30, configured to determine a high-order accuracy smooth scaling factor and a travel time field of the target work area according to the first-order accuracy smooth scaling factor by using a high-order accuracy fast sweeping solution of the factored eikonal equation.

[0121] In one embodiment, see Figure 15 , the time field acquisition unit 10 includes:

[0122] A velocity model acquisition module 101, configured to acquire the velocity model of the target work area;

[0123] A time field generation module 102, configured to generate an analytical travel time field corresponding to the velocity model according to a pre-generated uniform velocity field and the velocity model.

[0124] In one embodiment, see Figure 16 , the seismic wave travel time calculation device further includes: a velocity field generation unit 40, configured to generate the uniform velocity field, see Figure 17 , the velocity field generation unit 40 includes:

[0125] A uniform velocity field determination module 401, configured to determine the uniform velocity field according to the background velocity of the target work area;

[0126] The error between the analytical travel time calculated according to the uniform velocity field and the true travel time of the velocity model is within a preset range.

[0127] In one embodiment, referring to Figure 18 , the seismic wave travel time calculation device further includes:

[0128] An undulating surface determination unit 50 for determining the undulating surface in the velocity model;

[0129] A velocity filling unit 60 for filling the velocity of the strata above the undulating surface in the velocity model;

[0130] A velocity smoothing unit 70 for smoothing the velocity of the strata below the undulating surface in the velocity model.

[0131] In one embodiment, the first-order factor generation unit is specifically configured to use the first-order accuracy fast sweeping method for the factored eikonal equation to calculate the first arrival travel time according to the initial smooth scaling factor, the source position of the target block, and the position of the geophone, so as to generate the first-order accuracy smooth scaling factor.

[0132] In one embodiment, referring to Figure 19 , the travel time field determination unit 30 includes:

[0133] A scaling factor setting module 301 for setting the first-order accuracy smooth scaling factor as the initial smooth scaling factor of the high-order accuracy factored eikonal equation;

[0134] A travel time field determination single module 302 for inputting the initial smooth scaling factor into the high-order accuracy factored eikonal equation to determine the high-order accuracy smooth scaling factor and travel time field of the target work area.

[0135] In one embodiment, the initial smooth scaling factor, the first-order accuracy smooth scaling factor, and the high-order accuracy smooth scaling factor are scale coefficient fields.

[0136] As can be seen from the above description, the seismic wave travel time calculation device provided by the embodiment of the present invention first obtains the analytical travel time field of the homogeneous medium in the target work area, and then generates the first-order accuracy smooth scaling factor according to the analytical travel time field by using the first-order solution of the factored eikonal equation; finally, uses the high-order accuracy fast sweeping method of the factored eikonal equation to determine the high-order accuracy smooth scaling factor and travel time field of the target work area according to the first-order accuracy smooth scaling factor. This application directly solves the eikonal equation to solve the problem of the shadow area; and solves the problem of source singularity by adopting the factored eikonal equation; adopts a method combining high and low order solutions to improve the accuracy while accelerating the calculation speed.

[0137] Embodiments of the present application also provide a specific implementation manner of an electronic device that can implement all steps in the seismic wave travel time calculation method in the above embodiments. Refer to Figure 20 , the electronic device specifically includes the following content:

[0138] A processor 1201, a memory 1202, a communication interface 1203, and a bus 1204;

[0139] Among them, the processor 1201, the memory 1202, and the communication interface 1203 communicate with each other through the bus 1204; the communication interface 1203 is used to implement information transmission between related devices such as server-side devices, power measurement devices, and user-side devices.

[0140] The processor 1201 is used to call the computer program in the memory 1202. When the processor executes the computer program, all steps in the seismic wave travel time calculation method in the above embodiments are implemented. For example, when the processor executes the computer program, the following steps are implemented:

[0141] Step 100: Obtain the analytical travel time field of the homogeneous medium in the target work area,

[0142] Step 200: Use the first-order solution of the factored eikonal equation to generate a first-order accurate smooth scaling factor according to the analytical travel time field;

[0143] Step 300: Use the high-order accurate fast marching solution of the factored eikonal equation to determine the high-order accurate smooth scaling factor and travel time field of the target work area according to the first-order accurate smooth scaling factor.

[0144] Embodiments of the present application also provide a computer-readable storage medium that can implement all steps in the seismic wave travel time calculation method in the above embodiments. A computer program is stored on the computer-readable storage medium. When the computer program is executed by a processor, all steps in the seismic wave travel time calculation method in the above embodiments are implemented. For example, when the processor executes the computer program, the following steps are implemented:

[0145] Step 100: Obtain the analytical travel time field of the homogeneous medium in the target work area,

[0146] Step 200: Use the first-order solution of the factored eikonal equation to generate a first-order accurate smooth scaling factor according to the analytical travel time field;

[0147] Step 300: Use the high-order accurate fast marching solution of the factored eikonal equation to determine the high-order accurate smooth scaling factor and travel time field of the target work area according to the first-order accurate smooth scaling factor.

[0148] In summary, the computer-readable storage medium provided by the embodiments of the present invention can support the service provider to perform self-adaptive offline and online of services according to the availability of its own software and hardware resources, realize the self-isolation ability of the service provider, and ensure the success rate of the service provider's response to service requests.

[0149] Each embodiment in this specification is described in a progressive manner. The same or similar parts among the embodiments can be referred to each other, and the key points of each embodiment are the differences from other embodiments. In particular, for the hardware + program type embodiments, since they are basically similar to the method embodiments, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiments.

[0150] The specific embodiments of this specification have been described above. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recited in the claims can be performed in a different order than in the embodiments and still achieve the desired results. Additionally, the processes depicted in the figures do not necessarily require the specific order or sequential order shown to achieve the desired results. In certain embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0151] Although this application provides method operation steps such as in the embodiments or flowcharts, it may include more or fewer operation steps based on routine or non-creative labor. The step order listed in the embodiments is only one way among the execution orders of numerous steps and does not represent the only execution order. When the actual device or client product is executed, it can be executed in the order shown in the embodiments or the figures or in parallel (such as in an environment of parallel processors or multi-threaded processing).

[0152] Those skilled in the art should understand that the embodiments of the present invention can be provided as a method, a system, or a computer program product. Therefore, the present invention can take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present invention can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk memories, CD-ROMs, optical memories, etc.) containing computer-usable program code.

[0153] It should also be noted that the term "comprising", "including" or any other variant thereof is intended to cover non-exclusive inclusion, such that a process, method, commodity or device comprising a series of elements not only includes those elements but also includes other elements not expressly listed, or elements inherent to such process, method, commodity or device. Without further limitation, an element defined by the statement "comprising an..." does not exclude the presence of additional identical elements in the process, method, commodity or device comprising said element.

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

[0155] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, such that the instructions stored in the computer-readable memory produce a manufactured article including an instruction means that implements the functions specified in Figure 1 one or more of the flows Figure 1 or blocks or a combination of blocks.

[0156] These computer program instructions can also be loaded onto a computer or other programmable data processing device, such that a series of operation steps are executed on the computer or other programmable device to produce a computer-implemented process, and thus the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in Figure 1 one or more of the flows Figure 1 or blocks or a combination of blocks.

[0157] Specific embodiments are applied in the present invention to elaborate on the principles and implementation manners of the present invention. The description of the above embodiments is only used to help understand the method of the present invention and its core idea; at the same time, for those of ordinary skill in the art, according to the idea of the present invention, there will be changes in the specific implementation manners and application scopes. In summary, the content of this specification should not be construed as a limitation to the present invention.

Claims

1. A method for calculating seismic wave travel time, characterized in that, it includes: Obtain the analytical travel time field of the homogeneous medium in the target work area, Using the first-order solution of the factored eikonal equation, generate a first-order accurate smooth scaling factor according to the analytical travel time field; Using the high-order accurate fast marching method of the factored eikonal equation, determine the high-order accurate smooth scaling factor and travel time field of the target work area according to the first-order accurate smooth scaling factor; The step of using the first-order solution of the factored eikonal equation to generate a first-order accurate smooth scaling factor according to the analytical travel time field includes: Using the first-order accurate fast marching method of the factored eikonal equation, calculate the first arrival travel time according to the initial smooth scaling factor, the source position in the target work area, and the position of the geophone, so as to generate the first-order accurate smooth scaling factor; The step of using the high-order accurate fast marching method of the factored eikonal equation to determine the high-order accurate smooth scaling factor and travel time field of the target work area according to the first-order accurate smooth scaling factor includes: Set the first-order accurate smooth scaling factor as the initial smooth scaling factor of the high-order accurate factored eikonal equation; Input the initial smooth scaling factor into the high-order accurate factored eikonal equation to determine the high-order accurate smooth scaling factor and travel time field of the target work area.

2. The seismic wave travel time calculation method according to claim 1, characterized in that, The step of obtaining the analytical travel time field of the homogeneous medium in the target work area includes: Obtain the velocity model of the target work area; Generate the analytical travel time field corresponding to the velocity model according to the pre-generated homogeneous velocity field and the velocity model.

3. The seismic wave travel time calculation method according to claim 2, characterized in that, The step of generating the homogeneous velocity field includes: Determine the homogeneous velocity field according to the background velocity of the target work area; The error between the analytical travel time calculated according to the homogeneous velocity field and the true travel time of the velocity model is within a preset range.

4. The seismic wave travel time calculation method according to claim 2, characterized in that, Before generating the analytical travel time field corresponding to the velocity model according to the pre-generated homogeneous velocity field and the velocity model, it further includes: Determine the undulating surface in the velocity model; In the velocity model, fill the velocity of the strata above the undulating surface; In the velocity model, smooth the velocity of the strata below the undulating surface.

5. The seismic wave travel time calculation method according to claim 1, characterized in that, The initial smooth scaling factor, the first-order accurate smooth scaling factor, and the high-order accurate smooth scaling factor are scale coefficient fields.

6. A seismic wave travel time calculation device, characterized in that, it includes: A travel time field acquisition unit for acquiring the analytical travel time field of the homogeneous medium in the target work area, A first-order factor generation unit for using the first-order solution of the factored eikonal equation to generate a first-order accurate smooth scaling factor according to the analytical travel time field; A travel-time field determination unit, which uses a high-order accuracy fast marching method for the factorized eikonal equation to determine the high-order accuracy smooth scaling factor and the travel-time field of the target work area according to the first-order accuracy smooth scaling factor; The first-order factor generation unit is specifically configured to use the first-order accuracy fast marching method for the factorized eikonal equation to calculate the first arrival travel time according to the initial smooth scaling factor, the source position of the target work area, and the position of the geophone, so as to generate the first-order accuracy smooth scaling factor; The travel-time field determination unit includes: A scaling factor setting module, which is configured to set the first-order accuracy smooth scaling factor as the initial smooth scaling factor of the high-order accuracy factorized eikonal equation; A travel-time field determination single module, which is configured to input the initial smooth scaling factor into the high-order accuracy factorized eikonal equation to determine the high-order accuracy smooth scaling factor and the travel-time field of the target work area.

7. The seismic wave travel-time calculation device according to claim 6, wherein, The time field acquisition unit includes: A velocity model acquisition module, which is configured to acquire the velocity model of the target work area; A time field generation module, which is configured to generate an analytical travel-time field corresponding to the velocity model according to a pre-generated uniform velocity field and the velocity model.

8. The seismic wave travel-time calculation device according to claim 7, wherein, It further includes: A velocity field generation unit, which is configured to generate the uniform velocity field, and the velocity field generation unit includes: A uniform velocity field determination module, which is configured to determine the uniform velocity field according to the background velocity of the target work area; The error between the analytical travel time calculated according to the uniform velocity field and the true travel time of the velocity model is within a preset range.

9. The seismic wave travel-time calculation device according to claim 7, wherein, It further includes: An undulating surface determination unit, which is configured to determine the undulating surface in the velocity model; A velocity filling unit, which is configured to fill the velocity of the formation above the undulating surface in the velocity model; A velocity smoothing unit, which is configured to smooth the velocity of the formation below the undulating surface in the velocity model.

10. The seismic wave travel-time calculation device according to claim 6, wherein, The initial smooth scaling factor, the first-order accuracy smooth scaling factor, and the high-order accuracy smooth scaling factor are scale coefficient fields.

11. An electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein, When the processor executes the program, it implements the steps of the seismic wave travel-time calculation method according to any one of claims 1 to 5.

12. A computer-readable storage medium, on which a computer program is stored, wherein, When the computer program is executed by a processor, it implements the steps of the seismic wave travel-time calculation method according to any one of claims 1 to 5.

Images