Efficient collaborative construction method, device and equipment for simultaneous operation of multiple geophysical prospecting teams

Abstract

The invention provides an efficient cooperative construction method, device and equipment for simultaneous operation of multiple geophysical prospecting teams, relates to the technical field of seismic prospecting, and aims to solve the problem that the construction efficiency is low under the condition that the distance between shot points constructed by the multiple geophysical prospecting teams in the same time window in adjacent work areas is smaller than the minimum distance of synchronous excitation in the prior art. The method comprises the following steps: integrating adjacent work areas for analysis, calculating the lower limit value of the root-mean-square energy of the upper cannon interference of a conventional production cannon, and drawing a first time distance curve meeting the condition that the upper cannon interference is not greater than that of the production cannon by utilizing a single-cannon data single-cannon energy attenuation time-space domain distribution plane graph of a long-record long-arrangement-point test; according to the seismic data of the work area, obtaining the lowest refracted wave speed, and calculating a second time distance relation curve when the last shot refracted wave does not enter the single shot data; and an area clamped by the first gun loading interference time distance curve and the second gun loading interference time distance curve is used as a distance range of staggered excitation of two guns in adjacent work areas, a blasting time interval is dynamically calculated, and staggered excitation is carried out.

Description

TECHNICAL FIELD

[0001] Embodiments of the present application relate to the technical field of seismic exploration, in particular, relate to a multi-geophysical team simultaneous operation efficient collaborative construction method, device and equipment. BACKGROUND

[0002] With the acceleration of the seismic exploration process, the situation of adjacent work areas and simultaneous window data collection is becoming more and more common in seismic data acquisition work. The construction team in the adjacent work area avoids affecting the quality of seismic data and the construction efficiency by using synchronous excitation. Synchronous excitation needs to set the minimum excitation distance to meet the mutual non-interference of seismic data.

[0003] However, in the actual construction process, there are two cases where the shotpoint distance of the two work areas is less than the minimum excitation distance of synchronous excitation. At this time, the time-sharing excitation method is usually used. Although this method can avoid mutual influence in the seismic data acquisition process, it will greatly reduce the acquisition efficiency and increase the acquisition cost. SUMMARY

[0004] Embodiments of the present application provide a multi-geophysical team simultaneous operation efficient collaborative construction method, device and equipment, aiming to solve the problem of low construction efficiency of related technologies in the case of multi-geophysical team simultaneous operation in adjacent work areas with the same window construction shotpoint distance less than the minimum distance of synchronous excitation.

[0005] The first aspect of the embodiments of the present application provides a multi-geophysical team simultaneous operation efficient collaborative construction method, comprising:

[0006] integrating adjacent work areas to obtain a target work area;

[0007] selecting a target historical seismic data acquisition work area according to the target work area, using a statistical method to analyze the root mean square value of the background noise energy of the original single shot data of the target historical seismic data acquisition work area, and obtaining the upper shot interference energy lower limit value by removing outliers;

[0008] performing a long-distance long-array point test in the target work area to obtain seismic data, and drawing a single shot energy decay space-time domain distribution plan according to the target single shot data in the seismic data;

[0009] According to the single shot energy decay space-time domain distribution plan, the contour line corresponding to the upper shot interference energy lower limit value is extracted as a first upper shot interference time-distance curve; the first upper shot interference time-distance curve represents the minimum distance of the two shots of the time-delay excitation;

[0010] According to the lowest refraction wave velocity obtained from the seismic data, a second upper shot interference time-distance curve of the upper shot refraction wave not entering the lower shot data body is calculated; the second upper shot interference time-distance curve represents the maximum distance of the two shots of the time-delay excitation;

[0011] The area enclosed by the first upper shot interference time-distance curve and the second upper shot interference time-distance curve is defined as the distance range between the two shots that are staggered in the adjacent work area.

[0012] Based on the distance range between the two shots during staggered firing in the adjacent work areas, a dynamic staggered firing interval is set according to the divided shot point areas of the adjacent work areas to perform staggered firing.

[0013] In one alternative implementation, the target historical seismic data acquisition area refers to a region with the same surface type, near-surface conditions, and excitation parameters as the target acquisition area.

[0014] In one optional implementation, the long-distance, long-arrangement point test conducted in the target work area satisfies the following conditions:

[0015] The single-shot record length of the long-distance, long-arrangement point test is greater than the firing interval of the target historical seismic data acquisition area.

[0016] The longitudinal receiving arrangement length of the long-distance long-arrangement point test is the longitudinal full arrangement of the target work area.

[0017] In one optional implementation, the step of calculating the second upper-shot interference time-distance curve, based on the lowest refracted wave velocity obtained from the seismic data, to prevent the upper-shot refracted wave from entering the lower-shot data volume, includes:

[0018] The second upper gun interference time-distance curve is calculated using the following formula:

[0019] S = T × VX max ;

[0020] Where S represents the maximum distance between the two guns; T represents the staggered firing time; V represents the minimum refracted wave velocity; X max This indicates the maximum shot-receiver distance of the three-dimensional observation system;

[0021] The minimum refracted wave velocity and the maximum shot-receiver distance of the three-dimensional observation system are known. Based on the linear relationship of the formula, the second upper shot interference time-distance curve is obtained.

[0022] In one optional implementation, before setting the dynamic time-shifting interval and performing time-shifting excitation, the following steps are included:

[0023] The computer is timed using the Global Positioning System and the BeiDou Navigation Satellite System. Based on the original single time slot for synchronous excitation, multiple time slots are added to meet the requirements of the staggered excitation.

[0024] In an alternative embodiment, the method further comprises:

[0025] According to the two-shot distance range of the adjacent work area, a time-lapse interval table is set, in which the time-lapse interval time corresponds to the two-shot distance range one by one;

[0026] According to the distance between the corresponding shot points in the shot point area, the shortest time-lapse interval time corresponding to the two-shot distance range matched in the time-lapse interval table is selected as the dynamic time-lapse interval, and the time-lapse shooting is performed.

[0027] In an alternative embodiment, the method further comprises:

[0028] According to the two-shot distance range of the adjacent work area, a time-lapse interval table is set, in which the time-lapse interval time corresponds to the two-shot distance range one by one;

[0029] According to the two-shot distance range of the adjacent work area, a time-lapse interval table is set, in which the time-lapse interval time corresponds to the two-shot distance range one by one;

[0030] The corresponding two-shot distance range of each segment of the fixed time-lapse interval is the two-shot distance range corresponding to the longest time-lapse interval time in each segment of the fixed time-lapse interval;

[0031] According to the longest time-lapse interval time in each segment of the fixed time-lapse interval, the multiple segments of the fixed time-lapse interval are arranged in ascending order;

[0032] According to the distance between the corresponding shot points in the shot point area, the shortest time-lapse interval time corresponding to the two-shot distance range matched in the time-lapse interval table is selected as the dynamic time-lapse interval, and the time-lapse shooting is performed.

[0033] In an alternative embodiment, the time-lapse shooting is performed when the shot point distance of the multiple geophysical exploration teams in the adjacent work area is less than the minimum distance of the synchronous shooting.

[0034] The second aspect of the embodiments of the present application provides a multi-geophysical exploration team simultaneous operation efficient collaborative construction device, comprising:

[0035] The integration module is configured to integrate the adjacent work areas to obtain a target work area;

[0036] The analysis module is configured to select a target historical seismic data acquisition work area according to the target work area, analyze the root mean square value of the background noise energy of the original single-shot data of the target historical seismic data acquisition work area using a statistical method, and obtain a lower limit value of the upper shot interference energy by eliminating outliers.

[0037] a test module configured to perform a long-distance long-array point test in the target work area to obtain seismic data, and draw a single-shot energy attenuation space-time distribution plan according to target single-shot data in the seismic data;

[0038] an extraction module configured to extract an isogram corresponding to the upper shot interference energy lower limit value as a first upper shot interference time-distance curve according to the single-shot energy attenuation space-time distribution plan; the first upper shot interference time-distance curve represents a minimum distance between two shots in a staggered time excitation mode;

[0039] a calculation module configured to calculate a second upper shot interference time-distance curve of a refracted wave of an upper shot not entering a lower shot data volume according to a lowest refracted wave velocity obtained from the seismic data; the second upper shot interference time-distance curve represents a maximum distance between two shots in a staggered time excitation mode;

[0040] a region definition module configured to take a region sandwiched by the first upper shot interference time-distance curve and the second upper shot interference time-distance curve as a distance range between two shots in a staggered time excitation mode in the adjacent work area;

[0041] a staggered time excitation module configured to set a dynamic staggered time interval according to a divided shot point region of the adjacent work area based on the distance range between two shots in a staggered time excitation mode in the adjacent work area, and perform staggered time excitation.

[0042] The third aspect of the embodiment of the present application provides an electronic device, which comprises a memory, a processor, and a computer program stored in the memory, and the processor executes the computer program to implement the steps of the multi-geophysical exploration team simultaneous operation efficient collaborative construction method of the first aspect.

[0043] In the embodiment, the distance range between two shots in a staggered time excitation mode in the adjacent work area is obtained by calculating a first time-distance relationship curve of an upper shot interference not greater than a production shot and a second time-distance relationship curve of an upper shot refracted wave not entering single-shot data, a dynamic staggered time interval is calculated, and staggered time excitation is performed. The distance range between two shots in a staggered time excitation mode that meets the qualitative and quantitative evaluation criteria of seismic data is scientifically calculated, the dynamic staggered time interval is calculated in the construction process, staggered time excitation is performed, the technical advantages of staggered time excitation are fully utilized, the quality of seismic data is ensured, the exploration efficiency is improved, and the exploration cost is reduced. BRIEF DESCRIPTION OF DRAWINGS

[0044] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the following will briefly introduce the drawings needed to be used in the description of the embodiments of the present application. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor.

[0045] Figure 1is a step flow chart of the efficient collaborative construction method of simultaneous operation of multiple geophysical exploration teams according to an embodiment of the present application;

[0046] Figure 2 is a schematic diagram of the lower limit value of the energy of the shot interference in statistics of the efficient collaborative construction method of simultaneous operation of multiple geophysical exploration teams according to an embodiment of the present application;

[0047] Figure 3 is a planar diagram of the spatial and temporal domain distribution of the single-shot energy attenuation of the efficient collaborative construction method of simultaneous operation of multiple geophysical exploration teams according to an embodiment of the present application;

[0048] Figure 4 is a schematic diagram of the distance range of the two shots of the staggered time excitation of the adjacent work areas of the efficient collaborative construction method of simultaneous operation of multiple geophysical exploration teams according to an embodiment of the present application;

[0049] Figure 5 is a schematic diagram of the setting of the dynamic staggered time interval of the efficient collaborative construction method of simultaneous operation of multiple geophysical exploration teams according to an embodiment of the present application;

[0050] Figure 6 is a schematic diagram of the increase of a time slot on the basis of the synchronous excitation of the efficient collaborative construction method of simultaneous operation of multiple geophysical exploration teams according to an embodiment of the present application;

[0051] Figure 7 is a schematic diagram of the three-dimensional planar position of the adjacent work areas of the efficient collaborative construction method of simultaneous operation of multiple geophysical exploration teams according to an embodiment of the present application;

[0052] Figure 8 is a schematic diagram of an electronic device according to an embodiment of the present application. DETAILED DESCRIPTION

[0053] The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are some embodiments of the present application, rather than all the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative work fall within the scope of protection of the present application.

[0054] In the drawings, sometimes the size of the constituent elements, the thickness of the layers, or the area is exaggerated for the sake of clarity, and therefore, any one of the implementations of the present disclosure is not necessarily limited to the size shown in the drawings. The shapes and sizes of the components in the drawings do not reflect the actual proportions. Furthermore, the drawings schematically show ideal examples, and any one of the implementations of the present disclosure is not limited to the shapes or values shown in the drawings.

[0055] REFERENCE Figure 1 , Figure 1is a step flow chart of a high-efficiency collaborative construction method of multiple geophysical exploration teams simultaneously operating, which is proposed by an embodiment of the present application, and includes the following steps S11-S17.

[0056] Step S11: Integrating adjacent work areas as a whole to obtain a target work area.

[0057] In a seismic data acquisition project, multiple geophysical exploration teams usually simultaneously operate in adjacent work areas. Integrating these originally independent seismic acquisition projects into a whole seismic acquisition project, i.e., integrating adjacent work areas as a whole. This involves unified planning, coordinating the operation time, firing parameters, data acquisition standards, etc. of each geophysical exploration team, so as to facilitate subsequent analysis and calculation.

[0058] Step S12: Selecting a target historical seismic data acquisition work area according to the target work area, using statistical methods to analyze the root mean square value of background noise energy of the original single-shot data of the target historical seismic data acquisition work area, and obtaining the upper-shot interference energy lower limit value by removing outliers.

[0059] Background noise energy is the sum of all other seismic wave energy before or after each shot. Upper-shot interference is the interference of the seismic wave generated by the previous shot on the subsequent shot record in the continuous seismic data acquisition process.

[0060] Select a previous seismic data acquisition work area with the same surface type, similar near-surface conditions, and the same firing parameters as the target work area, obtain the original single-shot data from the seismic data storage library of the target historical seismic data acquisition work area, calculate the root mean square value of the energy of each segment of background noise data. Use statistical methods (such as box plot, Z-score method, etc.) to identify and remove outliers for all background noise energy root mean square values to ensure the accuracy and reliability of subsequent analysis. Perform statistical analysis on the background noise energy root mean square values after removing outliers, and select a suitable statistical quantity as the upper-shot interference energy lower limit value according to the distribution of the data, such as shown in the formula. Figure 2 The upper-shot interference energy lower limit value represents the lowest limit value of the energy level of the interference signal in most cases, and provides an important basis for subsequent data analysis and processing.

[0061] Further, the target historical seismic data acquisition work area refers to a region with the same surface type, near-surface conditions, and firing parameters as the target work area.

[0062] The previous seismic data acquisition area near the target work area can be selected as a reference, which has similar surface types (such as soil types, rock types, vegetation coverage, etc.), near-surface conditions (such as terrain undulations, groundwater levels, weathering layer thicknesses, etc.), and the same excitation parameters (such as source types, excitation energy, excitation depth, etc.) as the target work area. It can provide high-quality reference data for the determination of the subsequent upper gun interference energy lower limit value, and improve the accuracy and reliability of the data analysis results.

[0063] Step S13: performing a long-distance long-array point test in the target work area to obtain seismic data, and drawing a single-shot energy attenuation space-time distribution plan according to target single-shot data in the seismic data.

[0064] The long-distance long-array point test is performed in the target work area to collect seismic data. The "long distance" refers to a relatively long path of seismic wave propagation, so as to observe how the energy attenuates with the increase of distance. The "long array" refers to that the geophones for receiving seismic waves are arranged very long in the longitudinal direction, so as to record the information of the arrival of seismic waves at different positions at different times.

[0065] The energy extraction of the target single-shot data in the collected seismic data usually involves calculating the energy (such as the square sum of amplitude) or energy density (such as energy per unit time) of each channel data, and can also use instantaneous energy, spectral energy, etc. The extracted energy data is arranged in time and space order, and a single-shot energy attenuation space-time distribution plan is drawn by using a seismic data processing software or a drawing tool, as shown in FIG. 2, so that the changes of energy with time and space after well shooting can be directly observed. Figure 3

[0066] Further, the long-distance long-array point test in the target work area meets the following conditions:

[0067] The single-shot recording length of the long-distance long-array point test is greater than the shooting time interval of the target historical seismic data acquisition area.

[0068] The longitudinal receiving array length of the long-distance long-array point test is the longitudinal full array of the target work area.

[0069] ​The long-distance long-array point test on the target work area typically surface needs to meet: single shot record length is greater than the target historical seismic data acquisition work area shooting time interval; longitudinal receiving array length is the longitudinal full array of the target work area of this construction. The length of single shot record is set to be greater than the previous shooting time interval, to ensure that the seismic wave generated by the previous excitation is completely recorded after complete attenuation, and then the next excitation is carried out, avoiding the influence of the previous excitation on the record of the next excitation. The length of longitudinal receiving array is the longitudinal full array length of this construction, that is, the geophones are arranged according to the maximum length of the construction design in the longitudinal direction, so as to obtain more comprehensive information.

[0070] Step S14: According to the single shot energy attenuation space-time distribution plan, the contour line corresponding to the upper shot interference energy lower limit value is extracted as the first upper shot interference time-distance curve; the first upper shot interference time-distance curve represents the minimum distance of two shots of staggered excitation.

[0071] On the single shot energy attenuation space-time distribution plan, according to the upper shot interference energy lower limit value, the energy contour line of the shooting time interval corresponding to different distances is extracted, such as Figure 3 The arrow points to the contour line, which is the first upper shot interference time-distance curve that meets the quantitative evaluation. The curve directly shows that under different shooting time intervals, in order to ensure that the shooting interference does not exceed the lower limit of background noise energy, the minimum distance of two shots is required, and the first upper shot interference time-distance curve can represent the minimum distance of two shots of staggered excitation. It provides a scientific basis for the shooting operation in actual production to ensure the quality of seismic data acquisition.

[0072] Step S15: According to the lowest refraction wave velocity obtained from the seismic data, calculate the second upper shot interference time-distance curve of the upper shot refraction wave not entering the lower shot data body; the second upper shot interference time-distance curve represents the maximum distance of two shots of staggered excitation.

[0073] The lowest refraction wave velocity is the apparent velocity of the refraction wave in the underground medium due to the change of medium velocity. It can be determined by the refraction wave analysis technology in seismic data processing, and it is necessary to ensure that the selected refraction wave is stable and its velocity can represent the lowest velocity layer in the underground medium.

[0074] According to the determined minimum refraction wave velocity, in combination with the distance between the shot points and the propagation characteristics of the seismic wave, the shooting time interval required for the upper shot refraction wave not to enter the lower shot data body is calculated. The relationship between the distance between the shot points and the shooting time interval is the second upper shot interference time-distance curve that meets the qualitative evaluation. This curve intuitively shows the maximum distance between the two shots required to ensure that the upper shot refraction wave does not enter the lower shot data body under different shooting time intervals. The second upper shot interference time-distance curve can represent the maximum distance between the two shots of the staggered shooting. It provides a scientific basis for the shooting operation in actual production to ensure the quality and reliability of the seismic acquisition data.

[0075] Further, the second upper shot interference time-distance curve for the upper shot refraction wave not to enter the lower shot data body is calculated according to the minimum refraction wave velocity obtained from the seismic data, comprising:

[0076] The second upper shot interference time-distance curve is calculated according to the following formula:

[0077] S=T×V-X max ;

[0078] Wherein, S represents the maximum distance between the two shots; T represents the staggered shooting time; V represents the minimum refraction wave velocity; X max represents the maximum offset distance of the three-dimensional observation system.

[0079] The minimum refraction wave velocity and the maximum offset distance of the three-dimensional observation system are known, and the second upper shot interference time-distance curve is obtained according to the linear relationship of the formula.

[0080] The stable minimum velocity refraction wave velocity information, i.e. the minimum refraction wave velocity, is extracted from the seismic data, and the maximum offset distance of the three-dimensional observation system (the maximum straight-line distance between the shot points and the receiver points in the seismic observation system) is directly obtained in the seismic observation system file. In the case where the minimum refraction wave velocity and the maximum offset distance of the three-dimensional observation system are known, the second upper shot interference time-distance curve is obtained according to the linear relationship between the maximum distance between the two shots and the staggered shooting time in the following formula (1).

[0081] S=T×V-X max (1)

[0082] Wherein, S represents the maximum distance between the two shots; T represents the staggered shooting time; V represents the minimum refraction wave velocity; X max represents the maximum offset distance of the three-dimensional observation system.

[0083] Step S16: The area between the first upper shot interference time-distance curve and the second upper shot interference time-distance curve is taken as the range of the distance between the two shots of the staggered shooting in the adjacent work area.

[0084] The first upshot interference time-distance curve is used to determine that the interference strength of the upshot will not be higher than the lower limit of background noise energy, so as not to affect the data quality. The second upshot interference time-distance curve is used to ensure that the refracted wave of the upshot will not enter the data record of the downshot, so as to avoid interference between the data. The area between the two curves drawn in the same coordinate system is the time-lapse shooting distance range meeting the energy quantitative evaluation and single-shot qualitative evaluation requirements, as shown in FIG. 8. Figure 4 The above provides a scientific basis and technical support for seismic data acquisition, which is helpful to guide the time-lapse shooting in actual production, and improves the exploration efficiency while ensuring the data quality.

[0085] Step S17: Based on the two-shot distance range of the time-lapse shooting of the adjacent work area, a dynamic time-lapse interval is set according to the divided shot point area of the adjacent work area, and the time-lapse shooting is performed.

[0086] According to the actual situation of the engineering site, the acquisition shot points of the adjacent work area on the same day are finely divided, and a dynamic time-lapse interval is set for the shot points of the adjacent work area according to the divided shot point area and the determined two-shot distance range of the time-lapse shooting, as shown in FIG. 9. Figure 5 According to the set dynamic time-lapse interval, the shot points of the adjacent work area are sequentially shot. In the shooting process, the shooting time of each shot point should be strictly controlled to ensure the smooth progress of the time-lapse shooting. Giving full play to the technical advantages of the time-lapse shooting can effectively reduce the signal interference between adjacent work areas, improve the quality of seismic data, and realize the simultaneous acquisition of multiple work areas, thereby improving the efficiency and progress of seismic exploration.

[0087] Further, before the setting of the dynamic time-lapse interval and the time-lapse shooting, the method further includes:

[0088] The global positioning system and the Beidou satellite navigation system are used to time the computer, and based on the original single time slot of the synchronous shooting, multiple time slots are added to meet the requirements of the time-lapse shooting.

[0089] The global positioning system and the Beidou satellite navigation system are used to time the computer, and based on the original single time slot of the synchronous shooting, multiple time slots are added to meet the requirements of the time-lapse shooting. Figure 6

[0090] ​Furthermore, the step of setting a dynamic time-staggered interval based on the distance range between the two shots in the adjacent work areas and according to the divided shot point areas of the target work area, and performing time-staggered firing, includes:

[0091] Based on the distance range between the two shots in the adjacent work areas, a staggered time interval table is set up, and the staggered time excitation time in the staggered time interval table corresponds one-to-one with the distance range between the two shots.

[0092] Based on the distance between corresponding shot points in the shot point area, the shortest staggered firing time corresponding to the two shot distance range that matches the staggered firing interval table is selected as the dynamic staggered firing interval, and staggered firing is performed.

[0093] Based on the distance range between two shots fired at staggered intervals in adjacent work areas, the staggered firing time corresponding to the two shot distance ranges is calculated, and a staggered firing interval table is set. The results are shown in Table 1. Based on the distance between corresponding shot points in the shot point areas, the matching two shot distance ranges are found in the set dynamic staggered firing interval table, and the shortest staggered firing time within that range is selected as the dynamic staggered firing interval for staggered firing. The set dynamic staggered firing interval is as follows: Figure 5 As shown, selecting the shortest time can maximize the efficiency of seismic exploration.

[0094] Table 1. Time Interval Table

[0095]

[0096] Furthermore, the method also includes:

[0097] Based on the distance range between the two shots during staggered firing in the adjacent work areas, and according to the defined shot point areas of the target work area, a fixed staggered firing interval is set for staggered firing, including:

[0098] Based on the distance range between the two shots fired at staggered times in adjacent work areas, multiple fixed staggered time intervals are divided;

[0099] The range of two-shot distances corresponding to each fixed staggered time interval is taken as the range of two-shot distances corresponding to the longest staggered firing time in each fixed staggered time interval.

[0100] Based on the longest time-out excitation time in each fixed time-out interval, the multiple fixed time-out intervals are arranged in ascending order;

[0101] Based on the distance between corresponding shot points in the shot point area, select the minimum fixed time interval corresponding to the two shot distance range that matches the fixed time interval, and perform time-staggered excitation.

[0102] Based on the distance range between two shots during staggered excitation in adjacent work areas, the entire distance range is divided into multiple sub-intervals, i.e., multiple fixed staggered excitation intervals. Each sub-interval corresponds to a fixed staggered excitation time interval. Within each sub-interval, the distance range between two shots corresponding to the longest staggered excitation time within that interval is selected as the distance range between two shots in that sub-interval. The multiple fixed staggered excitation intervals are sorted in ascending order according to their corresponding longest staggered excitation times, as shown in Table 1. Based on the distance between corresponding shot points in the shot point areas, the smallest fixed staggered excitation interval corresponding to the distance range between two shots in the fixed staggered excitation interval is selected for staggered excitation. Fixed staggered excitation is simple to operate, easy to manage and standardize, but compared with dynamic staggered excitation, the efficiency of seismic exploration is reduced. In actual construction, it can be selected according to different needs.

[0103] Furthermore, when multiple geophysical exploration teams are working in adjacent work areas at the same time window, and the distance between their blast points is less than the minimum distance for synchronous excitation, the staggered excitation is carried out.

[0104] When multiple geophysical exploration teams are working in adjacent work areas within the same time window, and the distance between shot points is less than the minimum distance for synchronous excitation, staggered excitation can be implemented. This fully leverages the technical advantages of staggered excitation, achieving the goal of improving acquisition efficiency while ensuring seismic data quality.

[0105] To facilitate understanding, the embodiments of this disclosure will be described through specific application scenarios.

[0106] In the Tazhong region of the Tarim Basin, where the surface is entirely desert, all drilling operations were conducted using well-fired seismic instruments. The Manshen 17-20 West 3D seismic survey and the Shunbei 6 3D seismic survey are adjacent. The Manshen 17-20 West 3D seismic survey was conducted by BGP (Broadway Geophysical Exploration) while the Shunbei 6 3D seismic survey was conducted by Shengli Geophysical Exploration, both using wired seismic instruments. The minimum distance between the shot points in the two areas is only 11 km, far less than the minimum distance for synchronous seismic acquisition in the region (42.5 km). Traditionally, a time-sharing acquisition method was used, which generally has two forms: time-segmented acquisition and alternating acquisition. Time-segmented acquisition divides a 24-hour day into several equal periods, with one acquisition item acquired separately within each period. This allows for a maximum of 12 hours of acquisition per day, but due to other processes, the acquisition time is generally around 8 hours, significantly reducing the acquisition time. The second type of time-sharing acquisition, derived from synchronous acquisition, uses alternating acquisition. For the same acquisition item, the acquisition time increases from 37 seconds in synchronous acquisition to 74 seconds, significantly increasing the time interval between shot blasts. The previous practice resulted in a significant decrease in data collection efficiency. It was estimated that by adopting time-sharing data collection, the efficiency would decrease from 1358 shots to 612 shots. According to the data collection plan, time-sharing data collection accounted for 17.6% of the data collection points, and 7.7% of the data collection points exceeded the required time window (data collection was to be completed before March 30th, with April being the windy season).

[0107] Analysis showed that increasing the daily efficiency from the planned 600 shots to 960 shots based on time-sharing data acquisition could achieve the target (a 60% increase in daily efficiency). Therefore, staggered-time excitation technology was adopted in this project. The background noise per shot was lower than that of normal production shots, and refracted waves did not enter the data volume of a single shot. Compared to the designed time-sharing excitation, the actual average daily efficiency increased by 72%, reaching an average of 1055 shots, exceeding the designed 960 shots. The data acquisition task was completed on March 28th, exceeding the set target. This technology saved 6 days of data acquisition time and directly reduced costs by 12.3 million yuan. Specific implementation details are as follows:

[0108] Following step S11, clarify the positional relationship between the 2022 full-depth 17-20 western three-dimensional geostationary network of PetroChina Tarim Oilfield Company and the Shunbei 6 three-dimensional geostationary network of Sinopec Northwest Bureau, such as... Figure 7 As shown, adjacent work areas are analyzed as a whole as the target work area, and the distance between the shot points in the two work areas is accurately calculated. The closest distance between the shot points in the two work areas is 11km, and some shot points do not meet the minimum distance requirement for synchronous firing.

[0109] According to step S12, near the target work area, the 3D excitation, reception, and instrument parameters of the 2021 Middle East 21 Block development are the same as those of the full-depth 17-full-depth 20 West 3D, which can be used as the original single-shot background noise statistical value for calculation. Based on statistical principles, the root mean square value of the background noise energy of the original single-shot data is analyzed. After removing outliers, 20uv is the lower limit of the upper-shot interference energy. Figure 2 As shown.

[0110] Following step S13, a long-distance, long-arrangement point test was conducted within a three-dimensional area of ​​full depth 17-20 m. The single-shot recording length was 60 s (compared to the previous firing interval of 35 s), and the longitudinal receiving array length was 26237.5 m. Based on the single-shot data obtained from this data, a planar distribution map of single-shot energy attenuation was plotted, as shown below. Figure 3 As shown.

[0111] According to step S14, the threshold value (20uv) is set as the upper gun interference intensity not exceeding the lower limit of the background noise energy of the production gun (lower limit of upper gun interference energy). Based on the spatiotemporal distribution map of single gun energy attenuation obtained in step S13, the contour lines corresponding to the lower limit of upper gun interference energy are extracted, such as... Figure 3 The middle arrow points to the contour line, which is the time-distance curve of the upper gun interference that meets the qualitative evaluation.

[0112] According to step S15, the maximum shot-receiver distance of the three-dimensional observation system is 8600m. The minimum refracted wave velocity information obtained through the seismic data of the target work area is 2000m / s. The specific data of the second upper shot interference time distance curve calculated according to formula (1) is shown in Table 2. The linear relationship between the maximum distance between the two shots and the staggered excitation time is the second upper shot interference time distance curve.

[0113] Table 2. Data Table

[0114]

[0115] According to step S16, the region enclosed by the first and second upper-shot interference time-distance curves satisfies both the quantitative evaluation requirement that the interference intensity of the upper shot will not exceed the lower limit of the background noise energy, and the qualitative evaluation requirement that the refracted wave of the upper shot will not enter the data record of the lower shot. This region is the distance range between the two shots in adjacent work areas that are staggered in time. Figure 4 As shown in Table 3, based on the distance range between two shots fired at different staggered times in adjacent work areas, a table of distance ranges corresponding to different staggered firing times can be obtained.

[0116] Table 3. Distance Ranges Corresponding to Different Staggered Time Periods

[0117] Maximum distance 23.4 25.4 27.4 29.4 31.4 33.4 35.4 37.4 39.4 41.4 43.4 45.4 kilometres Minimum distance 21 20 19 18 17 16 15 14 13 12 11 10 kilometres Time offset 16 17 18 19 20 21 22 23 24 25 26 27 seconds

[0118] Following step S17, based on the actual conditions of the engineering site, the blasting points in adjacent work areas are finely divided for the day. Based on the divided blasting point areas and the determined distance range between the two blasts for staggered firing, a dynamic staggered firing interval is set for the blasting points in adjacent work areas, such as... Figure 5 As shown. The firing points in adjacent work areas are sequentially fired according to the set dynamic time intervals.

[0119] The efficient collaborative construction method for multiple geophysical exploration teams operating simultaneously, as disclosed in this embodiment, obtains the distance range between two shots in adjacent work areas for staggered excitation by calculating a first time-distance relationship curve where the interference from the upper shot is no greater than that of the production shot and a second time-distance relationship curve where the refracted wave from the upper shot does not enter the single-shot data. The staggered time interval is dynamically calculated, and staggered excitation is performed. This scientifically calculates the distance range between two shots for staggered excitation that meets the qualitative and quantitative evaluation standards of seismic data. By dynamically calculating the staggered time interval during construction and performing staggered excitation, the technical advantages of staggered excitation are fully utilized, improving exploration efficiency and reducing exploration costs while ensuring the quality of seismic data.

[0120] It should be noted that, for the sake of simplicity, the method embodiments are all described as a series of actions. However, those skilled in the art should understand that the embodiments of the present invention are not limited to the described order of actions, because according to the embodiments of the present invention, some steps can be performed in other orders or simultaneously. Furthermore, those skilled in the art should also understand that the embodiments described in the specification are preferred embodiments, and the actions involved are not necessarily essential to the embodiments of the present invention.

[0121] Based on the same inventive concept, this application discloses a device for efficient collaborative construction of multiple geophysical exploration teams operating simultaneously. The device includes:

[0122] The integration module is used to integrate adjacent work areas to obtain the target work area;

[0123] The analysis module is used to select the target historical seismic data acquisition area based on the target work area, use statistical methods to analyze the root mean square value of the background noise energy of the original single shot data of the target historical seismic data acquisition area, and obtain the lower limit value of the upper shot interference energy by removing outliers.

[0124] The test module is used to conduct long-distance, long-arrange point tests in the target work area to obtain seismic data, and to draw a planar map of the spatiotemporal distribution of single-shot energy attenuation based on the target single-shot data in the seismic data.

[0125] The extraction module is used to extract the contour line corresponding to the lower limit of the upper gun interference energy as the first upper gun interference time-distance curve based on the spatiotemporal distribution plan of the single gun energy attenuation; the first upper gun interference time-distance curve represents the minimum distance between the two guns that are triggered at different times.

[0126] The calculation module is used to calculate the second upper shot interference time-distance curve, which indicates that the upper shot refracted wave does not enter the lower shot data volume, based on the lowest refracted wave velocity obtained from the seismic data; the second upper shot interference time-distance curve represents the maximum distance between the two shots that are triggered at different times.

[0127] The region definition module is used to define the region enclosed by the first upper shot interference time-distance curve and the second upper shot interference time-distance curve as the distance range between the two shots that are staggered in the adjacent work area.

[0128] The staggered firing module is used to set a dynamic staggered firing interval based on the distance range between the two shots in the adjacent work areas and the divided shot point areas of the adjacent work areas.

[0129] In one optional implementation, the analysis module is specifically used for:

[0130] The target historical earthquake data acquisition area refers to the region with the same surface type, near-surface conditions, and excitation parameters as the target acquisition area.

[0131] In one optional implementation, the test module is specifically used for:

[0132] The single-shot record length of the long-distance, long-arrangement point test is greater than the firing interval of the target historical seismic data acquisition area.

[0133] The longitudinal receiving arrangement length of the long-distance long-arrangement point test is the longitudinal full arrangement of the target work area.

[0134] In one optional implementation, the computing module is specifically used for:

[0135] The second upper gun interference time-distance curve is calculated using the following formula:

[0136] S = T × VX max ;

[0137] Where S represents the maximum distance between the two guns; T represents the staggered firing time; V represents the minimum refracted wave velocity; X max This indicates the maximum shot-receiver distance of the three-dimensional observation system;

[0138] The minimum refracted wave velocity and the maximum shot-receiver distance of the three-dimensional observation system are known. Based on the linear relationship of the formula, the second upper shot interference time-distance curve is obtained.

[0139] In one optional implementation, the time-staggered excitation module is specifically used for:

[0140] The computer is timed using the Global Positioning System and the BeiDou Navigation Satellite System. Based on the original single time slot for synchronous excitation, multiple time slots are added to meet the requirements of the staggered excitation.

[0141] In one optional implementation, the time-staggered excitation module is specifically used for:

[0142] Based on the distance range between the two shots in the adjacent work areas, a staggered time interval table is set up, and the staggered time excitation time in the staggered time interval table corresponds one-to-one with the distance range between the two shots.

[0143] Based on the distance between corresponding shot points in the shot point area, the shortest staggered firing time corresponding to the two shot distance range that matches the staggered firing interval table is selected as the dynamic staggered firing interval, and staggered firing is performed.

[0144] In one alternative embodiment, the device further includes:

[0145] A fixed-time staggered firing module is used to set a fixed-time staggered interval for staggered firing based on the distance range between two shots in adjacent work areas and the defined shot point areas of the adjacent work areas. This includes: dividing the distance range between two shots in adjacent work areas into multiple fixed-time staggered intervals; taking the distance range between two shots corresponding to the longest staggered firing time in each fixed-time staggered interval; arranging the multiple fixed-time staggered intervals in ascending order based on the longest staggered firing time in each fixed-time staggered interval; and selecting the smallest fixed-time staggered interval corresponding to the distance range between two shots that matches the distance between the shot point areas, and then performing staggered firing.

[0146] In one optional implementation, the time-staggered excitation module is specifically used for:

[0147] When multiple geophysical exploration teams are working in adjacent work areas at the same time window, and the distance between their blast points is less than the minimum distance for synchronous excitation, the staggered excitation is carried out.

[0148] This disclosure also provides an electronic device, with reference to... Figure 8 , Figure 8 This is a schematic diagram of an electronic device illustrated in an embodiment of this disclosure. For example... Figure 8 As shown, the electronic device 100 includes a memory 110 and a processor 120. The memory 110 and the processor 120 are connected via a bus for communication. The memory 110 stores a computer program, which can run on the processor 120 to realize the steps in the efficient collaborative construction method for multiple geophysical exploration teams operating simultaneously as disclosed in this embodiment.

[0149] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on the differences from other embodiments. The same or similar parts between the various embodiments can be referred to each other.

[0150] This application describes embodiments with reference to flowchart illustrations and / or block diagrams of methods, apparatuses, electronic devices, and computer program products according to embodiments of this application. It should be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing terminal device to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing terminal device, generate instructions for implementing the flowchart... Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0151] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing terminal device to operate in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0152] These computer program instructions can also be loaded onto a computer or other programmable data processing terminal equipment, causing a series of operational steps to be performed on the computer or other programmable terminal equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable terminal equipment for implementing the process. Figure 1 One or more processes and / or boxesFigure 1 The steps of the function specified in one or more boxes.

[0153] Although preferred embodiments of the present application have been described, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments as well as all changes and modifications falling within the scope of the embodiments of the present application.

[0154] Finally, it should be noted that in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or terminal device that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or terminal device. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or terminal device that includes said element.

[0155] The above provides a detailed description of a method, apparatus, and equipment for efficient collaborative construction involving multiple geophysical exploration teams operating simultaneously, as provided in this application. Specific examples have been used to illustrate the principles and implementation methods of this application. The descriptions of the above embodiments are only for the purpose of helping to understand the method and its core ideas. At the same time, those skilled in the art will recognize that there will be changes in the specific implementation methods and application scope based on the ideas of this application. Therefore, the content of this specification should not be construed as a limitation of this application.

Claims

1. A method for efficient collaborative construction involving multiple geophysical exploration teams operating simultaneously, characterized in that, include: By integrating adjacent work areas, the target work area is obtained; Based on the target work area, select the target historical earthquake data acquisition work area, use statistical methods to analyze the root mean square value of background noise energy of the original single shot data of the target historical earthquake data acquisition work area, and obtain the lower limit value of the upper shot interference energy by removing outliers. Long-distance, long-arrange point tests were conducted in the target work area to obtain seismic data. Based on the target single-shot data in the seismic data, a spatiotemporal distribution map of single-shot energy attenuation was drawn. Based on the spatiotemporal distribution plan of the single-shot energy attenuation, the contour line corresponding to the lower limit of the upper shot interference energy is extracted as the first upper shot interference time-distance curve; the first upper shot interference time-distance curve represents the minimum distance between the two shots that are triggered at different times. Based on the lowest refracted wave velocity obtained from the earthquake data, the second upper shot interference time-distance curve is calculated so that the upper shot refracted wave does not enter the lower shot data volume; the second upper shot interference time-distance curve represents the maximum distance between the two shots triggered at different times. The area enclosed by the first upper shot interference time-distance curve and the second upper shot interference time-distance curve is defined as the distance range between the two shots that are staggered in the adjacent work area. Based on the distance range between the two shots during staggered firing in the adjacent work areas, a dynamic staggered firing interval is set according to the divided shot point areas of the adjacent work areas to perform staggered firing.

2. The efficient collaborative construction method for multiple geophysical exploration teams operating simultaneously according to claim 1, characterized in that, The target historical earthquake data acquisition area refers to the region with the same surface type, near-surface conditions, and excitation parameters as the target acquisition area.

3. The method for efficient collaborative construction of multiple geophysical exploration teams operating simultaneously according to claim 1, characterized in that, The long-distance, long-arrangement point test conducted in the target work area shall meet the following conditions: The single-shot record length of the long-distance, long-arrangement point test is greater than the firing interval of the target historical seismic data acquisition area. The longitudinal receiving arrangement length of the long-distance long-arrangement point test is the longitudinal full arrangement of the target work area.

4. The efficient collaborative construction method for multiple geophysical exploration teams operating simultaneously according to claim 1, characterized in that, The step of calculating the second upper shot interference time-distance curve, based on the lowest refracted wave velocity obtained from the seismic data, whereby the upper shot refracted wave does not enter the lower shot data volume, includes: The second upper gun interference time-distance curve is calculated using the following formula: S=T×V-X max ; Where S represents the maximum distance between the two guns; T represents the staggered firing time; V represents the minimum refracted wave velocity; X max This indicates the maximum shot-receiver distance of the three-dimensional observation system; The minimum refracted wave velocity and the maximum shot-receiver distance of the three-dimensional observation system are known. Based on the linear relationship of the formula, the second upper shot interference time-distance curve is obtained.

5. The efficient collaborative construction method for multiple geophysical exploration teams operating simultaneously according to claim 1, characterized in that, Before setting the dynamic time-off interval and performing time-off triggering, the following are included: The computer is timed using the Global Positioning System and the BeiDou Navigation Satellite System. Based on the original single time slot for synchronous excitation, multiple time slots are added to meet the requirements of the staggered excitation.

6. The method for efficient collaborative construction of multiple geophysical exploration teams operating simultaneously according to claim 1, characterized in that, The method of staggered firing of two shots based on the distance range between adjacent work areas, and setting a dynamic staggered firing interval according to the divided shot point area of ​​the target work area, includes: Based on the distance range between the two shots in the adjacent work areas, a staggered time interval table is set up, and the staggered time excitation time in the staggered time interval table corresponds one-to-one with the distance range between the two shots. Based on the distance between corresponding shot points in the shot point area, the shortest staggered firing time corresponding to the two shot distance range that matches the staggered firing interval table is selected as the dynamic staggered firing interval, and staggered firing is performed.

7. The efficient collaborative construction method for multiple geophysical exploration teams operating simultaneously according to claim 1, characterized in that, The method further includes: Based on the distance range between the two shots during staggered firing in the adjacent work areas, and according to the defined shot point areas of the adjacent work areas, a fixed staggered firing interval is set for staggered firing, including: Based on the distance range between the two shots fired at staggered times in adjacent work areas, multiple fixed staggered time intervals are divided; The range of two-shot distances corresponding to each fixed staggered time interval is taken as the range of two-shot distances corresponding to the longest staggered firing time in each fixed staggered time interval. Based on the longest time-out excitation time in each fixed time-out interval, the multiple fixed time-out intervals are arranged in ascending order; Based on the distance between corresponding shot points in the shot point area, select the minimum fixed time interval corresponding to the two shot distance range that matches the fixed time interval, and perform time-staggered excitation.

8. The efficient collaborative construction method for multiple geophysical exploration teams operating simultaneously according to claim 1, characterized in that, When multiple geophysical exploration teams are working in adjacent work areas at the same time window, and the distance between their blast points is less than the minimum distance for synchronous excitation, the staggered excitation is carried out.

9. A highly efficient collaborative construction device for multiple geophysical exploration teams to work simultaneously, characterized in that, include: The integration module is used to integrate adjacent work areas to obtain the target work area; The analysis module is used to select the target historical seismic data acquisition area based on the target work area, use statistical methods to analyze the root mean square value of the background noise energy of the original single shot data of the target historical seismic data acquisition area, and obtain the lower limit value of the upper shot interference energy by removing outliers. The test module is used to conduct long-distance, long-arrange point tests in the target work area to obtain seismic data, and to draw a planar map of the spatiotemporal distribution of single-shot energy attenuation based on the target single-shot data in the seismic data. The extraction module is used to extract the contour line corresponding to the lower limit of the upper gun interference energy as the first upper gun interference time-distance curve based on the spatiotemporal distribution plan of the single gun energy attenuation; the first upper gun interference time-distance curve represents the minimum distance between the two guns that are triggered at different times. The calculation module is used to calculate the second upper shot interference time-distance curve, which indicates that the upper shot refracted wave does not enter the lower shot data volume, based on the lowest refracted wave velocity obtained from the seismic data; the second upper shot interference time-distance curve represents the maximum distance between the two shots that are triggered at different times. The region definition module is used to define the region enclosed by the first upper shot interference time-distance curve and the second upper shot interference time-distance curve as the distance range between the two shots that are staggered in the adjacent work area. The staggered firing module is used to set a dynamic staggered firing interval based on the distance range between the two shots in the adjacent work areas and the divided shot point areas of the adjacent work areas.

10. An electronic device, characterized in that, include: A processor, a memory, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the efficient collaborative construction method for simultaneous operation of multiple geophysical exploration teams as described in any one of claims 1 to 8.

Images