阅读说明：本技术 积分法叠前深度偏移方法及装置 (Integral method prestack depth migration method and device ) 是由 岳玉波 张建磊 王成祥 金守利 张巍毅 孟阳 于 2020-03-31 设计创作，主要内容包括：本发明提供了一种积分法叠前深度偏移方法及装置,该方法包括：根据地震数据和计算参数,确定每一地层深度对应的旅行时表步长；在地层深度方向上进行采样,得到目标区域的旅行时表；根据每一地层深度对应的旅行时表步长和旅行时表,对采样点进行线性插值,进行积分法处理,得到目标区域地震数据的偏移结果；相较于现有技术中较小且固定的旅行时表的步长,确定每一地层深度对应的旅行时表步长,即旅行时表步长随地层深度变化而变化,在不同的地层深度处,选取不同的步长,即在不同的地层深度处,旅行时表的网格深度间隔不同,使得在保证成像精度的前提下,减少偏移成像的计算量和执行计算的硬盘的I/O吞吐量,提高了偏移效率。(The invention provides a method and a device for prestack depth migration by an integral method, wherein the method comprises the following steps: determining the travel time table step length corresponding to each stratum depth according to the seismic data and the calculation parameters; sampling in the depth direction of the stratum to obtain a travel time table of a target area; according to the travel time table step length and the travel time table corresponding to each stratum depth, carrying out linear interpolation on the sampling points, and carrying out integral processing to obtain the migration result of the seismic data of the target area; compared with the step length of a small and fixed travel time table in the prior art, the step length of the travel time table corresponding to each stratum depth is determined, namely the step length of the travel time table changes along with the change of the stratum depth, different step lengths are selected at different stratum depths, namely grid depth intervals of the travel time table are different at different stratum depths, so that the calculation amount of offset imaging and the I/O throughput of a hard disk for executing calculation are reduced on the premise of ensuring the imaging precision, and the offset efficiency is improved.)

1. An integration prestack depth migration method, comprising:

acquiring seismic data of a target area;

determining the travel time table step length corresponding to each stratum depth according to the seismic data and the calculation parameters;

sampling in the depth direction of the stratum according to the step length of the travel time table corresponding to each stratum depth to obtain the travel time table of the target area;

and performing linear interpolation on sampling points according to the travel time table step length corresponding to each stratum depth and the travel time table, and performing integral processing to obtain the migration result of the seismic data of the target area.

2. The method of claim 1, wherein the calculating parameters comprises at least: the initial step length, the stratum depth and the effective frequency band of the travel time table;

determining the travel time table step length corresponding to each stratum depth according to the seismic data and the calculation parameters, wherein the steps comprise:

obtaining the maximum effective frequency of the seismic data of the target area according to the effective frequency band;

determining a minimum value of the velocity model at each formation depth according to the velocity model and the formation depth;

determining an adjustment parameter at each formation depth according to the formation depth;

and determining the travel time table step length corresponding to each stratum depth according to the maximum effective frequency of the seismic data of the target area, the minimum value of the velocity model at each stratum depth, the depth interval of the imaging points, the adjusting parameters at each stratum depth and the initial step length of the travel time table.

3. The method of claim 2, wherein the travel time table step size corresponding to each formation depth is determined based on the maximum effective frequency of the seismic data for the target area, the minimum value of the velocity model at each formation depth, the imaging point depth separation, the adjustment parameters at each formation depth, and the travel time table initial step size according to the following formula:

wherein L represents the formation depth; TMINC [ L ]]Representing the step length of a travel time table corresponding to the stratum depth L; INT represents an integer function; v. of_min[L]Represents the minimum of the velocity model at the formation depth L; para [ L]Represents a conditioning parameter at the formation depth L; f. of_maxRepresenting a maximum effective frequency of seismic data for the target area; dz represents the depth interval of the imaging points and is a preset value; TMINC0 represents the initial step size of the travel time table, which is a preset value.

4. A method according to claim 3, characterized in that the adjustment parameter at each formation depth is determined in dependence on the formation depth according to the following formula:

wherein para [ L ] represents the conditioning parameter at the formation depth L; l represents the formation depth.

5. The method as claimed in claim 3, wherein the step length of the travel time table corresponding to each stratum depth and the travel time table are used for carrying out linear interpolation on sampling points and carrying out integral processing to obtain the migration result of the seismic data of the target area, and the method comprises the following steps:

when linear interpolation is carried out on the sampling points, the number of interpolation sampling points between any two sampling points is TMINC [ L ] -2.

6. An integration prestack depth migration apparatus, comprising:

the data acquisition module is used for acquiring seismic data of a target area;

the step length determining module is used for determining the step length of the travel time table corresponding to each stratum depth according to the seismic data and the calculation parameters;

the sampling module is used for sampling in the depth direction of the stratum according to the step length of the travel time table corresponding to each stratum depth to obtain the travel time table of the target area;

and the data migration module is used for performing linear interpolation on sampling points according to the travel time table step length corresponding to each stratum depth and the travel time table, and performing integral processing to obtain a migration result of the seismic data of the target area.

7. The apparatus of claim 6, wherein the calculating parameters comprises at least: the initial step length, the stratum depth and the effective frequency band of the travel time table;

the step length determining module is specifically configured to:

obtaining the maximum effective frequency of the seismic data of the target area according to the effective frequency band;

determining a minimum value of the velocity model at each formation depth according to the velocity model and the formation depth;

determining an adjustment parameter at each formation depth according to the formation depth;

and determining the travel time table step length corresponding to each stratum depth according to the maximum effective frequency of the seismic data of the target area, the minimum value of the velocity model at each stratum depth, the depth interval of the imaging points, the adjusting parameters at each stratum depth and the initial step length of the travel time table.

8. The apparatus of claim 7, wherein the traveller table step size for each formation depth is determined based on the maximum effective frequency of the target zone seismic data, the minimum value of the velocity model at each formation depth, the imaging point depth interval, the adjustment parameters at each formation depth, and the traveller table initial step size according to the following equation:

wherein L represents the formation depth; TMINC [ L ]]Representing the step length of a travel time table corresponding to the stratum depth L; INT represents an integer function; v. of_min[L]Represents the minimum of the velocity model at the formation depth L; para [ L]Represents a conditioning parameter at the formation depth L; f. of_maxRepresenting a maximum effective frequency of seismic data for the target area; dz represents the depth interval of the imaging points and is a preset value; TMINC0 represents the initial step size of the travel time table, which is a preset value.

9. The apparatus of claim 8, wherein the adjustment parameter at each formation depth is determined based on the formation depth according to the following formula:

wherein para [ L ] represents the conditioning parameter at the formation depth L; l represents the formation depth.

10. The apparatus of claim 8, wherein the data migration module is specifically configured to:

when linear interpolation is carried out on the sampling points, the number of interpolation sampling points between any two sampling points is TMINC [ L ] -2.

11. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any one of claims 1 to 5 when executing the computer program.

12. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method of any one of claims 1 to 5.

Technical Field

The invention relates to the technical field of geophysical exploration of petroleum, in particular to a method and a device for prestack depth migration by an integral method.

Background

The integral method of prestack depth migration is a practical migration method in industrial production, and the implementation process mainly comprises two parts, namely traveling calculation and superposition imaging processing along the traveling calculation. Therefore, the accuracy of offset imaging depends mainly on the calculation accuracy when traveling. However, in the integral prestack depth migration process, if the grid depth interval of the travel time table is selected too much, the deviation between interpolation travel time and real travel time is easy to be too large, and the imaging precision is affected; in order to meet imaging accuracy in the prior art, the step length of the travel time table is set to be the same value which is small enough, that is, the grid depth interval of the travel time table is small enough, so that the computation amount of offset imaging and the I/O throughput of a hard disk for performing computation are huge, and the offset efficiency is low.

Disclosure of Invention

The embodiment of the invention provides an integral prestack depth migration method, which is used for reducing the calculation amount of migration imaging and the I/O throughput of a hard disk executing calculation on the premise of ensuring the imaging precision and improving the migration efficiency and comprises the following steps:

acquiring seismic data of a target area;

determining the travel time table step length corresponding to each stratum depth according to the seismic data and the calculation parameters;

sampling in the depth direction of the stratum according to the step length of the travel time table corresponding to each stratum depth to obtain the travel time table of the target area;

and performing linear interpolation on sampling points according to the travel time table step length corresponding to each stratum depth and the travel time table, and performing integral processing to obtain the migration result of the seismic data of the target area.

An embodiment of the present invention further provides an integral prestack depth migration apparatus, configured to reduce computation amount of migration imaging and I/O throughput of a hard disk that performs computation on the premise of ensuring imaging accuracy, and improve migration efficiency, where the apparatus includes:

the data acquisition module is used for acquiring seismic data of a target area;

the step length determining module is used for determining the step length of the travel time table corresponding to each stratum depth according to the seismic data and the calculation parameters;

the sampling module is used for sampling in the depth direction of the stratum according to the step length of the travel time table corresponding to each stratum depth to obtain the travel time table of the target area;

and the data migration module is used for performing linear interpolation on sampling points according to the travel time table step length corresponding to each stratum depth and the travel time table, and performing integral processing to obtain a migration result of the seismic data of the target area.

The embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and when the processor executes the computer program, the integration prestack depth migration method is implemented.

Embodiments of the present invention also provide a computer-readable storage medium storing a computer program for executing the integration prestack depth migration method.

In the embodiment of the invention, the seismic data of the target area are obtained; determining the travel time table step length corresponding to each stratum depth according to the seismic data and the calculation parameters; sampling in the depth direction of the stratum according to the step length of the travel time table corresponding to each stratum depth to obtain the travel time table of the target area; according to the travel time table step length and the travel time table corresponding to each stratum depth, carrying out linear interpolation on the sampling points, and carrying out integral processing to obtain the migration result of the seismic data of the target area; compared with the step length of a small and fixed travel time table in the prior art, the step length of the travel time table corresponding to each stratum depth is determined, namely the step length of the travel time table changes along with the change of the stratum depth, different step lengths are selected at different stratum depths, namely grid depth intervals of the travel time table are different at different stratum depths, so that the calculation amount of offset imaging and the I/O throughput of a hard disk for executing calculation are reduced on the premise of ensuring the imaging precision, and the offset efficiency is improved.

Drawings

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

FIG. 1 is a diagram illustrating an integration prestack depth migration method according to an embodiment of the present disclosure.

FIG. 2 is a schematic diagram illustrating a method for determining a step size of a travel time table corresponding to each formation depth according to an embodiment of the present invention.

FIG. 3(a) is a diagram of an offset common reflection point gather in a fixed step travel time table in an embodiment of the present invention.

FIG. 3(b) is a schematic diagram of an offset common reflection point gather of a variable step travel time table in an embodiment of the present invention.

FIG. 3(c) is a cross-sectional view of the residual error between FIG. 3(a) and FIG. 3(b) in an embodiment of the present invention.

FIG. 4(a) is a schematic diagram of an offset common offset gather for a fixed step travel time table in an embodiment of the present invention.

FIG. 4(b) is a schematic diagram of an offset common offset gather for a variable stride travel time table in an embodiment of the present invention.

FIG. 4(c) is a cross-sectional view of the residual error between FIG. 4(a) and FIG. 4(b) in an embodiment of the present invention.

FIG. 5 is a block diagram of an integrated prestack depth migration apparatus according to an embodiment of the present invention.

Detailed Description

The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.

An embodiment of the present invention provides an integral prestack depth migration method, which is used to reduce the calculation amount of migration imaging and the I/O throughput of a hard disk executing calculation on the premise of ensuring imaging accuracy, and improve migration efficiency, as shown in fig. 1, and includes:

step 101: acquiring seismic data of a target area;

step 102: determining the travel time table step length corresponding to each stratum depth according to the seismic data and the calculation parameters;

step 103: sampling in the depth direction of the stratum according to the step length of the travel time table corresponding to each stratum depth to obtain the travel time table of the target area;

step 104: and performing linear interpolation on the sampling points according to the travel time table step length and the travel time table corresponding to each stratum depth, and performing integral processing to obtain the migration result of the seismic data of the target area.

As can be known from the process shown in FIG. 1, in the embodiment of the present invention, the seismic data of the target area is obtained; determining the travel time table step length corresponding to each stratum depth according to the seismic data and the calculation parameters; sampling in the depth direction of the stratum according to the step length of the travel time table corresponding to each stratum depth to obtain the travel time table of the target area; according to the travel time table step length and the travel time table corresponding to each stratum depth, carrying out linear interpolation on the sampling points, and carrying out integral processing to obtain the migration result of the seismic data of the target area; compared with the step length of a small and fixed travel time table in the prior art, the step length of the travel time table corresponding to each stratum depth is determined, namely the step length of the travel time table changes along with the change of the stratum depth, different step lengths are selected at different stratum depths, namely grid depth intervals of the travel time table are different at different stratum depths, so that the calculation amount of offset imaging and the I/O throughput of a hard disk for executing calculation are reduced on the premise of ensuring the imaging precision, and the offset efficiency is improved.

In specific implementation, firstly, the seismic data of a target area is acquired.

And after the seismic data are obtained, determining the travel time table step length corresponding to each stratum depth according to the seismic data and the calculation parameters. Wherein, calculating the parameters at least comprises: travel time table initial step size, formation depth and effective frequency band. The specific process is shown in fig. 2, and includes:

step 201: obtaining the maximum effective frequency of the seismic data of the target area according to the effective frequency band;

step 202: determining the minimum value of the velocity model at each stratum depth according to the velocity model and the stratum depth;

step 203: determining an adjustment parameter at each formation depth according to the formation depth;

step 204: and determining the travel time table step length corresponding to each stratum depth according to the maximum effective frequency of the seismic data of the target area, the minimum value of the velocity model at each stratum depth, the depth interval of the imaging points, the adjusting parameters at each stratum depth and the initial step length of the travel time table.

The above calculation parameters and the speed model may be calculated according to actual engineering, or may be initially given by a technician.

In this embodiment of the present invention, step 204 may determine the step size of the travel time table corresponding to each formation depth according to the following formula:

wherein L represents the formation depth; TMINC [ L ]]Pair at depth L of stratumThe step length of the corresponding travel time table; INT represents an integer function; v. of_min[L]Represents the minimum of the velocity model at the formation depth L; para [ L]Represents a conditioning parameter at the formation depth L; f. of_maxRepresenting a maximum effective frequency of seismic data for the target area; dz represents the depth interval of the imaging points and is a preset value; TMINC0 represents the initial step length of the travel time table, which is a preset value, and the engineering generally takes 4-8 times of dz value.

In the embodiment, para [ L ] is generally greater than or equal to 2, the larger the adjusting parameter is, the more severe the condition for increasing the step length is, and according to the trial and error of the applicant on a large amount of actual seismic data, the following formula is obtained, and the adjusting parameter para [ L ] at each stratum depth is determined according to the stratum depth L:

different adjusting parameter values are obtained according to different stratum depths, when the depth is shallow, the larger adjusting parameter value is selected to ensure the imaging precision of the shallow layer, and when the depth is larger, the smaller adjusting parameter value is selected to ensure the imaging efficiency of the deep layer.

And after the step length of the travel time table corresponding to each stratum depth is determined, sampling in the stratum depth direction according to the step length of the travel time table corresponding to each stratum depth, and obtaining the travel time table of the target area. The specific implementation process may be, for example: according to the initial step length of the travel time table, the travel time table of the target area corresponding to the initial step length is obtained through ray tracing, resampling is carried out according to the travel time table step length corresponding to each stratum depth in the travel time table of the target area corresponding to the initial step length, at the moment, the number of sample points of the travel time table after resampling in the depth direction is greatly reduced, the depth distance between the sample points is not equal any more and is changed into dz multiplied by TMINC [ L ], the travel time table step length TMINC [ L ] is changed along with the change of the stratum depth L, and the depth distance between the sample points is changed along with the change of the stratum depth, so that the travel time table with the variable step length is obtained.

And then, according to the step length of the travel time table corresponding to each stratum depth and the variable step length travel time table, carrying out linear interpolation on sampling points, and carrying out integral processing to obtain a prestack depth migration result of the seismic data of the target area. In the specific embodiment, when linear interpolation is carried out on the sampling points, the number of interpolation sampling points between any two sampling points is TMINC [ L ] -2. In a specific embodiment, a Kirchhoff integral shift algorithm is used when the integration method is performed.

It should be understood that the above calculation formula is only an example, and the above formula may be modified or other formulas or methods may be adopted in the implementation, and those skilled in the art may understand that the above formula is only an example, and the modification of the formula or other solving methods all fall into the protection scope of the present invention, and are not described in detail in the embodiments.

A specific example is given below to illustrate how the integration method prestack depth migration is performed according to an embodiment of the present invention. The method is applied to three-dimensional actual seismic data of a certain area to prove that the method improves the calculation efficiency on the basis of ensuring the migration effect.

The applicant finds that the requirement on the sampling rate is gradually relaxed due to the gradual reduction of the main frequency of the seismic signals along with the increase of the depth of the stratum, and meanwhile, the travel time difference of the same depth interval is gradually reduced due to the continuous increase of the speed of the underground medium. Therefore, the depth spacing of the travel time table at shallow layers needs to be small enough to ensure imaging accuracy; but high precision imaging can be obtained at deep layers even with large step sizes. Therefore, the embodiment of the invention adopts the step length-variable travel time table, avoids redundant calculation caused by adopting the same small grid, and can obtain an imaging result with similar precision.

Performing Kirchhoff integration method prestack depth migration twice on three-dimensional seismic data of a certain exploration area, and respectively setting the travel time table depth interval as a fixed step length (20m) and a variable step length (the depth interval is gradually increased from 20m to 120m along with the increase of the depth) in the application to obtain a common reflection point gather as shown in fig. 3(a) -3 (b) and a common offset gather as shown in fig. 4(a) -4 (b). Fig. 3(a) and 4(a) show the fixed-step integral prestack depth migration as a reference profile, fig. 3(b) and 4(b) show the variable-step migration results, and the two are subtracted from each other to obtain a residual profile, as shown in fig. 3(c) and 4 (c).

As can be seen from the offset results and the residual profiles shown in FIGS. 3(c) and 4(c), the step-variable travel time table can obtain imaging results with precision equivalent to that of a fixed step, and the error is small and almost negligible.

Compare the computational efficiency of two shifts:

in this test case, the throughput of the variable step travel time table is reduced by 73% compared to the fixed step travel time table. The offset time of the fixed step travel table is 8896s, the offset time of the variable step offset is 4302s, and the offset efficiency is improved by 107%.

When the integral method of the variable-step travel time table is used for prestack depth migration, the number of stored travel time table sample points is small, the I/O pressure and migration calculation amount of the travel time table in the migration process are reduced, and the error of the sample point amplitude value obtained by interpolation is small, so that the migration precision is guaranteed, and the migration efficiency can be greatly improved.

The implementation of the above specific application is only an example, and the rest of the embodiments are not described in detail.

Based on the same inventive concept, embodiments of the present invention further provide an integral method prestack depth migration apparatus, and since the principle of the problem solved by the integral method prestack depth migration apparatus is similar to that of the integral method prestack depth migration method, the implementation of the apparatus may refer to the implementation of the integral method prestack depth migration method, repeated details are not repeated, and the specific structure is shown in fig. 5:

a data acquisition module 501, configured to acquire seismic data of a target area;

a step length determining module 502, configured to determine a travel time table step length corresponding to each formation depth according to the seismic data and the calculation parameters;

the sampling module 503 is configured to sample in the depth direction of the stratum according to the step length of the travel time table corresponding to each depth of the stratum to obtain the travel time table of the target area;

and the data migration module 504 is configured to perform linear interpolation on the sampling points according to the travel time table step length and the travel time table corresponding to each stratum depth, and perform integral processing to obtain a migration result of the seismic data of the target area.

In a specific embodiment, calculating the parameters includes at least: travel time table initial step size, formation depth and effective frequency band. The step length determining module 502 is specifically configured to:

obtaining the maximum effective frequency of the seismic data of the target area according to the effective frequency band;

determining the minimum value of the velocity model at each stratum depth according to the velocity model and the stratum depth;

determining an adjustment parameter at each formation depth according to the formation depth;

and determining the travel time table step length corresponding to each stratum depth according to the maximum effective frequency of the seismic data of the target area, the minimum value of the velocity model at each stratum depth, the depth interval of the imaging points, the adjusting parameters at each stratum depth and the initial step length of the travel time table.

In specific implementation, according to the following formula, determining the travel time table step length corresponding to each stratum depth according to the maximum effective frequency of the seismic data of the target area, the minimum value of the velocity model at each stratum depth, the imaging point depth interval, the adjusting parameter at each stratum depth and the travel time table initial step length:

wherein L represents the formation depth; TMINC [ L ]]Representing the step length of a travel time table corresponding to the stratum depth L; INT represents an integer function; v. of_min[L]Represents the minimum of the velocity model at the formation depth L; para [ L]Represents a conditioning parameter at the formation depth L; f. of_maxRepresenting a maximum effective frequency of seismic data for the target area; dz represents the depth interval of the imaging points and is a preset value; TMINC0 represents the initial step size of the travel time table, which is a preset value.

Determining the adjustment parameter at each formation depth according to the following formula:

wherein para [ L ] represents the conditioning parameter at the formation depth L; l represents the formation depth.

In a specific embodiment, the data shifting module 504 is specifically configured to: when linear interpolation is carried out on the sampling points, the number of interpolation sampling points between any two sampling points is TMINC [ L ] -2.

The embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and when the processor executes the computer program, the integration prestack depth migration method is implemented.

An embodiment of the present invention further provides a computer-readable storage medium storing a computer program for executing the integration prestack depth migration method.

In summary, the method and the device for prestack depth migration by the integration method provided by the embodiments of the present invention have the following advantages:

acquiring seismic data of a target area; determining the travel time table step length corresponding to each stratum depth according to the seismic data and the calculation parameters; sampling in the depth direction of the stratum according to the step length of the travel time table corresponding to each stratum depth to obtain the travel time table of the target area; according to the travel time table step length and the travel time table corresponding to each stratum depth, carrying out linear interpolation on the sampling points, and carrying out integral processing to obtain the migration result of the seismic data of the target area; compared with the step length of a small and fixed travel time table in the prior art, the step length of the travel time table corresponding to each stratum depth is determined, namely the step length of the travel time table changes along with the change of the stratum depth, different step lengths are selected at different stratum depths, namely grid depth intervals of the travel time table are different at different stratum depths, so that the calculation amount of offset imaging and the I/O throughput of a hard disk for executing calculation are reduced on the premise of ensuring the imaging precision, and the offset efficiency is improved.

As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, apparatus, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.

The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

The above description is only a preferred embodiment of the present invention, and is not intended to limit the present invention, and various modifications and changes may be made to the embodiment of the present invention by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.