Comprehensive interpretation method and system for low-measurement data based on tomographic inversion

Abstract

The invention belongs to the field of seismic exploration of petroleum and natural gas, and relates to a comprehensive interpretation method and system for low-measurement data based on tomographic inversion. According to the method, a micro-logging routine interpretation result and a micro-logging tomography inversion result are obtained through the micro-logging first arrival time. And obtaining a small-refraction conventional interpretation result and a small-refraction tomography inversion result through the small-refraction first arrival time. And obtaining a comprehensive interpretation result by taking a micro-logging tomography inversion result and a small refraction tomography inversion result from a shallow layer and taking a micro-logging routine interpretation result and a small refraction routine interpretation result from a deep layer. According to the method, the advantages of high shallow layer precision of tomography inversion and stability of a deep layer conventional interpretation method are integrated, personal errors are eliminated, the detection capability of small refraction data is fully excavated, the precision of comprehensive interpretation of low measurement data is improved, the precision and stability of the low measurement data are improved, and through constraint tomography inversion used for multiple work areas, the accuracy and the stability of the low measurement data are improved. Therefore, the extremely shallow layer better conforms to the actual geological condition, and the imaging quality is improved.

Description

Technical Field

[0001] The invention belongs to the field of oil and natural gas seismic exploration, and relates to a method, system, equipment and medium for comprehensive interpretation of low-measurement data based on tomographic inversion. Background Art

[0002] Small refraction and micro-logging observations are two technical means of investigating the distribution of near-surface velocity in onshore oil and gas exploration. The survey results are mainly used for field well depth design. In recent years, with the development of first-arrival wave tomographic inversion technology, the extremely shallow near-surface velocity information obtained by small refraction and micro-logging is used to make up for the shortcomings of large-scale first-arrival tomographic inversion in extremely shallow layer modeling. By improving the accuracy of the near-surface model, the purpose of improving the integrated modeling effect of static correction and pre-stack depth migration is achieved. Therefore, higher requirements are put forward for the interpretation of small refraction and micro-logging.

[0003] At present, the interpretation results of micro-well logging and small refraction generally only give information on three layers: low-velocity layer, deceleration layer and high-velocity layer. When the low-velocity layer is thin, only two layers of information are given, namely low-velocity layer and high-velocity layer. This interpretation result basically meets the requirements of well depth design. However, it is used to constrain the first-arrival tomographic inversion of the cannon, and its accuracy is low. First, this "layered" interpretation result is an approximation of the near-surface velocity distribution, because most of the near-surface is Quaternary sediments, which are severely weathered and have a continuous increase in velocity rather than a layered distribution. Secondly, the interpretation methods of both data rely on human-computer interaction to determine the apparent velocity, which is subject to human errors. For small refraction, conventional interpretation uses the delayed time method. In theory, the thickness of the last layer cannot be obtained, and the detection depth of small refraction data cannot be fully explored. Summary of the invention

[0004] In order to overcome the shortcomings of the above-mentioned prior art, the purpose of the present invention is to provide a comprehensive interpretation method, system, equipment and medium for low-measurement data based on tomographic inversion. The present invention combines the advantages of high shallow-layer accuracy of tomographic inversion and stable deep-layer conventional interpretation methods, eliminates human errors, can obtain the maximum detection depth of small refractions, fully tap the detection capability of small refraction data, and improve the accuracy of comprehensive interpretation of low-measurement data.

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

[0006] In a first aspect, the present invention discloses a method for comprehensive interpretation of low-measurement data based on tomographic inversion, comprising the following steps:

[0007] Obtain the micro-logging first arrival time and well depth, and obtain the average zero-bias first arrival time based on the micro-logging first arrival time;

[0008] According to the average zero deviation first arrival time and well depth, the thickness and apparent velocity of the corresponding formation are determined to obtain the conventional interpretation results of micro-well logging;

[0009] Perform tomographic inversion on the micro-well logging first arrival time to obtain the micro-well logging tomographic inversion result;

[0010] Obtain the first arrival time of small refraction, interpret the first arrival time of small refraction according to the delayed time method, and obtain the conventional interpretation result of small refraction;

[0011] Perform first arrival wave tomographic inversion on the first arrival time of small refraction to obtain the tomographic inversion result of small refraction;

[0012] The micro-logging tomography inversion results and small refraction tomography inversion results are taken for the shallow layer, and the micro-logging conventional interpretation results and small refraction conventional interpretation results are taken for the deep layer to obtain the comprehensive interpretation results.

[0013] In a second aspect, the present invention discloses a low-measurement data comprehensive interpretation system based on tomographic inversion, comprising:

[0014] Data acquisition module: used to obtain the micro-logging first arrival time and well depth, and obtain the average zero deviation first arrival time based on the micro-logging first arrival time;

[0015] Micro-logging conventional interpretation module: used to determine the thickness and apparent velocity of the corresponding formation according to the average zero deviation first arrival time and well depth, and obtain the conventional interpretation results of micro-logging;

[0016] Micro-well logging tomographic inversion module: used to perform tomographic inversion on the micro-well logging first arrival time and obtain micro-well logging tomographic inversion results;

[0017] Small refraction conventional interpretation module: used to obtain the first arrival time of small refraction, interpret the first arrival time of small refraction according to the delayed time method, and obtain the conventional interpretation result of small refraction;

[0018] Small refraction tomographic inversion module: used to perform first arrival wave tomographic inversion on the first arrival time of small refraction to obtain the small refraction tomographic inversion result;

[0019] Comprehensive interpretation module: used to obtain micro-logging tomography inversion results and small refraction tomography inversion results in shallow layers, and micro-logging conventional interpretation results and small refraction conventional interpretation results in deep layers, to obtain comprehensive interpretation results.

[0020] In a third aspect, the present invention provides an electronic device, comprising: a processor; a memory for storing computer program instructions; and steps for implementing a comprehensive interpretation method for low-measurement data based on tomographic inversion when executing the computer program.

[0021] In a fourth aspect, the present invention provides a storage medium storing computer program instructions. When the computer program instructions are loaded and executed by a processor, the processor executes a comprehensive interpretation method for low-measurement data based on tomographic inversion.

[0022] Compared with the prior art, the present invention has the following beneficial effects:

[0023] The method of the present invention obtains the micro-well logging first arrival time and the well depth, and obtains the average zero-bias first arrival time according to the micro-well logging first arrival time. The thickness and apparent velocity of the corresponding formation are determined according to the average zero-bias first arrival time and the well depth, and the conventional interpretation result of the micro-well logging is obtained. The micro-well logging first arrival time is subjected to tomographic inversion to obtain the micro-well logging tomographic inversion result. The tomographic inversion is a data-driven algorithm, and there is no human error in the conventional interpretation method, which is conducive to improving the accuracy of the comprehensive interpretation of low-measurement data. The small refraction first arrival time is obtained, and the small refraction first arrival time is interpreted according to the delayed time method to obtain the conventional interpretation result of the small refraction. The small refraction first arrival time is subjected to first arrival wave tomographic inversion to obtain the small refraction tomographic inversion result, thereby fully exploring the detection capability of the small refraction data. The micro-well logging tomographic inversion result and the small refraction tomographic inversion result are obtained in the shallow layer, and the micro-well logging conventional interpretation result and the small refraction conventional interpretation result are obtained in the deep layer to obtain the comprehensive interpretation result. The present invention combines the advantages of high shallow-layer precision of tomographic inversion and stability of conventional deep-layer interpretation methods, eliminates human errors, can obtain the maximum detection depth of small refraction data, fully taps the detection capability of small refraction data, improves the precision of comprehensive interpretation of low-refraction data, improves the precision and stability of low-refraction data, and through constrained tomographic inversion used in multiple work areas, makes the extremely shallow layer more in line with the actual geological conditions and improves the imaging quality.

[0024] The system of the present invention includes a data acquisition module, a micro-well logging conventional interpretation module, a micro-well logging tomographic inversion module, a small refraction conventional interpretation module, a small refraction tomographic inversion module and a comprehensive interpretation module. The data acquisition module is used to obtain the micro-well logging first arrival time and well depth, and obtain the average zero deviation first arrival time according to the micro-well logging first arrival time. The micro-well logging conventional interpretation module is used to determine the thickness and apparent velocity of the corresponding formation according to the average zero deviation first arrival time and well depth, and obtain the micro-well logging conventional interpretation result. The micro-well logging tomographic inversion module is used to perform tomographic inversion on the micro-well logging first arrival time to obtain the micro-well logging tomographic inversion result. The small refraction conventional interpretation module is used to obtain the small refraction first arrival time, and interpret the small refraction first arrival time according to the delayed time method to obtain the small refraction conventional interpretation result. The small refraction tomographic inversion module is used to perform first arrival wave tomographic inversion on the small refraction first arrival time to obtain the small refraction tomographic inversion result. The comprehensive interpretation module is used to obtain the micro-logging tomography inversion results and small refraction tomography inversion results in the shallow layer, and the micro-logging conventional interpretation results and small refraction conventional interpretation results in the deep layer to obtain the comprehensive interpretation results. The modules cooperate with each other, combining the advantages of high shallow layer accuracy of tomography inversion and stable deep layer conventional interpretation methods, eliminating human errors, and being able to obtain the maximum detection depth of small refraction, fully tap the detection capability of small refraction data, and improve the accuracy of comprehensive interpretation of low-level measurement data. BRIEF DESCRIPTION OF THE DRAWINGS

[0025] Figure 1is a flow chart of the method of the present invention;

[0026] Figure 2 This is a schematic diagram of the micro-well logging wellhead time correction of the present invention;

[0027] Figure 3 This is an example diagram of conventional interpretation of a micro-logging data of the present invention;

[0028] Figure 4 This is an example diagram of the micro-well logging tomography inversion result of the present invention;

[0029] Figure 5 It is an example diagram of the conventional interpretation result of micro-well logging and the tomography inversion result of the present invention;

[0030] Figure 6 It is an example diagram of conventional interpretation results and comprehensive interpretation results of micro-well logging of the present invention;

[0031] Figure 7 This is an example diagram of the small refraction conventional explanation of the present invention;

[0032] Figure 8 for Figure 7 Example of small and medium refraction tomography inversion results;

[0033] Fig. 9 for Figure 8 Small and medium refraction tomography inversion ray density distribution map;

[0034] Fig.10 for Figure 7 Example diagrams of conventional interpretation results of small and medium refraction and tomographic inversion results;

[0035] Fig.11 for Figure 7 Example diagrams of conventional interpretation results and comprehensive interpretation results of small and medium refractions;

[0036] Fig.12 It is a system structure diagram of the present invention. DETAILED DESCRIPTION

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

[0038] It should be noted that the terms "first", "second", etc. in the specification and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units that are clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0039] The present invention is further described in detail below in conjunction with the accompanying drawings:

[0040] See also Figure 1 The present invention discloses a method for comprehensive interpretation of low-measurement data based on tomographic inversion, comprising the following steps:

[0041] S1. Obtain the micro-logging first arrival time and well depth, and obtain the average zero deviation first arrival time according to the micro-logging first arrival time;

[0042] S2. Determine the thickness and apparent velocity of the corresponding formation according to the average zero deviation first arrival time and the well depth, and obtain the conventional interpretation results of micro-well logging;

[0043] S3. Perform tomographic inversion on the micro-well logging first arrival time to obtain a micro-well logging tomographic inversion result;

[0044] S4. Obtain the first arrival time of the small refraction, interpret the first arrival time of the small refraction according to the delayed time method, and obtain the conventional interpretation result of the small refraction;

[0045] S5. performing first arrival wave tomographic inversion on the first arrival time of the small refraction to obtain a small refraction tomographic inversion result;

[0046] S6. Take the micro-logging tomography inversion results and small refraction tomography inversion results for the shallow layer, and the micro-logging conventional interpretation results and small refraction conventional interpretation results for the deep layer to obtain a comprehensive interpretation result.

[0047] See also Figure 1In another feasible embodiment of the present invention, the following is adaptively modified according to the situation. It includes the following steps: obtaining the micro-logging first arrival time and the well depth, and obtaining the average zero-bias first arrival time according to the micro-logging first arrival time. Determine the thickness and apparent velocity of the corresponding formation according to the average zero-bias first arrival time and the well depth, and obtain the conventional interpretation result of the micro-logging. Perform tomographic inversion on the micro-logging first arrival time to obtain the micro-logging tomographic inversion result. Tomographic inversion is a data-driven algorithm. There is no human error in the conventional interpretation method, which is conducive to improving the accuracy of the comprehensive interpretation of low-measurement data. Obtain the small refraction first arrival time, interpret the small refraction first arrival time according to the delayed time method, and obtain the small refraction conventional interpretation result. Perform the first arrival wave tomographic inversion on the small refraction first arrival time to obtain the small refraction tomographic inversion result. For small refraction data, the maximum detection depth can be obtained based on the ray density, and the detection capability of small refraction data can be fully exploited. Take the micro-logging tomographic inversion results and the small refraction tomographic inversion results for the shallow layer, and take the micro-logging conventional interpretation results and the small refraction conventional interpretation results for the deep layer to obtain the comprehensive interpretation result. The present invention performs tomographic inversion on micro-well logging and small refraction data respectively. First, a "grid" modeling is adopted for the near-surface model, which has higher accuracy than the conventional "layered" interpretation result. The present invention combines the advantages of high shallow layer accuracy of tomographic inversion and stability of deep layer conventional interpretation methods, eliminates human errors, can obtain the maximum detection depth of small refraction, fully taps the detection capability of small refraction data, improves the accuracy of comprehensive interpretation of low-reflection data, improves the accuracy and stability of low-reflection data, and makes the extremely shallow layer more in line with the actual geological conditions through constrained tomographic inversion used in multiple work areas, improving the imaging quality.

[0048] Embodiment 1:

[0049] See also Figure 1 This embodiment discloses a method for comprehensive interpretation of low-measurement data based on tomographic inversion, comprising the following steps:

[0050] S1. Obtain the micro-logging first arrival time and well depth, and obtain the average zero deviation first arrival time according to the micro-logging first arrival time;

[0051] The details are as follows:

[0052] The micro-logging first arrival time is projected onto the zero deviation according to the trigonometric function to obtain the zero deviation first arrival time. The average of all zero deviation first arrival times at the same depth is taken as the average zero deviation first arrival time at that depth.

[0053] S2. Determine the thickness and apparent velocity of the corresponding formation based on the average zero deviation first arrival time and well depth, and obtain the conventional interpretation results of micro-well logging, as follows:

[0054] The average zero deviation first arrival time and well depth are displayed as a distance-time curve graph. The slope of the distance-time curve graph is the apparent velocity at the corresponding depth. The thickness and apparent velocity of each layer are interactively determined. According to the apparent velocity trend, the velocity and depth distribution range of the low-velocity layer, the decreasing-velocity layer and the high-velocity layer are interactively determined to obtain the conventional interpretation results of micro-well logging.

[0055] S3. Perform tomographic inversion on the micro-well logging first arrival time to obtain a micro-well logging tomographic inversion result;

[0056] The process of obtaining micro-well logging tomography inversion results is as follows:

[0057] An initial velocity model is established according to the micro-logging first arrival time, and tomographic inversion is performed to obtain the inverted velocity model. The velocity at the wellhead is extracted to obtain the continuously distributed micro-logging tomographic inversion results.

[0058] S4. Obtain the first arrival time of the small refraction, interpret the first arrival time of the small refraction according to the delayed time method, and obtain the conventional interpretation result of the small refraction;

[0059] The process of obtaining the conventional interpretation results of small refraction is as follows:

[0060] The left and right branches of the first arrival time of the small refraction are displayed as time-distance scatter plots respectively. The slope of the scatter plot is the apparent velocity. Based on the delayed time method, the stratification and inflection point are determined interactively. The left and right interpretation results are averaged to obtain the conventional interpretation results of the small refraction, that is, the velocity and thickness of the low-speed layer and the deceleration layer, as well as the velocity of the high-speed layer. When interpreted in two layers, the low-speed layer and the deceleration layer are merged into one layer.

[0061] S5. performing first arrival wave tomographic inversion on the first arrival time of the small refraction to obtain a small refraction tomographic inversion result;

[0062] The process of obtaining the small refraction tomography inversion results is as follows:

[0063] Each small refraction observation is regarded as a two-dimensional survey line, and tomographic inversion is performed to obtain the velocity model and ray density within the range of the arrangement. The velocity curve at the midpoint of the arrangement is extracted as the tomographic inversion result at the small refraction observation point, and the detection depth is taken as the bottom limit of the ray density.

[0064] The conventional interpretation results of small refraction are the thickness and apparent velocity of each layer. The last layer has only velocity but no thickness.

[0065] S6. Take the micro-logging tomography inversion results and small refraction tomography inversion results for the shallow layer, and the micro-logging conventional interpretation results and small refraction conventional interpretation results for the deep layer to obtain a comprehensive interpretation result.

[0066] Embodiment 2:

[0067] See also Figure 1This embodiment discloses a method for comprehensive interpretation of low-measurement data based on tomographic inversion, comprising the following steps:

[0068] Step 1: Obtain the average zero deviation first arrival time, as follows:

[0069] Micro-logging observations usually involve burying several detectors or excitation shot points a few meters away from the wellhead. Conventional interpretation results require zero-bias time, that is, the detector points or shot points that deviate from the wellhead position need to be corrected according to trigonometric functions, and the micro-logging first arrival time is multiplied by the cosine of the angle with the well wall. Since the same depth often corresponds to several detector points or shot points, the zero-bias first arrival times of these traces are averaged to obtain the average zero-bias first arrival time at each depth. Figure 2 As shown in the figure, it is a schematic diagram of the wellhead time correction for non-zero deviation of micro-well logging. The dotted box represents the micro-well logging, the five-pointed star represents the location of a certain earthquake source, the well depth is h, and the solid triangle represents the actual position of a certain detection point, which has a certain deviation D from the wellhead. The actual first arrival propagation path is shown as the solid line in the figure. The wellhead correction is equivalent to moving the detection point to the position of the dotted triangle, and the ray propagation path becomes as shown in the dotted line in the figure.

[0070] Step 2: Conventional interpretation of micro-well logging, as follows:

[0071] After the trigonometric function correction, it is equivalent to that the excitation point and the receiving point are distributed along the micro-logging well wall and wellhead, and the lateral offset between the excitation point and the receiving point is zero. After the average zero offset first arrival time is displayed as a distance-time scatter plot, the slope is the apparent velocity. According to the apparent velocity trend, the processing personnel interactively determine the velocity and depth distribution range of the low-velocity layer, the velocity reduction layer and the high-velocity layer, and obtain the conventional interpretation results of the micro-logging. Figure 3 The figure shows an example of the zero deviation first arrival time and conventional interpretation of a micro-well logging, where the horizontal axis represents time, the vertical axis represents depth, the asterisks represent the propagation time and depth of each trace, and the broken lines in the figure represent two horizons of conventional interpretation.

[0072] Step 3: Perform tomographic inversion on micro-well logging data, as follows:

[0073] The initial velocity model is established based on the micro-logging first arrival time or the average zero-bias first arrival time after trigonometric function projection correction, and tomographic inversion is performed to obtain the inverted velocity model, such as Figure 4 As shown in the figure, the velocity at the wellhead is extracted to obtain the continuously distributed micro-logging tomography inversion results, such as Figure 5 As shown. Among them, Figure 4 It represents the velocity model obtained by tomographic inversion of the micro-well logging. The horizontal direction is distance and the vertical direction is depth. The velocity information at the midpoint is extracted to obtain the tomographic inversion result. Figure 5 Yes Figure 3The superposition of the conventional interpretation results of medium and micro-logging and the tomographic inversion results of micro-logging shows that the tomographic inversion results of micro-logging are disturbed at deep depths and the results are unstable.

[0074] Step 4: The general explanation of the method for delaying the first arrival time of the small refraction is as follows:

[0075] The left and right branches of the small refraction first arrival time are displayed as time-distance scatter plots, and their slopes are the apparent velocities. Based on the delayed time method, the layers and inflection points are determined interactively. The interpretation results of the left and right branches are averaged to obtain the conventional interpretation results, i.e., the velocity and thickness of the low-velocity layer and the deceleration layer, as well as the velocity of the high-velocity layer. When the interpretation is divided into two layers, the low-velocity layer and the deceleration layer are merged into one layer. Figure 7 The figure shows an example of conventional interpretation of small refraction data. The horizontal axis is distance and the vertical axis is time. Two layers are obtained along the left and right branches. The velocity and thickness of the low-speed layer are 910 m / s and 3.3 m respectively, and the velocity of the high-speed layer is 2049 m / s. Note that the conventional interpretation method cannot obtain the thickness of the high-speed layer.

[0076] Step 5: Perform tomographic inversion on the first arrival time of small refraction, as follows:

[0077] The first arrival time of each small refraction observation point is taken as a two-dimensional survey line, and the tomographic inversion is performed to obtain the velocity model and ray density within the range of the arrangement slice, respectively as follows: Figure 8 and Fig. 9 As shown, the velocity curve at the midpoint of the arrangement is extracted as the tomographic inversion result at the small refraction observation point, and the detection depth is taken as the bottom limit of the ray density. Figure 8 Yes Figure 7 Velocity model after tomographic inversion of the first arrival time of small and medium refractions. The horizontal axis represents distance and the vertical axis represents depth. The velocity at the midpoint of the black dotted arrow is extracted to obtain the tomographic inversion result. Fig. 9 Yes Figure 7 The ray density for tomographic inversion of small and medium refractions has distance as the horizontal coordinate and depth as the vertical coordinate. The maximum depth corresponding to its envelope is the detection depth of the small refraction.

[0078] Step 6: Comprehensive interpretation of small refraction results and micro-well logging results, as follows:

[0079] After the micro-logging and small refraction are processed by conventional interpretation and tomographic inversion respectively, each micro-logging or small refraction point has the conventional interpretation results and tomographic inversion results respectively. The common characteristics of them are that the tomographic inversion results are high in shallow layer precision. For micro-logging, due to the low number of deep layer coverage and small data redundancy, the inversion results are unstable and the reliability is low; for small refraction, due to the single incident / emergence angle of deep layer rays, the inversion results are unstable and the reliability is low. Therefore, the present invention takes the tomographic inversion results in the shallow layer and the conventional interpretation results in the deep layer, which combines the advantages of high shallow layer precision of tomographic inversion and stable deep layer of conventional interpretation results, as shown in the following examples. Figure 6 Micro-logging comprehensive interpretation, Fig.10 for Figure 7 Example diagram of conventional interpretation results of small and medium refraction and tomographic inversion results. Fig.11 An example diagram of a comprehensive explanation of small refraction. Figure 6 This is an example of comprehensive interpretation of micro-well logging, in which the tomographic inversion results are taken for the shallow layer and the conventional interpretation results are taken for the deep layer to form the comprehensive interpretation results.

[0080] The present invention performs tomographic inversion on micro-logging and small refraction data respectively. First, "grid" modeling is adopted for the near-surface model, which has higher accuracy than the conventional "layered" interpretation result. Tomographic inversion is a data-driven algorithm, and there is no human error in conventional interpretation methods. For micro-logging data, the error introduced by projection to zero bias is eliminated. For small refraction data, the maximum detection depth can be obtained based on ray density. In view of the low redundancy of deep observation data of micro-logging data and the single deep ray angle of small refraction data, and the unstable inversion results, the velocity information of conventional interpretation is adopted. The new invention combines the advantages of high shallow layer accuracy of tomographic inversion and stable deep layer conventional interpretation methods, improves the accuracy and stability of low-measurement data, and makes the extremely shallow layer more in line with the actual geological conditions through constrained tomographic inversion used in multiple work areas, improves the imaging quality, and compared with similar technologies, this intellectual property is at the leading level.

[0081] Based on the above method, the present invention discloses a low-measurement data comprehensive interpretation system based on tomographic inversion, comprising:

[0082] Data acquisition module: used to obtain the micro-logging first arrival time and well depth, and obtain the average zero deviation first arrival time based on the micro-logging first arrival time;

[0083] Micro-logging conventional interpretation module: used to determine the thickness and apparent velocity of the corresponding formation according to the average zero deviation first arrival time and well depth, and obtain the conventional interpretation results of micro-logging;

[0084] Micro-well logging tomographic inversion module: used to perform tomographic inversion on the micro-well logging first arrival time and obtain micro-well logging tomographic inversion results;

[0085] Small refraction conventional interpretation module: used to obtain the first arrival time of small refraction, interpret the first arrival time of small refraction according to the delayed time method, and obtain the conventional interpretation result of small refraction;

[0086] Small refraction tomographic inversion module: used to perform first arrival wave tomographic inversion on the first arrival time of small refraction to obtain the small refraction tomographic inversion result;

[0087] Comprehensive interpretation module: used to obtain micro-logging tomography inversion results and small refraction tomography inversion results in shallow layers, and micro-logging conventional interpretation results and small refraction conventional interpretation results in deep layers, to obtain comprehensive interpretation results.

[0088] See also Fig.12 In another feasible embodiment of the present invention, the following is adaptively modified according to the situation. It includes a data acquisition module, a micro-well logging conventional interpretation module, a micro-well logging tomographic inversion module, a small refraction conventional interpretation module, a small refraction tomographic inversion module and a comprehensive interpretation module. The data acquisition module is used to obtain the micro-well logging first arrival time and well depth, and obtain the average zero deviation first arrival time according to the micro-well logging first arrival time. Micro-well logging conventional interpretation module: used to determine the thickness and apparent velocity of the corresponding formation according to the average zero deviation first arrival time and well depth, and obtain the micro-well logging conventional interpretation result. The micro-well logging tomographic inversion module is used to perform tomographic inversion on the micro-well logging first arrival time to obtain the micro-well logging tomographic inversion result. The small refraction conventional interpretation module is used to obtain the small refraction first arrival time, and interpret the small refraction first arrival time according to the delayed time method to obtain the small refraction conventional interpretation result. The small refraction tomographic inversion module is used to perform first arrival wave tomographic inversion on the small refraction first arrival time to obtain the small refraction tomographic inversion result. The comprehensive interpretation module is used to obtain the micro-logging tomography inversion results and small refraction tomography inversion results in the shallow layer, and the micro-logging conventional interpretation results and small refraction conventional interpretation results in the deep layer to obtain the comprehensive interpretation results. The modules cooperate with each other, combining the advantages of high shallow layer accuracy of tomography inversion and stable deep layer conventional interpretation methods, eliminating human errors, and being able to obtain the maximum detection depth of small refraction, fully tap the detection capability of small refraction data, and improve the accuracy of comprehensive interpretation of low-level measurement data.

[0089] An electronic device comprises: a processor; a memory for storing computer program instructions; and steps for implementing a comprehensive interpretation method of low-measurement data based on tomographic inversion when executing the computer program.

[0090] A storage medium stores computer program instructions. When the computer program instructions are loaded and run by a processor, the processor executes a comprehensive interpretation method for low-measurement data based on tomographic inversion.

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

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

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

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

[0095] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the above embodiments, ordinary technicians in the relevant field should understand that the specific implementation methods of the present invention can still be modified or replaced by equivalents. Any modification or equivalent replacement that does not depart from the spirit and scope of the present invention should be covered within the scope of protection of the claims of the present invention.

[0096] The above contents are only for explaining the technical idea of ​​the present invention and cannot be used to limit the protection scope of the present invention. Any changes made on the basis of the technical solution in accordance with the technical idea proposed by the present invention shall fall within the protection scope of the claims of the present invention.

Claims

1. A comprehensive interpretation method for low-measurement data based on tomographic inversion, characterized in that: The following steps are involved: Obtain the micro-logging first arrival time and well depth, and obtain the average zero-bias first arrival time based on the micro-logging first arrival time; According to the average zero deviation first arrival time and well depth, the thickness and apparent velocity of the corresponding formation are determined to obtain the conventional interpretation results of micro-well logging; Perform tomographic inversion on the micro-well logging first arrival time to obtain the micro-well logging tomographic inversion result; Obtain the first arrival time of small refraction, interpret the first arrival time of small refraction according to the delayed time method, and obtain the conventional interpretation result of small refraction; Perform first arrival wave tomographic inversion on the first arrival time of small refraction to obtain the tomographic inversion result of small refraction; The micro-logging tomography inversion results and small refraction tomography inversion results are taken for the shallow layer, and the micro-logging conventional interpretation results and small refraction conventional interpretation results are taken for the deep layer to obtain the comprehensive interpretation results.

2. The method for comprehensive interpretation of low-measurement data based on tomographic inversion according to claim 1, characterized in that: The specific method of obtaining the average zero deviation first arrival time according to the micro-logging first arrival time and deviation distance is as follows: The micro-logging first arrival time is projected onto the zero deviation according to the trigonometric function to obtain the zero deviation first arrival time. The average of all zero deviation first arrival times at the same depth is taken as the average zero deviation first arrival time at that depth.

3. The method for comprehensive interpretation of low-measurement data based on tomographic inversion according to claim 1, characterized in that: The conventional interpretation results of micro-well logging are specifically as follows: The average zero deviation first arrival time and well depth are displayed as a distance-time curve graph. The slope of the distance-time curve graph is the apparent velocity at the corresponding depth. The thickness and apparent velocity of each layer are interactively determined. According to the apparent velocity trend, the velocity and depth distribution range of the low-velocity layer, the decreasing-velocity layer and the high-velocity layer are interactively determined to obtain the conventional interpretation results of micro-well logging.

4. The method for comprehensive interpretation of low-measurement data based on tomographic inversion according to claim 1, characterized in that: The process of obtaining the micro-logging tomography inversion results is as follows: An initial velocity model is established according to the micro-logging first arrival time, and tomographic inversion is performed to obtain the inverted velocity model. The velocity at the wellhead is extracted to obtain the continuously distributed micro-logging tomographic inversion results.

5. The method for comprehensive interpretation of low-measurement data based on tomographic inversion according to claim 1, characterized in that: The process of obtaining the conventional interpretation result of small refraction is as follows: The left and right branches of the first arrival time of the small refraction are displayed as time-distance scatter plots respectively. The slope of the scatter plot is the apparent velocity. Based on the delayed time method, the stratification and inflection point are determined interactively. The interpretation results of the left and right branches are averaged to obtain the conventional interpretation results of the small refraction.

6. The method for comprehensive interpretation of low-measurement data based on tomographic inversion according to claim 1, characterized in that: The process of obtaining the small refraction tomography inversion results is as follows: Each small refraction observation is regarded as a two-dimensional survey line, and tomographic inversion is performed to obtain the velocity model and ray density within the range of the arrangement. The velocity curve at the midpoint of the arrangement is extracted as the tomographic inversion result at the small refraction observation point, and the detection depth is taken as the bottom limit of the ray density.

7. A method for comprehensive interpretation of low-measurement data based on tomographic inversion as claimed in claim 6, characterized in that: The conventional interpretation results of small refraction are the thickness and apparent velocity of each layer, wherein the last layer has only velocity but no thickness.

8. A comprehensive interpretation system for low-level data based on tomographic inversion, characterized in that: include: Data acquisition module: used to obtain the micro-logging first arrival time and well depth, and obtain the average zero deviation first arrival time based on the micro-logging first arrival time; Micro-logging conventional interpretation module: used to determine the thickness and apparent velocity of the corresponding formation according to the average zero deviation first arrival time and well depth, and obtain the conventional interpretation results of micro-logging; Micro-well logging tomographic inversion module: used to perform tomographic inversion on the micro-well logging first arrival time and obtain micro-well logging tomographic inversion results; Small refraction conventional interpretation module: used to obtain the first arrival time of small refraction, interpret the first arrival time of small refraction according to the delayed time method, and obtain the conventional interpretation result of small refraction; Small refraction tomographic inversion module: used to perform first arrival wave tomographic inversion on the first arrival time of small refraction to obtain the small refraction tomographic inversion result; Comprehensive interpretation module: used to obtain micro-logging tomography inversion results and small refraction tomography inversion results in shallow layers, and micro-logging conventional interpretation results and small refraction conventional interpretation results in deep layers, to obtain comprehensive interpretation results.

9. An electronic device, comprising: A processor; a memory for storing computer program instructions; characterized in that it is used to implement the steps of the comprehensive interpretation method of low-measurement data based on tomographic inversion as described in any one of claims 1 to 7 when executing the computer program.

10. A storage medium storing computer program instructions, characterized in that: When the computer program instructions are loaded and executed by the processor, the processor executes the comprehensive interpretation method of low-measurement data based on tomographic inversion as described in any one of claims 1 to 7.

Images