Seismic data acquisition method, device, equipment, storage medium and product
By optimizing the distribution grid of seabed nodes and encrypting the shot point grid, the high cost problem caused by the insufficient number of seabed nodes was solved, high-precision wide-azimuth or full-angle seismic data acquisition was achieved, and the acquisition cost was reduced.
Patent Information
- Application Number
- CN202111405283.4
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2021-11-24
- Publication Date
- 2025-09-09
- Estimated Expiration
- 2041-11-24
AI Technical Summary
In the prior art, in order to achieve wide-azimuth or full-azimuth seismic data acquisition, it is necessary to increase the number of seabed nodes, which leads to excessively high seismic data acquisition costs.
By determining the distribution grid area and distribution grid type of the seabed nodes, the layout of the seabed nodes is optimized, the number of nodes is reduced while achieving wide azimuth or omnidirectional seismic data acquisition, and an encrypted shot point grid is used to improve imaging accuracy.
Under the condition of a limited number of seabed nodes, the cost of seismic data acquisition is reduced, while the imaging accuracy and feasibility of seismic data are improved.
Smart Images

Figure CN116165699B_ABST
Abstract
Description
Technical Field
[0001] The present application relates to the field of geophysical exploration technology, and more particularly to a seismic data acquisition method, device, equipment, storage medium, and product. Background Art
[0002] The seafloor node seismic data acquisition method involves deploying seafloor nodes at regular intervals onto the seafloor and then collecting seismic data through shot-point firing. The seafloor nodes are the seismic instruments used to collect seismic data. This method can collect wide-azimuth (aspect ratio greater than 0.5 and less than 1) or even omnidirectional (aspect ratio of 1) seismic data. Processing and imaging this seismic data reveals subsurface structural features. Wide-azimuth or omnidirectional seismic data is particularly useful for imaging complex geological structures within seismic work areas.
[0003] In the related art, seismic data is collected by setting smaller receiving point spacing and receiving line spacing, where the receiving point spacing is the distance between two adjacent seabed nodes on a receiving line, and the receiving line spacing is the distance between two adjacent receiving lines.
[0004] However, for a seismic work area, the smaller the spacing between receiving points and receiving lines, and the greater the number of seabed nodes deployed when using a wide-azimuth or omnidirectional observation system to collect seismic data, the higher the cost of collecting seismic data. Summary of the Invention
[0005] The embodiments of the present application provide a seismic data acquisition method, apparatus, device, storage medium, and product that can reduce the cost of seismic data acquisition. The specific technical solution is as follows:
[0006] In one aspect, an embodiment of the present application provides a method for collecting seismic data, the method comprising:
[0007] determining a first number of seafloor nodes available for seismic data acquisition in a seismic work area;
[0008] Determining a distribution grid area of the seabed nodes based on the first number, the aspect ratio, the maximum longitudinal offset, the length of the receiving line, and the type of the observation system, where the distribution grid area is the product of a first distance between two adjacent seabed nodes on a receiving line and a second distance between two adjacent receiving lines, where the receiving line is a straight line formed by multiple seabed nodes along a seismic data acquisition direction, the aspect ratio is the ratio of the maximum lateral offset to the maximum longitudinal offset, and the aspect ratio is greater than a preset threshold, and both the maximum longitudinal offset and the maximum lateral offset are parameters representing the distance between a shot point and a seabed node;
[0009] Determining a distribution grid type of a seabed node, and determining the first distance and the second distance based on the distribution grid area and the distribution grid type;
[0010] determining a second number of receive lines based on the longitudinal maximum offset, the aspect ratio, the second distance, and the observation system type;
[0011] Determining a third distance and a fourth distance based on a bin size parameter, wherein the third distance is used to represent the distance between two adjacent shot points on a shot line, and the fourth distance is used to represent the distance between two adjacent shot lines. The bin size parameter is a parameter representing imaging accuracy when imaging seismic data;
[0012] Based on the first distance, the second distance and the second number, the first positions of multiple seabed nodes are determined, and based on the third distance and the fourth distance, the second positions of multiple shot points are determined; and seismic data are collected through an observation system consisting of each seabed node at the first position and each shot point at the second position.
[0013] In a possible implementation, determining the distribution grid area of the seabed nodes based on the first number, the aspect ratio, the maximum longitudinal offset, the length of the receiving line, and the type of the observation system includes:
[0014] Determining, based on the type of the observation system, a width parameter of the first number of seabed nodes;
[0015] Determining a third number of seabed nodes per unit area based on the first number, the aspect ratio, the maximum longitudinal offset, the length of the receiving line, and the width parameter;
[0016] The reciprocal of the third quantity is determined to obtain the distribution grid area.
[0017] In another possible implementation, determining a third number of seabed nodes per unit area based on the first number, the aspect ratio, the maximum longitudinal offset, the length of the receiving line, and the width parameter includes:
[0018] Determine the product of the aspect ratio, the maximum longitudinal offset, the length of the receiving line, and the width parameter to obtain the area of a first working area of the seismic working area;
[0019] The ratio of the first quantity to the area of the first work area is determined to obtain the third quantity.
[0020] In another possible implementation, determining the first distance and the second distance based on the distribution grid area and the distribution grid type includes:
[0021] If the distribution grid type is square, determining the square root of the area of the distribution grid to obtain the first distance and the second distance, and the first distance and the second distance are the same;
[0022] If the distribution grid type is rectangular, the first distance and the second distance are assigned values so that the product of the first distance and the second distance is equal to the distribution grid area.
[0023] In another possible implementation, the method further includes:
[0024] Determining a second work area of an infilled seismic work area and a third work area of a non-infilled seismic work area in the seismic work area, wherein the distribution density of the seabed nodes in the infilled seismic work area is greater than the distribution density of the seabed nodes in the non-infilled seismic work area;
[0025] Based on the area of the second work area, the area of the third work area, and the first number, a fifth distance, a sixth distance, a seventh distance, and an eighth distance are determined, wherein the fifth distance is used to represent the distance between two adjacent seabed nodes on a receiving line in the encrypted seismic work area, the sixth distance is used to represent the distance between two adjacent receiving lines in the encrypted seismic work area, the seventh distance is used to represent the distance between two adjacent seabed nodes on a receiving line in the non-encrypted seismic work area, and the eighth distance is used to represent the distance between two adjacent receiving lines in the non-encrypted seismic work area;
[0026] Determine a third position of a plurality of seabed nodes in the infilled seismic work area based on the fifth distance and the sixth distance; determine a fourth position of a plurality of seabed nodes in the non-infilled seismic work area based on the seventh distance and the eighth distance;
[0027] The acquisition of seismic data by an observation system composed of each seafloor node at the first position and each shot point at the second position includes:
[0028] Seismic data is collected through an observation system consisting of an ocean floor node at each fourth position in the non-encrypted seismic work area, an ocean floor node at each third position in the encrypted seismic work area, and a shot point at each second position.
[0029] In another possible implementation, determining the fifth distance, the sixth distance, the seventh distance, and the eighth distance based on the area of the second work area, the area of the third work area, and the first number includes:
[0030] Determine an encrypted grid type of the encrypted seismic work area and a non-encrypted grid type of the non-encrypted seismic work area;
[0031] Assigning a first encryption parameter and a second encryption parameter based on the encryption grid type;
[0032] Determine the non-encrypted grid area of the non-encrypted seismic work area based on the assigned first encryption parameter, the assigned second encryption parameter, the area of the second work area, the area of the third work area, and the first number, so that the sum of the fourth number of seabed nodes in the non-encrypted seismic work area and the fifth number of seabed nodes in the encrypted seismic work area is the first number, the fourth number is the ratio of the area of the third work area to the non-encrypted grid area, the fifth number is the ratio of the area of the second work area to the encrypted grid area, and the encrypted grid area is the product of the non-encrypted grid area, the assigned first encryption parameter, and the assigned second encryption parameter;
[0033] Determining the seventh distance and the eighth distance based on the non-encrypted grid type and the non-encrypted grid area;
[0034] determining a product of the seventh distance and the assigned first encryption parameter to obtain the fifth distance;
[0035] A product of the eighth distance and the assigned second encryption parameter is determined to obtain the sixth distance.
[0036] In another possible implementation, assigning the first encryption parameter and the second encryption parameter based on the encrypted grid type includes:
[0037] If the encrypted grid type is a square, assigning the same value to the first encryption parameter and the second encryption parameter;
[0038] If the encrypted grid type is rectangular, different values are assigned to the first encryption parameter and the second encryption parameter.
[0039] On the other hand, an embodiment of the present application provides a seismic data acquisition device, comprising:
[0040] A first determining module is used to determine a first number of seafloor nodes that can be used for seismic data acquisition in a seismic work area;
[0041] a second determining module, which determines a distribution grid area of the seafloor nodes based on the first number, the aspect ratio, the maximum longitudinal offset, the length of the receiving line, and the type of the observation system, wherein the distribution grid area is the product of a first distance between two adjacent seafloor nodes on a receiving line and a second distance between two adjacent receiving lines, wherein the receiving line is a straight line formed by multiple seafloor nodes along the seismic data acquisition direction, the aspect ratio is the ratio of the maximum lateral offset to the maximum longitudinal offset, and the aspect ratio is greater than a preset threshold, and the maximum longitudinal offset and the maximum lateral offset are both parameters representing the distance between the shot point and the seafloor node;
[0042] a third determining module, configured to determine a distribution grid type of the seabed node, and determine the first distance and the second distance based on the distribution grid area and the distribution grid type;
[0043] a fourth determining module, configured to determine a second number of receiving lines based on the maximum longitudinal offset, the aspect ratio, the second distance, and the observation system type;
[0044] a fifth determining module, configured to determine a third distance and a fourth distance based on a bin size parameter, wherein the third distance represents the distance between two adjacent shot points on a shot line, and the fourth distance represents the distance between two adjacent shot lines, and the bin size parameter represents the imaging accuracy of seismic data imaging;
[0045] An acquisition module is configured to determine a first position of a plurality of seafloor nodes based on the first distance, the second distance, and the second number, and to determine a second position of a plurality of shot points based on the third distance and the fourth distance; and to acquire seismic data through an observation system consisting of each seafloor node at the first position and each shot point at the second position.
[0046] In one possible implementation, the second determination module is used to determine a width parameter for laying the first number of seabed nodes based on the observation system type; determine a third number of seabed nodes per unit area based on the first number, the aspect ratio, the maximum longitudinal offset, the length of the receiving line, and the width parameter; and determine the reciprocal of the third number to obtain the distribution grid area.
[0047] In another possible implementation, the second determination module is used to determine the product of the aspect ratio, the maximum longitudinal offset, the length of the receiving line and the width parameter to obtain the first work area of the seismic work area; and determine the ratio of the first number and the first work area to obtain the third number.
[0048] In another possible implementation, the third determination module is used to determine the square root of the distribution grid area if the distribution grid type is a square, to obtain the first distance and the second distance, and the first distance and the second distance are the same; if the distribution grid type is a rectangle, to assign values to the first distance and the second distance so that the product of the first distance and the second distance is equal to the distribution grid area.
[0049] In another possible implementation, the apparatus further includes:
[0050] A seventh determining module is configured to determine a second work area of the infilled seismic work area and a third work area of the non-infilled seismic work area in the seismic work area, wherein the distribution density of the seabed nodes in the infilled seismic work area is greater than the distribution density of the seabed nodes in the non-infilled seismic work area;
[0051] an eighth determination module, configured to determine, based on the area of the second work area, the area of the third work area, and the first number, a fifth distance, a sixth distance, a seventh distance, and an eighth distance, wherein the fifth distance is used to represent the distance between two adjacent seabed nodes on a receiving line in the infilled seismic work area, the sixth distance is used to represent the distance between two adjacent receiving lines in the infilled seismic work area, the seventh distance is used to represent the distance between two adjacent seabed nodes on a receiving line in the non-infilled seismic work area, and the eighth distance is used to represent the distance between two adjacent receiving lines in the non-infilled seismic work area;
[0052] a ninth determining module, configured to determine a third position of the plurality of seafloor nodes in the infilled seismic work area based on the fifth distance and the sixth distance; and to determine a fourth position of the plurality of seafloor nodes in the non-infilled seismic work area based on the seventh distance and the eighth distance;
[0053] The acquisition module is used to collect seismic data through an observation system consisting of an ocean floor node at each fourth position in the non-encrypted seismic work area, an ocean floor node at each third position in the encrypted seismic work area, and a shot point at each second position.
[0054] In another possible implementation, the eighth determination module is used to determine the encrypted grid type of the encrypted seismic work area and the non-encrypted grid type of the non-encrypted seismic work area; based on the encrypted grid type, assign a first encryption parameter and a second encryption parameter; based on the assigned first encryption parameter, the assigned second encryption parameter, the second work area area, the third work area area and the first number, determine the non-encrypted grid area of the non-encrypted seismic work area, so that the sum of the fourth number of seabed nodes in the non-encrypted seismic work area and the fifth number of seabed nodes in the encrypted seismic work area is the first number , the fourth quantity is the ratio of the third work area area to the non-encrypted grid area, the fifth quantity is the ratio of the second work area area to the encrypted grid area, the encrypted grid area is the product of the non-encrypted grid area, the assigned first encryption parameter and the assigned second encryption parameter; based on the non-encrypted grid type and the non-encrypted grid area, the seventh distance and the eighth distance are determined; the fifth distance is obtained by determining the product of the seventh distance and the assigned first encryption parameter; the sixth distance is obtained by determining the product of the eighth distance and the assigned second encryption parameter.
[0055] In another possible implementation, the eighth determination module is configured to assign the same value to the first encryption parameter and the second encryption parameter if the encrypted grid type is a square; and to assign different values to the first encryption parameter and the second encryption parameter if the encrypted grid type is a rectangle.
[0056] On the other hand, an embodiment of the present application provides a computer device, which includes a processor and a memory, wherein the memory stores at least one program code, and the at least one program code is loaded and executed by the processor to implement the operations performed in the seismic data acquisition method described in the embodiment of the present application.
[0057] On the other hand, an embodiment of the present application provides a computer-readable storage medium, which stores at least one program code, and the at least one program code is loaded and executed by a processor to implement the operations performed in the seismic data acquisition method described in the embodiment of the present application.
[0058] On the other hand, an embodiment of the present application provides a computer program product, which stores at least one program code, and the at least one program code is loaded and executed by a processor to implement the operations performed in the seismic data acquisition method described in the embodiment of the present application.
[0059] The beneficial effects of the technical solution provided by the embodiments of the present application are:
[0060] An embodiment of the present application provides a seismic data acquisition method. Under the condition that the number of seabed nodes is limited, the method realizes the acquisition of wide-azimuth or full-azimuth seismic data by determining a first distance between two adjacent seabed nodes on a receiving line, a second distance between two adjacent receiving lines, and a second number of receiving lines, thereby reducing the cost of seismic data acquisition. BRIEF DESCRIPTION OF THE DRAWINGS
[0061] Figure 1 This is a flow chart of a seismic data acquisition method provided by an embodiment of the present application;
[0062] Figure 2 This is a flow chart for determining the distribution of seabed nodes in an infill seismic work area provided by an embodiment of the present application;
[0063] Figure 3 This is a schematic diagram showing that the grid size of the seabed nodes in the infilled seismic work area is inconsistent with the grid size of the seabed nodes in the non-infilled seismic work area provided in an embodiment of the present application;
[0064] Figure 4 This is a schematic diagram of an embodiment of the present application providing a method of collecting omnidirectional seismic data using a bilateral blasting beam-shaped observation system;
[0065] Figure 5 This embodiment of the present application provides a Figure 4 A partial enlarged view of
[0066] Figure 6 This embodiment of the present application provides a Figure 4 The rose diagram corresponding to the longitudinal offset limit of 5 km for the observing system in ;
[0067] Figure 7 This embodiment of the present application provides a Figure 4 Schematic diagram of the corresponding full coverage times after the longitudinal offset of the observation system is limited to 5 km;
[0068] Figure 8 This is a schematic diagram of another embodiment of the present application for collecting omnidirectional seismic data using a bilateral blasting beam observation system;
[0069] Figure 9 This embodiment of the present application provides a Figure 8 A partial enlarged view of
[0070] Figure 10 This embodiment of the present application provides a Figure 8 The rose diagram corresponding to the longitudinal offset limit of 5 km for the observing system in ;
[0071] Figure 11 This embodiment of the present application provides a Figure 8 Schematic diagram of the corresponding full coverage times after the longitudinal offset of the observation system is limited to 5 km;
[0072] Figure 12 This is a structural diagram of a seismic data acquisition device provided in an embodiment of the present application;
[0073] Figure 13 This is a structural block diagram of a computer device provided in an embodiment of the present application. DETAILED DESCRIPTION
[0074] In order to make the technical solutions and advantages of the present application clearer, the implementation methods of the present application are described in further detail below.
[0075] This application proposes a wide (full) azimuth seabed node seismic data acquisition method that matches the number of seabed nodes based on the characteristics of single-point acquisition of seabed node seismic data. The distance between two adjacent seabed nodes on a receiving line changes flexibly. This application mainly solves the imaging problem of complex underground geological bodies by appropriately increasing the distribution grid of seabed nodes and encrypting the shot point grid when necessary to achieve wide (full) azimuth and high-precision seismic data acquisition of seabed nodes. It saves cost investment by reducing the amount of seabed node investment and improves the feasibility of wide (full) azimuth seismic data acquisition of seabed nodes. The specific process is as follows:
[0076] The present invention provides a method for collecting seismic data, which is executed by a computer device. Figure 1 , the method comprising:
[0077] Step 101: A computer device determines a first number of seafloor nodes available for seismic data acquisition in a seismic work area.
[0078] The computer device may first determine the total number of submarine nodes expected to be deployed in the earthquake work area and the spare number of spare submarine nodes, and then determine the difference between the total number and the spare number to obtain the first number.
[0079] The number of backup submarine nodes can be determined based on the backup ratio. The process is as follows: a computer determines the product of the total number of submarine nodes and the backup ratio to obtain the backup number. The backup ratio is the percentage of backup submarine nodes to the total number of submarine nodes. For example, if N1 represents the first number, c represents the backup ratio, and N represents the total number, then N1 = N × (1-c).
[0080] It should be noted that the seabed node in the embodiment of the present application can be a seabed node that has been deployed on the seabed, or a seabed node that has not been deployed on the seabed. If the seabed node is a seabed node that has been deployed on the seabed, the method provided in the embodiment of the present application can adjust the position of the seabed node that has been deployed in the seismic work area, thereby realizing the collection of wide azimuth or full-angle seismic data. If the seabed node is a seabed node that has not been deployed on the seabed, the method provided in the embodiment of the present application can determine the position of each seabed node in the seismic work area, thereby realizing the collection of wide azimuth or full-angle seismic data.
[0081] Step 102: The computer device determines width parameters of the first number of seabed nodes based on the type of the observation system.
[0082] The width parameter is used to determine the maximum area corresponding to the first number of submarine nodes. The computer device can determine the width parameter corresponding to the first number of submarine nodes corresponding to the observation system type based on a pre-established correspondence between the observation system type and the width parameter of the submarine node layout.
[0083] Different observation system types correspond to different width parameters. When the observation system type is a linear intermediate shot mode, the width parameter is 2. When the observation system type is a linear bilateral shot mode, the width parameter is 1. When the observation system type is an orthogonal large cross observation system, the width parameter is the inverse of the shot repetition rate. Other special observation system types are determined based on their characteristics. The width parameter can be an integer or a decimal, and this is not specifically limited in the embodiments of this application.
[0084] Step 103: The computer device determines a third number of seabed nodes within a unit area based on the first number, the aspect ratio, the maximum longitudinal offset, and the length and width parameters of the receiving line.
[0085] The receiving line is a straight line composed of multiple seabed nodes along the direction of seismic data acquisition. The aspect ratio is the ratio of the maximum horizontal offset to the maximum vertical offset, and the aspect ratio is greater than or equal to the preset threshold. The maximum horizontal offset and the maximum vertical offset are both parameters representing the distance between the shot point and the seabed node.
[0086] If the seismic work area is regular and the corresponding area is rectangular, the computer device determines the aspect ratio, the maximum longitudinal offset, the product of the length of the receiving line and the width parameter to obtain the first work area area of the seismic work area. Among them, the product of the aspect ratio, the maximum longitudinal offset and the width parameter is used to represent the width of the first number of seabed nodes laid, and the length of the receiving line is used to represent the length of the first number of seabed nodes laid. The product of the width and length of the seabed nodes is the first work area area of the seismic work area. If the seismic work area is irregular, the first work area area can be obtained by measurement. The computer device then determines the ratio of the first number to the first work area area to obtain the third number.
[0087] For example, r represents the aspect ratio, L1 represents the maximum longitudinal offset, L2 represents the length of the receiving line, n represents the width parameter, S1 represents the area of the first working area, and N2 represents the third number. Then the area of the first working area S1 = L2 × (r × L1 × n), and the third number N2 = N1 / S1.
[0088] It should be noted that prior to step 103, the computer device first determines the aspect ratio, maximum longitudinal offset, and length of the receiving line. Then, in this step, the third quantity is determined based on the determined aspect ratio, maximum longitudinal offset, and length of the receiving line. The aspect ratio, maximum longitudinal offset, and length of the receiving line are all adjustable. For example, when the aspect ratio is greater than 0.5 and less than 1, seismic data is acquired over a wide azimuth angle. When the aspect ratio is equal to 1, seismic data is acquired over an omnidirectional angle. In other words, if wide azimuth seismic data acquisition is desired, the aspect ratio is determined to be any value between 0.5 and 1. If omnidirectional seismic data acquisition is desired, the aspect ratio is determined to be 1.
[0089] Step 104: The computer device determines the inverse of the third quantity to obtain the distribution grid area.
[0090] The distribution grid area is the product of the first distance between two adjacent seabed nodes on a receiving line and the second distance between two adjacent receiving lines. For example, if S2 represents the distribution grid area, then the distribution grid area S2 = 1 / N2.
[0091] Step 105: The computer device determines the first distance and the second distance based on the distribution grid area and the distribution grid type.
[0092] According to step 104 , the first distance is the distance between two adjacent seabed nodes on a receiving line, and the second distance is the distance between two adjacent receiving lines.
[0093] If the distribution grid type is square, the first distance and the second distance are the same, and the computer device determines the square root of the distribution grid area to obtain the first distance and the second distance. For example, if X1 represents the first distance and X2 represents the second distance, then
[0094] If the distribution grid type is rectangular, the computer device may assign values to the first distance and the second distance so that the product of the first distance and the second distance equals the distribution grid area. For example, if the value assigned to the first distance X1 is x, then the second distance X2 = S2 / x.
[0095] Alternatively, if the distribution grid type is rectangular, the computer device may determine multiple distance combinations in which the product of the first distance and the second distance is the area of the distribution grid, wherein each distance combination is a combination of the first distance and the second distance, and the product of the first distance and the second distance in each distance combination is the area of the distribution grid. The computer device determines a selected distance combination from the multiple distance combinations and obtains the first distance and the second distance in the selected distance combination.
[0096] It should be noted that the first and second distances determined using the above method may be integers or decimals. If the first and second distances are decimals, the computer device may round them down or up to adjust them to integers. Furthermore, rounding down may involve adding 1 to the original distance or increasing it by a value greater than 1, provided the increased value does not exceed a preset value. Similarly, rounding up may involve subtracting 1 from the original distance or decreasing it by a value greater than 1, provided the decreased value does not exceed a preset value. For example, if the first distance is 197m, rounding down may adjust the first distance to 200m. For another example, if the first distance is 202m, rounding up may adjust the first distance to 200m. Furthermore, the adjusted first distance does not necessarily represent an integer multiple of the bin size parameter along a shot line, and the adjusted second distance does not necessarily represent an integer multiple of the bin size parameter along a perpendicular shot line.
[0097] Step 106: The computer device determines the second number of receiving lines based on the maximum longitudinal offset, the aspect ratio, the second distance, and the width parameter.
[0098] The computer determines the product of the maximum longitudinal offset and the aspect ratio to obtain the maximum transverse offset. It then determines the ratio of the maximum transverse offset to the second distance, which represents the number of receiving lines within the span width of the maximum transverse offset. The computer then multiplies this ratio by the width parameter to obtain the second number. For example, if L3 represents the maximum transverse offset and N3 represents the second number, then L3 = L1 × r, and N3 = (L3 / X2) × n.
[0099] If the product of the ratio and the width parameter determined by the computer device is an integer, the product is directly determined as the second number of receiving lines. If the product determined by the computer device is not an integer, the product may be rounded up or down, and the rounded-up or rounded-down value is used as the second number of receiving lines.
[0100] In the related art, the distance between two adjacent seabed nodes on a receiving line is 25m or 50m. According to the principle of no block construction or less block construction according to the size of the work area, to achieve wide (full) azimuth seismic data acquisition, it is often necessary to invest a large number of seabed nodes, resulting in high cost of seismic data acquisition. The method provided in the embodiment of the present application can solve the problem that the wide (full) azimuth seismic data acquisition of seabed nodes cannot be achieved by adopting the traditional observation system design method due to the constraints of equipment investment and cost factors. This method can achieve wide (full) azimuth seismic data acquisition under the condition of limited seabed nodes. The number of seabed nodes invested is determined according to the size of the work area, and the geological requirements and total cost investment are comprehensively considered to achieve matching with the number of seabed nodes and optimize cost efficiency, thereby improving the feasibility of wide (full) azimuth seismic data acquisition of seabed nodes and reducing operating costs.
[0101] Step 107: The computer device determines the third distance and the fourth distance based on the bin size parameter.
[0102] The third distance is used to represent the distance between two adjacent shot points on a shot line, the fourth distance is used to represent the distance between two adjacent shot lines, and the bin size parameter is a parameter representing the imaging accuracy when imaging seismic data.
[0103] The commonly used bin size parameters in high-precision imaging are 12.5 meters * 12.5 meters or 25 meters * 25 meters. The grid size of the shot point is twice the bin size parameter. That is, if the bin size parameter is 25 meters * 25 meters, the grid size of the shot point is 50 meters * 50 meters, that is, the third distance X3 and the fourth distance X4 are both 50 meters.
[0104] Step 108: The computer device determines first positions of the plurality of seabed nodes based on the first distance, the second distance, and the second quantity, and determines second positions of the plurality of shot points based on the third distance and the fourth distance.
[0105] The computer device determines the positions of the second number of receiving lines based on the second distance. The computer device determines the ratio of the length of the receiving line to the first distance to obtain a sixth number of submarine nodes on each receiving line. The computer device then determines the first positions of the sixth number of submarine nodes on each receiving line based on the first distance.
[0106] In this embodiment of the present application, the product of the sixth number of seafloor nodes on each receiving line and the second number of receiving lines is the number of seafloor nodes acquiring seismic data in the seismic work area, and the difference between this number and the first number is within a preset difference range. This number can be greater than, less than, or equal to the first number, as long as the difference between the two is within the preset difference range. This preset difference range can be set and modified as needed and is not specifically limited in this embodiment of the present application.
[0107] The computer device determines the position of each gun line according to the fourth distance, and then determines the second position of each gun point on each gun line according to the third distance.
[0108] In an embodiment of the present application, after the computer device determines the first position of each seafloor node and the second position of each shot point, step 109 can be directly executed. Alternatively, the computer device can also determine a coverage count based on the first position of each seafloor node and the second position of each shot point. The coverage count represents the number of repeated observations of the same location on the subsurface interface. If the coverage count is greater than a preset coverage count, step 109 is executed. If the coverage count is not greater than the preset coverage count, the length of the receiving line is re-determined. The re-determined receiving line length is less than the previously determined receiving line length. Based on the re-determined receiving line length, steps 103 to 108 are re-executed, ultimately achieving a uniform distribution grid that matches the number of seafloor nodes. Alternatively, the computer device can appropriately increase the shot point density to increase the coverage count. Furthermore, if the computer device determines that the coverage count is greater than the preset coverage count, it can also perform attribute analysis based on the first position of each seafloor node and the second position of each shot point, and then conduct a comprehensive assessment based on geological requirements and economic investment. If the geological requirements are met and the economic investment is within the budget, step 109 is executed.
[0109] Step 109: The computer device collects seismic data through an observation system composed of each seabed node at the first position and each shot point at the second position.
[0110] The operator can carry out on-site construction based on the observation system composed of the first position of each seabed node and the second position of each shot point determined above, and then collect seismic data through the seabed nodes through shot point excitation, and the computer equipment obtains the seismic data collected by the seabed nodes. The seismic data is wide-azimuth or full-angle seismic data.
[0111] An embodiment of the present application provides a seismic data acquisition method. Under the condition that the number of seabed nodes is limited, the method realizes the acquisition of wide-azimuth or full-azimuth seismic data by determining a first distance between two adjacent seabed nodes on a receiving line, a second distance between two adjacent receiving lines, and a second number of receiving lines, thereby reducing the cost of seismic data acquisition.
[0112] In the embodiment of the present application, the distribution grid determined by the above steps 101 to 106 is consistent, that is, the number of adjacent two seabed nodes on each receiving line in the seismic work area is the same, and the distance between each two adjacent receiving lines is the same. However, in the actual exploration process, there may be areas that need to be explored in a focused manner. Therefore, it is necessary to increase the distribution density of seabed nodes in the focused exploration area. Accordingly, the computer equipment can determine the position of each seabed node in the focused exploration area and the position of each seabed node in the non-focus exploration area through the following steps 201 to 207 under the condition that the number of seabed nodes is limited, see Figure 2 , the method comprising:
[0113] Step 201: The computer device determines the second work area of the infilled seismic work area and the third work area of the non-infilled seismic work area in the seismic work area.
[0114] The distribution density of seabed nodes in the intensified seismic area is greater than that in the non-intensified seismic area.
[0115] The computer device first determines the intensified seismic work area and the non-intensified seismic work area within the seismic work area, and then determines the area of the intensified seismic work area and the area of the non-intensified seismic work area. If the intensified seismic work area is regular and the corresponding area is rectangular, the length and width of the intensified seismic work area can be determined, and then the product of the length and width of the intensified seismic work area can be determined to obtain the area of the second intensified seismic work area. If the intensified seismic work area is regular and the corresponding area is circular, the exploration radius of the intensified seismic work area can be determined, and the area of the second intensified seismic work area can be determined based on the exploration radius. If the intensified seismic work area is irregular, the area of the second intensified seismic work area can be determined through measurement.
[0116] The computer determines the difference between the area of the first seismic work area and the area of the second intensified seismic work area to obtain the area of the third seismic work area. For example, if S3 represents the area of the second seismic work area and S4 represents the area of the third seismic work area, then after S3 is determined, S4 = S1 - S3.
[0117] The process of determining the area of the first work area is the same as the process of determining the area of the first work area in step 103, and will not be repeated here.
[0118] Step 202: The computer device determines a fifth distance, a sixth distance, a seventh distance, and an eighth distance based on the second work area area, the third work area area, and the first quantity.
[0119] This step can be achieved by following the steps (1) to (6), including:
[0120] (1) Computer equipment determines the encrypted grid type of the encrypted seismic work area and the non-encrypted grid type of the non-encrypted seismic work area.
[0121] The encrypted grid type of the encrypted seismic work area can be square or rectangular. The non-encrypted grid type of the non-encrypted seismic work area can be square or rectangular. In other words, the encrypted grid type and the non-encrypted grid type can be the same or different. This is not specifically limited in the embodiments of the present application.
[0122] (2) The computer device assigns a first encryption parameter and a second encryption parameter based on the encryption grid type.
[0123] If the encrypted grid type is a square, the computer device assigns the same value to the first encryption parameter and the second encryption parameter; if the encrypted grid type is a rectangle, the computer device assigns different values to the first encryption parameter and the second encryption parameter. The first encryption parameter and the second encryption parameter are both values less than 1.
[0124] (3) The computer device determines the non-encrypted grid area of the non-encrypted seismic work area based on the assigned first encryption parameter, the assigned second encryption parameter, the second work area area, the third work area area and the first number, so that the sum of the fourth number of seabed nodes in the non-encrypted seismic work area and the fifth number of seabed nodes in the encrypted seismic work area is the first number.
[0125] Among them, the fourth quantity is the ratio of the area of the third work area to the area of the non-encrypted grid, the fifth quantity is the ratio of the area of the second work area to the area of the encrypted grid, and the encrypted grid area is the product of the non-encrypted grid area, the assigned first encryption parameter, and the assigned second encryption parameter.
[0126] In step (3), the area of the non-refined grid can be determined by the following formula:
[0127] S4 / S5+S3 / (a×b×S5)=N1.
[0128] Among them, S5 represents the non-encrypted grid area, a represents the first encryption parameter, b represents the second encryption parameter, S4 / S5 is the fourth quantity, and S3 / (a×b×S5) is the fifth quantity.
[0129] In this formula, only S5 is a variable, and the rest are constants. Therefore, S5 can be determined by other constants.
[0130] (4) The computer device determines the seventh distance and the eighth distance based on the non-encrypted grid type and the non-encrypted grid area.
[0131] If the non-encrypted grid type is square, the seventh distance is equal to the eighth distance, and the computer device determines the square root of the area of the non-encrypted grid to obtain the seventh and eighth distances. If the non-encrypted grid type is rectangular, the computer device assigns values to the seventh and eighth distances such that the product of the seventh and eighth distances equals the area of the non-encrypted grid.
[0132] The process of the computer device assigning values to the seventh distance and the eighth distance is similar to the process of assigning values to the first distance and the second distance, and will not be repeated here.
[0133] (5) The computer device determines the product of the seventh distance and the assigned first encryption parameter to obtain the fifth distance.
[0134] For example, X5 represents the fifth distance and X7 represents the seventh distance, then X5=X7×a.
[0135] (6) The computer device determines the product of the eighth distance and the assigned second encryption parameter to obtain the sixth distance.
[0136] For example, X6 represents the sixth distance and X8 represents the eighth distance, then X6=X8×b.
[0137] See also Figure 3 , Figure 3 This is a schematic diagram showing that the grid size of the seabed nodes in the dense seismic area is different from that in the non-densified seismic area. Figure 3 It can be seen from the figure that the grid size of the seabed nodes in the non-encrypted seismic area is 600m×600m, that is, the distance between two adjacent seabed nodes on a receiving line and the distance between two adjacent receiving lines in the non-encrypted seismic area are both 600m. The grid size of the seabed nodes in the encrypted seismic area is 300m×300m, that is, the distance between two adjacent seabed nodes on a receiving line and the distance between two adjacent receiving lines in the encrypted seismic area are both 300m.
[0138] Step 203: The computer device determines the third positions of the plurality of seabed nodes in the dense seismic work area based on the fifth distance and the sixth distance.
[0139] The computer equipment determines the position of each receiving line in the encrypted seismic work area according to the sixth distance, and determines the third position of each seabed node at the position of each receiving line in the encrypted seismic work area according to the fifth distance.
[0140] Step 204: The computer device determines the fourth positions of the plurality of seabed nodes in the non-reinforced seismic work area based on the seventh distance and the eighth distance.
[0141] The computer equipment determines the position of each receiving line in the non-encrypted seismic work area according to the eighth distance, and determines the fourth position of each seabed node on the position of each receiving line in the non-encrypted seismic work area according to the seventh distance.
[0142] Step 205: The computer device determines the third distance and the fourth distance based on the bin size parameter.
[0143] This step is the same as step 107 and will not be repeated here.
[0144] Step 206: The computer device determines second positions of the plurality of shot points based on the third distance and the fourth distance.
[0145] This step is the same as the process of the computer device determining the second positions of the plurality of shot points in step 108, and will not be described in detail here.
[0146] Step 207: The computer device collects seismic data through an observation system consisting of an ocean floor node at each fourth position in the non-encrypted seismic work area, an ocean floor node at each third position in the encrypted seismic work area, and a shot point at each second position.
[0147] Operators can carry out on-site construction based on the observation system composed of the third position of each seabed node in the encrypted seismic work area, the fourth position of each seabed node in the non-encrypted seismic work area, and the second position of each shot point determined above, and then collect seismic data through the seabed nodes through the shot point excitation, and the computer equipment obtains the seismic data collected by the seabed nodes.
[0148] In an embodiment of the present application, under the condition that the number of seabed nodes in a seismic work area is limited, wide-azimuth or full-angle seismic data acquisition in the key exploration area is achieved by determining the fifth distance between two adjacent seabed nodes on a receiving line in the key exploration area and the sixth distance between two adjacent receiving lines. Wide-azimuth or full-angle seismic data acquisition in the non-key exploration area is achieved by determining the seventh distance between two adjacent seabed nodes on a receiving line in the non-key exploration area and the eighth distance between two adjacent receiving lines, thereby reducing the cost of seismic data acquisition.
[0149] Example 1
[0150] In this example, the total number of seafloor nodes expected to be deployed in the seismic work area is 6,000, the seafloor node backup rate is 20%, the maximum longitudinal offset of the seismic work area is 5 km, the bin size parameter is 25 meters * 25 meters, an omnidirectional design is adopted, the aspect ratio is 1, the maximum length of the receiving line is 38 km, the observation system type is a bilateral shot beam observation system, and the distribution grid type is a square. The specific steps are as follows:
[0151] (1) Determine a first number N1 of seafloor nodes that can be used for seismic data acquisition in the seismic work area.
[0152] N1=6000×(1-20%)=4800.
[0153] (2) Determine the width parameter n of the bilateral blasting beam observation system to be 1.
[0154] (3) Determine the product of the aspect ratio r, the maximum longitudinal offset L1, the length L2 of the receiving line, and the width parameter n to obtain the first work area S1 of the seismic work area.
[0155] S1=1×5×38×1=190.
[0156] (4) Determine the ratio of the first number N1 to the area S1 of the first work area to obtain the third number N2 of seabed nodes per unit area.
[0157] N2=4800 / 190≈25.26.
[0158] (5) Determine the reciprocal of the third number N2 to obtain the distribution grid area S2.
[0159] S2=1 / 25.26=0.0396.
[0160] (6) Determine the square root of the distribution grid area S2 to obtain the first distance X1. The first distance is the same as the second distance X2.
[0161]
[0162] That is, when a square distribution grid is used, the distance between two adjacent seabed nodes on a receiving line is equal to the distance between two adjacent receiving lines.
[0163] The computer device adjusts the first distance and the second distance, rounds down the first distance and the second distance, and the adjusted first distance and the second distance are 200 meters.
[0164] (7) Determine the product of the maximum longitudinal offset L1 and the aspect ratio r to obtain the maximum lateral offset L3.
[0165] L3=5×1=5km=5000m.
[0166] (8) Determine the ratio of the maximum lateral offset distance L3 to the second distance X2, and obtain the number N4 of receiving lines corresponding to the maximum lateral offset distance within the spread width.
[0167] N4=5000 / 200=25 items.
[0168] (9) Determine the product of the number N3 of receiving lines corresponding to the maximum lateral offset within the spread width and the width parameter n to obtain a second number N3.
[0169] N3=25×1=25 pieces.
[0170] (10) According to the bin size parameter 25m×25m, the third distance X3 and the fourth distance X4 are determined.
[0171] X3 = X4 = 2 × 25 = 50 m. That is, the distance between two adjacent gun points on a gun line and the distance between two adjacent gun lines are both 50 m.
[0172] (11) Perform attribute analysis on the square grid observation system determined according to the above steps.
[0173] See also Figure 4 , Figure 4 The double-sided shooting beam observation system is used to collect all-directional seismic data. Figure 4 There are 25 receiving lines in total, each with a length of 38 km. The number of submarine nodes in each receiving line is 190, for a total of 4,750 submarine nodes.
[0174] See also Figure 5 , Figure 5 for Figure 4 A partial enlarged view of Figure 5 It can be seen from the figure that the distribution grid of the seabed nodes is 200m×200 meters, that is, the distance between two adjacent seabed nodes on a receiving line and the distance between two adjacent receiving lines are both 200m, and the grid of the shot points is 50m×50m, that is, the distance between two adjacent shot points on a shot line and the distance between two adjacent shot lines are both 50m.
[0175] See also Figure 6 , Figure 6 for Figure 4 The rose diagram corresponding to the observation system longitudinal offset limit of 5km is shown in the figure. Figure 6 The offset in is the distance between the shot point and the seabed node, Figure 6 It can be seen from the figure that the rose diagram is circular, indicating that the observation system can realize omnidirectional seismic data acquisition.
[0176] See also Figure 7 , Figure 7 for Figure 4 The number of full coverage times corresponding to the longitudinal offset limit of 5 km for the observation system in Figure 7 It can be seen that the full coverage number is 625 times, which is relatively high and can meet the requirements.
[0177] Example 2
[0178] In this embodiment, the distribution grid type is rectangular, and parameters such as the total number of seabed nodes, spare rate, maximum longitudinal offset, bin size parameters, aspect ratio, maximum extent of receiving lines, and observation system type are the same as those in Example 1 and are not repeated here.
[0179] Since the above parameters are the same as those in Example 1, steps (1) to (5) in this embodiment are also the same as steps (1) to (5) in Example 1, which are all 0.0396, and will not be repeated here.
[0180] (6) Assign a value to the first distance, determine the ratio of the distribution grid area to the assigned first distance, and obtain the second distance.
[0181] For example, if the first distance is 100m, the second distance is rounded down to 400m. For another example, if the first distance is 133m, the second distance is rounded down to 300m. In this embodiment, only the first distance of 133m and the second distance of 300m are used as an example for description.
[0182] (7) Determine the product of the maximum longitudinal offset L1 and the aspect ratio r to obtain the maximum lateral offset L3.
[0183] Step (7) is the same as step (7) in Example 1 and will not be repeated here.
[0184] (8) Determine the ratio of the maximum lateral offset distance L3 to the second distance X2, and obtain the number N3 of receiving lines corresponding to the maximum lateral offset distance within the spread width.
[0185] N3=5000÷300≈17 pieces.
[0186] (9) Determine the product of the number N3 of receiving lines corresponding to the maximum lateral offset within the spread width and the width parameter n to obtain a second number N4.
[0187] Step (9) is the same as step (9) in Example 1 and will not be repeated here.
[0188] (10) According to the bin size parameter 25m×25m, the third distance X3 and the fourth distance X4 are determined.
[0189] Step (10) is the same as step (10) in Example 1 and will not be repeated here.
[0190] (11) Perform attribute analysis on the rectangular grid observation system determined according to the above steps.
[0191] See also Figure 8 , Figure 8 The double-sided shooting beam observation system is used to collect all-directional seismic data. Figure 4There are 17 receiving lines in total, each of which is 38 km long. There are 286 submarine nodes in each receiving line, for a total of 4,862 submarine nodes.
[0192] See also Figure 9 , Figure 9 for Figure 8 A partial enlarged view of Figure 9 It can be seen from the figure that the distribution grid of the seabed nodes is 133m×300m, that is, the distance between two adjacent seabed nodes on a receiving line is 133m, and the distance between two adjacent receiving lines is 300m. The grid of the shot points is 50m×50m, that is, the distance between two adjacent shot points on a shot line and the distance between two adjacent shot lines are both 50m.
[0193] See also Figure 10 , Figure 10 for Figure 8 The corresponding rose diagram after the longitudinal offset of the observation system is limited to 5km is shown in the figure. Figure 10 It can be seen from the figure that the rose diagram is circular, indicating that the observation system can realize omnidirectional seismic data acquisition.
[0194] See also Figure 11 , Figure 11 for Figure 8 The number of full coverage times corresponding to the longitudinal offset limit of 5 km for the observation system in Figure 11 It can be seen that the full coverage times are 629 and 646, which are relatively high and distributed in strips, which can meet the requirements.
[0195] The present application provides a seismic data acquisition device. Figure 12 , the device comprises:
[0196] A first determining module 1201 is configured to determine a first number of seafloor nodes capable of acquiring seismic data in a seismic work area;
[0197] A second determining module 1202 is configured to determine a distribution grid area of the seafloor nodes based on the first number, the aspect ratio, the maximum longitudinal offset, the length of the receiving line, and the type of the observation system, where the distribution grid area is the product of a first distance between two adjacent seafloor nodes on a receiving line and a second distance between two adjacent receiving lines. A receiving line is a straight line formed by multiple seafloor nodes along the seismic data acquisition direction. The aspect ratio is the ratio of the maximum transverse offset to the maximum longitudinal offset, and the aspect ratio is greater than a preset threshold. The maximum longitudinal offset and the maximum transverse offset are both parameters representing the distance between a shot point and a seafloor node.
[0198] The third determining module 1203 is configured to determine a distribution grid type of the seabed node, and determine a first distance and a second distance based on the distribution grid area and the distribution grid type;
[0199] A fourth determining module 1204 is configured to determine a second number of receiving lines based on the maximum longitudinal offset, the aspect ratio, the second distance, and the observation system type;
[0200] A fifth determining module 1205 is configured to determine a third distance and a fourth distance based on a bin size parameter, wherein the third distance represents the distance between two adjacent shot points on a shot line, and the fourth distance represents the distance between two adjacent shot lines. The bin size parameter represents the imaging accuracy of seismic data imaging.
[0201] Acquisition module 1206 is configured to determine a first position of a plurality of seafloor nodes based on the first distance, the second distance, and the second quantity, and to determine a second position of a plurality of shot points based on the third distance and the fourth distance; and to acquire seismic data through an observation system consisting of each seafloor node at the first position and each shot point at the second position.
[0202] In one possible implementation, the second determination module 1202 is used to determine the width parameter of the first number of seabed nodes based on the observation system type; determine the third number of seabed nodes per unit area based on the first number, the aspect ratio, the maximum longitudinal offset, the length of the receiving line and the width parameter; and determine the reciprocal of the third number to obtain the distribution grid area.
[0203] In another possible implementation, the second determination module 1202 is used to determine the product of the aspect ratio, the maximum longitudinal offset, and the length and width parameters of the receiving line to obtain the first work area of the seismic work area; and determine the ratio of the first quantity and the first work area to obtain the third quantity.
[0204] In another possible implementation, the third determination module 1203 is used to determine the square root of the distribution grid area if the distribution grid type is a square, to obtain the first distance and the second distance, and the first distance and the second distance are the same; if the distribution grid type is a rectangle, then assign values to the first distance and the second distance so that the product of the first distance and the second distance is equal to the distribution grid area.
[0205] In another possible implementation, the apparatus further includes:
[0206] A seventh determination module is configured to determine a second work area of the infilled seismic work area and a third work area of the non-infilled seismic work area in the seismic work area, wherein the distribution density of the seabed nodes in the infilled seismic work area is greater than the distribution density of the seabed nodes in the non-infilled seismic work area;
[0207] an eighth determination module, configured to determine a fifth distance, a sixth distance, a seventh distance, and an eighth distance based on the area of the second work area, the area of the third work area, and the first quantity, wherein the fifth distance is used to represent the distance between two adjacent seabed nodes on a receiving line in the intensified seismic work area, the sixth distance is used to represent the distance between two adjacent receiving lines in the intensified seismic work area, the seventh distance is used to represent the distance between two adjacent seabed nodes on a receiving line in the non-intensified seismic work area, and the eighth distance is used to represent the distance between two adjacent receiving lines in the non-intensified seismic work area;
[0208] a ninth determination module configured to determine a third position of the plurality of seafloor nodes in the infilled seismic work area based on the fifth distance and the sixth distance; and to determine a fourth position of the plurality of seafloor nodes in the non-infilled seismic work area based on the seventh distance and the eighth distance;
[0209] The acquisition module 1206 is used to acquire seismic data through an observation system consisting of an ocean floor node at every fourth position in the non-encrypted seismic work area, an ocean floor node at every third position in the encrypted seismic work area, and a shot point at every second position.
[0210] In another possible implementation, the eighth determination module is used to determine the encrypted grid type of the encrypted seismic work area and the non-encrypted grid type of the non-encrypted seismic work area; based on the encrypted grid type, the first encryption parameter and the second encryption parameter are assigned; based on the assigned first encryption parameter, the assigned second encryption parameter, the second work area area, the third work area area and the first quantity, the non-encrypted grid area of the non-encrypted seismic work area is determined, so that the sum of the fourth number of seabed nodes in the non-encrypted seismic work area and the fifth number of seabed nodes in the encrypted seismic work area is the first quantity, the fourth quantity is the ratio of the third work area area to the non-encrypted grid area, the fifth quantity is the ratio of the second work area area to the encrypted grid area, and the encrypted grid area is the product of the non-encrypted grid area, the assigned first encryption parameter and the assigned second encryption parameter; based on the non-encrypted grid type and the non-encrypted grid area, the seventh distance and the eighth distance are determined; the product of the seventh distance and the assigned first encryption parameter is determined to obtain the fifth distance; the product of the eighth distance and the assigned second encryption parameter is determined to obtain the sixth distance.
[0211] In another possible implementation, the eighth determination module is configured to assign the same value to the first encryption parameter and the second encryption parameter if the encrypted grid type is a square; and to assign different values to the first encryption parameter and the second encryption parameter if the encrypted grid type is a rectangle.
[0212] An embodiment of the present application provides a seismic data acquisition device, which, under the condition that the number of seabed nodes is limited, realizes the acquisition of wide-azimuth or full-azimuth seismic data by determining a first distance between two adjacent seabed nodes on a receiving line, a second distance between two adjacent receiving lines, and a second number of receiving lines, thereby reducing the cost of seismic data acquisition.
[0213] It should be noted that the seismic data acquisition device provided in the above embodiment is merely illustrated by the division of the aforementioned functional modules when collecting seismic data. In actual applications, the aforementioned functions can be assigned to different functional modules as needed, i.e., the internal structure of the device can be divided into different functional modules to perform all or part of the functions described above. Furthermore, the seismic data acquisition device provided in the above embodiment and the seismic data acquisition method embodiment are based on the same concept. The specific implementation process is detailed in the method embodiment and will not be repeated here.
[0214] Figure 13 The following is a block diagram of a computer device 1300 according to an exemplary embodiment of the present application. Computer device 1300 may be a portable mobile computer device, such as a smartphone, a tablet computer, an MP3 player (Moving Picture Experts Group Audio Layer III), an MP4 player (Moving Picture Experts Group Audio Layer IV), a laptop computer, or a desktop computer. Computer device 1300 may also be referred to as a user device, a portable computer device, a laptop computer device, a desktop computer device, or other similar names.
[0215] Typically, the computer device 1300 includes a processor 1301 and a memory 1302 .
[0216] The processor 1301 may include one or more processing cores, such as a 4-core processor, an 8-core processor, etc. The processor 1301 may be implemented in at least one hardware form of DSP (Digital Signal Processing), FPGA (Field-Programmable Gate Array), or PLA (Programmable Logic Array). The processor 1301 may also include a main processor and a coprocessor. The main processor is a processor for processing data in the awake state, also known as a CPU (Central Processing Unit); the coprocessor is a low-power processor for processing data in the standby state. In some embodiments, the processor 1301 may be integrated with a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content to be displayed on the display screen. In some embodiments, the processor 1301 may also include an AI (Artificial Intelligence) processor, which is used to process computing operations related to machine learning.
[0217] Memory 1302 may include one or more computer-readable storage media, which may be non-transitory. Memory 1302 may also include high-speed random access memory and non-volatile memory, such as one or more disk storage devices and flash memory storage devices. In some embodiments, the non-transitory computer-readable storage medium in memory 1302 is used to store at least one program code, which is executed by processor 1301 to implement the seismic data acquisition method provided in the method embodiment of the present application.
[0218] In some embodiments, computer device 1300 may optionally include a peripheral device interface 1303 and at least one peripheral device. Processor 1301, memory 1302, and peripheral device interface 1303 may be connected via a bus or signal lines. Each peripheral device may be connected to peripheral device interface 1303 via a bus, signal lines, or circuit boards. Specifically, the peripheral device may include at least one of a radio frequency circuit 1304, a display screen 1305, a camera assembly 1306, an audio circuit 1307, a positioning assembly 1308, and a power supply 1309.
[0219] The peripheral device interface 1303 can be used to connect at least one I / O (Input / Output)-related peripheral device to the processor 1301 and the memory 1302. In some embodiments, the processor 1301, the memory 1302, and the peripheral device interface 1303 are integrated on the same chip or circuit board; in some other embodiments, any one or two of the processor 1301, the memory 1302, and the peripheral device interface 1303 can be implemented on separate chips or circuit boards, which is not limited in this embodiment.
[0220] The RF circuit 1304 is used to receive and transmit RF (Radio Frequency) signals, also known as electromagnetic signals. The RF circuit 1304 communicates with communication networks and other communication devices via electromagnetic signals. The RF circuit 1304 converts electrical signals into electromagnetic signals for transmission, or converts received electromagnetic signals into electrical signals. Optionally, the RF circuit 1304 includes an antenna system, an RF transceiver, one or more amplifiers, a tuner, an oscillator, a digital signal processor, a codec chipset, a user identity module card, and the like. The RF circuit 1304 can communicate with other computer devices via at least one wireless communication protocol. Such wireless communication protocols include, but are not limited to, the World Wide Web, metropolitan area networks, intranets, various generations of mobile communication networks (2G, 3G, 4G, and 5G), wireless local area networks, and / or WiFi (Wireless Fidelity) networks. In some embodiments, the RF circuit 1304 may also include circuitry related to Near Field Communication (NFC), although this application does not limit this.
[0221] The display screen 1305 is used to display a UI (User Interface). The UI may include graphics, text, icons, videos, and any combination thereof. When the display screen 1305 is a touch screen display, the display screen 1305 also has the ability to collect touch signals on the surface or above the surface of the display screen 1305. The touch signal can be input as a control signal to the processor 1301 for processing. In this case, the display screen 1305 can also be used to provide virtual buttons and / or virtual keyboards, also known as soft buttons and / or soft keyboards. In some embodiments, there can be one display screen 1305, which is set on the front panel of the computer device 1300; in other embodiments, there can be at least two display screens 1305, which are respectively set on different surfaces of the computer device 1300 or in a folding design; in other embodiments, the display screen 1305 can be a flexible display screen, which is set on the curved surface or folding surface of the computer device 1300. Even more, the display screen 1305 can be set to a non-rectangular irregular shape, that is, a special-shaped screen. The display screen 1305 can be made of materials such as LCD (Liquid Crystal Display) and OLED (Organic Light-Emitting Diode).
[0222] The camera assembly 1306 is used to capture images or videos. Optionally, the camera assembly 1306 includes a front camera and a rear camera. Typically, the front camera is set on the front panel of the computer device, and the rear camera is set on the back of the computer device. In some embodiments, there are at least two rear cameras, which are any one of a main camera, a depth of field camera, a wide-angle camera, and a telephoto camera, so as to realize the fusion of the main camera and the depth of field camera to realize the background blur function, the fusion of the main camera and the wide-angle camera to realize panoramic shooting and VR (Virtual Reality) shooting function or other fusion shooting functions. In some embodiments, the camera assembly 1306 may also include a flash. The flash can be a single-color temperature flash or a dual-color temperature flash. A dual-color temperature flash refers to a combination of a warm light flash and a cold light flash, which can be used for light compensation under different color temperatures.
[0223] The audio circuit 1307 may include a microphone and a speaker. The microphone is used to collect sound waves from the user and the environment, and convert the sound waves into electrical signals that are input into the processor 1301 for processing, or input into the radio frequency circuit 1304 to achieve voice communication. For the purpose of stereo sound collection or noise reduction, there may be multiple microphones, each located in different parts of the computer device 1300. The microphone may also be an array microphone or an omnidirectional collection microphone. The speaker is used to convert electrical signals from the processor 1301 or the radio frequency circuit 1304 into sound waves. The speaker may be a traditional thin film speaker or a piezoelectric ceramic speaker. When the speaker is a piezoelectric ceramic speaker, it can not only convert electrical signals into sound waves audible to humans, but also convert electrical signals into sound waves inaudible to humans for purposes such as distance measurement. In some embodiments, the audio circuit 1307 may also include a headphone jack.
[0224] The positioning component 1308 is used to locate the current geographic location of the computer device 1300 to implement navigation or LBS (Location Based Service). The positioning component 1308 can be a positioning component based on the US GPS (Global Positioning System), China's Beidou system, or Russia's Galileo system.
[0225] Power supply 1309 is used to power various components in computer device 1300. Power supply 1309 can be AC power, DC power, a disposable battery, or a rechargeable battery. When power supply 1309 includes a rechargeable battery, the rechargeable battery can be a wired rechargeable battery or a wireless rechargeable battery. A wired rechargeable battery is a battery that is charged via a wired line, while a wireless rechargeable battery is a battery that is charged via a wireless coil. The rechargeable battery can also be used to support fast charging technology.
[0226] In some embodiments, the computer device 1300 further includes one or more sensors 1310 , including but not limited to an acceleration sensor 1311 , a gyroscope sensor 1312 , a pressure sensor 1313 , a fingerprint sensor 1314 , an optical sensor 1315 , and a proximity sensor 1316 .
[0227] The accelerometer 1311 can detect the magnitude of acceleration along the three coordinate axes of the coordinate system established by the computer device 1300. For example, the accelerometer 1311 can be used to detect the components of gravity acceleration along the three coordinate axes. The processor 1301 can control the display screen 1305 to display the user interface in a landscape or portrait view based on the gravity acceleration signal collected by the accelerometer 1311. The accelerometer 1311 can also be used to collect game or user motion data.
[0228] The gyroscope sensor 1312 can detect the orientation and rotation angle of the computer device 1300. It can also work with the accelerometer 1311 to collect 3D motions of the user on the computer device 1300. Based on the data collected by the gyroscope sensor 1312, the processor 1301 can implement the following functions: motion sensing (e.g., changing the UI based on the user's tilt), image stabilization during shooting, game control, and inertial navigation.
[0229] The pressure sensor 1313 can be installed on the side frame of the computer device 1300 and / or below the display screen 1305. When the pressure sensor 1313 is installed on the side frame of the computer device 1300, it can detect the user's grip signal on the computer device 1300. The processor 1301 can perform left and right hand recognition or shortcut operations based on the grip signal collected by the pressure sensor 1313. When the pressure sensor 1313 is installed below the display screen 1305, the processor 1301 controls the operational controls on the UI interface based on the user's pressure operation on the display screen 1305. The operational controls include at least one of a button control, a scroll bar control, an icon control, and a menu control.
[0230] The fingerprint sensor 1314 is used to collect the user's fingerprint. The processor 1301 identifies the user's identity based on the fingerprint collected by the fingerprint sensor 1314, or the fingerprint sensor 1314 identifies the user's identity based on the collected fingerprint. When the user's identity is recognized as a trusted identity, the processor 1301 authorizes the user to perform relevant sensitive operations, such as unlocking the screen, viewing encrypted information, downloading software, making payments, and changing settings. The fingerprint sensor 1314 can be set on the front, back, or side of the computer device 1300. When a physical button or manufacturer logo is set on the computer device 1300, the fingerprint sensor 1314 can be integrated with the physical button or manufacturer logo.
[0231] Optical sensor 1315 is used to detect ambient light intensity. In one embodiment, processor 1301 can control the display brightness of display screen 1305 based on the ambient light intensity detected by optical sensor 1315. Specifically, when the ambient light intensity is high, the display brightness of display screen 1305 is increased; when the ambient light intensity is low, the display brightness of display screen 1305 is decreased. In another embodiment, processor 1301 can also dynamically adjust the shooting parameters of camera assembly 1306 based on the ambient light intensity detected by optical sensor 1315.
[0232] Proximity sensor 1316, also known as a distance sensor, is typically located on the front panel of computer device 1300. Proximity sensor 1316 is used to detect the distance between the user and the front of computer device 1300. In one embodiment, when proximity sensor 1316 detects that the distance between the user and the front of computer device 1300 is gradually decreasing, processor 1301 controls display screen 1305 to switch from the screen-on state to the screen-off state. When proximity sensor 1316 detects that the distance between the user and the front of computer device 1300 is gradually increasing, processor 1301 controls display screen 1305 to switch from the screen-off state to the screen-on state.
[0233] Those skilled in the art will understand that Figure 13 The structure shown in the figure does not constitute a limitation on the computer device 1300, and the computer device 1300 may include more or fewer components than shown in the figure, or combine some components, or adopt a different component arrangement.
[0234] An embodiment of the present application further provides a computer-readable storage medium, in which at least one program code is stored. The at least one program code is loaded and executed by a processor to implement the operations performed in the seismic data acquisition method in the embodiment of the present application.
[0235] An embodiment of the present application further provides a computer program product, which stores at least one program code, and the at least one program code is loaded and executed by a processor to implement the operations performed in the seismic data acquisition method described in the embodiment of the present application.
[0236] In some embodiments, the computer program involved in the embodiments of the present application may be deployed and executed on a computer device, or on multiple computer devices located at one location, or on multiple computer devices distributed at multiple locations and interconnected through a communication network. Multiple computer devices distributed at multiple locations and interconnected through a communication network may constitute a blockchain system.
[0237] The above description is only for the purpose of facilitating those skilled in the art to understand the technical solution of this application and is not intended to limit this application. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of this application shall be included in the scope of protection of this application.
Claims
1. A seismic data acquisition method, characterized in that: The method comprises: determining a first number of seafloor nodes available for seismic data acquisition in a seismic work area; Determining a distribution grid area of the seabed nodes based on the first number, the aspect ratio, the maximum longitudinal offset, the length of the receiving line, and the type of the observation system, where the distribution grid area is the product of a first distance between two adjacent seabed nodes on a receiving line and a second distance between two adjacent receiving lines, where the receiving line is a straight line formed by multiple seabed nodes along a seismic data acquisition direction, the aspect ratio is the ratio of the maximum lateral offset to the maximum longitudinal offset, and the aspect ratio is greater than a preset threshold, and both the maximum longitudinal offset and the maximum lateral offset are parameters representing the distance between a shot point and a seabed node; Determining a distribution grid type of a seabed node, and determining the first distance and the second distance based on the distribution grid area and the distribution grid type; determining a second number of receive lines based on the longitudinal maximum offset, the aspect ratio, the second distance, and the observation system type; Determining a third distance and a fourth distance based on a bin size parameter, wherein the third distance is used to represent the distance between two adjacent shot points on a shot line, and the fourth distance is used to represent the distance between two adjacent shot lines. The bin size parameter is a parameter representing imaging accuracy when imaging seismic data; Based on the first distance, the second distance and the second number, the first positions of multiple seabed nodes are determined, and based on the third distance and the fourth distance, the second positions of multiple shot points are determined; and seismic data are collected through an observation system consisting of each seabed node at the first position and each shot point at the second position.
2. The method according to claim 1, characterized in that The determining of the distribution grid area of the seabed nodes based on the first number, the aspect ratio, the maximum longitudinal offset, the length of the receiving line, and the type of the observation system includes: Determining, based on the type of the observation system, a width parameter of the first number of seabed nodes; Determining a third number of seabed nodes per unit area based on the first number, the aspect ratio, the maximum longitudinal offset, the length of the receiving line, and the width parameter; The reciprocal of the third quantity is determined to obtain the distribution grid area.
3. The method according to claim 2, characterized in that The determining, based on the first number, the aspect ratio, the maximum longitudinal offset, the length of the receiving line, and the width parameter, of a third number of seabed nodes per unit area includes: Determine the product of the aspect ratio, the maximum longitudinal offset, the length of the receiving line, and the width parameter to obtain the area of a first working area of the seismic working area; The ratio of the first quantity to the area of the first work area is determined to obtain the third quantity.
4. The method according to claim 1, wherein The determining the first distance and the second distance based on the distribution grid area and the distribution grid type includes: If the distribution grid type is square, determining the square root of the area of the distribution grid to obtain the first distance and the second distance, and the first distance and the second distance are the same; If the distribution grid type is rectangular, the first distance and the second distance are assigned values so that the product of the first distance and the second distance is equal to the distribution grid area.
5. The method according to claim 1, wherein The method further comprises: Determining a second work area of an infilled seismic work area and a third work area of a non-infilled seismic work area in the seismic work area, wherein the distribution density of the seabed nodes in the infilled seismic work area is greater than the distribution density of the seabed nodes in the non-infilled seismic work area; Based on the area of the second work area, the area of the third work area, and the first number, a fifth distance, a sixth distance, a seventh distance, and an eighth distance are determined, wherein the fifth distance is used to represent the distance between two adjacent seabed nodes on a receiving line in the encrypted seismic work area, the sixth distance is used to represent the distance between two adjacent receiving lines in the encrypted seismic work area, the seventh distance is used to represent the distance between two adjacent seabed nodes on a receiving line in the non-encrypted seismic work area, and the eighth distance is used to represent the distance between two adjacent receiving lines in the non-encrypted seismic work area; Determine a third position of a plurality of seafloor nodes in the infilled seismic work area based on the fifth distance and the sixth distance; determine a fourth position of a plurality of seafloor nodes in the non-infilled seismic work area based on the seventh distance and the eighth distance; The acquisition of seismic data by an observation system composed of each seafloor node at the first position and each shot point at the second position includes: Seismic data is collected through an observation system consisting of an ocean floor node at each fourth position in the non-encrypted seismic work area, an ocean floor node at each third position in the encrypted seismic work area, and a shot point at each second position.
6. The method according to claim 5, characterized in that The determining of the fifth distance, the sixth distance, the seventh distance, and the eighth distance based on the area of the second work area, the area of the third work area, and the first number includes: Determine an encrypted grid type of the encrypted seismic work area and a non-encrypted grid type of the non-encrypted seismic work area; Assigning a first encryption parameter and a second encryption parameter based on the encryption grid type; Determine the non-encrypted grid area of the non-encrypted seismic work area based on the assigned first encryption parameter, the assigned second encryption parameter, the area of the second work area, the area of the third work area, and the first number, so that the sum of the fourth number of seabed nodes in the non-encrypted seismic work area and the fifth number of seabed nodes in the encrypted seismic work area is the first number, the fourth number is the ratio of the area of the third work area to the non-encrypted grid area, the fifth number is the ratio of the area of the second work area to the encrypted grid area, and the encrypted grid area is the product of the non-encrypted grid area, the assigned first encryption parameter, and the assigned second encryption parameter; Determining the seventh distance and the eighth distance based on the non-encrypted grid type and the non-encrypted grid area; determining a product of the seventh distance and the assigned first encryption parameter to obtain the fifth distance; A product of the eighth distance and the assigned second encryption parameter is determined to obtain the sixth distance.
7. A seismic data acquisition device, characterized in that: The device comprises: A first determining module is used to determine a first number of seafloor nodes that can be used for seismic data acquisition in a seismic work area; a second determining module, which determines a distribution grid area of the seafloor nodes based on the first number, the aspect ratio, the maximum longitudinal offset, the length of the receiving line, and the type of the observation system, wherein the distribution grid area is the product of a first distance between two adjacent seafloor nodes on a receiving line and a second distance between two adjacent receiving lines, wherein the receiving line is a straight line formed by multiple seafloor nodes along the seismic data acquisition direction, the aspect ratio is the ratio of the maximum lateral offset to the maximum longitudinal offset, and the aspect ratio is greater than a preset threshold, and the maximum longitudinal offset and the maximum lateral offset are both parameters representing the distance between the shot point and the seafloor node; a third determining module, configured to determine a distribution grid type of the seabed node, and determine the first distance and the second distance based on the distribution grid area and the distribution grid type; a fourth determining module, configured to determine a second number of receiving lines based on the maximum longitudinal offset, the aspect ratio, the second distance, and the observation system type; a fifth determining module, configured to determine a third distance and a fourth distance based on a bin size parameter, wherein the third distance represents the distance between two adjacent shot points on a shot line, and the fourth distance represents the distance between two adjacent shot lines, and the bin size parameter represents the imaging accuracy of seismic data imaging; An acquisition module is configured to determine a first position of a plurality of seafloor nodes based on the first distance, the second distance, and the second number, and to determine a second position of a plurality of shot points based on the third distance and the fourth distance; and to acquire seismic data through an observation system consisting of each seafloor node at the first position and each shot point at the second position.
8. A computer device, characterized in that: The computer device includes a processor and a memory, wherein the memory stores at least one program code, and the at least one program code is loaded and executed by the processor to implement the seismic data acquisition method according to any one of claims 1 to 6.
9. A computer-readable storage medium, characterized in that At least one program code is stored in the storage medium, and the at least one program code is loaded and executed by the processor to implement the seismic data acquisition method according to any one of claims 1 to 6.
10. A computer program product, characterized in that The computer program product stores at least one program code, and the at least one program code is loaded and executed by a processor to implement the seismic data acquisition method according to any one of claims 1 to 6.
Citation Information
Patent Citations
Simulation full three-dimensional earthquake acquisition method for improving exploration precision of reservoir layer
CN103823238A
Seismic data acquisition method and device, computer equipment and storage medium
CN113009573A