CN108761523B - Method and device for optimizing combined excitation among seismic sources - Google Patents
Method and device for optimizing combined excitation among seismic sources Download PDFInfo
- Publication number
- CN108761523B CN108761523B CN201810371633.1A CN201810371633A CN108761523B CN 108761523 B CN108761523 B CN 108761523B CN 201810371633 A CN201810371633 A CN 201810371633A CN 108761523 B CN108761523 B CN 108761523B
- Authority
- CN
- China
- Prior art keywords
- excitation
- seismic
- source
- sources
- combination
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/02—Generating seismic energy
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The embodiment of the application provides a method and a device for optimizing combined excitation among seismic sources, wherein the method comprises the following steps: determining the positions of a plurality of seismic sources meeting preset excitation conditions in a target area when a time point is appointed; determining all possible source excitation combinations of the plurality of sources at the specified time point, determining distances between the sources in each source excitation combination according to the positions, and selecting one source excitation combination from all the possible source excitation combinations as a preferred source excitation combination according to the distances; determining an order of excitation of the preferred source excitation combination. The embodiment of the application can improve the operation efficiency of seismic data acquisition.
Description
Technical Field
The application relates to the technical field of multi-seismic-source seismic data acquisition, in particular to a method and a device for optimizing combined excitation among seismic sources.
Background
Seismic data acquisition is the primary working link of oil and gas exploration, and with the development of seismic acquisition technology, high-efficiency acquisition technology is more and more widely adopted by oil companies at home and abroad. The application of the high-efficiency acquisition method has more and more requirements on the number of acquisition excitation devices (seismic sources), which not only brings more and more challenges to field management work of field acquisition receiving devices, but also provides higher requirements on how to more scientifically and effectively perform acquisition excitation management.
In general, regardless of the manner in which the seismic sources are fired, the time-distance rule may need to be satisfied. For example, assuming that there are two seismic sources a, B and the distance between the two seismic sources is d, the corresponding time-distance rule may be set as:
1) if d <1km, B must be excited after A (if A is excited first), 16 seconds later;
2) if 1< d <4km, B must be excited after a (if a is excited first), 12 seconds later;
3) if d <12km, A and B must be satisfied by A; the sliding time is higher than the time determined by the oblique lines in the figure;
4) if 12km < ═ d. AB may be excited simultaneously.
In a conventional seismic acquisition process, an instrument (for example, a Secel instrument) needs to manage a plurality of seismic sources on site, generally, in an acquisition and production process, production excitation sequence sequencing is determined in the instrument according to the distance from each seismic source to be vibrated to an activated seismic source, when a new seismic source is added into a sequencing state, the distance between the new seismic source and the activated seismic source and the distance between the new seismic source and the seismic source to be vibrated are compared, a time-distance rule is followed in the instrument, the distance is selected to be the largest, and the seismic sources meeting synchronous excitation conditions are preferentially excited.
However, the inventors of the present application have studied and found that: under the condition of high-efficiency acquisition, partial seismic sources in a plurality of seismic sources are often redefined and sequenced in the instrument, and because a sequencing excitation principle of synchronous excitation priority may exist, the principle is favorable for grouping the seismic sources meeting the synchronous excitation to the most preferential excitation position. This will result in that shots at both ends of the spread (because the seismic sources at both ends of the spread are easier to combine with other seismic sources to form a seismic source group satisfying synchronous excitation) will have a large probability of being selected as a synchronous excitation mode, while the excitation point at the middle position is often selected as a sliding scan excitation mode or an alternate scan excitation mode in terms of excitation mode. Therefore, the excitation acquisition modes of the whole line of the shot points are distributed unevenly, and the probability that the shot points at two ends are arranged in a synchronous excitation mode is higher; and the shot point positioned in the middle of the beam has higher probability of adopting a sliding scanning excitation mode or an alternate scanning excitation mode. Especially in the case of super-arrays, the recovery of the entire array is dependent on the acquisition completion time of the last shot using the array. Therefore, the collection efficiency is inevitably affected under the unbalanced collection mode condition.
Therefore, how to perform optimal selection among seismic source combinations under the condition of satisfying the time-distance rule for a great number of seismic source excitation devices so as to further improve the field acquisition operation efficiency is a technical problem to be solved urgently at present.
Disclosure of Invention
The embodiment of the application aims to provide a method and a device for optimizing combined excitation among seismic sources so as to improve the working efficiency of seismic data acquisition.
In order to achieve the above object, in one aspect, an embodiment of the present application provides a method for optimizing combined excitation among seismic sources, including:
determining the positions of a plurality of seismic sources meeting preset excitation conditions in a target area when a time point is appointed;
determining all possible source excitation combinations of the plurality of sources at the specified time point, determining distances between the sources in each source excitation combination according to the positions, and selecting one source excitation combination from all the possible source excitation combinations as a preferred source excitation combination according to the distances;
determining an order of excitation of the preferred source excitation combination.
In the inter-seismic source combined excitation optimization method according to the embodiment of the application, the specified time points are time points with equal time intervals.
In the inter-seismic source combined excitation optimization method in the embodiment of the application, the duration of the equal time interval is positively correlated with the number of seismic sources which can be input into the target area and the average time required for each seismic source to complete one-time excitation acquisition.
In the inter-seismic source combined excitation optimization method according to the embodiment of the application, the determining all possible combinations of seismic source excitations of the plurality of seismic sources at the specified time point includes:
determining all possible combinations of source shots for the plurality of sources at the specified point in time by a traversal algorithm.
In the inter-seismic source combined excitation optimization method according to the embodiment of the application, the selecting one seismic source excitation combination from all the possible seismic source excitation combinations as a preferred seismic source excitation combination according to the distance includes:
for each seismic source excitation combination, determining the maximum value of the number of seismic sources meeting the synchronous excitation condition in the seismic source excitation combination according to the distance between the seismic sources in the seismic source excitation combination, and determining the excitation batch required by completing one-time excitation of all the seismic sources in the seismic source excitation combination;
and selecting the seismic source excitation combination with the maximum seismic source number and the minimum excitation batch as the preferred seismic source excitation combination from all possible seismic source excitation combinations at the specified time point.
In the inter-seismic source combination excitation optimization method according to the embodiment of the present application, the excitation sequence of the preferred seismic source excitation combination sequentially includes: synchronous excitation, sliding scan excitation, and alternating scan excitation.
On the other hand, the embodiment of the present application further provides an inter-seismic source combined excitation optimization apparatus, including:
the seismic source position determining module is used for determining the positions of a plurality of seismic sources meeting preset excitation conditions in a target area when a time point is appointed;
a firing combination determination module to determine all possible combinations of firing of the seismic sources at the specified time point for the plurality of seismic sources; determining the distance between the seismic sources in each seismic source excitation combination according to the positions, and selecting one seismic source excitation combination from all the possible seismic source excitation combinations as a preferred seismic source excitation combination according to the distance;
and the excitation sequence determination module is used for determining the excitation sequence of the preferred seismic source excitation combination.
In the inter-seismic source combined excitation optimization apparatus according to the embodiment of the present application, the specified time points are time points at equal time intervals.
In the inter-seismic source combined excitation optimization device in the embodiment of the application, the duration of the equal time interval is positively correlated with the number of seismic sources which can be input into a target area and the average time required for completing one-time excitation acquisition of each seismic source.
In the inter-seismic source combined excitation optimization apparatus according to the embodiment of the present application, the determining all possible combinations of seismic source excitations of the plurality of seismic sources at the specified time point includes:
determining all possible combinations of source shots for the plurality of sources at the specified point in time by a traversal algorithm.
In the inter-seismic source combination excitation optimization apparatus according to the embodiment of the present application, the selecting one seismic source excitation combination from the all possible seismic source excitation combinations as the preferred seismic source excitation combination according to the distance includes:
for each seismic source excitation combination, determining the maximum value of the number of seismic sources meeting the synchronous excitation condition in the seismic source excitation combination according to the distance between the seismic sources in the seismic source excitation combination, and determining the excitation batch required by completing one-time excitation of all the seismic sources in the seismic source excitation combination;
and selecting the seismic source excitation combination with the maximum seismic source number and the minimum excitation batch as the preferred seismic source excitation combination from all possible seismic source excitation combinations at the specified time point.
In the inter-seismic source combined excitation optimization apparatus according to the embodiment of the present application, the excitation sequence of the preferred seismic source excitation combination sequentially includes: synchronous excitation, sliding scan excitation, and alternating scan excitation.
In another aspect, an embodiment of the present application further provides another inter-seismic source combined excitation optimization apparatus, including a memory, a processor, and a computer program stored on the memory, where the computer program, when executed by the processor, performs the following steps:
determining the positions of a plurality of seismic sources meeting preset excitation conditions in a target area when a time point is appointed;
determining all possible source excitation combinations of the plurality of sources at the specified time point, determining distances between the sources in each source excitation combination according to the positions, and selecting one source excitation combination from all the possible source excitation combinations as a preferred source excitation combination according to the distances;
determining an order of excitation of the preferred source excitation combination.
According to the technical scheme provided by the embodiment of the application, the positions of a plurality of seismic sources meeting the preset excitation condition in the target area are determined when the time point is appointed; secondly, determining all possible seismic source excitation combinations of a plurality of seismic sources at a specified time point, then determining the distance between the seismic sources in each seismic source excitation combination according to the positions, and selecting one seismic source excitation combination from all the possible seismic source excitation combinations as a preferred seismic source excitation combination according to the distance; and finally determining the excitation sequence of the preferred source excitation combination. Therefore, under the condition of a time-distance rule, the number of synchronously excited seismic sources is found out to the maximum, and the excitation sequence among seismic source groups is redefined, so that the working efficiency of seismic data acquisition is improved.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, it is obvious that the drawings in the following description are only some embodiments described in the present application, and for those skilled in the art, other drawings can be obtained according to the drawings without any creative effort. In the drawings:
FIG. 1 is a flow chart of a method for optimizing the combined excitation among seismic sources according to an embodiment of the present application;
FIG. 2a is a schematic diagram of a traversal algorithm based on depth search according to an embodiment of the present application;
FIG. 2b is a schematic diagram of a traversal algorithm based on breadth search according to an embodiment of the present application;
FIG. 3 is a schematic representation of a seismic source distribution in an embodiment of the present application;
FIG. 4 is a schematic diagram illustrating the combined spacing of four secondary seismic sources according to an embodiment of the present application;
FIG. 5 is a schematic diagram illustrating the combined spacing of three secondary seismic sources according to an embodiment of the present application;
FIG. 6 is a schematic diagram illustrating a combined spacing between two secondary seismic sources according to an embodiment of the present application;
FIG. 7 is a block diagram of an inter-source combined excitation optimization apparatus according to an embodiment of the present application;
fig. 8 is a block diagram of an inter-source combined excitation optimization apparatus according to another embodiment of the present application.
Detailed Description
In order to make those skilled in the art better understand the technical solutions in the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application. For example, in the following description, forming the second component over the first component may include embodiments in which the first and second components are formed in direct contact, embodiments in which the first and second components are formed in non-direct contact (i.e., additional components may be included between the first and second components), and so on.
Also, for ease of description, some embodiments of the present application may use spatially relative terms such as "above …," "below …," "top," "below," etc., to describe the relationship of one element or component to another (or other) element or component as illustrated in the various figures of the embodiments. It will be understood that the spatially relative terms are intended to encompass different orientations of the device in use or operation in addition to the orientation depicted in the figures. For example, if the device in the figures is turned over, elements or components described as "below" or "beneath" other elements or components would then be oriented "above" or "over" the other elements or components.
Referring to fig. 1, the method for optimizing the combined excitation among the seismic sources according to the embodiment of the present application may include the following steps:
s101, when a time point is appointed, the positions of a plurality of seismic sources meeting preset excitation conditions in a target area are determined.
In some embodiments of the present application, the designated time point may be a time point with equal time intervals. For example, starting from a certain time point, every 30 minutes is a designated time point. The time length of the equal time interval is an empirical value, and is generally in positive correlation with the number of the seismic sources which can be input in the target area and the average time required by each seismic source to complete one-time excitation acquisition, that is, the more the number of the input seismic sources and/or the longer the average time required by each seismic source to complete one-time excitation acquisition, the longer the corresponding time length of the equal time interval is set, so as to ensure that all or most of the seismic sources are ready for the next excitation (that is, all or most of the seismic sources are distributed in the whole seismic array in a ready-to-excite state) when the next specified time point arrives. Wherein, each source is independent, and the distribution mode between the sources can be random layout in the seismic arrangement.
In some embodiments of the present application, the predetermined excitation condition is satisfied, i.e., the seismic source is ready for the next excitation. In some embodiments of the present application, each seismic source has its activation readiness flag, such as ready (OK), not ready (NO), etc., and when the seismic source is in which activation readiness state, it can be marked with the corresponding activation readiness flag, so that it can be confirmed whether it satisfies the activation condition by looking at the activation readiness flag of each seismic source.
In some embodiments of the present application, each source, regardless of its state, has a corresponding coordinate location at a given point in time that can be used for a preferred calculation of a firing pattern between subsequent sources. In some embodiments of the present application, the coordinate location may be determined by the coordinates of the shots in the spread, or by a locating device onboard the seismic source, or the like.
S102, determining all possible seismic source excitation combinations of the plurality of seismic sources at the specified time point; and determining the distance between the seismic sources in each seismic source excitation combination according to the positions, and selecting one seismic source excitation combination from all the possible seismic source excitation combinations as a preferred seismic source excitation combination according to the distance.
In some embodiments of the present application, the selecting one of the possible source excitation combinations as the preferred source excitation combination according to the distance may be, for example:
firstly, for each seismic source excitation combination, determining the maximum value of the number of seismic sources meeting synchronous excitation conditions in the seismic source excitation combination according to the distance between the seismic sources in the seismic source excitation combination, and determining an excitation batch required by one-time excitation of all the seismic sources in the seismic source excitation combination;
then, from all possible combinations of the seismic source excitations at the specified time point, the combination of the seismic source excitations with the maximum seismic source number and the minimum excitation batch is selected as the preferred combination of the seismic source excitations.
In some embodiments of the present application, all possible combinations of source shots at the specified time point for the plurality of sources may be determined by a traversal algorithm. The traversal algorithm is also called exhaustion method, which refers to enumerating all possible states one by one to obtain results or checking until all solutions are found or the checking is completed. While the optimal solution or result in this state can be selected from the obtained solution and result, the traversal algorithm can generally perform a deep search (for example, as shown in fig. 2 a) and a wide search (for example, as shown in fig. 2 b), and when the number of the seismic source samples to be analyzed and the target data volume to be studied is relatively small, the traversal method can be considered to perform the optimal solution of the result.
Referring to FIG. 3, in an exemplary embodiment of the present application, there are a total of 5 sources in the spread, each labeled A, B, C, D, E. And each seismic source adopts a single acquisition mode to acquire data seismic. The source locations may be arbitrarily distributed over the entire spread.
Considering the principle of synchronous excitation priority, firstly judging whether 5 seismic source synchronous excitations are met:
at a specified time point, when the distances between every two of the 5 seismic sources meet the synchronous excitation distanceWhen the seismic sources are separated, 5 seismic sources can be synchronously excited; the corresponding seismic source excitation combination number is C 5 51, i.e. only one source excitation combination: ABCDE. Since 5 seismic sources can be fired simultaneously, the corresponding firing batch is only 1 time in this case.
In an embodiment of the present application, if the distances between every two of the 5 seismic sources do not satisfy the synchronous excitation distance, 4 seismic source excitation combinations synchronously excited by the seismic sources are found, and the number of the corresponding seismic source excitation combinations is C 5 45, there are 5 combinations of source excitations: ABCD, BCDE, CDEA, DEAB and EABC. For the 5 seismic source excitation combinations, respectively calculating the distance between every two 4 seismic sources in each seismic source excitation combination; for each seismic source excitation combination, the corresponding seismic source excitation combination can be synchronously excited by 4 seismic sources only when the distances between every two 4 seismic sources in the combination meet the synchronous excitation condition. The corresponding calculation and judgment can be seen in fig. 4. In FIG. 4, dsIs the synchronous excitation distance, d is the separation between two seismic sources, and λ is an empirical coefficient greater than 1.
In an embodiment of the present application, if there are more than 5 seismic source excitation combinations of ABCD, BCDE, CDEA, DEAB, and EABC that satisfy the synchronous excitation condition, one of the combinations may be selected and combined with the remaining seismic sources to form a preferred seismic source excitation combination. For example, ABCD and BCDE meet the synchronous firing conditions, ABCD + E or BCDE + a may optionally be used as a preferred source firing combination.
Correspondingly, since 4 of the 5 seismic sources can be synchronously excited, the rest of the 5 seismic sources need to be separately excited, and the corresponding excitation batch is 2 times in this case. For example, when the preferred source excitation combination is ABCD + E, one excitation field is required for ABCD synchronous excitation; after the synchronous excitation of the ABCD, E needs to be excited separately, and therefore, there is also one excitation field, and correspondingly, 2 excitation batches are needed for completing one excitation of each of the 5 seismic sources in the ABCD + E.
In an embodiment of the present application, if none of the 5 source excitation combinations ABCD, BCDE, CDEA, DEAB, and EABC satisfies 4 source synchronous excitations, 3 source synchronous excitations can be further foundThe number of the corresponding seismic source excitation combinations is C 5 310, there are 10 combinations of source excitations: ABC, ABD, ABE, ACD, ACE, ADE, BCD, BCE, BDE, and CDE. Similarly, for the 10 seismic source excitation combinations, respectively calculating the distance between every two 3 seismic sources in each seismic source excitation combination; for each seismic source excitation combination, only when the distance between every two 3 seismic sources in each seismic source excitation combination meets the synchronous excitation condition, the corresponding seismic source excitation combination can synchronously excite the 3 seismic sources. The corresponding calculation and judgment can be seen in fig. 5. In FIG. 5, dsIs the synchronous excitation distance, d is the separation between two seismic sources, and λ is an empirical coefficient greater than 1.
In one embodiment of the present application, if the 10 source shots are combined: when a plurality of ABC, ABD, ABE, ACD, ACE, ADE, BCD, BCE, BDE and CDE meet the synchronous excitation conditions, for the combination meeting the synchronous excitation conditions, whether the rest seismic sources can be synchronously excited or not needs to be correspondingly judged. For example, assuming that ABC and ABD satisfy the synchronous excitation condition, for ABC, it is necessary to determine whether D and E can be synchronously excited; for ABD, it is necessary to determine whether C and E can be fired simultaneously. The method comprises the following specific steps:
if for ABC, D and E can be fired simultaneously; and for ABD, C and E may not be fired simultaneously; then ABC + DE is the preferred source excitation combination. The corresponding challenge batch at this time was 2 times. If for ABC, D and E cannot be fired synchronously; and for ABD, C and E may be excited simultaneously; then ABD + CE is the preferred source excitation combination. The corresponding challenge batch at this time was 2 times.
If for ABC, D and E can be fired simultaneously; and for ABD, C and E may be excited simultaneously; then either ABC + DE or ABD + CE may be arbitrarily selected as the preferred source excitation combination. The corresponding challenge batch at this time was 2 times.
If for ABC, D and E cannot be fired synchronously; and for ABD, C and E can not be synchronously excited; then ABC + D + E or ABD + C + E may be arbitrarily selected as the preferred source excitation combination. The corresponding excitation batches were 3 times.
In one embodiment of the present application, if the 10 source shots are combined: when ABC, ABD, ABE, ACD, ACE, ADE, BCD, BCE, BDE and CDE do not meet the synchronous excitation condition, 2 seismic source excitation combinations synchronously excited by seismic sources can be further searched, and the number of the corresponding seismic source excitation combinations is C 5 210, there are 10 combinations of source excitations: AB. AC, AD, AE, BC, BD, BE, CD, CE, and DE. For the 10 seismic source excitation combinations, respectively calculating the spacing between 2 seismic sources in each seismic source excitation combination; for each seismic source excitation combination, the corresponding seismic source excitation combination can be synchronously excited by 2 seismic sources only when the spacing between 2 seismic sources in the combination meets the synchronous excitation condition. The corresponding calculation and judgment can be seen in fig. 6. In FIG. 6, dsIs the synchronous excitation distance, d is the separation between two seismic sources, and λ is an empirical coefficient greater than 1.
In one embodiment of the present application, if the 10 source shots are combined: AB. When a plurality of AC, AD, AE, BC, BD, BE, CD, CE, and DE satisfy the synchronous excitation conditions, one of the combinations satisfying the synchronous excitation conditions may BE selected and combined with the remaining sources to form a preferred source excitation combination. For example, AB, AC, and AD all satisfy the simultaneous excitation condition, any one of AB + C + D + E, AC + B + D + E and AD + B + C + E may be the preferred source excitation combination. At this time, the corresponding challenge batch was 4 times.
In one embodiment of the present application, if the 10 source shots are combined: AB. When AC, AD, AE, BC, BD, BE, CD, CE and DE do not satisfy the synchronous excitation condition, only one possible combination of the 5 seismic sources is A + B + C + D + E. In this case, the 5 seismic source excitation batches are 5 times.
In view of this, through the above calculation, the source excitation combination with the largest maximum value of the number of sources and the smallest excitation batch can be selected as the preferred source excitation combination from all possible source excitation combinations.
And S103, determining the excitation sequence of the preferred seismic source excitation combination.
In some embodiments of the present application, the inventors have found that when the excitation batch of the preferred source excitation combination is multiple times, the excitation priority order is: the synchronous excitation is preferred, the sliding scan is excited to the next time, and the operation efficiency is highest when the alternate scan is last.
Therefore, when the excitation batch of the preferred seismic source excitation combination is multiple times, if synchronous excitation is available, synchronous excitation operation is preferentially carried out; when there are a plurality of simultaneous excitation operations, the largest number of seismic sources in the simultaneous excitation is preferentially performed, for example, in the above combination of ABD + CE in step S102, ABD may be simultaneously excited, CE may also be simultaneously excited, since ABD is three seismic sources and CE is two seismic sources, the simultaneous excitation of ABD is preferentially performed, and then the simultaneous excitation of CE is performed. Accordingly, the ABD + CE combination corresponds to an excitation sequence ABD → CE.
And for the preferred seismic source excitation combination, when seismic sources needing to be excited individually exist in the preferred seismic source excitation combination, if excitation of the latest field in the preferred seismic source excitation combination is synchronous excitation, judging whether the distance between the seismic source needing to be excited individually and the seismic source closest to the nearest field in the group of the synchronous excitation seismic sources of the latest field meets a sliding scanning excitation condition, if so, determining the excitation mode of the seismic sources needing to be excited individually as sliding scanning excitation, otherwise, determining the excitation mode of the seismic sources needing to be excited individually as alternate scanning excitation. Taking ABC + D + E in step S102 as an example, when the distance between the D seismic source and the seismic source closest to the D seismic source in ABC satisfies the sliding scan excitation condition, and the distance between the E seismic source and the seismic source closest to the E seismic source in ABC does not satisfy the sliding scan excitation condition, the excitation sequence corresponding to the ABC + D + E combination is ABC → D → E. Conversely, if the distance between the D seismic source and the closest seismic source in ABC does not satisfy the sliding scan excitation condition, and the distance between the E seismic source and the closest seismic source in ABC satisfies the sliding scan excitation condition, the excitation sequence corresponding to the combination of ABC + D + E is ABC → E → D.
In addition, for the preferred source firing combination, the firing order between the sources may be randomly selected when there are multiple sources that can be fired in sliding sweep, and likewise, the firing order between the sources may be randomly selected when there are multiple sources that can be fired in alternating sweep.
Referring to fig. 7, an inter-source combined excitation optimization apparatus according to an embodiment of the present application may include:
a seismic source position determination module 71, configured to determine positions of a plurality of seismic sources within the target area that satisfy a preset excitation condition at a specified time point;
a firing combination determination module 72 operable to determine all possible combinations of source firing at the specified time points for the plurality of seismic sources; determining the distance between the seismic sources in each seismic source excitation combination according to the positions, and selecting one seismic source excitation combination from all the possible seismic source excitation combinations as a preferred seismic source excitation combination according to the distance;
an activation sequence determination module 73 may be used to determine the activation sequence of the preferred source activation combination.
Referring to fig. 8, another inter-source combined excitation optimization apparatus according to an embodiment of the present application may include a memory, a processor, and a computer program stored on the memory, the computer program when executed by the processor performing the steps of:
determining the positions of a plurality of seismic sources meeting preset excitation conditions in a target area when a time point is appointed;
determining all possible combinations of source shots for the plurality of sources at the specified point in time;
determining the distance between the seismic sources in each seismic source excitation combination according to the positions, and selecting one seismic source excitation combination from all the possible seismic source excitation combinations as a preferred seismic source excitation combination according to the distance;
determining an order of excitation of the preferred source excitation combination.
While the process flows described above include operations that occur in a particular order, it should be appreciated that the processes may include more or less operations that are performed sequentially or in parallel (e.g., using parallel processors or a multi-threaded environment).
For convenience of description, the above devices are described as being divided into various units by function, and are described separately. Of course, the functionality of the units may be implemented in one or more software and/or hardware when implementing the present application.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In a typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include forms of volatile memory in a computer readable medium, Random Access Memory (RAM) and/or non-volatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). Memory is an example of a computer-readable medium.
Computer-readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), Static Random Access Memory (SRAM), Dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), Read Only Memory (ROM), Electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), Digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, a computer readable medium does not include a transitory computer readable medium such as a modulated data signal and a carrier wave.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other like elements in a process, method or apparatus that comprises the element.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The application may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. The application may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the system embodiment, since it is substantially similar to the method embodiment, the description is simple, and for the relevant points, reference may be made to the partial description of the method embodiment.
The above description is only an example of the present application and is not intended to limit the present application. Various modifications and changes may occur to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the scope of the claims of the present application.
Claims (9)
1. A method for optimizing combined excitation among seismic sources is characterized by comprising the following steps:
determining the positions of a plurality of seismic sources meeting preset excitation conditions in a target area when a time point is appointed;
determining all possible source excitation combinations of the plurality of sources at the specified time point, determining distances between the sources in each source excitation combination according to the positions, and selecting one source excitation combination from all the possible source excitation combinations as a preferred source excitation combination according to the distances;
determining an order of excitation of the preferred source excitation combination; wherein the selecting one of the all possible source excitation combinations as a preferred source excitation combination according to the distance comprises:
for each seismic source excitation combination, determining the maximum value of the number of seismic sources meeting the synchronous excitation condition in the seismic source excitation combination according to the distance between the seismic sources in the seismic source excitation combination, and determining the excitation batch required by completing one-time excitation of all the seismic sources in the seismic source excitation combination;
selecting the seismic source excitation combination with the maximum value of the number of seismic sources and the minimum excitation batch as a preferred seismic source excitation combination from all possible seismic source excitation combinations at the specified time point;
wherein the determining all possible combinations of source shots for the plurality of sources at the specified point in time comprises: determining all possible combinations of source shots for the plurality of sources at the specified time point by a traversal algorithm; the traversal algorithm is as follows: and listing and obtaining results of all possible combinations of the source excitations of the plurality of sources at the specified time point one by one until all solutions are found so as to select an optimal solution from the solutions.
2. The method for inter-source compound shot optimization of claim 1, wherein the designated time points are equally spaced time points.
3. The method for optimizing excitation combining between seismic sources of claim 2, wherein the duration of the time intervals is positively correlated to the number of seismic sources available in the target area and the average time required for each seismic source to complete one excitation acquisition.
4. The method for optimizing excitation of a combination of seismic sources of claim 1, wherein the preferred combination of seismic source excitations is activated in the order: synchronous excitation, sliding scan excitation, and alternating scan excitation.
5. An optimization apparatus for combined excitation between seismic sources, comprising:
the seismic source position determining module is used for determining the positions of a plurality of seismic sources meeting preset excitation conditions in a target area when a time point is appointed;
a firing combination determination module to determine all possible combinations of firing of the seismic sources at the specified time point for the plurality of seismic sources; determining the distance between the seismic sources in each seismic source excitation combination according to the positions, and selecting one seismic source excitation combination from all the possible seismic source excitation combinations as a preferred seismic source excitation combination according to the distance;
an excitation order determination module for determining an excitation order of the preferred source excitation combination; wherein the selecting one of the all possible source excitation combinations as a preferred source excitation combination according to the distance comprises:
for each seismic source excitation combination, determining the maximum value of the number of seismic sources meeting the synchronous excitation condition in the seismic source excitation combination according to the distance between the seismic sources in the seismic source excitation combination, and determining the excitation batch required by completing one-time excitation of all the seismic sources in the seismic source excitation combination;
selecting the seismic source excitation combination with the maximum value of the number of seismic sources and the minimum excitation batch as a preferred seismic source excitation combination from all possible seismic source excitation combinations at the specified time point;
wherein the determining all possible combinations of source shots for the plurality of sources at the specified point in time comprises: determining all possible combinations of source shots for the plurality of sources at the specified time point by a traversal algorithm; the traversal algorithm is as follows: and listing and obtaining results of all possible combinations of the source excitations of the plurality of sources at the specified time point one by one until all solutions are found so as to select an optimal solution from the solutions.
6. The inter-source compound shot optimization apparatus of claim 5, wherein the designated time points are equally spaced time points.
7. The inter-source combined excitation optimization apparatus of claim 6, wherein the duration of the equal time intervals is positively correlated to the number of sources that can be placed in the target area and the average time required for each source to complete one excitation acquisition.
8. The inter-source combination excitation optimization apparatus of claim 5, wherein the preferred source excitation combinations have an excitation sequence that is, in order: synchronous excitation, sliding scan excitation, and alternating scan excitation.
9. An apparatus for optimizing the combined excitation between seismic sources, comprising a memory, a processor, and a computer program stored on the memory, wherein the computer program when executed by the processor performs the steps of:
determining the positions of a plurality of seismic sources meeting preset excitation conditions in a target area when a time point is appointed;
determining all possible source excitation combinations of the plurality of sources at the specified time point, determining distances between the sources in each source excitation combination according to the positions, and selecting one source excitation combination from all the possible source excitation combinations as a preferred source excitation combination according to the distances;
determining an order of excitation of the preferred source excitation combination; wherein the selecting one of the all possible source excitation combinations as a preferred source excitation combination according to the distance comprises:
for each seismic source excitation combination, determining the maximum value of the number of seismic sources meeting the synchronous excitation condition in the seismic source excitation combination according to the distance between the seismic sources in the seismic source excitation combination, and determining the excitation batch required by completing one-time excitation of all the seismic sources in the seismic source excitation combination;
selecting the seismic source excitation combination with the maximum value of the number of seismic sources and the minimum excitation batch as a preferred seismic source excitation combination from all possible seismic source excitation combinations at the specified time point;
wherein the determining all possible combinations of source shots for the plurality of sources at the specified point in time comprises: determining all possible combinations of source shots for the plurality of sources at the specified time point by a traversal algorithm; the traversal algorithm is as follows: and listing and obtaining results of all possible combinations of the source excitations of the plurality of sources at the specified time point one by one until all solutions are found so as to select an optimal solution from the solutions.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810371633.1A CN108761523B (en) | 2018-04-24 | 2018-04-24 | Method and device for optimizing combined excitation among seismic sources |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810371633.1A CN108761523B (en) | 2018-04-24 | 2018-04-24 | Method and device for optimizing combined excitation among seismic sources |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108761523A CN108761523A (en) | 2018-11-06 |
| CN108761523B true CN108761523B (en) | 2020-08-11 |
Family
ID=64011435
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810371633.1A Active CN108761523B (en) | 2018-04-24 | 2018-04-24 | Method and device for optimizing combined excitation among seismic sources |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108761523B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112764102B (en) * | 2019-11-06 | 2024-06-25 | 中国石油天然气集团有限公司 | Seismic data acquisition method and device |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101825724B (en) * | 2010-04-13 | 2011-12-21 | 中国石化集团华北石油局地球物理勘探公司 | Method for improving seismic prospecting distinguishability of loess tableland by combined earthquake focus excitation |
| US20120075955A1 (en) * | 2010-09-28 | 2012-03-29 | Timothy Dean | Efficient seismic source operation in connection with a seismic survey |
| CN104375165B (en) * | 2013-08-15 | 2017-02-15 | 中国石油天然气集团公司 | Controlled earthquake source partitioning type simultaneous scanning motivating method |
| CN104422951B (en) * | 2013-09-04 | 2017-02-15 | 中国石油天然气集团公司 | Controllable dynamic sliding scanning excitation method for seismic sources |
| CN104536034B (en) * | 2015-01-08 | 2017-02-22 | 吉林大学 | Multi-seismic-resource parallel motivation collecting and hybrid seismic record separating method |
-
2018
- 2018-04-24 CN CN201810371633.1A patent/CN108761523B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108761523A (en) | 2018-11-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20180004810A1 (en) | Splitting of a join operation to allow parallelization | |
| CN110460460B (en) | Service link fault positioning method, device and equipment | |
| CN104915251B (en) | Method for scheduling task and device | |
| RU2017127028A (en) | METHOD AND DEVICE FOR MONITORING AUTOMATED DRILLING PROCESSES | |
| CN106598822B (en) | A kind of abnormal deviation data examination method and device for Capacity Assessment | |
| CN111400868B (en) | Distributed workshop scheduling optimization method and system with order and robot carrying functions | |
| CN108761523B (en) | Method and device for optimizing combined excitation among seismic sources | |
| US11125898B2 (en) | Method and system for positioning seismic source in microseism monitoring | |
| CN103942209A (en) | Data processing method | |
| US9921951B2 (en) | Optimizations for regression tracking and triaging in software testing | |
| CN112905441A (en) | Test case generation method, test method, device and equipment | |
| JP2007264908A (en) | Business analysis system | |
| CN107422377B (en) | Utilize the method and device of seismic data prediction hydrocarbon source rock total content of organic carbon | |
| CN106648839A (en) | Method and device for processing data | |
| CN112764102B (en) | Seismic data acquisition method and device | |
| CN110647459A (en) | Application testing method and device | |
| CN117196850A (en) | Method and device for determining financial product combination, electronic equipment and storage medium | |
| CN111637897B (en) | Map updating method, map updating device, storage medium, and processor | |
| US10422911B2 (en) | Method and device for processing well data | |
| US20200201299A1 (en) | Method and apparatus for computer aided optimization of an occupancy of magazine slots by tools within at least one tool magazine for a machine tool | |
| CN113761746B (en) | Algebraic ring identification method and device | |
| CN106059691B (en) | Sensitivity decaying (desense) test method and device | |
| CN111435171B (en) | Clock drift correction method and device for ocean bottom node seismic data and storage medium | |
| CN106155898B (en) | Method and device for acquiring path of flow chart | |
| CN114003034A (en) | Method, equipment and storage medium for making routing inspection plan of routing inspection robot |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |