EA005450B1 - Использование скоростей кусочков породы для прогнозирования в реальном времени порового давления и градиента давления гидравлического разрыва пласта - Google Patents
Использование скоростей кусочков породы для прогнозирования в реальном времени порового давления и градиента давления гидравлического разрыва пласта Download PDFInfo
- Publication number
- EA005450B1 EA005450B1 EA200400568A EA200400568A EA005450B1 EA 005450 B1 EA005450 B1 EA 005450B1 EA 200400568 A EA200400568 A EA 200400568A EA 200400568 A EA200400568 A EA 200400568A EA 005450 B1 EA005450 B1 EA 005450B1
- Authority
- EA
- Eurasian Patent Office
- Prior art keywords
- rock
- pore pressure
- stress
- pieces
- equation
- Prior art date
Links
- 239000011148 porous material Substances 0.000 title claims abstract description 62
- 238000005520 cutting process Methods 0.000 title abstract description 5
- 238000000034 method Methods 0.000 claims abstract description 92
- 239000011435 rock Substances 0.000 claims description 68
- 238000005259 measurement Methods 0.000 claims description 43
- 238000005553 drilling Methods 0.000 claims description 42
- 230000010355 oscillation Effects 0.000 claims description 28
- 230000015572 biosynthetic process Effects 0.000 claims description 25
- 238000010586 diagram Methods 0.000 claims description 8
- 230000005484 gravity Effects 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000001419 dependent effect Effects 0.000 claims description 4
- 239000004576 sand Substances 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000004590 computer program Methods 0.000 claims description 2
- 238000001514 detection method Methods 0.000 claims 2
- 238000010276 construction Methods 0.000 claims 1
- 229920006395 saturated elastomer Polymers 0.000 claims 1
- 238000006467 substitution reaction Methods 0.000 claims 1
- 238000005755 formation reaction Methods 0.000 description 19
- 239000000523 sample Substances 0.000 description 17
- 238000004458 analytical method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 239000012530 fluid Substances 0.000 description 3
- 150000003839 salts Chemical class 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 235000015076 Shorea robusta Nutrition 0.000 description 2
- 244000166071 Shorea robusta Species 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 2
- 239000004927 clay Substances 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 239000012535 impurity Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 229920005479 Lucite® Polymers 0.000 description 1
- 241000364021 Tulsa Species 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 238000011088 calibration curve Methods 0.000 description 1
- 150000001768 cations Chemical class 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 231100001261 hazardous Toxicity 0.000 description 1
- 238000007654 immersion Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000011545 laboratory measurement Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000004058 oil shale Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000004926 polymethyl methacrylate Substances 0.000 description 1
- 239000010454 slate Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000008961 swelling Effects 0.000 description 1
- 238000004448 titration Methods 0.000 description 1
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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
- E21B49/005—Testing the nature of borehole walls or the formation by using drilling mud or cutting data
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Environmental & Geological Engineering (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Geophysics (AREA)
- Fluid Mechanics (AREA)
- Acoustics & Sound (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Geophysics And Detection Of Objects (AREA)
- Perforating, Stamping-Out Or Severing By Means Other Than Cutting (AREA)
Abstract
Способ предназначен для точного прогнозирования в реальном времени порового давления и градиента давления гидравлического разрыва пласта на месте нахождения буровой установки путем определения скорости волн по кусочкам породы посредством портативного прибора для измерений методом незатухающих гармонических колебаний, с помощью которого проводят измерения кусочков породы на высокой резонансной частоте, и затем полученную скорость используют в сочетании с новым способом настоящего изобретения для получения точных прогнозов значений порового давления и градиента давления гидравлического разрыва пласта. При использовании способа обеспечивается прогнозирование в реальном времени порового давления на месте нахождения буровой установки с небольшим пределом ошибки, который нельзя получить, используя сейсмическую разведку, вертикальное сейсмическое профилирование или контроль скорости по взрыву.
Description
Область техники
Настоящее изобретение относится к анализу и интерпретации скоростей акустических волн в образцах породы, полученных при проходке. Более конкретно, изобретение относится к новому способу очень точного прогнозирования в реальном времени порового давления, градиента давления гидравлического разрыва и прочности пласта на месте нахождения буровой установки путем определения скорости кусочка породы посредством портативного прибора для измерений методом незатухающих гармонических колебаний (НГК) на высокой резонансной частоте и используя измерение скорости в сочетании со способом настоящего изобретения для получения точных прогнозов значений порового давления и градиента давления гидравлического разрыва пласта. При использовании нового способа обеспечивается прогнозирование в реальном времени порового давления на месте нахождения буровой установки с небольшим полем ошибки, который нельзя получить, используя сейсмическую разведку, вертикальное сейсмическое профилирование или контроль скорости по взрыву.
Предшествующий уровень техники
Решающее значение во время бурения скважины имеет знание точных значений порового давления, градиента давления гидравлического разрыва и прочности пласта. Кроме того, поровое давление и градиент давления гидравлического разрыва пласта являются управляющими входными параметрами при моделировании устойчивости ствола скважины, проектировании и расчете скважины и оптимизации траектории ствола скважины. Хотя отсутствуют промышленные приборы для измерения порового давления перед буровой коронкой, но разработаны методики для вычисления порового давления на интервалах проведения каротажа по данным удельного сопротивления и/или акустического каротажа.
Известно несколько способов для получения информации во время бурения, касающейся других характеристик. Например, известен способ определения пористости пласта по характеристикам бурения. В патенте США № 4064749 раскрыт способ определения пористости пласта по характеристикам бурения, в котором при вращении буровой коронки, прикрепленной к нижнему концу бурильной колонны, управляют направленной вниз силой, действующей на буровую коронку. Способ включает в себя этапы измерения частоты вращения буровой коронки, измерения глубины погружения буровой коронки в стволе скважины, измерения массы на указанной буровой коронке, определяющей затупление зубьев указанной буровой коронки, измерения вращающего момента, приложенного к бурильной колонне, определения эмпирически эталонного вращающего момента и определения пористости путем объединения указанных результатов измерений и определений и использования уравнения.
В патенте США № 4949575 раскрыт способ анализа пласта, в котором результаты измерений характеристик бурения используют в качестве входных данных, зависящих от пористости. Способ включает в себя получение сигнала бурения, характеризующего сопротивление бурению пласта, пробуриваемого буровой коронкой, получение большого количества дополнительных сигналов, отражающих свойства пласта, и выполнение на основании указанного сигнала бурения и указанных дополнительных сигналов волюметрического анализа приповерхностного пласта.
В патенте США № 4876512 раскрыт способ определения на буровой площадке содержания разбухающей глины в сланцах и в сланцеватых песчаниках путем измерений площади проводящей поверхности. Образцы промывают жидкостью, имеющей существенно меньшую активность по сравнению с активностью воды, и которая может содержать растворимый катион, и на заранее выбранной частоте (1 МГц) проводят измерения диэлектрической проницаемости образцов для последующего сравнения с калибровочными кривыми, тем самым получая данные по содержанию разбухающей глины в пласте.
В патенте США № 5282384 (Холбрук) раскрыт усовершенствованный способ, основанный на принципах механического зондирования, с использованием данных каротажа скважины для вычисления давления жидкости, содержащейся в осадочной породе, которая естественным образом уплотнилась под действием силы тяжести. Часть способа, относящаяся к определению эффективного напряжения, включает в себя внутренние и внешние измерения деформации зернистой матрицы породы. Поэтому одна и та же калибровка эффективного напряжения может быть равным образом применена к данным о толщине породы, полученным путем измерения, и к петрофизическим данным о пористости породы. Степенная зависимость между эффективным напряжением и деформацией для любой осадочной породы может быть определена из взвешенного среднего степенных функций для минералов, которые образуют осадочную породу. В настоящем изобретении для прогнозирования порового давления и градиента давления гидравлического разрыва пласта способ Холбрука модифицирован как в части методики, так и в части типа данных.
Проблема, часто встречающаяся при бурении скважин во многих частях мира, заключается в небольших полях ошибок, что требует большой точности при прогнозировании как порового давления, так и градиента давления гидравлического разрыва пласта, чтобы предотвратить нестабильность сланцев, создающую опасность прекращения циркуляции и/или выброса газа/фонтанирования. Например, при большой глубине Мексиканского залива допуск по бурению, выраженный через поровое давление и через градиент давления гидравлического разрыва пласта, может быть меньше ±0,5 фунта/галлон. Поэтому точность скоростей, получаемых на основании данных сейсмических исследований, каротажа во время бурения и/или кабельного каротажа, является очень важной для бурения в дополнение к другим известным применениям скорости для петрофизических или пластовых исследований.
Кусочки породы, образующиеся при бурении, представляют собой источник информации в квазиреальном времени, который можно получать на месте нахождения буровой установки. В прошлом использование кусочков породы ограничивалось из-за трудностей, возникавших при выполнении измерений на очень малых образцах.
В отчете №в е! а1., Кщ-вбс апб 1аЬога£огу иве о£ ΟΨΤ асоивйс νοίοοίίν тсазигстсШз оп сиШпдв. 8ос1е£у о£ Рс£го1сит Епдшссгв Рарсг № 36854, 1996, показана возможность использования методики незатухающих гармонических колебаний для измерения скорости акустических волн в кусочках породы при использовании портативного оборудования, которое пригодно для контроля небольших образцов, и тем самым предложен новый источник данных, который может быть использован в квазиреальном времени на месте нахождения буровой установки.
Для калибровки скоростей сейсмических волн, данных каротажа во время бурения и/или результатов кабельного акустического каротажа и, что еще более важно, для получения в реальном времени точных результатов измерений скорости волн в пласте кусочки породы обеспечивают потенциально бесценный источник информации.
В настоящее время необходим способ, использование которого позволяет точно прогнозировать поровое давление и градиент давления гидравлического разрыва пласта путем выполнения в реальном времени измерений на месте нахождения буровой установки. В случае получения таких данных они также будут полезными для распознавания опасности мелководных зон, оптимизации массы бурового раствора, обнаружения опасных пустот, обнаружения зон аномального давления, определения прочности пласта для оптимизации траектории скважины и в общем случае для осуществления наиболее безопасного, экономически эффективного бурения.
Краткое изложение существа изобретения
Задачей настоящего изобретения является создание способа прогнозирования в реальном времени порового давления и градиента давления гидравлического разрыва пласта на месте нахождения буровой установки путем использования метода незатухающих гармонических колебаний на резонансных часто тах для получения в реальном времени значения скорости кусочка породы и путем использования данных в сочетании со способом, раскрытым в настоящем изобретении для точного прогнозирования порового давления и градиента давления гидравлического разрыва пласта. Кроме того, значения скоростей также используют в реальном времени для определения механических характеристик и прочности пла ста.
Способ, согласно изобретению, заключается в том, что
1) получают результаты измерений скорости Ур кусочка породы, для чего
a) подготавливают образец с двумя противоположными параллельными поверхностями, разделенными расстоянием Ь,
b) помещают образец между преобразователями Р- и 8-волн,
c) непрерывно пропускают синусоидальный сигнал в пределах надлежащего частотного диапазона, б) извлекают и регистрируют огибающую принятого амплитудно-модулированного сигнала,
с) измеряют А£, ширину кривой на уровне половинной мощности для каждого максимума амплитуды и осуществляют подстановку значения в уравнение Ур=2Ь-А£,
2) используют полученные результаты измерения скорости для определения пористости по скорости с помощью акустического метода пластового коэффициента при использовании уравнения
где β=0,625 и Ур в матрице =18000 футов/с для песка,
3) прогнозируют поровое давление на основании зарегистрированной скорости, полученной методом незатухающих гармонических колебаний, для чего
а) вычисляют суммарное вертикальное напряжение (напряжение в покрывающем слое) путем ин тегрирования плотности
8ν=ίρ(ζ)·§·6ζ, где ζ - глубина; |}(ζ) - объемная плотность породы пласта как функция глубины; б - плотность; д постоянная силы тяжести,
b) получают эффективное вертикальное напряжение при использовании корреляции Холбрука (Холбрук, 1996) σν ^тах(1-ф) , где атах и α представляют собой постоянные параметры, зависящие от литологии, определенные Холбруком,
c) определяют поровое давление при использовании закона ТсгеадЫ
-2005450
Ρρ=(8ν-σν)/γ, где Рр - поровое давление; 8ν - суммарное напряжение в покрывающем пласте, вычисленное на этапе 3а; σν - эффективное вертикальное напряжение, вычисленное на этапе 3Ь; γ - коэффициент Байота, равный 1-(Кр/Ку), где Кр - объемная сжимаемость, К.. - сжимаемость зерен,
4) прогнозируют градиент давления гидравлического разрыва пласта с помощью модифицированного способа Холбрука при использовании скорости кусочка породы, полученной методом незатухающих гармонических колебаний, из уравнения
8|, ωΐη=δ[Ρρ+σν(1-φ)], где 8|, - минимальное горизонтальное напряжение, которое представляет собой нижний предел градиента давления гидравлического разрыва пласта; φ - пористость, полученная на этапе 2; Рр - поровое давление, полученное на этапе 3(с); δν - эффективное напряжение, вычисленное на этапе 3(Ь); δ показатель локальной прочности (для района разведочных работ используют 0,85 по умолчанию),
5) или прогнозируют градиент давления гидравлического разрыва пласта с помощью модифицированного способа Сат С1ау при использовании диэлектрических измерений кусочков породы (ДИКП), для чего
a) измеряют диэлектрическую проницаемость кусочка породы,
b) вычисляют площадь поверхности (ПП) на основании диэлектрических измерений кусочков по роды при использовании уравнения
ПП=0,0274-(ДИКП)1,9467,
с) определяют безразмерный коэффициент М, характеризующий механическую прочность пласта, при использовании площади поверхности, полученной на основании диэлектрических измерений кусочков породы на этапе 5(Ь), из уравнения
31П1ОЛ
М = ---------- , + 3ϊη10 где А=1,54-1,36Е-3-(площадь поверхности),
б) получают значения минимального горизонтального напряжения при использовании уравнения
35у - (5У - ЗМ)Рр
2М + 3
Краткое описание чертежей
В дальнейшем изобретение поясняется описанием предпочтительного варианта выполнения со ссылками на сопровождающие чертежи, на которых фиг. 1 изображает схему экспериментальной установки, согласно изобретению;
фиг. 2 - резонансный спектр частот для кусочка сланца толщиной 2,5 мм, согласно изобретению;
фиг. 3 - диаграмму зависимости скоростей, полученных методом незатухающих гармонических колебаний, от частоты для контрольных материалов - люцита и соли, согласно изобретению;
фиг. 4 - диаграмму зависимости скоростей, полученных в реальном времени методом незатухающих гармонических колебаний на скважине в Мексиканском заливе, согласно изобретению;
фиг. 5 - спрогнозированные в реальном времени значения порового давления и градиента давления гидравлического разрыва пласта при использовании данных полученных кусочков породы методом незатухающих гармонических колебаний и диэлектрических измерений, в зависимости от глубины на выбранных интервалах скважины в Мексиканском заливе, согласно изобретению;
фиг. 6 - диаграмму зависимости значений скоростей, полученных методом незатухающих гармонических колебаний и кабельным акустическим каротажем, от глубины на выбранных интервалах скважины в Мексиканском заливе, согласно изобретению;
фиг. 7 - диаграмму зависимости спрогнозированных значений порового давления и градиента давления гидравлического разрыва пласта от глубины для скважины в Мексиканском заливе после бурения интервала и сбора данных кабельного акустического каротажа; прогнозы по методике незатухающих гармонических колебаний и на основе диэлектрических измерений кусочков породы сделаны в реальном времени во время бурения на том же самом интервале и приведены для сравнения, согласно изобретению;
фиг. 8 - диаграмму зависимости спрогнозированных значений порового давления и градиента давления гидравлического разрыва пласта от глубины для другой скважины в Мексиканском заливе, полученных при использовании данных метода незатухающих гармонических колебаний и диэлектрических измерений кусочков породы.
Подробное описание предпочтительных вариантов воплощения изобретения
Использование заявленного способа согласно изобретению обеспечивает ряд преимуществ по сравнению с любым способом, известным в настоящее время. Прежде всего, обеспечивается возможность точного прогнозирования в реальном времени порового давления и градиента гидравлического разрыва пласта при использовании результатов измерений на образцах реальной породы на месте нахождения буровой установки. Способ также обеспечивает возможность прогнозирования прочности пласта в буро
-3005450 вой скважине без проведения акустического каротажа, что предоставляет операторам возможность получать все входные данные, необходимые для анализа устойчивости ствола скважины, и проводить анализ устойчивости ствола скважины, а для успешного бурения скважины давать рекомендации относительно оптимальных давлений в скважине. Потенциально эти особенности обеспечивают значительное снижение затрат на бурение, связанных с точной оптимизацией траектории скважины.
Кроме того, результаты измерений могут быть использованы для сравнения с получаемыми позднее данными каротажа при бурении для анализа соответствия, при этом выявление различий может способствовать определению в реальном времени опасных зон, смягчению возникающих проблем, предотвращению неустойчивости ствола скважины и минимизации времени простоя при эксплуатации.
Прибор для измерений методом незатухающих гармонических колебаний, используемый в изобретении для измерения характеристик небольших образцов кусочков породы на месте нахождения буровой установки, представляет собой портативное устройство, содержащее держатель образца с двумя пьезоэлектрическими ультразвуковыми преобразователями, расположенными навстречу друг другу. Один преобразователь используют для возбуждения акустических волн, а второй преобразователь используют для их обнаружения. Один из преобразователей закреплен на микрометрическом винте микрометра, что позволяет одновременно измерять длину образца. Образец располагают между двумя преобразователями. При подаче частоты возбуждения в пределах диапазона, соответствующего нескольким резонансам стоячих волн, в составном резонаторе создаются резонансы ультразвуковых стоячих волн. Огибающую принятого сигнала усиливают с помощью широкополосного приемника с регулируемым коэффициентом усиления, который подключают к переносному компактному персональному компьютеру. При поступлении нескольких входных параметров с помощью вспомогательной программы вычисляются скорости. Пример результатов измерений, выполненных для образца сланца из бурового шлама толщиной 2,5 мм (0,5 дюйма), показан на фиг. 2. Портативное устройство для измерений методом незатухающих гармонических колебаний, использованное в настоящем изобретении (фиг. 1), можно получить от ТЕМСО Сотрапу, Талса, штат Оклахома.
Ниже приведены определения терминов, используемых в описании изобретения:
α - затухание звуковых колебаний, Нп/см; НГК - незатухающие гармонические колебания; МНГК метод незатухающих гармонических колебаний; £ - частота, Гц, с-1; Ό - диаметр, м; Ь, ζ - толщина образца, м; λ - длина волны, м; ИГ СБР - измеренная глубина стола бурового ротора, м; ФВГ - фактическая вертикальная глубина; п - число половин длин волн, целое число; Ур - скорость звука, м/с, мкс/фут; Р первичная волна; 8 -поперечная/вторичная волна; σν - вертикальное эффективное напряжение; φ - пористость; птах и α - постоянные параметры, зависящие от литологии; 8ν - суммарное напряжение в покрывающем пласте; ρ(ζ) - объемная плотность породы как функция глубины; ζ - глубина; д - постоянная силы тяжести; б - плотность; Рр - прогнозируемое поровое давление; уРр - поровое давление, где γ - коэффициент Байота=1-(Кр/Ку), Κβ - объемная сжимаемость, Κγ -сжимаемость зерен.
Нижние индексы:
п - число половин длин волн; р - первичный.
Устройство для измерения методом незатухающих гармонических колебаний (МНГК) основано на создании резонансов стоячих волн в составном резонаторе при подаче частоты £ возбуждения в пределах диапазона, соответствующего нескольким резонансам стоячих волн. Посредством второго преобразователя принимается амплитудно-модулированный сигнал, включающий резонансы, показанные на фиг. 2. Каждый из п резонансов стоячих волн характеризуется резонансной частотой £п. На практике при наличии плоских волн предполагается, что отношение λ/Ό является достаточно малым, где λ - длина акустической волны, Ό - диаметр. Резонанс будет наблюдаться, по существу, на каждом интервале λ=ν/£, чему соответствует Τ=η·λ/2, где Ь - толщина образца или расстояние между поверхностями двух преобразователей, п - целое число. Фазовая скорость акустических волн для Р- или 8-волн может быть выражена как Ур = 2ΚΔΓ. где Δ£ = £п+1-£п представляет собой разность частот двух последовательных резонансов. Например, для образца сланца с Ур=6550 футов/с и Ь=0,06 дюйма Δ£ будет составлять »0,67 МГц, а практическая частота будет находиться в пределах 1-10 МГц.
При подготовке к измерениям оператор получает образец кусочка породы с двумя противоположными параллельными поверхностями, отстоящими друг от друга на расстоянии Ь (расстояние Ь соответствует длине образца кусочка породы), прикрепляет к каждой из двух поверхностей преобразователи Рили 8-волн, в соответствующем частотном диапазоне непрерывно изменяет частоту синусоидального сигнала, выделяет и регистрирует огибающую принятого амплитудно-модулированного сигнала, измеряет Δ£ и подставляет значения в уравнение Ур=2^-Δ£.
Подготовка образца включает в себя формирование двух гладких и параллельных поверхностей и может быть выполнена шкуркой. Кусочки породы следует отбирать внимательно, чтобы уменьшить число ошибок, обусловленных неопределенностью глубины, неоднородностями и механическими повреждениями кусочков буровой породы и/или обрушениями. Можно подготавливать и измерять образцы с мелкими зернами толщиной в долю дюйма меньше, чем через 10 мин после помещения бурового шлама на вибрационное сито.
-4005450
Для проведения измерений образец зажимают между двумя преобразователями. Микрометр прикрепляют к металлической раме и используют для позиционирования преобразователя, тем самым обеспечивают возможность измерения длины образца под нагрузкой. Чтобы гарантировать достаточное качество сигнала без влияния на характеристики образца, на обе поверхности подготовленного образца наносят соответствующим образом выбранный гель для улучшения акустической связи.
Требуемый диапазон подаваемой частоты обеспечивают посредством генератора частот. Огибающую принятого сигнала усиливают с помощью широкополосного приемника с регулируемым коэффициентом усиления, который подключают к переносному компактному персональному компьютеру. Скорости вычисляют с помощью вспомогательной программы компьютера при поступлении нескольких входных параметров.
Использование способа согласно настоящему изобретению обеспечивает получение точных данных в реальном времени на образцах меньше дюйма. Однако было замечено, что, хотя значительные изменения точности отсутствуют, для образцов от 2,5 до 3 мм имеется возможность получать наиболее точные данные. Для образцов меньше 2,5 мм или образцов, которые толще 4-5 мм, иногда получают неверные данные.
После определения значений скоростей для подготовленных образцов, имеющих предпочтительную длину, каждую точку данных скорости преобразуют в прогнозируемое поровое давление, используя модифицированный способ Холбрука.
Модифицированный способ Холбрука представляет собой один из наиболее надежных способов для прогнозирования порового давления. Способ Холбрука выгодно отличается независимостью от выбора ориентации сланца, что требуется для большинства других известных способов. При использовании этого способа непосредственно учитываются литологические изменения при прогнозировании порового давления. Модифицированный способ Холбрука основан на трех основных уравнениях. Первое уравнение представляет собой определение вертикального эффективного напряжения, то есть закон эффективного напряжения Тег/ац1ч σν=5ν-λΡρ <=> Ρρ=(3ν-σν)/γ(1) где 8ν - суммарное напряжение в покрывающем пласте; Рр -поровое давление; γ - коэффициент Байота, равный 1-(Κβ/Κγ), Κβ - объемная сжимаемость, Κγ - сжимаемость зерен. 8ν обычно получают путем интегрирования плотности по данным каротажа, образуя второе уравнение
3ν=Ιρ (ζ) · д · άζ(2)
Третье уравнение представляет собой эмпирическое уравнение, полученное на основании расширенного множества данных каротажа скважины и измерений пластового давления. Этим уравнением связь эффективного напряжения с пористостью пласта представлена в виде σν=σιη3χ(1-φ) ,(3) где Зтах и α представляют собой постоянные параметры, зависящие от литологии (Холбрук, 1996).
Для уравнения (3) в качестве исходных данных требуется пористость. Пористость может быть измерена непосредственно в лаборатории при использовании образцов керна или может быть получена из разнообразных каротажных данных (сопротивления, проводимости или акустических), сейсмических данных или из данных лабораторных измерений кернов, или из скорости кусочков породы.
Предпочтительно определять пористость по скорости, чтобы минимизировать влияние изменений температуры и солености на величины каротажных данных. Кроме того, небольшие изменения пористости, связанные с пластами, находящимися под аномальными давлениями, в значительной степени отражаются на акустических характеристиках. Применив метод незатухающих гармонических колебаний, скорости, полученные на основании измерений, можно использовать для определения пористости и объемной сжимаемости.
Пористость может быть определена по скорости с помощью акустического пластового коэффициента согласно уравнению ( ν р измеренная
V у р в матрице где β=0,625 и νρ в матрице=18000 футов/с для песка. После определения пористости уравнение (3) используют для получения вертикального эффективного напряжения σν. После нахождения вертикального эффективного напряжения уравнение (1) используют для прогнозирования порового давления.
Во втором варианте осуществления изобретения можно спрогнозировать градиент давления гидравлического разрыва пласта путем выполнения двух различных вычислений. В первом случае градиент давления гидравлического разрыва пласта прогнозируют с помощью модифицированного способа Холбрука, используя скорость кусочков породы, определенную методом незатухающих гармонических колебаний
8Ь ,^δίΡρ+σνΟ-φ)] где 8Ь т;п - минимальное горизонтальное напряжение, которое характеризует нижний предел градиента гидравлического разрыва пласта; φ - пористость; Рр - поровое давление; σν -эффективное напряже
(4) ние, вычисление которого рассматривалось выше; δ - коэффициент локального напряжения (0,85 по умолчанию для района разведочных работ).
Альтернативно градиент давления гидравлического разрыва пласта может быть спрогнозирован с помощью модифицированного метода Сат С1ау при использовании результатов диэлектрических измерений кусочков породы (ДИКП) следующим образом:
a) измеряют диэлектрическую проницаемость кусочка породы,
b) на основании результата диэлектрических измерений вычисляют площадь поверхности (ПП) ПП=0,0274- (ДИКП)1,9467,
c) определяют безразмерный коэффициент М, характеризующий механическую прочность пласта, при использовании площади поверхности, полученной на основании диэлектрических измерений кусочка породы
63ίη10Α м = ---------Г / + δίηΙΟ где А=1,54-1,36Е-3-(площадь поверхности),
б) получают значения минимального горизонтального напряжения согласно уравнению = 35ν - (5ν - ЗМ)Рр Ига1П 2М + 3 *
После вычисления прогнозируемых значений порового давления и градиента давления гидравлического разрыва пласта значения могут быть нанесены на график в зависимости от глубины. Дополнительно это рассмотрено в примерах 2 и 3 и отображено на фиг. 5, 7 и 8.
В настоящем изобретении метод незатухающих гармонических колебаний для измерения скоростей сочетается с новыми способами, разработанными заявителем, что обеспечивает возможность точного прогнозирования порового давления и градиента давления гидравлического разрыва пласта по скорости продольных волн. Способ согласно настоящему изобретению может быть использован для проверки прогноза порового давления на основании скоростей сейсмических волн в районе разведочных работ, а также для калибровки сейсмических данных до того, как станут доступны данные акустического каротажа. Настоящее изобретение будет также весьма полезным для обнаружения слабых зон, которые предрасположены к быстрому увеличению и утрате циркуляции, и для решения проблем устойчивости ствола скважины, например заклинивания труб и налипания породы на долото.
Способ согласно настоящему изобретению является удобным, относительно недорогим и простым в реализации. Используя оборудование, которое может быть доставлено на буровую установку в портфеле, и используя образцы размером меньше дюйма, можно в течение нескольких минут определить скорости и преобразовать их в предсказанные значения порового давления и градиента давления гидравлического разрыва пласта. Опыт работы с прибором для измерений методом незатухающих гармонических колебаний на разведочно-эксплуатационной скважине во время работы буровой установки показал, что существует значительное отклонение получаемых результатов от интервала скоростей сейсмических волн, однако было обнаружено, что они намного лучше согласуются с данными кабельного акустического каротажа, когда указанные данные кабельного акустического каротажа собирались на более поздней стадии.
При использовании изобретения на месторождении было замечено, что при интервалах глубин около ±30 футов в большинстве случаев точность была достаточной, но для повышения точности кусочки породы можно собирать с более частыми интервалами, например в пределах глубин ±10 футов. Кроме того, измерения выполняются при напряжении по одной оси, а не по трем осям, однако было обнаружено, особенно для сланцев, небольшое различие между результатами, полученными при напряжении по одной оси и по трем осям, не приводит к значительному отклонению от реального давления.
Нижеследующие примеры служат иллюстрацией изобретения, раскрытого в настоящей заявке.
Эксперимент
Пример 1.
Контрольные материалы исследовались в лаборатории и на месте нахождения буровой установки для проверки точности методики измерения скорости (методом незатухающих гармонических колебаний), которая является источником исходных данных для способа, рассмотренного в настоящем изобретении. Результаты контрольных исследований представлены на фиг. 3. Образцы соли, использованные при измерениях, имели примесь. Поэтому были получены два различных значения Αί, соответствующие чистой поваренной соли и примеси. На базе этих двух различных значений Δί зависимости скорости от частоты получены методом незатухающих гармонических колебаний для двух различных скоростей.
Пример 2.
Прибор проверяли в полевых условиях на месте нахождения буровой установки в Мексиканском заливе, а измерения в реальном времени методом незатухающих гармонических колебаний выполняли для отдельных интервалов глубин буровой скважины. Из-за больших размеров скважины акустический каротаж во время бурения не проводили. Поэтому скорости, определяемые в реальном времени методом незатухающих гармонических колебаний, были единственной основой во время бурения для калибровки порового давления, предсказанного до бурения на основании сейсмических данных. Полученные зави
-6005450 симости скорости от глубины для отдельных интервалов представлены на фиг. 4. Соответствующие диаграммы порового давления и градиента давления гидравлического разрыва пласта, полученные на основании измерений скоростей методом незатухающих гармонических колебаний, представлены на фиг. 5, а полученные на основании измерений диэлектрической проницаемости при использовании кусочков породы, показаны на фиг. 6. По окончании бурения интервала собирали данные кабельного акустического каротажа. Сравнение скоростей, полученных с помощью кабельного зонда, и скоростей, измеренных методом незатухающих гармонических колебаний, представлено на фиг. 6.
Сравнение итоговых значений порового давления и градиента давления гидравлического разрыва пласта, полученных на основании сейсмических данных, данных кабельного акустического каротажа, скоростей, измеренных методом незатухающих гармонических колебаний, и градиента давления гидравлического разрыва пласта, полученного на основании диэлектрических измерений кусочков породы, приведено на фиг. 7. На фиг. 7 видно, что спрогнозированные в реальном времени значения порового давления и градиента гидравлического разрыва пласта на основании измерений методом незатухающих гармонических колебаний и диэлектрических измерений кусочков породы являются точными и обеспечивают чрезвычайно полезную информацию для подтверждения глубины нахождения уступа в скважине для опоры башмака обсадной колонны, массы бурового раствора, для анализа устойчивости ствола скважины и другой буровой активности с целью безопасного бурения.
Пример 3.
Кусочки породы из другой скважины в Мексиканском заливе использовали для прогнозирования порового давления и градиента давления гидравлического разрыва пласта при использовании скорости, измеренной методом незатухающих гармонических колебаний, и данных диэлектрических измерений кусочков породы. Спрогнозированные значения порового давления и градиента давления гидравлического разрыва пласта, полученные на основании скоростей и данных диэлектрических измерений кусочков породы, представлены на фиг. 8, как и полученные из каротажа профили порового давления и градиента давления гидравлического разрыва пласта, а также реальные данные о давлении жидкости, измеренные при периодических испытаниях пласта, и данные испытаний на просачивание, уже проведенных при различных глубинах башмака обсадной колонны на соседних скважинах месторождения. Видно хорошее соответствие спрогнозированных значений порового давления, полученных на основании измерений методом незатухающих гармонических колебаний, данным, измеренным при периодических испытаниях пласта.
Claims (15)
- ФОРМУЛА ИЗОБРЕТЕНИЯ1. Способ прогнозирования порового давления во время бурения по кусочкам породы, заключающийся в том, чтоa) получают кусочки породы при бурении,b) вычисляют скорость по кусочкам породы, используя метод незатухающих гармонических колебаний,c) вычисляют пористость φ пласта,б) используют пористость φ для вычисления вертикального эффективного напряжения σν согласно уравнению σν=σοαχ(1-φ) , где σ^ и α - постоянные параметры, зависящие от литологии,е) вычисляют суммарное напряжение 8ν покрывающего пласта путем интегрирования плотности по данным каротажа8ν=ίρ(ζ)·§·6ζ, где ζ - глубина; ρ(ζ) - объемная плотность породы пласта как функция глубины; б - плотность; д - постоянная силы тяжести,ί) прогнозируют поровое давление согласно уравнениюРр=(8у-а„)/у, где 8ν - суммарное напряжение в покрывающем пласте; Рр - поровое давление; γ - коэффициент Байота, равный 1-(Κβ/Κγ), Κβ - объемная сжимаемость, Κγ - сжимаемость зерен.
- 2. Способ по п.1, отличающийся тем, что осуществляют построение диаграммы зависимости спрогнозированного порового давления от глубины.
- 3. Способ по п.1, отличающийся тем, что используют метод незатухающих гармонических колебаний, который реализуют при использовании прибора, содержащего портативное устройство, состоящее из держателя образца с двумя пьезоэлектрическими ультразвуковыми преобразователями, обращенными поверхностями друг к другу, причем один преобразователь используют для возбуждения акустических волн, а второй преобразователь используют для обнаружения.
- 4. Способ по п.3, отличающийся тем, что дополнительноа) подготавливают кусочки породы для помещения в держатель образца с двумя противоположными параллельными поверхностями, разделенными расстоянием Ь, где Ь - длина образца кусочка породы,-7005450b) размещают образец между преобразователями Р- и 8-волн,c) непрерывно подают синусоидальный сигнал в пределах заданного частотного диапазона,ά) извлекают и регистрируют огибающую принятого амплитудно-модулированного сигнала,е) измеряют ширину Δ£ полосы и осуществляют подстановку этого значения в уравнение у=2Б-А£ для получения результата измерения скорости для каждого максимума амплитуды, где А£ = £η+ι-£η представляет собой разность частот между двумя последовательными резонансными значениями.
- 5. Способ по п.2, отличающийся тем, что дополнительно осуществляют сравнение указанного спрогнозированного порового давления с прогнозами на основании скоростей сейсмических волн для калибровки данных до получения данных акустического каротажа.
- 6. Способ определения прогнозируемого порового давления при бурении по субдюймовым кусочкам породы, заключающийся в том, чтоa) получают результат измерения скорости кусочка породы путем подготовки образца с двумя противоположными параллельными поверхностями, отделенными расстоянием Б (Б - длина образца),b) помещают образец между преобразователями Р- и 8-волн,c) непрерывно подают синусоидальный сигнал в пределах заданного частотного диапазона,ά) определяют и регистрируют огибающую принятого амплитудно-модулированного сигнала, ие) измеряют Δ£ для каждого максимума амплитуды, где Δ£ = £η+ι-£η представляет собой разность частот между двумя последовательными резонансными значениями, и путем подстановки Δ£ в уравнениеУр=2БА£ вычисляют скорость Ур для каждого максимума амплитуды, где Б - расстояние между парал лельными поверхностями или длина кусочка породы,£) преобразуют значения Ур скорости в прогнозируемое поровое давление, для чего вычисляют пористость пласта, используют пористость и вычисляют вертикальное эффективное напряжение σν согласно уравне нию σν σтаx(1-Φ)α, где стах и α представляют собой постоянные параметры, зависящие от литологии,д) вычисляют суммарное напряжение 8ν покрывающего пласта путем интегрирования плотности по данным каротажа §ν=ίρ(ζ)·§·άζ, где ζ - глубина, ρ(ζ) - объемная плотность породы пласта как функция глубины, ά - плотность, д - посто янная силы тяжести,й) прогнозируют поровое давление согласно уравнениюΡρ=(8ν-σν)/γ, где 8ν - суммарное напряжение в покрывающем пласте, Рр - поровое давление, γ - коэффициент Байота, равный 1-(Κβ/Κγ), Κβ - объемная сжимаемость, Κγ - сжимаемость зерен,1) осуществляют построение кривой порового давления в зависимости от глубины.
- 7. Способ прогнозирования градиента давления гидравлического разрыва пласта по кусочкам породы при бурении, заключающийся в том, чтоa) получают результат измерения скорости Ур в кусочке породы, для чегоb) изготавливают образец с двумя противоположными параллельными поверхностями, разделенными расстоянием Б (Б - длина образца),c) помещают образец между преобразователями Р- и 8-волн,ά) непрерывно подают синусоидальный сигнал в пределах заданного частотного диапазона,е) определяют и регистрируют огибающую принятого амплитудно-модулированного сигнала,£) измеряют Δ£ - разность частот между двумя последовательными резонансными значениями для каждого максимума амплитуды и осуществляют подстановку Δ£ в уравнение νρ=2Ε·Δ£,д) используют полученный результат измерения скорости для определения пористости по скорости с помощью акустического метода пластового коэффициента при использовании уравнения где β=0,625 и Ур в ма1рице=18000 футов/с для песка,й) прогнозируют поровое давление на основании скорости, полученной методом незатухающих гармонических колебаний, для чего1) вычисляют суммарное вертикальное напряжение (напряжение в покрывающем пласте) путем ин тегрирования плотности
- 8ν=ίρ(ζ)·§·άζ, где ζ - глубина; ρ(ζ) - объемная плотность породы пласта как функция глубины; ά - плотность; д - посто янная силы тяжести,Ί) получают эффективное вертикальное напряжение при использовании корреляции Холбрука-8005450 σν σιηίΐχ(1-φ) , где атах и α представляют собой постоянные параметры, зависящие от литологии, определенные Холбруком,k) определяют поровое давление согласно закону Ί οΐ'ζαμίιίΡρ=(8ν-σν)/γ, где Рр - поровое давление; 8ν - суммарное напряжение в покрывающем пласте, вычисленное на этапе ί; σν - эффективное вертикальное напряжение, вычисленное на этапе _|; γ - коэффициент Байота, равный 1-(Κβ/Κγ), Κβ объемная сжимаемость, Κγ - сжимаемость зерен,l) прогнозируют градиент давления гидравлического разрыва пласта с помощью модифицированного способа Холбрука при использовании скорости кусочка породы, полученной методом незатухающих гармонических колебаний8|, тт=б[Рр+^<1-ф)], где 8Ь т;п - минимальное горизонтальное напряжение, которое представляет собой нижний предел градиента давления гидравлического разрыва пласта; ф - пористость, полученная на этапе д; Рр - поровое давление, полученное на этапе к; δν - эффективное напряжение, вычисленное на этапе _|; δ - показатель локальной прочности (для района разведочных работ используют 0,85 по умолчанию).8. Способ по п.7, отличающийся тем, что дополнительно осуществляют построение диаграммы зависимости спрогнозированного градиента давления гидравлического разрыва пласта от глубины.
- 9. Способ по п.7, отличающийся тем, что осуществляют частотные и амплитудные измерения кусочков породы с использованием прибора незатухающих гармонических колебаний, содержащего портативное устройство, состоящее из держателя образца с двумя пьезоэлектрическими ультразвуковыми преобразователями, расположенными поверхностями друг к другу, причем один преобразователь используют для возбуждения акустических волн, а второй преобразователь используют для обнаружения.
- 10. Способ прогнозирования в реальном времени градиента давления гидравлического разрыва пласта при бурении скважины при использовании субдюймовых кусочков породы, заключающийся в том, чтоa) получают результат измерения скорости νρ в кусочке породы, для чегоb) изготавливают образец с двумя противоположными параллельными поверхностями, разделенными расстоянием Б,c) помещают образец между преобразователями Р- и 8-волн,ά) непрерывно подают синусоидальный сигнал в пределах заданного частотного диапазона,е) определяют и регистрируют огибающую принятого амплитудно-модулированного сигнала,ί) измеряют Δί - разность частот между двумя последовательными резонансами для каждого максимума амплитуды и осуществляют подстановку Δί в уравнение νρ=2Ε·Δ£,д) используют полученный результат измерения скорости для определения пористости по скорости с помощью акустического метода пластового коэффициента согласно уравнению ίν Υ ф = 1_ Т----* р в ма трице у ггде β=0,625 и νρ в матрице=18000 футов/с для песка,
- 11) прогнозируют поровое давление на основании скорости, полученной методом незатухающих гармонических колебаний, для чего1) вычисляют суммарное вертикальное напряжение (напряжение в покрывающем пласте) путем интегрирования плотности8ν=/ρ(ζ)·§·άζ, где ζ - глубина; ρ(ζ) - объемная плотность породы пласта как функция глубины; ά - плотность; д - постоянная силы тяжести,|) получают эффективное вертикальное напряжение при использовании корреляции Холбрука σν=σтаx(1-φ) , где σ№Χ и α представляют собой постоянные параметры, зависящие от литологии, определенные Холбруком,k) определяют поровое давление согласно закону ТегаадЫРр=(8^)/у, где Рр, - поровое давление; 8ν - суммарное напряжение в покрывающем пласте, вычисленное на этапе ί; σν - эффективное вертикальное напряжение, вычисленное на этапе _р γ - коэффициент Байота, равный 1-(Κβ/Κγ), Κβ объемная сжимаемость; Κγ - сжимаемость зерен,l) прогнозируют градиент давления гидравлического разрыва пласта с помощью модифицированного способа Сат С1ау при использовании диэлектрических измерений кусочков породы (ДИКП), для чегот) измеряют диэлектрическую проницаемость кусочка породы,-9005450η) вычисляют площадь поверхности (ПП) на основании диэлектрических измерений кусочков породы согласно уравнениюПП=0,0274 - (ДИКП)1’9467;o) определяют безразмерный коэффициент М, характеризующий механическую прочность пласта, при использовании площади поверхности, полученной на основании диэлектрического измерения кусочка породы на этапе η,6з1п10ЛМ = ----------- ,3 + 3Ϊη10Α где А=1,54-1,36Е-3-(площадь поверхности),p) получают значения для минимального горизонтального напряжения пласта согласно уравнению _ 38 ν - (5ν - ЗМ)Рр ‘т1Г1 2М + 3 *11. Способ по любому из пп.1, 6, 7, 10, отличающийся тем, что осуществляют сбор кусочков породы на интервалах, при которых обеспечивается точность от ±5 до ±50 футов (от ±1,524 до ±15,24 м).
- 12. Способ по п.11, отличающийся тем, что кусочки породы собирают на интервалах, где обеспечивается точность ±30 футов (±9,144 м).
- 13. Способ по любому из пп.4, 6, 7, 10, отличающийся тем, что длина I. образца находится в пределах от около 1,5 до 4 мм.
- 14. Способ по любому из пп.1, 7, 10, отличающийся тем, что кусочки породы насыщены.
- 15. Способ по любому из пп.1, 6, 7, 10, отличающийся тем, что дополнительно выполняют вычисления согласно компьютерной программы.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US33634801P | 2001-10-24 | 2001-10-24 | |
PCT/US2002/034200 WO2003036044A1 (en) | 2001-10-24 | 2002-10-24 | Use of cutting velocities for real time pore pressure and fracture gradient prediction |
Publications (2)
Publication Number | Publication Date |
---|---|
EA200400568A1 EA200400568A1 (ru) | 2004-10-28 |
EA005450B1 true EA005450B1 (ru) | 2005-02-24 |
Family
ID=23315675
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EA200400568A EA005450B1 (ru) | 2001-10-24 | 2002-10-24 | Использование скоростей кусочков породы для прогнозирования в реальном времени порового давления и градиента давления гидравлического разрыва пласта |
Country Status (8)
Country | Link |
---|---|
US (1) | US6968274B2 (ru) |
EP (1) | EP1444416B1 (ru) |
AU (1) | AU2002335137B2 (ru) |
BR (1) | BR0213516A (ru) |
CA (1) | CA2462973C (ru) |
EA (1) | EA005450B1 (ru) |
OA (1) | OA13137A (ru) |
WO (1) | WO2003036044A1 (ru) |
Families Citing this family (38)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6751558B2 (en) * | 2001-03-13 | 2004-06-15 | Conoco Inc. | Method and process for prediction of subsurface fluid and rock pressures in the earth |
WO2003036288A1 (en) * | 2001-10-24 | 2003-05-01 | Shell Internationale Research Maatschappij B.V. | Use of cuttings for real time attenuation prediction |
WO2004095077A1 (en) | 2003-04-23 | 2004-11-04 | Commonwealth Scientific And Industrial Research Organisation | Method for predicting pore pressure |
AU2004232863B2 (en) * | 2003-04-23 | 2010-08-26 | Commonwealth Scientific And Industrial Research Organisation | Method for predicting pore pressure |
US7299132B2 (en) * | 2005-08-08 | 2007-11-20 | Schlumberger Technology Corp. | Method and system for pre-drill pore pressure prediction |
US7844405B2 (en) * | 2005-09-08 | 2010-11-30 | Africk Steven A | Static and ultrasonic methods for measuring compressibility and density of small particles |
US8627889B2 (en) | 2007-09-27 | 2014-01-14 | Schlumberger Technology Corporation | Drilling and fracturing fluid |
US8499830B2 (en) * | 2008-07-07 | 2013-08-06 | Bp Corporation North America Inc. | Method to detect casing point in a well from resistivity ahead of the bit |
US7861801B2 (en) * | 2008-07-07 | 2011-01-04 | Bp Corporation North America Inc. | Method to detect coring point from resistivity measurements |
US8061442B2 (en) * | 2008-07-07 | 2011-11-22 | Bp Corporation North America Inc. | Method to detect formation pore pressure from resistivity measurements ahead of the bit during drilling of a well |
GB201001833D0 (en) | 2010-02-04 | 2010-03-24 | Statoil Asa | Method |
US8229671B2 (en) * | 2009-08-13 | 2012-07-24 | Pritchard David M | Method and system for riserless casing seat optimization |
US8494827B2 (en) * | 2009-09-25 | 2013-07-23 | Exxonmobil Upstream Research Company | Method of predicting natural fractures and damage in a subsurface region |
US9051815B2 (en) * | 2009-09-28 | 2015-06-09 | Baker Hughes Incorporated | Apparatus and method for predicting vertical stress fields |
US8214152B2 (en) * | 2009-09-28 | 2012-07-03 | Baker Hughes Incorporated | Apparatus and method for predicting vertical stress fields |
EP2591381A4 (en) * | 2010-07-08 | 2015-12-02 | Geco Technology Bv | METHODS AND DEVICES FOR TRANSFORMING COLLECTED DATA TO IMPROVE VISUALIZATION CAPACITY |
US9238962B2 (en) * | 2010-12-21 | 2016-01-19 | Schlumberger Technology Corporation | Pore pressure from spectroscopy and sonic data |
US9222350B2 (en) | 2011-06-21 | 2015-12-29 | Diamond Innovations, Inc. | Cutter tool insert having sensing device |
US20140321240A1 (en) * | 2013-04-26 | 2014-10-30 | Siemens Medical Solutions Usa, Inc. | Elastography for cement integrity inspection |
RU2505675C1 (ru) * | 2012-09-03 | 2014-01-27 | Шлюмберже Текнолоджи Б.В. | Способ определения свойств углеводного пласта и добываемых флюидов в процессе добычи |
CN103089253B (zh) * | 2013-01-22 | 2015-06-03 | 中国石油大学(北京) | 一种利用小波变换计算地层孔隙压力的方法 |
US10145216B2 (en) | 2014-12-02 | 2018-12-04 | Landmark Graphics Corporation | Determining dominant scenarios for slowing down trip speeds |
US10190406B2 (en) * | 2014-12-23 | 2019-01-29 | Baker Hughes, A Ge Company, Llc | Formation fracturing potential using surrounding pore pressures |
US10400591B2 (en) | 2016-05-24 | 2019-09-03 | Saudi Arabian Oil Company | Systems and methods for acoustic testing of laminated rock to determine total organic carbon content |
CN106522919B (zh) * | 2016-11-07 | 2019-07-02 | 中国石油大学(北京) | 三探头密度测井的围岩校正方法及装置 |
CN106526697B (zh) * | 2016-12-07 | 2019-03-05 | 中国地质大学(北京) | 一种三角洲建设与破坏主要转换面识别方法 |
CN110320574B (zh) * | 2018-03-30 | 2021-06-15 | 中国石油化工股份有限公司 | 基于缓坡三角洲薄层砂体刻画的方法 |
CN109667573B (zh) * | 2018-12-12 | 2022-07-15 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 三维页岩储层孔隙压力预测方法、装置和电子设备 |
CN110485977A (zh) * | 2019-08-15 | 2019-11-22 | 中石化石油工程技术服务有限公司 | 快速预测页岩气层地层破裂压力梯度的测井方法 |
CN110410069B (zh) * | 2019-08-15 | 2022-08-19 | 中石化石油工程技术服务有限公司 | 页岩气水平井地层破裂压力梯度测井预测方法 |
CN110471104B (zh) * | 2019-08-26 | 2021-03-16 | 电子科技大学 | 基于智能特征学习的叠后地震反射模式识别方法 |
CN111255435B (zh) * | 2020-01-17 | 2021-09-07 | 西安石油大学 | 一种复杂储层泥质含量计算方法 |
CN111594101B (zh) * | 2020-03-19 | 2021-03-12 | 中国石油大学(华东) | 一种垂直井管外出砂亏空剖面预测及高饱和砾石充填施工方法 |
CN111502647B (zh) * | 2020-03-27 | 2021-02-02 | 中国石油化工股份有限公司石油工程技术研究院 | 一种钻井地质环境因素确定方法、装置及存储介质 |
CN112065361B (zh) * | 2020-09-16 | 2021-03-12 | 成都理工大学 | 一种基于声波衰减的致密储层含气饱和度确定方法 |
US11719094B2 (en) * | 2020-10-23 | 2023-08-08 | Saudi Arabian Oil Company | Reservoir characterization using rock geochemistry for lithostratigraphic interpretation of a subterranean formation |
CN114396257B (zh) * | 2021-12-31 | 2023-10-31 | 中国石油大学(北京) | 基于地震数据的三维地层压力预测方法及装置 |
CN114783531B (zh) * | 2022-04-02 | 2022-11-22 | 广州海洋地质调查局 | 计算砂岩型水合物中胶结型和孔隙充填型饱和度的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4876512A (en) * | 1988-03-30 | 1989-10-24 | Exxon Production Research Company | Method for rapidly determining the swelling-clay content in shales and shaly sandstone formations by high-frequency dielectric constant measurements |
US4981037A (en) * | 1986-05-28 | 1991-01-01 | Baroid Technology, Inc. | Method for determining pore pressure and horizontal effective stress from overburden and effective vertical stresses |
US5200929A (en) * | 1992-03-31 | 1993-04-06 | Exxon Production Research Company | Method for estimating pore fluid pressure |
US5282384A (en) * | 1992-10-05 | 1994-02-01 | Baroid Technology, Inc. | Method for calculating sedimentary rock pore pressure |
US5285692A (en) * | 1990-08-31 | 1994-02-15 | Exxon Production Research Company | Methods for measuring physical parameters of a low permeability rock formation in situ |
-
2002
- 2002-10-24 US US10/493,370 patent/US6968274B2/en not_active Expired - Fee Related
- 2002-10-24 EA EA200400568A patent/EA005450B1/ru not_active IP Right Cessation
- 2002-10-24 BR BR0213516-7A patent/BR0213516A/pt not_active Application Discontinuation
- 2002-10-24 CA CA2462973A patent/CA2462973C/en not_active Expired - Fee Related
- 2002-10-24 OA OA1200400118A patent/OA13137A/en unknown
- 2002-10-24 WO PCT/US2002/034200 patent/WO2003036044A1/en active IP Right Grant
- 2002-10-24 EP EP02802212A patent/EP1444416B1/en not_active Expired - Fee Related
- 2002-10-24 AU AU2002335137A patent/AU2002335137B2/en not_active Ceased
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4981037A (en) * | 1986-05-28 | 1991-01-01 | Baroid Technology, Inc. | Method for determining pore pressure and horizontal effective stress from overburden and effective vertical stresses |
US4876512A (en) * | 1988-03-30 | 1989-10-24 | Exxon Production Research Company | Method for rapidly determining the swelling-clay content in shales and shaly sandstone formations by high-frequency dielectric constant measurements |
US5285692A (en) * | 1990-08-31 | 1994-02-15 | Exxon Production Research Company | Methods for measuring physical parameters of a low permeability rock formation in situ |
US5200929A (en) * | 1992-03-31 | 1993-04-06 | Exxon Production Research Company | Method for estimating pore fluid pressure |
US5282384A (en) * | 1992-10-05 | 1994-02-01 | Baroid Technology, Inc. | Method for calculating sedimentary rock pore pressure |
Non-Patent Citations (3)
Title |
---|
HOLBROOK P.W. ET AL.: "A Petrophysical-Mechanical Math Model for Real-Time Wellsite Pore Pressure/Fracture Gradient Prediction", SPE ANNUAL TECHNICAL CONFERENCE AND EXHIBITION, no. 16666, 27 September 1987 (1987-09-27) - 30 September 1997 (1997-09-30), pages 163-171, XP002229035, Dallas page 164, right-hand column * |
NES. O.M. ET AL: "Rig-site and Laboratory use of CWT acoustic Velocity Measurements on Cuttings", SPE ANNUAL TECHNICAL CONFERENCE AND EXHIBITION, no. 36854, 22 - 24 October 1996, pages 349-356, XP002229036, Milan cited in the application page 349, right-hand column * |
RAIGA-CLEMENCEAU J. ET AL.: "The concept of acoustic formation factor for more accurate porosity determination from sonic transit time data", SPWLA TWENTY-SEVENTH ANNUAL LOGGING SYMPOSIUM, no. Paper G, 9-13 June 1986, XP002179122, abstract * |
Also Published As
Publication number | Publication date |
---|---|
CA2462973C (en) | 2012-09-11 |
CA2462973A1 (en) | 2003-05-01 |
EP1444416A1 (en) | 2004-08-11 |
US20040236513A1 (en) | 2004-11-25 |
WO2003036044A1 (en) | 2003-05-01 |
AU2002335137B2 (en) | 2007-11-22 |
EP1444416B1 (en) | 2006-10-11 |
BR0213516A (pt) | 2004-10-19 |
EA200400568A1 (ru) | 2004-10-28 |
OA13137A (en) | 2006-12-13 |
US6968274B2 (en) | 2005-11-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EA005450B1 (ru) | Использование скоростей кусочков породы для прогнозирования в реальном времени порового давления и градиента давления гидравлического разрыва пласта | |
AU2002335137A1 (en) | Use of cutting velocities for real time pore pressure and fracture gradient prediction | |
US6681185B1 (en) | Method of seismic signal processing | |
US6098021A (en) | Estimating formation stress using borehole monopole and cross-dipole acoustic measurements: theory and method | |
US8004932B2 (en) | Identification of stress in formations using angles of fast and slow dipole waves in borehole acoustic logging | |
WO2009055196A2 (en) | In-situ determination of yield stress state of earth formations | |
US8553493B2 (en) | Method for permeable zone detection | |
US8681582B2 (en) | Method for sonic indication of formation porosity and lithology | |
WO2004095077A1 (en) | Method for predicting pore pressure | |
Nes et al. | Rig site and laboratory use of CWT acoustic velocity measurements on cuttings | |
US8077545B2 (en) | Method for detecting gas influx in wellbores and its application to identifying gas bearing formations | |
CA2837493C (en) | Determining the location of a material located behind an obstruction | |
Marsala et al. | Sonic while drilling: Have you thought about cuttings? | |
AU2002349944B2 (en) | Use of cutting for real time attenuation prediction | |
Donald et al. | Adaptation of a triaxial cell for ultrasonic P-wave attenuation, velocity and acoustic emission measurements | |
AU2002349944A1 (en) | Use of cutting for real time attenuation prediction | |
US11366049B2 (en) | Estimation of objective driven porous material mechanical properties | |
Sinha et al. | Applications of sonics and ultrasonics in geophysical prospecting | |
Brata et al. | Mechanical Earth Model with Uniaxial Core Test Validation from HPHT Exploration Field in North Sumatra | |
GB2174201A (en) | Method of measuring pore pressure in impermeable rocks | |
NO20191505A1 (en) | Method for estimating rock brittleness from well-log data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | Lapse of a eurasian patent due to non-payment of renewal fees within the time limit in the following designated state(s) |
Designated state(s): AM AZ BY KG MD TJ TM |
|
MM4A | Lapse of a eurasian patent due to non-payment of renewal fees within the time limit in the following designated state(s) |
Designated state(s): KZ RU |