A joint interactive interpretation method of P-wave and S-wave micro-logging data in surface structure investigation

Abstract

The present invention discloses a method for joint interactive interpretation of P-wave and S-wave micro-logging data in surface structure investigation. The method comprises four steps: establishing a joint interactive processing and interpretation framework for P-wave and S-wave micro-logging data, performing real-time joint interactive communication, and realizing interactive interpretation and selection and output of interactive interpretation schemes. The method improves the existing method in which P-wave and S-wave velocity stratification can only be interpreted separately, and establishes real-time joint interactive communication of P-wave and S-wave data information between sub-windows within the framework; adopts a human-computer interaction method to realize a human-computer interactive interpretation mode and an automatic interpretation mode completely controlled by interpretation parameters; realizes the optimization of interpretation schemes through interactive interpretation scheme selection, thereby improving the accuracy of P-wave and S-wave micro-logging data processing results.

Description

Technical Field

[0001] The invention belongs to the technical field of near-surface formation parameter survey, and relates to a data joint interactive interpretation method in surface structure survey, and specifically to a P- and S-wave micro-logging data joint interactive interpretation method in surface structure survey. Background Art

[0002] In industries such as oil and gas exploration, geological and mineral exploration, and civil engineering exploration, the combined P- and S-wave surface structure survey can more precisely characterize the characteristic parameters of the surface structure. It has always been a topic of continuous research by geophysicists engaged in surface structure surveys and P- and S-wave seismic data processing. However, further research and technological breakthroughs are still needed in the simultaneous observation and investigation of P- and S-wave surface structures.

[0003] The seismic P-wave surface structure survey characterizes the P-wave characteristics of the surface structure, which is a comprehensive response of the surface structure skeleton and the fluid contained therein, reflecting the P-wave characteristics of the surface structure of the low-velocity reduction zone and the high-velocity layer. The surface structure below the water table is rich in water, so its skeleton characteristics are difficult to be effectively characterized. The seismic S-wave surface structure survey can more accurately characterize the lithological characteristics and skeleton of the surface structure. The characterized surface structure is richer, the lithological characteristics are clearer and more reliable, and it will not be affected by the water content of the water table and the surface structure. It can be seen that the combination of P-wave and S-wave in the surface structure survey can more comprehensively and objectively characterize the surface structure, and the research is of great significance.

[0004] From the analysis of P-wave and S-wave micro-logging data in the existing surface structure survey, the variation law of P-wave and S-wave velocities near the surface has its own influencing factors, the velocity stratification interface correspondence is unclear, and it is impossible to give a unified conclusion in terms of geophysical characteristics and geological sedimentary characteristics. The existing P-wave and S-wave micro-logging data have no interactive integration ideas in the processing interpretation process, interpretation results, and application of results, which has a great impact on the interpretation of S-wave micro-logging data, the calculation of S-wave static correction and the application of surface S-wave data. At the same time, due to the multi-solution nature of the P-wave and S-wave micro-logging data velocity stratification interpretation, there is currently no pre-storage idea of ​​multiple interpretation results, and there is a lack of methods for optimizing interpretation results. Summary of the invention

[0005] The purpose of the present invention is to provide a method for the joint interactive interpretation of P-wave and S-wave micro-logging data in surface structure investigation, which achieves the purpose of improving the efficiency of the joint interactive interpretation of P-wave and S-wave logging data and the accuracy of the results in surface structure investigation by using the same-window and same-scale display technology and synchronously performing P-wave and S-wave joint manual and automatic interactive interpretation and other technical means.

[0006] To achieve the above purpose, the technical solution adopted by the present invention is as follows:

[0007] The combined interactive interpretation method of P-wave and S-wave micro-logging data in surface structure investigation includes the following steps:

[0008] S1. Establish a framework for joint interactive processing and interpretation of P-wave and S-wave micrologging data

[0009] Display the first arrival data of P-wave and S-wave micro-logging in the main window, and establish a joint interactive processing and interpretation framework of P-wave and S-wave micro-logging data, including manual interactive interpretation mode, automatic interactive interpretation mode and auxiliary function control module;

[0010] The P-wave and S-wave micro-logging first arrival data include P-wave first arrival data and S-wave first arrival data, and the two types of data are displayed in sub-windows under the main window with the same scale respectively;

[0011] Among them, the P-wave and S-wave micro-logging data processing window is divided into two sub-windows, left (P-wave first arrival) and right (S-wave first arrival), and the P-wave and S-wave first arrival data are loaded at the same time. Two interpretation modes, "manual interactive interpretation" and "automatic interactive interpretation", and an auxiliary function control module are set in the main window.

[0012] S2. Conduct real-time joint interactive communication

[0013] In the framework of joint interactive processing and interpretation of P-wave and S-wave micro-logging data, joint interactive interpretation is adopted for the P-wave and S-wave micro-logging first arrival data. By monitoring the information of P-wave and S-wave velocity layers in their respective sub-windows and calling the P-wave and S-wave interpretation process code, the intermediate results and final results are interpreted using structure and common variable pointers.

[0014] Use private data structure pointer variables to process private data in intermediate results and final results, dynamically allocate storage space and save private data;

[0015] The common data common body pointer variable is used to process the common data in the first arrival data of the P-wave and S-wave micro-logging, and the storage space is dynamically allocated to realize the real-time joint interactive communication of the velocity layered interpretation of the first arrival data of the P-wave and S-wave micro-logging.

[0016] S3. Implementing interactive explanation

[0017] The longitudinal wave first arrival data is processed by linear function fitting to obtain the longitudinal wave layered node results, which are displayed in a sub-window on one side;

[0018] The shear wave first arrival data is processed by nonlinear function fitting, and the shear wave layered node results are obtained by human-computer interaction and displayed in the sub-window on the other side;

[0019] Among them, in the process of transverse wave interpretation, by monitoring the corresponding events, calling the longitudinal wave sub-window interpretation process code, intermediate result data and final result data in the event process, using structure pointer variables to process the intermediate result data, and using common pointer variables to process the common data of the intermediate result data and the final result data, the processing results are obtained, the storage space is dynamically allocated and managed, and real-time joint interactive communication is realized;

[0020] Graphics are drawn using the VCL graphics component of C++Builder, and are drawn in the canvas of a graphic object (such as Form, Image, etc.) rather than directly in the object. The canvas is a property of the graphic object, and it is also an object itself. The program can draw directly on the control's canvas by calling the drawing method of the canvas object.

[0021] S4. Selection and output of interactive interpretation scheme

[0022] The selection criteria for interactive interpretation schemes include: consistency of water-free stratum interfaces in P- and S-wave stratification results, consistency of lithologic stratification interfaces, and consistency of interpretation results with environmental areas and uses;

[0023] The output of the interactive interpretation scheme is the interpretation result data and the joint interpretation result diagram of longitudinal and transverse waves.

[0024] Preferably, the interactive interpretation scheme selects a scheme with high consistency of the water-free stratum interface in the P- and S-wave stratification results, high consistency of the lithology stratification interface, high consistency of the interpretation results with other survey points in the environmental area, and high consistency of the interpretation results with the problem to be solved (i.e., purpose), specifically:

[0025] ① The interface of water-free strata is highly consistent, and the interface of lithologic stratification is highly consistent. To reduce the contradiction between the P- and S-wave interpretation results, the interface of water-free strata should be as consistent as possible, the top surface of the aquifer may not correspond, but the interface of lithologic stratification should be highly consistent. The S-wave velocity is more sensitive to compaction than the P-wave, and a nonlinear stratification method is used to improve the correspondence of the lithologic stratification interface.

[0026] ② The interpretation results are highly consistent with other survey points in the environmental area. The optimization principle is to take both vertical and horizontal directions into consideration, not the finer the stratification, the better. The plane interpolation accuracy is directly related to the coarseness of the stratification. For different surface geographical features and sedimentary environments, the corresponding relationship between the longitudinal and transverse wave velocity stratification of the near-surface structure is different. Not only should the interpretation parameters of individual survey points be optimized, but they should also be highly consistent with the interpretation schemes of various survey points within the same geographical features and sedimentary environment area.

[0027] ③ The interpretation results are highly consistent with the problem being solved (i.e. purpose). The interpretation results vary depending on their purpose. The emphasis of P- and S-wave velocity stratification is also different for different application requirements. For seismic exploration, the selection of the best lithology section for the design of the well depth for excitation is more focused on velocity and lithology, and a small S-wave velocity fitting time difference parameter should be selected; if the selection is for the extraction of high-precision static correction, the consistency of the interpretation accuracy of most measuring points should be emphasized, and the stratification fitting parameters should be regular and not change dramatically.

[0028] As a limitation of the present invention, the auxiliary function control module includes a parameter setting module, an automatic interactive refinement and comparison module, an interpretation scheme saving module and a pre-stored interpretation scheme retrieval module.

[0029] Wherein, the parameter setting module includes fixed parameters;

[0030] The function of the automatic interactive refinement comparison module is to apply fixed parameters to coordinate and control the velocity layering nodes of one wave to achieve the consistency of the longitudinal and transverse wave layering nodes;

[0031] The function of the interpretation scheme saving module is to save the interactive interpretation results in the form of independent files;

[0032] The function of the pre-stored interpretation scheme retrieval module is to call and view the interpretation scheme saving module to realize the interpretation scheme selection and output the interactive interpretation scheme.

[0033] As another limitation of the present invention, the processing results of the P-wave and S-wave micro-logging first arrival data are displayed in respective sub-windows through interpretation graphics.

[0034] As a further limitation of the present invention, the interpretation graphics are drawn by using the VCL graphics component of C++Builder to process the first arrival data of the longitudinal and shear wave micrologging.

[0035] As a third limitation of the present invention,

[0036] The linear function is H=V×T,

[0037] Among them, H represents depth (unit: meter); V represents velocity (unit: meter / second); T represents the corresponding time (unit: second).

[0038] The nonlinear function is Vs=a i +b i ×H,

[0039] Where Vs represents the shear wave layer velocity (unit: m / s); a i represents the intercept of the i-th segment fitting; b i represents the fitting shear wave velocity gradient of the ith segment; H represents the depth (unit: meter).

[0040] According to the above two stratification functions, the velocity stratification preprocessing of the longitudinal and shear wave micro-logging data is performed respectively to obtain the velocity stratification information of the longitudinal and shear wave micro-logging data, that is, the preprocessing result submits the stratification node information such as velocity and thickness of each velocity layer.

[0041] As a fourth limitation of the present invention,

[0042] The human-computer interaction method consists of manual interaction interpretation and automatic interaction interpretation;

[0043] The manual interactive interpretation is to manually move the hierarchical results of the sub-window on one side to control the corresponding update of the hierarchical results in the sub-window on the other side.

[0044] As a fifth limitation of the present invention, the automatic interactive interpretation is completed by inputting the fixed parameters after the manual interactive interpretation is completed, and the P-wave and S-wave micro-logging velocity layered interpretation result diagram is obtained.

[0045] As a further limitation, the fixed parameters include: a P-wave velocity fitting threshold value, a S-wave fitting time difference coefficient, an allowable range value of the P-wave and S-wave velocity stratification interface, a percentage parameter of the velocity difference between adjacent P-wave and S-wave layers, and a maximum distance parameter of the S-wave stratification node from the interface.

[0046] As a sixth limitation of the present invention, the method is implemented by software having modules for manual interactive interpretation, automatic interactive interpretation, detailed comparison, interactive interpretation control parameter setting and interpretation scheme saving.

[0047] As a further limitation, the software is PSWISS, which has applied for software copyright, copyright certificate number: Software Registration No. 7487061, software registration number: 2021SR0764435, and can be obtained and used by the public.

[0048] Due to the adoption of the above technical solution, the technical progress achieved by the present invention compared with the prior art is:

[0049] The method of the present invention can reasonably interpret the surface P-wave and S-wave micro-logging data by displaying the data at the same scale in a main window and performing P-wave and S-wave joint interactive interpretation simultaneously, thereby drawing an accurate velocity model of the surface P-wave and S-wave low-speed drop zone, and improving the accuracy of P-wave and S-wave surface investigation and the efficiency of data processing.

[0050] The present invention establishes a framework for the joint interactive processing and interpretation of P- and S-wave micro-logging data, which improves the existing method in which P- and S-wave velocity stratification can only be interpreted separately, establishes real-time joint interactive communication of P- and S-wave data information between sub-windows within the framework, improves the existing interpretation method in which P- and S-wave velocity stratification do not refer to each other, adopts a human-computer interaction method, realizes a human-computer interaction interpretation mode and an automatic interpretation mode completely controlled by interpretation parameters, realizes the optimization of interpretation schemes through interactive interpretation scheme selection, thereby improving the accuracy of P- and S-wave micro-logging data processing results.

[0051] The present invention can be applied to the field of near-surface formation parameter survey technology for petroleum seismic exploration. It can further meet the multiple requirements of high-precision surface survey, conventional seismic exploration and multi-component seismic exploration data processing and interpretation, provide high-precision longitudinal and shear wave surface static corrections for multi-wave seismic exploration data processing, and improve the effectiveness of oil and gas exploration.

[0052] The present invention is further described in detail below with reference to the accompanying drawings and specific embodiments. BRIEF DESCRIPTION OF THE DRAWINGS

[0053] Figure 1 This is a window diagram of the framework for the joint interactive processing and interpretation of P-wave and S-wave micro-logging data in Example 1;

[0054] Figure 2 It is a diagram for setting interactive interpretation parameters for automatic interactive interpretation in Example 1;

[0055] Figure 3 This is a schematic diagram of the division of surface geographical features of the Alaxin three-component block in Example 2;

[0056] Figure 4 This is the result diagram of the joint interpretation of the longitudinal and transverse waves in the Jiangwan area in Example 2, which is a 0.5ms fitting time difference;

[0057] Figure 5 This is the result diagram of the joint interpretation of longitudinal and transverse waves in Jiangwan District in Example 2 - fitting time difference 1.0ms;

[0058] Figure 6 This is the result diagram of the joint interpretation of the longitudinal and transverse waves in the Jiangwan area in Example 2 - fitting time difference 1.5ms;

[0059] Figure 7 This is the result diagram of the joint interpretation of the longitudinal and transverse waves in the Jiangwan area in Example 2, which is a 2.0ms fitting time difference;

[0060] Figure 8 This is the result diagram of the joint interpretation of the dryland area-fitting time difference 0.5 longitudinal and transverse waves in Example 2;

[0061] Fig. 9 This is the result diagram of the joint interpretation of the dryland area-fitting time difference 1.0 longitudinal and transverse waves in Example 2;

[0062] Fig.10 This is the result diagram of the joint interpretation of the dryland area-fitting time difference 1.5 longitudinal and transverse waves in Example 2;

[0063] Fig.11 This is the result diagram of the joint interpretation of the Gaogang area in Example 2 - fitting time difference 0.5 longitudinal and transverse waves;

[0064] Fig.12 This is the result diagram of the joint interpretation of the Gaogang area-fitting time difference 1.0 longitudinal and transverse waves in Example 2;

[0065] Fig.13 This is the result diagram of the joint interpretation of the fitting time difference 1.5 longitudinal and transverse waves in the Gaogang area of ​​Example 2;

[0066] Fig.14 This is the result diagram of the joint interpretation of the Gaogang area-fitting time difference 2.0 longitudinal and transverse waves in Example 2;

[0067] Fig.15 It is a plane diagram of the static correction of longitudinal wave detection points in the Alaxin work area in Example 2;

[0068] Fig.16 This is a plane diagram of the static correction values ​​of the shear wave detection points in the Alaxin work area in Example 2. DETAILED DESCRIPTION

[0069] Example 1 A method for joint interactive interpretation of P-wave and S-wave micro-logging data in surface structure investigation

[0070] The method of this specific embodiment is implemented by the software PSWISS having modules for manual interactive interpretation, automatic interactive interpretation, detailed comparison, interactive interpretation control parameter setting and interpretation scheme saving, and includes the following steps performed in sequence:

[0071] S1. Establish a framework for joint interactive processing and interpretation of P-wave and S-wave micrologging data

[0072] The P-wave first arrival data and S-wave first arrival data are displayed in the sub-windows under the main window at the same scale, and a P-wave and S-wave micro-logging data joint interactive processing and interpretation framework including manual interactive interpretation mode, automatic interactive interpretation mode and auxiliary function control module is established;

[0073] Among them, the auxiliary function control module includes a parameter setting module, an automatic interactive refinement comparison module, an interpretation scheme saving module, and a pre-stored interpretation scheme retrieval module. Figure 1 shown.

[0074] Among them, the parameter setting module includes fixed parameters; the fixed parameters include: P-wave velocity fitting threshold value, S-wave fitting time difference coefficient, P-wave and S-wave velocity stratification interface allowable range value, P-wave and S-wave adjacent layer velocity difference percentage parameter and S-wave stratification node maximum distance parameter from the interface.

[0075] The automatic interactive refinement comparison module controls the velocity stratification nodes of shear waves through the velocity stratification nodes of longitudinal waves, or controls the velocity stratification nodes of shear waves through the velocity stratification nodes of longitudinal waves. Finally, the corresponding relationship of the interactive interpretation stratification nodes is controlled by the corresponding allowable range error of the position of the velocity stratification interface of longitudinal and shear waves and the velocity difference percentage parameters between adjacent layers of longitudinal and shear waves, so as to better reflect the consistency of the stratification nodes of longitudinal and shear waves.

[0076] The interpretation scheme saving module saves the interpretation results of a set of interactive interpretation parameters into a separate file. This file contains the first arrival time of the longitudinal and transverse wave micro-logging data of the physical point and the corresponding velocity and thickness layer parameter results. Multiple interpretation schemes can be saved for easy selection.

[0077] The pre-stored interpretation scheme retrieval module is to call and view the pre-stored interpretation schemes for the convenience of optimization, and finally retain and output the best joint interactive interpretation results of P-wave and S-wave micro-logging data.

[0078] S2. Conduct real-time joint interactive communication

[0079] In the framework of joint interactive processing and interpretation of P-wave and S-wave micro-logging data, interactive interpretation is adopted for the first arrival data of P-wave and S-wave micro-logging. The information of each sub-module window is layered by monitoring the corresponding P-wave and S-wave velocity. In the monitoring process, the P-wave and S-wave interpretation process codes are called. The intermediate results and final results are interpreted by using structure and common variable pointers. The private data structure pointer variables are used to process the private data in the final result. The storage space is dynamically allocated and managed according to the needs. The private data is saved by dynamically allocating the storage space.

[0080] For the common data in the first arrival data of longitudinal and shear wave micro-logging, the common data common body pointer variable is used to dynamically allocate storage space, which can be called by the longitudinal and shear wave interpretation process. At the same time, the VCL graphic component of C++Builder is used to draw the corresponding interpretation graphics in the corresponding sub-window display area to realize the real-time joint interactive communication of the first arrival data of longitudinal and shear wave micro-logging; S3. Realize interactive interpretation

[0081] The P-wave first arrival data are processed by fitting the linear function H=V×T to obtain the P-wave stratification results, which are displayed in the subwindow;

[0082] Through the nonlinear function Vs = a i +b i ×H fitting processes the shear wave first arrival data, adopts a human-computer interaction method consisting of manual interactive interpretation and automatic interactive interpretation, obtains the shear wave stratification results, and displays them in the subwindow;

[0083] In the linear function, H represents the depth (unit: meter); V represents the velocity (unit: meter / second); and T represents the corresponding time (unit: second). In the nonlinear function, Vs represents the shear wave layer velocity; a i represents the intercept of the i-th segment fitting; b i represents the fitting shear wave velocity gradient of the i-th segment; H represents the depth;

[0084] According to the above two stratification functions, the velocity stratification preprocessing of the P-wave and S-wave micro-logging data is performed respectively. The obtained P-wave stratification results and S-wave stratification results are the velocity stratification information of the P-wave and S-wave micro-logging data, which are the stratification node information such as velocity and thickness of each velocity layer.

[0085] Manual interaction explanation is to manually move the hierarchical node information of the main window. By clicking and dragging the corresponding information to a new position and releasing it, the hierarchical node information in the sub-window can be controlled to be updated accordingly.

[0086] Automatic interactive interpretation is achieved by inputting fixed parameters after manual interactive interpretation is completed, such as Figure 2 shown.

[0087] S4. Selection and output of interactive interpretation scheme

[0088] The interactive interpretation scheme selects the scheme with high consistency of water-free stratum interface in the P- and S-wave stratification results, high consistency of lithology stratification interface, high consistency of interpretation results with other survey points in the environmental area, and high consistency of interpretation results with the problem to be solved (i.e. purpose), specifically:

[0089] ① The interface of water-free strata is highly consistent, and the interface of lithologic stratification is highly consistent. To reduce the contradiction between the P- and S-wave interpretation results, the interface of water-free strata should be as consistent as possible, the top surface of the aquifer may not correspond, but the interface of lithologic stratification should be highly consistent. Here, the S-wave velocity is more sensitive to compaction than the P-wave, and a nonlinear stratification method is used to improve the correspondence of the lithologic stratification interface.

[0090] ② The interpretation results are highly consistent with other survey points in the environmental area. The optimization principle is to take both vertical and horizontal directions into consideration, not the finer the stratification, the better. The plane interpolation accuracy is directly related to the coarseness of the stratification. For different surface geographical features and sedimentary environments, the corresponding relationship between the longitudinal and transverse wave velocity stratification of the near-surface structure is different. Not only should the interpretation parameters of individual survey points be optimized, but they should also be highly consistent with the interpretation schemes of various survey points within the same geographical features and sedimentary environment area.

[0091] ③ The interpretation results are highly consistent with the problem being solved (i.e. purpose). The interpretation results vary depending on their purpose. The emphasis of P- and S-wave velocity stratification is also different for different application requirements. For seismic exploration, the selection of the best lithology section for the design of the well depth for excitation is more focused on velocity and lithology, and a small S-wave velocity fitting time difference parameter should be selected; if the selection is for the extraction of high-precision static correction, the consistency of the interpretation accuracy of most measuring points should be emphasized, and the stratification fitting parameters should be regular and not change dramatically.

[0092] The output of the interactive interpretation scheme is the P-wave and S-wave joint interpretation result data and the P-wave and S-wave joint interpretation result map.

[0093] Example 2 Effect of the combined interactive interpretation method of P-wave and S-wave micro-logging data in surface structure investigation

[0094] In the Daqing exploration area, three three-component seismic exploration 3D acquisition projects were completed using the near-surface structure micro-logging longitudinal and shear wave joint synchronous investigation and processing and interpretation technology, and a total of 1,237 three-component micro-logging wells were processed.

[0095] This embodiment adopts the method of embodiment 1, and selects the longitudinal and transverse wave micro-logging data of the surface structure survey of the Alaxin work area in the western slope of the Songliao Basin in 2011, where the near-surface structure change law is relatively prominent. Figure 3 As shown, three points are selected for single-point P-wave and S-wave interactive interpretation. The five-star marks in the figure represent the surface structure investigation results of Jiangwan river channel area, dry field area and high hill area respectively. Through comparative analysis of various interactive interpretation parameter combinations, the process and effect of the joint interactive interpretation of P-wave and S-wave micro-logging data of the present invention are demonstrated, and the ideas of the interactive interpretation methods of P-wave and S-wave micro-logging data in areas with different surface sedimentary characteristics are summarized.

[0096] Analysis 1) Shear wave fitting time difference parameter analysis

[0097] Four fitting time difference parameters were analyzed for the three points in Jiangwan River Channel, Dry Field and Gaogang District. A single variable was maintained, that is, other interactive interpretation parameters and P-wave stratification results were fixed with the best results, and only the fitting time difference parameters were different. The results are as follows:

[0098] ①Jiangwan River Area

[0099] Table 1 Interactive interpretation results of Jiangwan District

[0100] Fitting time difference parameters Shear wave fitting layer accuracy Lithology interface Layering within the lithology section Interaction Explanation Effect 0.5ms high Match accurate very good 1.0ms accurate Basic match none good 1.5ms generally Basic match none generally 2.0ms generally Not consistent none not good

[0101] This physical point is located in the Nenjiang River Bay. This area is affected by the diversion of the river, seasonal water flow changes and the amount of water flowing, which makes the surface sediments in this area very complex and the surface sediments change rapidly. In addition, the water table in this area is shallowly buried. The longitudinal wave micro-logging can only observe two layers of velocity, and the lithology changes below the water table cannot be observed. However, the shear wave velocity stratification can clearly identify the lithology changes. By analyzing the four fitting time differences of this point, it can be seen that when the fitting time difference is small, that is, 0.5ms, Figure 4 As shown in the figure, the shear wave stratification is relatively fine, and it can distinguish the local small layers within the same lithology segment, which is conducive to finding the stable lithology medium layer, that is, finding the best exciting lithology layer. However, it is not conducive to distinguishing the interface of sedimentary lithology, and cannot clearly describe the changes in the shear wave surface structure. Figure 5 As shown in the figure, from the interactive interpretation results with a fitting time difference of 1.0 ms, it can be seen that the position of the stratified nodes is reasonable, the variation law of the shear wave velocity is clearly described, and the lithological interface below the P-wave high-speed layer can be well explained. Figure 6-7 The two shear wave fitting time differences of 1.5ms and 2.0ms cannot reflect the shear wave velocity stratification, and only the lithology interface position below the P-wave high-speed layer is observed. Therefore, after comprehensive analysis, it is believed that the fitting time difference of 0.5ms is Figure 4 The interactive explanation effect is better.

[0102] ② Dryland area

[0103] Table 2 Interaction interpretation results of dryland area

[0104] Fitting time difference parameters Shear wave fitting layer accuracy Lithology interface Layering within the lithology section Interaction Explanation Effect 0.5ms high Basic match exist generally 1.0ms accurate Not consistent none good 1.5ms generally Not consistent none generally 2.0ms generally Basic match none generally

[0105] This physical point is located in the dry field, which is a low-lying area with gentle terrain. The interpretation results of different fitting time difference parameters are shown in Figure 8-10 Since there are two optional locations for the top node of the longitudinal wave at this point, which are about 8 meters and 10 meters, and the node of the transverse wave at about 8 meters is relatively clear, you can see the details. Fig. 9 Therefore, it can be seen from this point that the stratification node position of the shear wave has a certain guiding role in the stratification of the longitudinal wave, and it is of great significance to improve the stratification accuracy of the longitudinal wave micro-logging.

[0106] ③ Gaogang area

[0107] This physical point is located on a high hill, which is a relatively flat area on the high hill. The water table is buried deep, and the surface sedimentary lithology is relatively stable. The P-wave micro-logging data of this physical point also has only two simple layers, and the high-speed layer velocity stratification interface of this physical point is the lithology stratification interface. The lithology interface above the high-speed layer interface is also not observed. Through the comparative analysis of four interactive interpretation parameters (see Table 3), it is concluded that the fitting time difference parameter 1 has the best interactive interpretation effect. The shear wave velocity stratification is more accurate in describing the interface of surface sedimentary lithology changes. The interpretation results of different fitting time difference parameters are shown in Table 3. Figure 11-14 .

[0108] Table 3 Interaction interpretation results of Gaogang District

[0109] Fitting time difference parameters Shear wave fitting layer accuracy Lithology interface Layering within the lithology section Interaction Explanation Effect 0.5ms high Basic match exist confusion 1.0ms accurate Match exist good 1.5ms generally Basic warming none generally 2.0ms generally Basic match exist generally

[0110] Through the comparative analysis of the four fitting time difference values ​​used in the joint interactive interpretation of the longitudinal and shear wave micro-logging data of the three physical points with different geographical conditions and sedimentary characteristics, the following two conclusions are drawn:

[0111] First, the fitting time difference used in different geographical environments should be different, that is, in areas where hydrodynamics or other dynamics occur frequently, a small fitting time difference parameter should be adopted, roughly in the range of 0.5-1.0, in order to distinguish the change of physical properties (lithology change) and the change of internal compaction degree of the same physical property conditions (change of shear wave velocity). In other areas where the surface structure is relatively stable, the fitting time difference parameter can be appropriately increased, because in this area, whether it is from the source, sedimentary environment and other factors, it is relatively stable, and the formation time is also relatively long, and the traceability and regularity are strong, so the fitting time difference of about 1.0-1.5ms is more suitable in this block, see Table 4 for details.

[0112] Table 4 Statistics of optimal fitting parameters in different regions

[0113] Test area type The best fitting time difference parameters The best fitting parameters for interaction Jiangwan River Area 0.5ms 0.5ms Dryland area 1.0ms 1.0ms Gaogang District 1.0ms 1.0ms

[0114] Second, the fitting time difference used should be different according to actual needs. The purpose of surface structure investigation is nothing more than two: to determine the excitation position and to build a fine and accurate longitudinal and transverse wave surface structure model. The determination of the precise excitation position requires that the transverse wave fitting time difference can distinguish the stable lithology medium of 2-3 meters, that is, the layer with a transverse wave velocity gradient of "0" at 2-3 meters. Through the analysis of the above points, a fitting time difference of about 0.5 can meet this requirement.

[0115] The conditions for constructing a fine and accurate P- and S-wave surface structure model are: it must be able to reflect the velocity change law of the S-wave in the vertical direction, and it must also reflect the change law and track the layers in the horizontal direction, thereby ensuring the traceability, regularity and rationality of the stratigraphic positions of the P- and S-wave models. Therefore, in the process of interactive interpretation of P- and S-wave micro-logging data, the lateral change law of sedimentary lithology is determined by analyzing the laws of the fitting time difference parameters in each area of ​​the exploration area, so as to complete the construction of an accurate and reasonable P- and S-wave near-surface model of the exploration area.

[0116] Analysis 2) Analysis of the implementation examples of longitudinal and transverse wave static correction results

[0117] Fig.15 is the plane diagram of static correction of longitudinal wave detection points, Fig.16 is a plane diagram of the static correction of shear wave detection points. The background of the two plane diagrams is the surface elevation color code change diagram of the work area. It is not difficult to see from the figure that in areas with high terrain and low diving surface, the relationship between the static correction of P-wave and S-wave and the change of elevation is consistent. It also shows that in this type of area, the consistency of P-wave and S-wave stratification and the consistency of velocity change law are good. In low-lying areas and river areas with high diving surface, the change regularity of P-wave static correction is not strong, because the high-speed top interface in this area is the diving surface position, and it does not respond to changes in other factors. The static correction of S-wave makes its regularity and the description of sedimentary lithology media clearer due to the application of velocity stratification method based on fitting parameter optimization, which also reflects the advantages of mutual restraint and complementarity in the joint interpretation process of near-surface P-wave and S-wave micro-logging data.

[0118] The above analysis shows that the method of the present invention can demonstrate its high efficiency and reasonable accuracy of results in the joint interpretation of P-wave and S-wave logging data in surface structure investigation, and can also draw accurate and reasonable conclusions by being applied to the construction of P-wave and S-wave surface structure models.

[0119] It should be noted that the above is only a preferred embodiment of the present invention and is not intended to limit the present invention. Although the present invention is described in detail with reference to the above embodiments, those skilled in the art can still modify the technical solutions described in the above embodiments or replace some of the technical features therein by equivalents. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the scope of protection of the claims of the present invention.

Claims

1. A joint interactive interpretation method of P-wave and S-wave micro-logging data in surface structure investigation. It is characterized in that The process includes the following steps: S1. Establishment of a joint interactive processing and interpretation framework for P-wave and S-wave micrologging data Display the first arrival data of P-wave and S-wave micro-logging in the main window, and establish a joint interactive processing and interpretation framework of P-wave and S-wave micro-logging data, including manual interactive interpretation mode, automatic interactive interpretation mode and auxiliary function control module; The P-wave and S-wave micro-logging first arrival data include P-wave first arrival data and S-wave first arrival data, and the two types of data are displayed in sub-windows under the main window with the same scale respectively; S2. Conduct real-time joint interactive communication In the framework of joint interactive processing and interpretation of P-wave and S-wave micro-logging data, joint interactive interpretation is adopted for the P-wave and S-wave micro-logging first arrival data. By monitoring the information of P-wave and S-wave velocity layers in their respective sub-windows and calling the P-wave and S-wave interpretation process code, the intermediate results and final results are interpreted using structure and common variable pointers. Use private data structure pointer variables to process private data in intermediate results and final results, dynamically allocate storage space and save private data; The common data common body pointer variable is used to process the common data in the first arrival data of the P-wave and S-wave micro-logging, and the storage space is dynamically allocated to realize the real-time joint interactive communication of the velocity layered interpretation of the first arrival data of the P-wave and S-wave micro-logging. S3. Implementing interactive explanation The longitudinal wave first arrival data is processed by linear function fitting to obtain the longitudinal wave layered node results, which are displayed in a sub-window on one side; The shear wave first arrival data is processed by nonlinear function fitting, and the shear wave layered node results are obtained by human-computer interaction and displayed in the sub-window on the other side; S4. Selection and output of interactive interpretation scheme The selection criteria for interactive interpretation schemes include: consistency of water-free stratum interfaces in P- and S-wave stratification results, consistency of lithologic stratification interfaces, and consistency of interactive interpretation results with environmental areas and uses; The output of the interactive interpretation scheme is the P-wave and S-wave joint interpretation result data and the P-wave and S-wave joint interpretation result map.

2. The method for combined interactive interpretation of P-wave and S-wave micro-logging data in surface structure investigation according to claim 1, It is characterized in that The auxiliary function control module includes a parameter setting module, an automatic interactive refinement comparison module, an interpretation scheme saving module and a pre-stored interpretation scheme retrieval module; Wherein, the parameter setting module includes fixed parameters; The function of the automatic interactive refinement comparison module is to apply fixed parameters to coordinate and control the velocity layering nodes of one wave to achieve the consistency of the longitudinal and transverse wave layering nodes; The function of the interpretation scheme saving module is to save the interactive interpretation results in the form of independent files; The function of the pre-stored interpretation scheme retrieval module is to call and view the interpretation scheme saving module to realize the interpretation scheme selection and output the interactive interpretation scheme.

3. The method for joint interactive interpretation of P-wave and S-wave micro-logging data in surface structure investigation according to claim 1, It is characterized in that The processing results of the first arrival data of the longitudinal and shear wave micrologging are displayed in respective sub-windows through interpretation graphics.

4. The combined interactive interpretation method of P-wave and S-wave micro-logging data in surface structure investigation according to claim 3, It is characterized in that The interpretation graphics are drawn by using the VCL graphics component of C++Builder to process the first arrival data of the longitudinal and shear wave micrologging.

5. The method for combined interactive interpretation of P-wave and S-wave micro-logging data in surface structure investigation according to claim 1, It is characterized in that The linear function is H=V×T, Among them, H represents depth; V represents speed; T represents the corresponding time; The nonlinear function is Vs=a i +b i ×H, Where Vs represents the shear wave layering velocity; a i represents the intercept of the i-th segment fitting curve; b i represents the fitting shear wave velocity gradient of the ith section; H represents the depth of the ith section.

6. The method for combined interactive interpretation of P-wave and S-wave micro-logging data in surface structure investigation according to claim 1, It is characterized in that The human-computer interaction method consists of manual interaction interpretation and automatic interaction interpretation; The manual interactive interpretation is to manually move the hierarchical results of the sub-window on one side to control the corresponding update of the hierarchical results in the sub-window on the other side.

7. The method for combined interactive interpretation of P-wave and S-wave micro-logging data in surface structure investigation according to claim 2, It is characterized in that The automatic interactive interpretation is completed by inputting the fixed parameters after the manual interactive interpretation is completed, and the P-wave and S-wave micro-logging velocity layered interpretation result diagram is obtained.

8. The method for combined interactive interpretation of P-wave and S-wave micro-logging data in surface structure investigation according to claim 7, It is characterized in that The fixed parameters include: a P-wave velocity fitting threshold value, a S-wave fitting time difference coefficient, an allowable range value of the P-wave and S-wave velocity stratification interface, a percentage parameter of the velocity difference between adjacent P-wave and S-wave layers, and a maximum distance parameter between a S-wave stratification node and the interface.

9. The method for joint interactive interpretation of P-wave and S-wave micro-logging data in surface structure survey according to any one of claims 1 to 8, It is characterized in that The method is implemented by software with modules for manual interactive interpretation, automatic interactive interpretation, detailed comparison, interactive interpretation control parameter setting and interpretation scheme saving.

Images