CN102121996B - Method and device for displaying seismic data collecting quality - Google Patents

Method and device for displaying seismic data collecting quality Download PDF

Info

Publication number
CN102121996B
CN102121996B CN 201010588367 CN201010588367A CN102121996B CN 102121996 B CN102121996 B CN 102121996B CN 201010588367 CN201010588367 CN 201010588367 CN 201010588367 A CN201010588367 A CN 201010588367A CN 102121996 B CN102121996 B CN 102121996B
Authority
CN
China
Prior art keywords
view
evaluation result
coordinate
dimensional image
resolution
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
Application number
CN 201010588367
Other languages
Chinese (zh)
Other versions
CN102121996A (en
Inventor
魏新建
杨午阳
徐永泽
李琳
禄娟
张建新
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Natural Gas Co Ltd
Original Assignee
China Petroleum and Natural Gas Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN 201010588367 priority Critical patent/CN102121996B/en
Publication of CN102121996A publication Critical patent/CN102121996A/en
Application granted granted Critical
Publication of CN102121996B publication Critical patent/CN102121996B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

The embodiment of the invention provides a method and a device for displaying seismic data collecting quality, wherein the method comprises the following steps of: loading satellite photo files and digital elevation files of a target work area, wherein the coordinates of the satellite photo files and the coordinates of the digital elevation files are geodetic coordinates; converting the coordinates of the satellite photo files; generating an initial three-dimensional image of the target work area according to the converted satellite photo files and digital elevation files; loading network measuring information, projecting the network measuring information on the initial three-dimensional image to obtain a three-dimensional image including the network measuring information; loading an evaluation result and projecting the evaluation result on the initial three-dimensional image to obtain the three-dimensional image including the evaluation result; and displaying the three-dimensional image including the network measuring information or the three-dimensional image including the evaluation result. The method provided by the invention can intuitively display the three-dimensional space distribution of the seismic data collecting quality to users.

Description

Seismic data collecting quality display packing and device
Technical field
The present invention relates to seismic prospecting, relate in particular to a kind of seismic data collecting quality display packing and device.
Background technology
Development along with seismic prospecting, seismic data acquisition work progressively to the earth's surface, the low signal-to-noise ratio area of inferior geological condition complex extend, the seismic exploration data amount increases, and for example, the reception channel number increases, degree of covering strengthens, bin diminishes and many components etc.Current seismic data collecting quality monitor procedure does not almost have much variations, gathers traditional state that supervision still rests on " with naked eyes, by experience ".
The inventor finds in realizing process of the present invention, and closely during the last ten years, along with the development of computer hardware technology and figure, image processing techniques, the three dimensions of geologic body and morphologic characteristics thereof reproduces and visualization technique has obtained fast development.And the seismic data collecting quality method for supervising of current main-stream still rests in the planar graph mode and shows on the evaluation result, can not provide more directly perceived and reliable scientific basis for the user.
Summary of the invention
The fundamental purpose of the embodiment of the invention is to provide a kind of seismic data collecting quality display packing and device, utilize high-resolution remote sensing image and digital elevation model (DEM), and in conjunction with the big gun evaluation result, carry out comprehensive composite diagram and the analysis of seismic data, think that the user provides more reliable scientific basis.
The above-mentioned purpose of the embodiment of the invention is achieved by the following technical solution:
A kind of seismic data collecting quality display packing, described method comprises:
The satellite photo file in loaded targets work area and digital elevation file, the coordinate of the coordinate of described satellite photo file and described digital elevation file is terrestrial coordinate;
Described satellite photo file is carried out coordinate conversion, according to the satellite photo file after the conversion and the initial three-dimensional image in the described target of described digital elevation file generated work area;
Load survey grid information, described survey grid information projection on described initial three-dimensional image, is obtained to comprise the 3-D view of survey grid information;
Load evaluation result, described evaluation result is projected on the described initial three-dimensional image, obtain to comprise the 3-D view of evaluation result;
Show described 3-D view or the described 3-D view that comprises evaluation result that comprises survey grid information, so that the acquisition quality of monitoring seismic data.
A kind of seismic data collecting quality display device, described device comprises:
The first loading unit, for satellite photo file and the digital elevation file in loaded targets work area, the coordinate of the coordinate of described satellite photo file and described digital elevation file is terrestrial coordinate;
The first processing unit is used for described satellite photo file is carried out coordinate conversion, according to the satellite photo file after the conversion and the initial three-dimensional image in the described target of described digital elevation file generated work area;
The second loading unit is used for loading survey grid information;
The second processing unit is used for described survey grid information projection obtaining to comprise the 3-D view of survey grid information to described initial three-dimensional image;
The 3rd loading unit is used for loading evaluation result;
The 3rd processing unit is used for described evaluation result is projected to described initial three-dimensional image, obtains to comprise the 3-D view of evaluation result;
Display unit is used for showing described 3-D view or the described 3-D view that comprises evaluation result that comprises survey grid information, so that the acquisition quality of monitoring seismic data.
The method and apparatus that the embodiment of the invention provides, utilize the true earth's surface information of remote sensing image and the relief information of digital elevation model (DEM), and in conjunction with the big gun evaluation result, three-dimensional spatial distribution situation from the earthquake data acquisition quality to the user that can show intuitively, the seismic data acquisition scheme of monitoring better the seismic data collecting quality in this work area and disposing next step for the user provides detailed scientific basis.
Description of drawings
Accompanying drawing described herein is used to provide a further understanding of the present invention, consists of the application's a part, does not consist of limitation of the invention.In the accompanying drawings:
Fig. 1 is the process flow diagram of the seismic data collecting quality display packing of the embodiment of the invention;
Fig. 2 a is the schematic diagram of first embodiment of the embodiment of the invention 3-D view that comprises the big gun energy;
Fig. 2 b is the schematic diagram of second embodiment of the embodiment of the invention 3-D view that comprises the big gun energy;
Fig. 2 c is the schematic diagram of the 3rd embodiment of the embodiment of the invention 3-D view that comprises the big gun energy;
Fig. 3 a is the schematic diagram that the embodiment of the invention comprises first embodiment of the 3-D view that contains the letter ratio;
Fig. 3 b is the schematic diagram that the embodiment of the invention comprises second embodiment of the 3-D view that contains the letter ratio;
Fig. 3 c is the schematic diagram that the embodiment of the invention comprises the 3rd embodiment of the 3-D view that contains the letter ratio;
Fig. 4 a is the schematic diagram of first embodiment of the embodiment of the invention 3-D view that comprises resolution;
Fig. 4 b is the schematic diagram of second embodiment of the embodiment of the invention 3-D view that comprises resolution;
Fig. 4 c is the schematic diagram of the 3rd embodiment of the embodiment of the invention 3-D view that comprises resolution;
Fig. 5 a be the embodiment of the invention comprise the big gun energy, contain letter than and the schematic diagram of first embodiment of the 3-D view of resolution;
Fig. 5 b be the embodiment of the invention comprise the big gun energy, contain letter than and the schematic diagram of second embodiment of the 3-D view of resolution;
Fig. 5 c be the embodiment of the invention comprise the big gun energy, contain letter than and the schematic diagram of the 3rd embodiment of the 3-D view of resolution;
Fig. 6 is the schematic diagram that the embodiment of the invention comprises the 3-D view of common midpoint coverage number of times;
Fig. 7 is the composition frame chart of the seismic data collecting quality display device of the embodiment of the invention.
Embodiment
For the purpose, technical scheme and the advantage that make the embodiment of the invention is clearer, below in conjunction with embodiment and accompanying drawing, the embodiment of the invention is described in further details.At this, illustrative examples of the present invention and explanation thereof are used for explanation the present invention, but not as a limitation of the invention.
The process flow diagram of a kind of seismic data collecting quality display packing that Fig. 1 provides for the embodiment of the invention please refer to Fig. 1, and the method comprises:
Step 101: the satellite photo file in loaded targets work area and digital elevation file, the coordinate of the coordinate of described satellite photo file and described digital elevation file is terrestrial coordinate;
Step 102: described satellite photo file is carried out coordinate conversion, according to the satellite photo file after the conversion and the initial three-dimensional image in the described target of described digital elevation file generated work area;
Step 103: load survey grid information, described survey grid information projection on described initial three-dimensional image, is obtained to comprise the 3-D view of survey grid information;
Step 104: load evaluation result, described evaluation result is projected on the described initial three-dimensional image, obtain to comprise the 3-D view of evaluation result;
Step 105: show described 3-D view or the described 3-D view that comprises evaluation result that comprises survey grid information, so that the acquisition quality of monitoring seismic data.
In the step 101 of the present embodiment, the satellite photo file including that loads the earth's surface information in target work area, the digital elevation file including that loads the relief information in target work area, these information provide foundation for the three-dimensional spatial distribution of user study earthquake data acquisition quality.
Wherein, the earth's surface packets of information contains lithology, environment etc.
In the step 102 of the present embodiment, because the coordinate of satellite photo file is terrestrial coordinate, numeral is larger, if press the actual coordinate drawing image, the scope of image is very large, easily causes into the low and slow-footed shortcoming of interactive operation of figure efficient.For this reason, the embodiment of the invention was carried out coordinate conversion to it before the satellite photo imaging.
In one embodiment, can carry out coordinate conversion according to following formula:
MPx = ( SL x - SL min x ) * ( MP max x - MP min x ) ( SL max x - SL min x ) ;
MPy = ( SL y - SL min y ) * ( MP max y - MP min y ) ( SL max y - SL min y ) ;
MPz = ( SL z - SL min z ) * ( MP max z - MP min z ) ( SL max z - SL min z ) ;
Wherein, MPx, MPy, MPz be for will draw coordinate x, y, the value of z, MPmax x, MPmax y, MPmax zFor drawing coordinate x, y, z maximal value, MPmin x, MPmin y, MPmin zFor drawing coordinate x, y, z minimum value, SL x, SL y, SL zBe the coordinate figure of satellite photo file, SLmax x, SLmax y, SLmax zBe the maximal value in the coordinate figure of satellite photo file, SLmin x, SLmin y, SLmin zBe the minimum value in the coordinate figure of satellite photo file.
In the step 103 of the present embodiment, the survey grid information projection that loads has been arrived on the 3-D view of step 102 generation.The survey grid information has here comprised the position of shot point, the position of geophone station etc., the survey grid information projection is here arrived on the 3-D view (initial three-dimensional image) that generates before, and physical location and the distribution in three dimensions provides foundation for user study shot point and geophone station.
In the step 104 of the present embodiment, the evaluation result that loads has been projected on the 3-D view of step 102 generation.The evaluation result has here comprised the big gun energy, has contained letter ratio and resolution etc., on the 3-D view that generates before the evaluation result here projected to, for the acquisition quality of these geological datas of user study and the earthquake-capturing scheme of disposing next step provide foundation.
In one embodiment, can with the big gun energy in this evaluation result, contain letter than and the 3-D view that generates before projecting to respectively of resolution on one-tenth figure, also this three can be merged one-tenth figure on the 3-D view that generates before projecting to.Becoming chart-pattern can be spot style, grid type or superposed type, has also namely comprised the figure that becomes of loose point and grid.
Fig. 2 a is for by spot style, with the big gun energy projection in the evaluation result to the 3-D view that generates before, the schematic diagram of the 3-D view that comprises the big gun energy of acquisition.
Fig. 2 b is for by grid type, with the big gun energy projection in the evaluation result to the 3-D view that generates before, the schematic diagram of the 3-D view that comprises the big gun energy of acquisition.
Fig. 2 c projects to the resolution in the evaluation result on the 3-D view that generates before the schematic diagram of the 3-D view that comprises resolution of acquisition for by superposed type.
Fig. 3 a is for by spot style, on the 3-D view that contains letter than generation before projecting in the evaluation result, and the schematic diagram that comprises the 3-D view that contains the letter ratio of acquisition.
Fig. 3 b is for by grid type, on the 3-D view that contains letter than generation before projecting in the evaluation result, and the schematic diagram that comprises the 3-D view that contains the letter ratio of acquisition.
Fig. 3 c is for by superposed type, on the 3-D view that contains letter than generation before projecting in the evaluation result, and the schematic diagram that comprises the 3-D view that contains the letter ratio of acquisition.
Fig. 4 a projects to the resolution in the evaluation result on the 3-D view that generates before the schematic diagram of the 3-D view that comprises resolution of acquisition for by spot style.
Fig. 4 b projects to the resolution in the evaluation result on the 3-D view that generates before the schematic diagram of the 3-D view that comprises resolution of acquisition for by grid type.
Fig. 4 c projects to the resolution in the evaluation result on the 3-D view that generates before the schematic diagram of the 3-D view that comprises resolution of acquisition for by superposed type.
Fig. 5 a is for by spot style, with the big gun energy in the evaluation result, contain letter than and the 3-D view that generates before projecting to of resolution on, acquisition comprise the big gun energy, contain letter than and the schematic diagram of the 3-D view of resolution.
Fig. 5 b is for by grid type, with the big gun energy in the evaluation result, contain letter than and the 3-D view that generates before projecting to of resolution on, acquisition comprise the big gun energy, contain letter than and the schematic diagram of the 3-D view of resolution.
Fig. 5 c is for by superposed type, with the big gun energy in the evaluation result, contain letter than and the 3-D view that generates before projecting to of resolution on, acquisition comprise the big gun energy, contain letter than and the schematic diagram of the 3-D view of resolution.
In the step 104 of the present embodiment, passing through spot style, grid type or superposed type, with the big gun energy in the evaluation result, contain letter than and the 3-D view that generates before projecting to of resolution on the time, at first will with this big gun energy, contain letter than and resolution merge; Then determine the minimum circumscribed rectangular region of all shot points that participation is estimated; Last again with the evaluation result data projection that merges to the 3-D view that generates before.
In the present embodiment, can be according to formula
Figure BDA0000038201140000051
With the big gun energy, contain letter than and resolution merge.Wherein, SN be shot point contain the letter ratio, The mean value that contains the letter ratio of described all shot points, λ SNIt is the fusion weight that contains the letter ratio; SE is the big gun energy of shot point,
Figure BDA0000038201140000053
The mean value of the big gun energy of described all shot points, λ SEIt is the fusion weight of big gun energy; RES is the resolution of shot point,
Figure BDA0000038201140000054
The resolution mean value of described all shot points, λ RESIt is the fusion weight of resolution; K is calibration factor.
In the present embodiment, because shot point is loose point, and its coordinate has certain angle with the terrestrial coordinate axle, so the present embodiment specifically is used for determining and the terrestrial coordinate axle has the minimum boundary rectangle of the loose point of certain angle.
In one embodiment, can be first according to formula
Figure BDA0000038201140000055
The coordinate of described all shot points is rotated, makes it be parallel to the terrestrial coordinate axle; Again with the minimum value of the transverse and longitudinal coordinate of described all shot point coordinates and the maximal value of transverse and longitudinal coordinate, the also i.e. minimum x of all shot point coordinates, y value and maximum x, the y value is as the four angular coordinate of the minimum circumscribed rectangular region at described all shot point places; Contrary formula by above formula rotates to original state with described four angular coordinate at last.Determine the minimum circumscribed rectangular region at all shot point places that participation is estimated with this.Wherein, Rotate x, Rotate yBe respectively postrotational coordinate, α is the angle of shot point and coordinate axis forward.
In the present embodiment, when the evaluation result data of this fusion are carried out grid type one-tenth figure, can participate in first the minimum circumscribed rectangular region of all shot points of evaluation according to the target work area, derive Continuous Rectangular zone; Then according to derivative Continuous Rectangular zone, with the evaluation result data projection that merges on described initial three-dimensional image.
In one embodiment, can pass through formula
Figure BDA0000038201140000061
Derive Continuous Rectangular zone according to minimum circumscribed rectangular region, wherein,
Figure BDA0000038201140000062
Be S 0The predicted value at place; N is the quantity of sampling point around the future position that will use in the prediction and calculation process; λ iBe the weight of each sampling point that will use in the prediction and calculation process, this value reduces along with the increase of distance between sampling point and the future position; Z (S i) be at S iThe measured value that the place obtains.Wherein, can be according to formula
Figure BDA0000038201140000063
Determine weight, here, d I0Future position S 0With known sampling point S iBetween distance, along with d I0Increase, sampling point also reduces simultaneously on the weight of future position impact.Wherein, N is for satisfying and S 0Distance less than the sampling point number of r, r is given search radius.
In the step 105 of the present embodiment, can be according to the 3-D view that comprises survey grid information that generates before and the 3-D view that comprises evaluation result, show to the user respectively, for example show that this comprises the 3-D view of survey grid information, perhaps show the 3-D view that comprises evaluation result, so that the user is according to the seismic data collecting quality in this target work area of this picture control and the earthquake-capturing scheme of disposing next step.
On the basis of seismic data collecting quality display packing shown in Figure 1, the embodiment of the invention can also further load the common midpoint coverage number of times, and this common midpoint coverage number of times is projected on the 3-D view of before generation, think that the user provides the reference of further seismic data collecting quality.
Fig. 6 is the schematic diagram that the present embodiment projects to the common midpoint coverage number of times 3-D view that generates before.In the present embodiment, step 105 can show this 3-D view that comprises the common midpoint coverage number of times and the 3-D view that comprises survey grid information simultaneously, and step 105 also can show simultaneously this 3-D view that comprises the common midpoint coverage number of times and comprise the 3-D view of evaluation result.
In the present embodiment, for obtaining of the obtaining of the obtaining of the obtaining of satellite photo file, digital elevation file, survey grid information, evaluation result etc., can realizing by existing means, is not the emphasis that the present invention studies, and does not repeat them here.
The seismic data collecting quality display packing of the embodiment of the invention, with earthquake-capturing information such as survey grid information, big gun evaluation result, common center degree of covering, organically combine with the true earth's surface information (photo files and elevation file obtain via satellite) of target work area remote sensing image, and show with 3-D view, the detailed information of the more closer to reality situation by the exploration targets work area is provided to the user, promote it further to improve the acquisition quality monitoring level, carry out the seismic exploration decision and deployment.
The composition frame chart of a kind of seismic data collecting quality display device that Fig. 7 provides for the embodiment of the invention please refer to Fig. 7, and this device comprises:
The first loading unit 71, for satellite photo file and the digital elevation file in loaded targets work area, the coordinate of the coordinate of satellite photo file and digital elevation file is terrestrial coordinate;
The first processing unit 72 is used for the satellite photo file is carried out coordinate conversion, according to the initial three-dimensional image in the digital elevation file generated target work area of the satellite photo file after the conversion and loading;
The second loading unit 73 is used for loading survey grid information;
The second processing unit 74 is used for the survey grid information projection obtaining to comprise the 3-D view of survey grid information to initial three-dimensional image;
The 3rd loading unit 75 is used for loading evaluation result;
The 3rd processing unit 76 is used for evaluation result is projected to initial three-dimensional image, obtains to comprise the 3-D view of evaluation result;
Display unit 77 is used for showing the 3-D view that comprises the 3-D view of survey grid information or comprise evaluation result, so that the acquisition quality of monitoring seismic data.
In one embodiment, the first processing unit 72 is concrete for according to following formula the satellite photo file being carried out coordinate conversion:
MPx = ( SL x - SL min x ) * ( MP max x - MP min x ) ( SL max x - SL min x ) ;
MPy = ( SL y - SL min y ) * ( MP max y - MP min y ) ( SL max y - SL min y ) ;
MPz = ( SL z - SL min z ) * ( MP max z - MP min z ) ( SL max z - SL min z ) ;
Wherein, MPx, MPy, MPz be for will draw coordinate x, y, the value of z, MPmax x, MPmax y, MPmax zFor drawing coordinate x, y, z maximal value, MPmin x, MPmin y, MPmin zFor drawing coordinate x, y, z minimum value, SL x, SL y, SL zBe the coordinate figure of satellite photo file, SLmax x, SLmax y, SLmax zBe the maximal value in the coordinate figure of satellite photo file, SLmin x, SLmin y, SLmin zBe the minimum value in the coordinate figure of satellite photo file.
In one embodiment, evaluation result comprises the big gun energy, contains letter ratio and resolution, and the 3rd processing unit 76 is concrete for by spot style, grid type or superposed type, and the big gun energy projection on initial three-dimensional image, is obtained to comprise the 3-D view of big gun energy; Perhaps, by spot style, grid type or superposed type, will contain the letter ratio and project on the initial three-dimensional image, obtain to comprise the 3-D view that contains the letter ratio; Perhaps, by spot style, grid type or superposed type, resolution is projected on the initial three-dimensional image, obtain to comprise the 3-D view of resolution; Perhaps, by spot style, grid type or superposed type, with the big gun energy, contain letter than and resolution project on the initial three-dimensional image, obtain to comprise the big gun energy, contain letter than and the 3-D view of resolution.
In one embodiment, the 3rd processing unit 76 comprises:
Fusion Module, be used for the big gun energy, contain letter than and resolution merge, obtain the evaluation result data that merge;
Determination module is used for definite minimum circumscribed rectangular region that participates in all shot point places of evaluation;
Projection module is used for according to the evaluation result data that merge and definite minimum circumscribed rectangular region, by grid type, the evaluation result data projection that merges on initial three-dimensional image, is obtained to comprise the 3-D view of the evaluation result data of fusion.
Wherein, Fusion Module specifically is used for: according to formula
Figure BDA0000038201140000081
With the big gun energy, contain letter than and resolution merge.Wherein, SN be shot point contain the letter ratio,
Figure BDA0000038201140000082
The mean value that contains the letter ratio of all shot points, λ SNIt is the fusion weight that contains the letter ratio; SE is the big gun energy of shot point,
Figure BDA0000038201140000083
The mean value of the big gun energy of all shot points, λ SEIt is the fusion weight of big gun energy; RES is the resolution of shot point, The resolution mean value of all shot points, λ RESIt is the fusion weight of resolution; K is calibration factor.
Wherein, determination module specifically is used for: according to formula The coordinate of all shot points is rotated, make it be parallel to the terrestrial coordinate axle, with the maximal value of the minimum value of the transverse and longitudinal coordinate of all shot point coordinates and the transverse and longitudinal coordinate four angular coordinate as the minimum circumscribed rectangular region at all shot point places, the contrary formula by above formula rotates to original state with four angular coordinate.
Wherein, projection module specifically is used for: derive Continuous Rectangular zone according to minimum circumscribed rectangular region, according to derivative Continuous Rectangular zone, by grid type with the evaluation result data projection that merges on initial three-dimensional image.
In one embodiment, projection module specifically is used for: pass through formula
Figure BDA0000038201140000086
Derive Continuous Rectangular zone according to minimum circumscribed rectangular region.Wherein, Be S 0The predicted value at place; λ iWeight for each sampling point that will use in the prediction and calculation process; Z (S i) be at S iThe measured value that the place obtains;
Figure BDA0000038201140000088
d I0Future position S 0With known sampling point S iBetween distance, d I0Future position S 0With known sampling point S jBetween distance; N is for satisfying and S 0Distance less than the sampling point number of r, r is given search radius.
In one embodiment, this device can also comprise:
The 4th loading unit 78 is used for loading the common midpoint coverage number of times;
The manages unit 79 everywhere, is used for this common midpoint coverage number of times is projected to initial three-dimensional image, obtains to comprise the 3-D view of common midpoint coverage number of times.
In the present embodiment, display unit 77 specifically is used for: show to comprise the 3-D view of common midpoint coverage number of times, and comprise the 3-D view of survey grid information or comprise the 3-D view of evaluation result.
Each ingredient of the device of the present embodiment is respectively applied to realize each step of preceding method embodiment, owing in preceding method embodiment, each step is had been described in detail, does not repeat them here.
The device of the embodiment of the invention by establishing target work area three-dimensional visualization image, reproduces the three dimensions feature in target work area, can realize on the whole more directly perceived, more fully comprehensive evaluation analysis research is carried out in the target work area.
The method of describing in conjunction with embodiment disclosed herein or the step of algorithm can directly use the software module of hardware, processor execution, and perhaps the combination of the two is implemented.Software module can place the storage medium of any other form known in random access memory (RAM), internal memory, ROM (read-only memory) (ROM), electrically programmable ROM, electrically erasable ROM, register, hard disk, moveable magnetic disc, CD-ROM or the technical field.
Above-described specific embodiment; purpose of the present invention, technical scheme and beneficial effect are further described; institute is understood that; the above only is specific embodiments of the invention; the protection domain that is not intended to limit the present invention; within the spirit and principles in the present invention all, any modification of making, be equal to replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (14)

1. a seismic data collecting quality display packing is characterized in that, described method comprises:
The satellite photo file in loaded targets work area and digital elevation file, the coordinate of the coordinate of described satellite photo file and described digital elevation file is terrestrial coordinate;
Described satellite photo file is carried out coordinate conversion, according to the satellite photo file after the conversion and the initial three-dimensional image in the described target of described digital elevation file generated work area;
Load survey grid information, described survey grid information projection on described initial three-dimensional image, is obtained to comprise the 3-D view of survey grid information;
Load evaluation result, described evaluation result is projected on the described initial three-dimensional image, obtain to comprise the 3-D view of evaluation result;
Show described 3-D view or the described 3-D view that comprises evaluation result that comprises survey grid information, so that the acquisition quality of monitoring seismic data;
Described satellite photo file is carried out coordinate conversion, comprising:
According to following formula the satellite photo file is carried out coordinate conversion:
MPx = ( SL x - SL min x ) * ( MP max x - MP min x ) ( SL max x - SL min x ) ;
MPy = ( SL y - SL min y ) * ( MP max y - MP min y ) ( SL max y - SL min y ) ;
MPz = ( SL z - SL min z ) * ( MP max z - MP min z ) ( SL max z - SL min z ) ;
Wherein, MPx, MPy, MPz be for will draw coordinate x, y, the value of z, MPmax x, MPmax y, MPmax zFor drawing coordinate x, y, z maximal value, MPmin x, MPmin y, MPmin zFor drawing coordinate x, y, z minimum value, SL x, SL y, SL zBe the coordinate figure of satellite photo file, SLmax x, SLmax y, SLmax zBe the maximal value in the coordinate figure of satellite photo file, SLmin x, SLmin y, SLmin zBe the minimum value in the coordinate figure of satellite photo file;
Described evaluation result comprises the big gun energy, contain letter than and resolution, described described evaluation result is projected on the described initial three-dimensional image, obtain to comprise the 3-D view of evaluation result, comprising:
By spot style, grid type or superposed type, described big gun energy projection on described initial three-dimensional image, is obtained to comprise the 3-D view of big gun energy; Perhaps
By spot style, grid type or superposed type, the described letter ratio that contains is projected on the described initial three-dimensional image, obtain to comprise the 3-D view that contains the letter ratio; Perhaps
By spot style, grid type or superposed type, described resolution is projected on the described initial three-dimensional image, obtain to comprise the 3-D view of resolution; Perhaps
By spot style, grid type or superposed type, with described big gun energy, contain letter than and resolution project on the described initial three-dimensional image, obtain to comprise the big gun energy, contain letter than and the 3-D view of resolution.
2. method according to claim 1 is characterized in that, the described grid type that passes through projects to described evaluation result on the described initial three-dimensional image, and acquisition comprises the big gun energy, contains the 3-D view of letter ratio and resolution, comprising:
With described big gun energy, contain letter than and resolution merge the evaluation result data that to merge;
Determine the minimum circumscribed rectangular region at all shot point places that participation is estimated;
According to evaluation result data and the described minimum circumscribed rectangular region of described fusion, by grid type, the evaluation result data projection of described fusion on described initial three-dimensional image, is obtained to comprise the 3-D view of the evaluation result data of described fusion.
3. method according to claim 2 is characterized in that, described with described big gun energy, contain letter than and resolution merge, be specially:
According to following formula with described big gun energy, contain letter than and resolution merge:
QC = ( SN SN ‾ × λ SN + SE SE ‾ × λ SE + RES RES ‾ × λ RES ) + k ;
Wherein, SN be shot point contain the letter ratio,
Figure FDA00001879457000022
The mean value that contains the letter ratio of described all shot points, λ SNIt is the fusion weight that contains the letter ratio; SE is the big gun energy of shot point,
Figure FDA00001879457000023
The mean value of the big gun energy of described all shot points, λ SEIt is the fusion weight of big gun energy; RES is the resolution of shot point, The resolution mean value of described all shot points, λ RESIt is the fusion weight of resolution; K is calibration factor.
4. method according to claim 2 is characterized in that, the minimum circumscribed rectangular region at all shot point places that described definite participation is estimated comprises:
Be rotated according to the coordinate of following formula with described all shot points, make it be parallel to the terrestrial coordinate axle:
Rotate x=x*cosa-y*sina
Rotate y=-x*sina+y*cosa
With the maximal value of the minimum value of the transverse and longitudinal coordinate of described all shot point coordinates and the transverse and longitudinal coordinate four angular coordinate as the minimum circumscribed rectangular region at described all shot point places;
Contrary formula by above formula rotates to original state with described four angular coordinate.
5. method according to claim 2 is characterized in that, according to evaluation result data and the described minimum circumscribed rectangular region of described fusion, by grid type, the evaluation result data projection of described fusion on described initial three-dimensional image, being comprised:
Derive Continuous Rectangular zone according to described minimum circumscribed rectangular region;
According to described derivative Continuous Rectangular zone, by grid type with the evaluation result data projection of described fusion on described initial three-dimensional image.
6. method according to claim 5 is characterized in that, derives Continuous Rectangular zone according to described minimum circumscribed rectangular region, is specially:
By following formula, derive Continuous Rectangular zone according to described minimum circumscribed rectangular region:
Z ^ ( S 0 ) = Σ i = 1 N λ i Z ( s i ) ;
Wherein,
Figure FDA00001879457000032
Be S 0The predicted value at place; λ iWeight for each sampling point that will use in the prediction and calculation process; Z (s i) be at S iThe measured value that the place obtains;
Wherein, d I0Future position S 0With known sampling point S iBetween distance, d J0Future position S 0With known sampling point S jBetween distance;
Wherein, N is for satisfying and S 0Distance less than the sampling point number of r, r is given search radius.
7. method according to claim 1 is characterized in that, described acquisition comprises after the 3-D view of survey grid information and evaluation result, also comprises:
Load the common midpoint coverage number of times, described common midpoint coverage number of times is projected on the described initial three-dimensional image, obtain to comprise the 3-D view of common midpoint coverage number of times;
Wherein, described 3-D view or the described 3-D view that comprises evaluation result that comprises survey grid information of described demonstration, be specially: show the described 3-D view that comprises the common midpoint coverage number of times, and the described 3-D view that comprises the 3-D view of survey grid information or comprise evaluation result.
8. a seismic data collecting quality display device is characterized in that, described device comprises:
The first loading unit, for satellite photo file and the digital elevation file in loaded targets work area, the coordinate of the coordinate of described satellite photo file and described digital elevation file is terrestrial coordinate;
The first processing unit is used for described satellite photo file is carried out coordinate conversion, according to the satellite photo file after the conversion and the initial three-dimensional image in the described target of described digital elevation file generated work area;
The second loading unit is used for loading survey grid information;
The second processing unit is used for described survey grid information projection obtaining to comprise the 3-D view of survey grid information to described initial three-dimensional image;
The 3rd loading unit is used for loading evaluation result;
The 3rd processing unit is used for described evaluation result is projected to described initial three-dimensional image, obtains to comprise the 3-D view of evaluation result;
Display unit is used for showing described 3-D view or the described 3-D view that comprises evaluation result that comprises survey grid information, so that the acquisition quality of monitoring seismic data;
Described the first processing unit specifically is used for according to following formula the satellite photo file being carried out coordinate conversion:
MPx = ( SL x - SL min x ) * ( MP max x - MP min x ) ( SL max x - SL min x ) ;
MPy = ( SL y - SL min y ) * ( MP max y - MP min y ) ( SL max y - SL min y ) ;
MPz = ( SL z - SL min z ) * ( MP max z - MP min z ) ( SL max z - SL min z ) ;
Wherein, MPx, MPy, MPz be for will draw coordinate x, y, the value of z, MPmax x, MPmax y, MPmax zFor drawing coordinate x, y, z maximal value, MPmin x, MPmin y, MPmin zFor drawing coordinate x, y, z minimum value, SL x, SL y, SL zBe the coordinate figure of satellite photo file, SLmax x, SLmax y, SLmax zBe the maximal value in the coordinate figure of satellite photo file, SLmin x, SLmin y, SLmin zBe the minimum value in the coordinate figure of satellite photo file;
Described evaluation result comprises the big gun energy, contains letter ratio and resolution, and described the first processing unit specifically is used for by spot style, grid type or superposed type, and described big gun energy projection on described initial three-dimensional image, is obtained to comprise the 3-D view of big gun energy; Perhaps, by spot style, grid type or superposed type, the described letter ratio that contains is projected on the described initial three-dimensional image, obtain to comprise the 3-D view that contains the letter ratio; Perhaps, by spot style, grid type or superposed type, described resolution is projected on the described initial three-dimensional image, obtain to comprise the 3-D view of resolution; Perhaps, by spot style, grid type or superposed type, with described big gun energy, contain letter than and resolution project on the described initial three-dimensional image, obtain to comprise the big gun energy, contain letter than and the 3-D view of resolution.
9. device according to claim 8 is characterized in that, it is characterized in that, described the first processing unit comprises:
Fusion Module, be used for described big gun energy, contain letter than and resolution merge, obtain the evaluation result data that merge;
Determination module is used for definite minimum circumscribed rectangular region that participates in all shot point places of evaluation;
Projection module, be used for evaluation result data and described minimum circumscribed rectangular region according to described fusion, by grid type, the evaluation result data projection of described fusion on described initial three-dimensional image, is obtained to comprise the 3-D view of the evaluation result data of described fusion.
10. device according to claim 9 is characterized in that, described Fusion Module specifically is used for:
According to formula QC = ( SN SN ‾ × λ SN + SE SE ‾ × λ SE + RES RES ‾ × λ RES ) + k With described big gun energy, contain letter than and resolution merge;
Wherein, SN be shot point contain the letter ratio,
Figure FDA00001879457000051
The mean value that contains the letter ratio of described all shot points, λ SNIt is the fusion weight that contains the letter ratio; SE is the big gun energy of shot point,
Figure FDA00001879457000052
The mean value of the big gun energy of described all shot points, λ SEIt is the fusion weight of big gun energy; RES is the resolution of shot point,
Figure FDA00001879457000053
The resolution mean value of described all shot points, λ RESIt is the fusion weight of resolution; K is calibration factor.
11. device according to claim 9 is characterized in that, described determination module specifically is used for:
According to formula Rotate x = x * cos a - y * sin a Rotate y = - x * sin a + y * cos a The coordinate of described all shot points is rotated, make it be parallel to the terrestrial coordinate axle, with the maximal value of the minimum value of the transverse and longitudinal coordinate of described all shot point coordinates and the transverse and longitudinal coordinate four angular coordinate as the minimum circumscribed rectangular region at described all shot point places, the contrary formula by above formula rotates to original state with described four angular coordinate.
12. device according to claim 9 is characterized in that, described projection module specifically is used for:
Derive Continuous Rectangular zone according to described minimum circumscribed rectangular region, according to described derivative Continuous Rectangular zone, by grid type with the evaluation result data projection of described fusion on described initial three-dimensional image.
13. device according to claim 12 is characterized in that, described projection module specifically is used for:
Pass through formula
Figure FDA00001879457000055
Derive Continuous Rectangular zone according to described minimum circumscribed rectangular region;
Wherein,
Figure FDA00001879457000056
Be S 0The predicted value at place; λ iWeight for each sampling point that will use in the prediction and calculation process; Z (s i) be at S iThe measured value that the place obtains;
Figure FDA00001879457000057
d I0Future position S 0With known sampling point S iBetween distance, d J0Future position S 0With known sampling point S jBetween distance; N is for satisfying and S 0Distance less than the sampling point number of r, r is given search radius.
14. device according to claim 8 is characterized in that, described device also comprises:
The 4th loading unit is used for loading the common midpoint coverage number of times;
The manages the unit everywhere, is used for described common midpoint coverage number of times is projected to described initial three-dimensional image, obtains to comprise the 3-D view of common midpoint coverage number of times;
Wherein, described display unit specifically is used for: show the described 3-D view that comprises the common midpoint coverage number of times, and described 3-D view or the described 3-D view that comprises evaluation result that comprises survey grid information.
CN 201010588367 2010-12-14 2010-12-14 Method and device for displaying seismic data collecting quality Active CN102121996B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201010588367 CN102121996B (en) 2010-12-14 2010-12-14 Method and device for displaying seismic data collecting quality

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201010588367 CN102121996B (en) 2010-12-14 2010-12-14 Method and device for displaying seismic data collecting quality

Publications (2)

Publication Number Publication Date
CN102121996A CN102121996A (en) 2011-07-13
CN102121996B true CN102121996B (en) 2013-05-01

Family

ID=44250589

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201010588367 Active CN102121996B (en) 2010-12-14 2010-12-14 Method and device for displaying seismic data collecting quality

Country Status (1)

Country Link
CN (1) CN102121996B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA3118864C (en) 2012-03-08 2022-11-29 Shell Internationale Research Maatschappij B.V. Integrated seismic monitoring system and method
CN104280764A (en) * 2013-07-01 2015-01-14 中国石油化工股份有限公司 Seismic exploration auxiliary data quality control method based on three-dimensional GIS
CN104318098A (en) * 2014-10-17 2015-01-28 中国石油天然气股份有限公司 Remote sensing based method for screening three-dimensional seismic receiving combination central point
CN107664778B (en) * 2017-08-24 2019-11-08 中国石油天然气股份有限公司 The determination method and apparatus of earthquake data quality
CN109031418A (en) * 2018-07-17 2018-12-18 中国石油天然气股份有限公司 A kind of earthquake collection quality analysis method and device
CN113945965A (en) * 2020-07-17 2022-01-18 中国石油天然气股份有限公司 Seismic data monitoring method and device
CN113989454B (en) * 2021-12-30 2022-03-18 山东省地质测绘院 Fusion method, device and system suitable for geological data and geographic information data

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU562787A1 (en) * 1975-07-10 1977-06-25 Всесоюзный Научно-Исследовательский Институт Геофизических Методов Разведки Seismic imaging device
CN101839998A (en) * 2009-03-18 2010-09-22 中国石油天然气集团公司 High precision prestack depth migration method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU562787A1 (en) * 1975-07-10 1977-06-25 Всесоюзный Научно-Исследовательский Институт Геофизических Методов Разведки Seismic imaging device
CN101839998A (en) * 2009-03-18 2010-09-22 中国石油天然气集团公司 High precision prestack depth migration method

Also Published As

Publication number Publication date
CN102121996A (en) 2011-07-13

Similar Documents

Publication Publication Date Title
CN102121996B (en) Method and device for displaying seismic data collecting quality
Liu et al. The community velocity model V. 1.0 of southwest China, constructed from joint body‐and surface‐wave travel‐time tomography
Colaninno et al. Quantitative comparison of methods for predicting the arrival of coronal mass ejections at Earth based on multiview imaging
Glauert et al. Three‐dimensional electron radiation belt simulations using the BAS Radiation Belt Model with new diffusion models for chorus, plasmaspheric hiss, and lightning‐generated whistlers
Lillis et al. Three‐dimensional structure of the Martian nightside ionosphere: Predicted rates of impact ionization from Mars Global Surveyor magnetometer and electron reflectometer measurements of precipitating electrons
Hassan et al. Monte Carlo performance studies for the site selection of the Cherenkov Telescope Array
Berat et al. Full simulation of space-based extensive air showers detectors with ESAF
Bai et al. 2-D/3-D irregular shortest-path ray tracing for multiple arrivals and its applications
Rahmati et al. Seasonal variability of neutral escape from Mars as derived from MAVEN pickup ion observations
Maguya et al. Adaptive algorithm for large scale dtm interpolation from lidar data for forestry applications in steep forested terrain
Vallée et al. Imaging coseismic rupture in far field by slip patches
Liang et al. 3‐D inversion of gravity data in spherical coordinates with application to the GRAIL data
Yang et al. Time‐dependent MHD modeling of the global solar corona for year 2007: Driven by daily‐updated magnetic field synoptic data
Gold et al. A terrestrial lidar-based workflow for determining three–dimensional slip vectors and associated uncertainties
Nordström et al. Venus ion outflow estimates at solar minimum: Influence of reference frames and disturbed solar wind conditions
Liu et al. Constraints on time-dependent volcanic source models at Long Valley Caldera from 1996 to 2009 using InSAR and geodetic measurements
Bergman et al. The influence of varying spacecraft potentials and Debye lengths on in situ low‐energy ion measurements
Twardzik et al. Robust features of the source process for the 2004 Parkfield, California, earthquake from strong-motion seismograms
Chaput et al. Imaging of Erebus volcano using body wave seismic interferometry of Strombolian eruption coda
Yuan et al. Seismic source tracking with six degree‐of‐freedom ground motion observations
Li et al. Source process featuring asymmetric rupture velocities of the 2021 Mw 7.4 Maduo, China, earthquake from teleseismic and geodetic data
Engel et al. Simulations of inner radiation belt proton loss during geomagnetic storms
Rodgers et al. WUS256: An adjoint waveform tomography model of the crust and upper mantle of the western United States for improved waveform simulations
Schoorlemmer et al. Characteristics of extensive air showers around the energy threshold for ground-particle-based γ-ray observatories
Kim et al. S‐velocity mantle structure of East Asia from teleseismic traveltime tomography: Inferred mechanisms for the Cenozoic intraplate volcanoes

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant