واهماميخت ناپايا با استفاده از بهينه سازي تصادفي تبلور شبيه سازي شده

Non-stationary deconvolution via stochastic simulated annealing optimization

واهماميخت يكي از مراحل مهم پردازش سيگنال هاي لرزه اي براي به دست آوردن تصاوير با تفكيك پذيري بالا از زير سطح زمين از جمله تصاوير مخازن هيدروكربن است. در اين مقاله فرآيند واهماميخت ناپايا به منظور تخمين دامنه و زمان ضرايب غير صفر سري بازتاب و همچنين تخمين فاكتور كيفيت لايه ها با فرض وجود محيط ميرا كننده براي انتشار موجك چشمه با استفاده از روشي آماري انجام شده است. بازيابي محل دقيق بازتاب ها و حذف اثر ميرايي زمين با استفاده از الگوريتم تبلور شبيه سازي شده صورت گرفته و فاكتور كيفيت براي هر لايه تخمين زده مي شود. همچنين دامنه سري ضرايب بازتاب نيز با روش كمترين مربعات به دقت به دست مي آيد. در روش پيشنهادي، هر سه پارامتر دامنه، زمان ضرايب و فاكتور كيفيت مربوط به هر لايه به صورت هم زمان تخمين زده مي شود كه ميتوان آن را مزيتي نسبت به ساير روش ها بيان كرد. از الگوريتم تبلور شبيه سازي شده كه الگوريتمي آماري و بر پايه ي تكرار است؛ براي حل مسئله ي مورد نظر استفاده شده است. نتايج حاصل از پياده سازي اين الگوريتم بر داده هاي مصنوعي و ميداني، دقت عملكرد آن را به خوبي نشان مي دهد.

Deconvolution is one of the significant steps in seismic signal processing for obtaining high resolution images of subsurface. In this paper, a new nonstationary deconvolution is developed to retrieve the impulse response of the earth and remove the earth attenuation effects from it via a simulated annealing (SA) optimization. The time location of each spike and the corresponding quality factor is estimated by the SA method and its amplitude is calculated using a least-squares method. All three parameters of the reflection coefficient, time lag, and the quality factor for each layer are estimated simultaneously. Available information from well-logs can also be incorporated as a priori information. Numerical examples from simulated and field seismic data show high performance of the proposed algorithm for estimating the reflectivity structure and the model of the quality factor. Numerical examples confirm promising applicability of the new method for processing and interpretation of seismic data. The goal of non-stationary deconvolution is to remove the effects of the source wavelet and the attenuation filter from the data to retrieve the earth impulse response with an enhanced temporal resolution.