[OpendTect_Developers] Get seismic in well position

[Arnaud Huck]

Dear Fernando,

An alternate way to do this would be to define a volume statistics 
attribute in your code, that will take care of extracting the relevant 
data in a given stepout, and return you the stacked trace. Put your 
track geometry in a datapointset, set the correct attribute definition 
and then call the function

Executor* Attrib::EngineMan::getTableExtractor( DataPointSet& datapointset,
                                         const Attrib::DescSet& descset,
                                         BufferString& errmsg, int firstcol,
                                         bool needprep )


that will compute the attribute values for all points (samples) of the 
datapointset. The attribute would need to be defined like in 
plugins/uiGapDecon/uigapdeconattrib.cc  (createVolStatsDesc function).

Note 1: There is a already a class that samples the track, such that you 
would get the positions along the track at a fixed sampling rate 
(wellextractdata.h I believe). That would handle deviated wells also, 
not only vertical wells.

Note 2: The synthetic-to-seismic tie module already does such wavelet 
extraction using a Wiener-Levinson deconvolution, although without the 
lateral trace averaging. In some way, it is equivalent to input the raw 
seismic data and perform the trace averaging in the code, or input a 
laterally filtered seismic to be extracted at the nearest trace 
location. Have a look at WellAttrib for the existing code.

Best regards,
Arnaud Huck.

-- Chief Geophysicist
-- dGB Earth Sciences
-- Nijverheidstraat 11-2, 7511 JM Enschede, The Netherlands
-- mailto: arnaud.huck at dgbes.com, http://www.dgbes.com
-- Tel: +31 53 4315155 , Fax: +31 53 4315104

Le 06/02/2014 18:00, Fernando Bordignon wrote :
> Hi, I'm developing a wavelet extraction plugin based on well logs and 
> seismic amplitude, so I need to collect the nearest trace to the well 
> track. My track is vertical.
>
> I was thinking in creating a SelData by computing a small square of 
> traces (say 3x3 or 4x4) that I expect to be at the well position. Then 
> I would get the traces by a SeisTrcReader limited by this SelData. 
> With the SeisTrcInfos in hand I cal get the coordinates and calculate 
> the one that is the nearest to the well coordinates. There are any 
> other ways to do it?
>
>
> Thanks,
> Fernando Bordignon
>
>
> _____________________________________________________________
> OpendTect Developers mailing list Developers at opendtect.org
> http://lists.opendtect.org/mailman/listinfo/developers
>
> You receive this mail because you are listed on developers at opendtect.org To unsubscribe please go to http://lists.opendtect.org/mailman/options/developers . If you encounter any problems, please contact support at opendtect.org .

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opendtect.org/pipermail/developers/attachments/20140212/4133f8fa/attachment.html>