CC BY
96
УДК 550.8
А. Н. Кремлев КОЛИЧЕСТВЕННЫЕ ХАРАКТЕРИСТИКИ МЕТОДА CSP
Предложены и исследованы количественные характеристики метода Common Scattering Point (CSP), позволяющие оценивать его эффективность при изучении геологических сред по рассеянным сейсмическим волнам, в частности, при поиске трещинных коллекторов углеводородов. Наиболее практически важный из них — коэффициент чувствительности, который показывает, насколько слабые, по сравнению с отражающими границами, рассеивающие объекты могут быть обнаружены методом CSP. Численно исследована зависимость этих характеристик от геометрии систем регистрации сейсмических волн и других параметров.
The quantitative characteristics of the CSP method, which poses to estimate its effectiveness for fracture hydrocarbon reservoirs imaging using scattered seismic waves, are defined and investigated. The most important is the sensitivity coefficient, which shows how small diffractors, compared with the shadowed reflectors, may be to be detecting by the CSP method. These dependencies from acquisition geometry and other parameters are numerically studied.
Ключевые слова: метод CSP, рассеянные волны, дифракторы, трещиннокавернозный коллектор, коэффициент чувствительности.
Key words: CSP method, scattered waves, fracture cavernous reservoir, sensitivity coefficient.
Введение
Общепризнанно, что от 25 до 50 % общемировых запасов углеводородов находятся в резервуарах трещинного и трещинно-кавернозного типа. Поиск и разведка таких резервуаров с помощью традиционных технологий сейсморазведки весьма затруднительны. Причина в том, что эти технологии ориентированы на резервуары с поровыми коллекторами и используют для поиска и разведки отраженные волны. В ловушках трещинно-кавернозного типа распределение нефти по разрезу контролируется каналами миграции флюидов, которые, в свою очередь, связаны с зонами трещиноватости, кавернозности и/или карсто-вания. Эти объекты не отражают, а рассеивают сейсмическую энергию. При этом амплитуда рассеянных волн на несколько порядков слабее отраженных [1], являющихся очень сильными волнами-помехами.
Все существующие сегодня подходы для изучения коллекторов трещинно-кавернозного типа стремятся максимально ослабить (вычесть) отраженные волны. Так, в методе МИРО (миграционное изображение рассеивающих объектов) [2] для этого используются оптико-геометрические представления о распространении энергии вдоль сейсмических лучей и их зеркального отражения. Аналогичные, по сути, идеи эксплуатируются и в фокусирующих преобразованиях [3; 4].
© Кремлев А. Н., 2013
Вестник Балтийского федерального университета им. И. Канта. 2013. Вып. 10. С. 121 — 126.
121
122
1. Метод CSP
Метод CSP (Common Scattering Point) [5 — 7] базируется на решении обратной задачи рассеяния акустических волн по данным МОГТ в приближении Борна. Рассеяние волн описывается волновым уравнением
я2 u
—г = с2[1 + а(г )]Ди + S(i)5(f - Гр), м|(<0 = 0, (1)
а обратная задача заключается в определении неоднородностей среды a(r) по рассеянному полю и0(р, р0, t) = u(r, r0, t)|z=z=0, зарегистрированному на дневной поверхности для различных положений приемника р и источника р0. Здесь r, Г0 є R3; р, р0 є R2; t — время; с = const — скорость акустических волн во вмещающей среде.
В результате решения задачи в явном виде построены два линейных
интегральных оператора Lrefl и Ldifr, действие которых на полное поле МОГТ и0(р, р0, t) приводит к его расщеплению на отраженную и рассеянную компоненты соответственно uiiffr = Liiffr • u0 и Uf = Lrefl • u0.
Затем с помощью престековой волновой миграции M раздельно вычисляются распределения по пространству рассеивающих и отражающих элементов соответственно adlffr (r) = M • udlffr и arefl (r) = M • urefl, причем a(r) = areil (r) + aiiffr (r). Рассеивающие объекты мы далее будем
называть дифракторами, а отражающие площадки — рефлекторами. Единственным параметром, который необходимо задать для их раздельного вычисления (кроме, конечно, геометрии профиля и скорости вмещающей среды), является наибольший угол наклона отражателей amax. Такие же операторы были построены также и для 2D случая, которым мы и ограничимся в настоящей работе.
Многие примеры тестирования метода CSP на синтетических данных и его применения для прогнозирования трещинно-кавернозных коллекторов, в том числе и по данным 3D МОГТ, описаны в [5 — 7].
2. Количественные характеристики метода CSP
Рассмотрим задачу количественной оценки возможностей метода CSP. Эти возможности прежде всего определяются тем, насколько хорошо нам удается вычесть (устранить) из полного поля МОГТ отраженные волны. Качество вычитания отраженных волн будем описывать с помощью коэффициента деструкции kde, равного отношению амплитуды рефлектора, полученного без процедуры вычитания, к максимальному значению остатков от рефлектора, оставшихся после применения этой процедуры.
Вычитание отраженных волн влияет на рассеянное поле. Для количественной характеристики этого влияния определим коэффициент восстановления (реконструкции) ке дифрактора как отношение амплитуды дифрактора, восстановленного после процедуры вычитания, к амплитуде дифрактора, вычисленного без применения этой процедуры. Очевидно, что этот коэффициент всегда меньше единицы и чем он ближе к единице, тем лучше.
Чувствительность метода — это возможность обнаружения с его помощью на временных разрезах слабых объектов. Будем считать, что дифрактор расположен на затеняющем его рефлекторе (наихудший случай) и что амплитуды рассеянной и отраженной волн (в апексе их годографов) равны. Определим коэффициент чувствительности ке как отношение амплитуды дифрактора, вычисленного после вычитания отраженных волн, к максимальной амплитуде остатков от рефлектора. Коэффициент показывает, во сколько раз рассеянное поле может быть слабее отраженного поля, чтобы этот дифрактор мог быть выявлен на временном разрезе с помощью метода СБР. Нужно понимать, что это предельное значение, так как при таком отношении амплитуда восстановленного дифрактора и остатков от рефлектора равны. Очевидно, что для надежной идентификации дифрактора амплитуда остатков должна быть в несколько раз меньше амплитуды целевого объекта.
Заданные коэффициенты зависят от геометрии профиля, точки расположения дифрактора на разрезе, угла наклона и положения рефлектора, скорости миграции, скорости волн во вмещающей среде и наибольшего угла наклона вычитаемых рефлекторов. Для исследования этих зависимостей использовано моделирование синтетических волновых полей, затем обработанных по методу СБР. Приближение Борна, в рамках которого решалась задача (1), это приближение однократного рассеяния/ отражения. Поэтому для моделирования рассеянного поля использовали однократное рассеяние от точечного дифрактора, а для моделирования отражения — зеркальное приближение.
Моделирование 20 волновых полей МОГТ проводилось на апертуре 9,6 км. Источники располагались на отрезке [0; 3,2] км с интервалом А = 25 м. База приема составляла 6,4 км, интервал между приемниками Дг = 25 м. Вынос первого источника в расстановке равнялся -3200 м, а вынос последнего — +3200 м. Скорость волн во вмещающей среде равнялась 2 км/ с, форма сигнала — импульс Риккера с доминантной частотой 20 Гц, шаг между временными отсчетами Дí = 2 мс. Обработка синтетических полей МОГТ проводилась по методу СБР на апертуре, занимающей отрезок [-3,2; 6,4] км для скорости миграции 2 км/с. Принималось, что наибольший угол наклона рефлекторов равен а тах = 5°.
На рисунке 1, а в логарифмическом масштабе дан временной разрез точечного дифрактора, восстановленного по синтетическим данным без применения процедуры вычитания отраженных волн, а на рисунке 1, Ь после применения этой процедуры. Эти разрезы нормированы на амплитуду дифрактора, вычисленного без применения вычитания.
123
124
^ --------------------------------------------------------------
На рисунке 1, с представлены остатки от рефлектора, наклоненного под углом 2°, после CSP-вычитания. Слева — для Ав = 25 м, а справа — для Ав = 100 м. Эти разрезы также были нормированы на амплитуду рефлектора, вычисленного без применения процедуры вычитания.
Из рисунка 1, а видно, что качество фокусировки, используемой в методе CSP, весьма высокое. Практически на всем временном разрезе, кроме непосредственной окрестности дифрактора, уровень шумов ниже -60 dB, то есть более чем в 1000 раз слабее амплитуды самого дифрактора. Процедура «вычитания» приводит к появления дополнительного шума (рис. 1, Ь), имеющего специфическую крестообразную форму с уровнем шума порядка -37 dB, то есть более чем в 70 раз слабее амплитуды дифрактора.
На рисунке 1, с при Ав = 25 м уровень шумов от остатков рефлектора практически везде ниже -47 dB (в 224 раза), а при Ав = 100 м — -41 dB, то есть в 112 раз слабее амплитуды исходного рефлектора. Увеличение уровня остатков рефлектора при росте интервала между приемниками связано с эффектом эляйсинга. Тем не менее в обоих случаях эти остатки существенно меньше амплитуды исходного рефлектора.
Рис. 1. Логарифмические временные разрезы точечного дифрактора и остатков от рефлектора: а, Ь — временной разрез (в логарифмическом масштабе) точечного дифрактора до (а) и после (Ь) применения метода CSP; с, й — остатки от рефлектора для Ав = 25 м и Ав = 100 м соответственно
На рисунке 2 приведен результат реконструкции точечного ди-фрактора после вычитания проходящего через него рефлектора. Угол наклона рефлектора 2°, амплитуда рассеянной волны в 100 раз меньше
амплитуды отраженной волны. Хорошо видно, что этот достаточно слабый дифрактор, который при стандартной обработке совершенно не виден на временном разрезе, очень отчетливо проявляет себя на разрезе СБР-дифракторов.
Рис. 2. Результат реконструкции дифрактора, затененного рефлектором. Амплитуда рассеянной волны в 100 раз меньше амплитуды отраженной волны: а — разрез в линейном; Ь — в логарифмическом масштабах
На рисунке 3 слева приведена зависимость коэффициента реконструкции дифрактора кге от времени £„ его расположения на временном разрезе. Видно, что его относительная величина порядка 80% и незначительно уменьшается при увеличении глубины расположения дифра-ктора. Эта зависимость соответствует шагу по источникам Ав = 25 м. Практически точно такие же графики получается и для Ав = 50, 100, 200 м.
На рисунке 3, Ь изображены графики зависимости коэффициента чувствительности метода СБР для Ав = 25, 50, 100, 200 м. Видно, что для всех этих Ав увеличение времени от 0,5 до 3 с приводит к росту коэффициента чувствительности, а затем — к его снижению. Увеличение шага между источниками монотонно уменьшает коэффициент чувствительности. Если в качестве целевого выбрать время 2 с, то кве = 1200 при Ав = 25 м и кве = 300 при Ав = 200 м. Эти значения для коэффициента чувствительности метода СБР, особенно когда шаг между приемниками и источниками равен 25 м, очень высоки.
Рис. 3. Зависимость коэффициента реконструкции точечного дифрактора кт от ^ (а) и зависимости коэффициента чувствительности кж от ^ для различных значении Ав (Ь)
126
Заключение
Введенные нами коэффициенты позволяют количественно оценить возможности метода CSP для поиска трещинно-кавернозных коллекторов. Эти коэффициенты можно назвать паспортными данными метода CSP. Они дают количественные характеристики применяемого метода обработки и повышают надежность интерпретации и вытекающих из нее интерпретационных выводов. Ограничения чувствительности метода CSP и других характеристик этого метода определяются фундаментальными свойствами волновых процессов и дискретных преобразований, используемых для обработки волновых полей МОГТ. Эти ограничения весьма высоки. Например, при As = Ar = 25 м коэффициент чувствительности порядка 1000. Это показывает очень высокую чувствительность метода CSP, которая обусловливает его эффективность при поиске трещинных и трещинно-кавернозных резервуаров углеводородов. Такие значения чувствительности близки к предельным и вряд ли могут быть превзойдены другими методами.
Список литературы
1. Кузнецов О. Л. и др. Экспериментальные исследования // Сейсмоакустика пористых и трещиноватых геологических сред : в 3 т. М., 2004. Т. 2.
2. Kozlov E. et al. Imagingscattering objects masked by specular reflections // 74th Annual International Meeting, SEG. Expanded Abstracts. 2004. P. 1131 — 1134.
3. Гольдин С. В. и др. Построение сейсмических изображений в рассеянных волнах как средство детализации сейсмического разреза // Геофизика, специальный выпуск. 2004. С. 23—29.
4. Поздняков В. А., Сафонов Д. В., Чеверда В. А. Оптимизация параметров фокусирующих преобразований с использованием численного моделирования // Геология и геофизика. 2000. Т. 41, № 6. С. 930—938.
5. Кремлев А. Н. и др. Прогноз коллекторов трещинно-кавернозного типа по рассеянным сейсмическим волнам // Технологии сейсморазведки. 2008. № 3. С. 36 — 39.
6. Kremlev A. N. et al. Fracture and cavernous reservoirs prospecting by the CSP prestack migration method // Extended Abstracts of 73th EAGE Conference & Exhibition. Vienna, Austria. 2011. URL: http://www.earthdoc. org/ detail. php?pubid=50142.
7. Erokhin G. N. et al. CSP-method rospecting of fracture-cavernous reservoirs in the Bazhen formation of the Salym oilfield / / Extended abstract. 74th EAGE Conference & Exhibition. Copenhagen, Denmark. 2012. URL: http://www.earthdoc. org/detail.php?pubid=59720.
Об авторе
Андрей Николаевич Кремлев — канд. физ.-мат. наук, доц., Балтийский федеральный университет им. И. Канта, Калининград.
E-mail: ankremlev@gmail.com
About the author
Dr Andrey Kremlev — Ass. Prof., I. Kant Baltic Federal University, Kaliningrad.
E-mail: ankremlev@gmail.com
