Использование скоростей кусочков породы для прогнозирования в реальном времени порового давления и градиента давления гидравлического разрыва пласта

005450 Область техники Настоящее изобретение относится к анализу и интерпретации скоростей акустических волн в образцах породы, полученных при проходке. Более конкретно, изобретение относится к новому способу очень точного прогнозирования в реальном времени порового давления, градиента давления гидравлического разрыва и прочности пласта на месте нахождения буровой установки путем определения скорости кусочка породы посредством портативного прибора для измерений методом незатухающих гармонических колебаний (НГК) на высокой резонансной частоте и используя измерение скорости в сочетании со способом настоящего изобретения для получения точных прогнозов значений порового давления и градиента давления гидравлического разрыва пласта. При использовании нового способа обеспечивается прогнозирование в реальном времени порового давления на месте нахождения буровой установки с небольшим полем ошибки, который нельзя получить, используя сейсмическую разведку, вертикальное сейсмическое профилирование или контроль скорости по взрыву. Предшествующий уровень техники Решающее значение во время бурения скважины имеет знание точных значений порового давления,градиента давления гидравлического разрыва и прочности пласта. Кроме того, поровое давление и градиент давления гидравлического разрыва пласта являются управляющими входными параметрами при моделировании устойчивости ствола скважины, проектировании и расчете скважины и оптимизации траектории ствола скважины. Хотя отсутствуют промышленные приборы для измерения порового давления перед буровой коронкой, но разработаны методики для вычисления порового давления на интервалах проведения каротажа по данным удельного сопротивления и/или акустического каротажа. Известно несколько способов для получения информации во время бурения, касающейся других характеристик. Например, известен способ определения пористости пласта по характеристикам бурения. В патенте США 4064749 раскрыт способ определения пористости пласта по характеристикам бурения,в котором при вращении буровой коронки, прикрепленной к нижнему концу бурильной колонны, управляют направленной вниз силой, действующей на буровую коронку. Способ включает в себя этапы измерения частоты вращения буровой коронки, измерения глубины погружения буровой коронки в стволе скважины, измерения массы на указанной буровой коронке, определяющей затупление зубьев указанной буровой коронки, измерения вращающего момента, приложенного к бурильной колонне, определения эмпирически эталонного вращающего момента и определения пористости путем объединения указанных результатов измерений и определений и использования уравнения. В патенте США 4949575 раскрыт способ анализа пласта, в котором результаты измерений характеристик бурения используют в качестве входных данных, зависящих от пористости. Способ включает в себя получение сигнала бурения, характеризующего сопротивление бурению пласта, пробуриваемого буровой коронкой, получение большого количества дополнительных сигналов, отражающих свойства пласта, и выполнение на основании указанного сигнала бурения и указанных дополнительных сигналов волюметрического анализа приповерхностного пласта. В патенте США 4876512 раскрыт способ определения на буровой площадке содержания разбухающей глины в сланцах и в сланцеватых песчаниках путем измерений площади проводящей поверхности. Образцы промывают жидкостью, имеющей существенно меньшую активность по сравнению с активностью воды, и которая может содержать растворимый катион, и на заранее выбранной частоте (1 МГц) проводят измерения диэлектрической проницаемости образцов для последующего сравнения с калибровочными кривыми,тем самым получая данные по содержанию разбухающей глины в пласте. В патенте США 5282384 (Холбрук) раскрыт усовершенствованный способ, основанный на принципах механического зондирования, с использованием данных каротажа скважины для вычисления давления жидкости, содержащейся в осадочной породе, которая естественным образом уплотнилась под действием силы тяжести. Часть способа, относящаяся к определению эффективного напряжения, включает в себя внутренние и внешние измерения деформации зернистой матрицы породы. Поэтому одна и та же калибровка эффективного напряжения может быть равным образом применена к данным о толщине породы, полученным путем измерения, и к петрофизическим данным о пористости породы. Степенная зависимость между эффективным напряжением и деформацией для любой осадочной породы может быть определена из взвешенного среднего степенных функций для минералов, которые образуют осадочную породу. В настоящем изобретении для прогнозирования порового давления и градиента давления гидравлического разрыва пласта способ Холбрука модифицирован как в части методики, так и в части типа данных. Проблема, часто встречающаяся при бурении скважин во многих частях мира, заключается в небольших полях ошибок, что требует большой точности при прогнозировании как порового давления, так и градиента давления гидравлического разрыва пласта, чтобы предотвратить нестабильность сланцев,создающую опасность прекращения циркуляции и/или выброса газа/фонтанирования. Например, при большой глубине Мексиканского залива допуск по бурению, выраженный через поровое давление и через градиент давления гидравлического разрыва пласта, может быть меньше 0,5 фунта/галлон. Поэтому точность скоростей, получаемых на основании данных сейсмических исследований, каротажа во время-1 005450 бурения и/или кабельного каротажа, является очень важной для бурения в дополнение к другим известным применениям скорости для петрофизических или пластовых исследований. Кусочки породы, образующиеся при бурении, представляют собой источник информации в квазиреальном времени, который можно получать на месте нахождения буровой установки. В прошлом использование кусочков породы ограничивалось из-за трудностей, возникавших при выполнении измерений на очень малых образцах. В отчете Nes et al., "Rig-site and laboratory use of CWT acoustic velocity measurements on cuttings". Society of Petroleum Engineers Paper36854, 1996, показана возможность использования методики незатухающих гармонических колебаний для измерения скорости акустических волн в кусочках породы при использовании портативного оборудования, которое пригодно для контроля небольших образцов, и тем самым предложен новый источник данных, который может быть использован в квазиреальном времени на месте нахождения буровой установки. Для калибровки скоростей сейсмических волн, данных каротажа во время бурения и/или результатов кабельного акустического каротажа и, что еще более важно, для получения в реальном времени точных результатов измерений скорости волн в пласте кусочки породы обеспечивают потенциально бесценный источник информации. В настоящее время необходим способ, использование которого позволяет точно прогнозировать поровое давление и градиент давления гидравлического разрыва пласта путем выполнения в реальном времени измерений на месте нахождения буровой установки. В случае получения таких данных они также будут полезными для распознавания опасности мелководных зон, оптимизации массы бурового раствора, обнаружения опасных пустот, обнаружения зон аномального давления, определения прочности пласта для оптимизации траектории скважины и в общем случае для осуществления наиболее безопасного, экономически эффективного бурения. Краткое изложение существа изобретения Задачей настоящего изобретения является создание способа прогнозирования в реальном времени порового давления и градиента давления гидравлического разрыва пласта на месте нахождения буровой установки путем использования метода незатухающих гармонических колебаний на резонансных частотах для получения в реальном времени значения скорости кусочка породы и путем использования данных в сочетании со способом, раскрытым в настоящем изобретении для точного прогнозирования порового давления и градиента давления гидравлического разрыва пласта. Кроме того, значения скоростей также используют в реальном времени для определения механических характеристик и прочности пласта. Способ, согласно изобретению, заключается в том, что 1) получают результаты измерений скорости Vp кусочка породы, для чегоa) подготавливают образец с двумя противоположными параллельными поверхностями, разделенными расстоянием L,b) помещают образец между преобразователями Р- и S-волн,c) непрерывно пропускают синусоидальный сигнал в пределах надлежащего частотного диапазона,d) извлекают и регистрируют огибающую принятого амплитудно-модулированного сигнала,e) измеряют f, ширину кривой на уровне половинной мощности для каждого максимума амплитуды и осуществляют подстановку значения в уравнение Vp=2Lf,2) используют полученные результаты измерения скорости для определения пористости по скорости с помощью акустического метода пластового коэффициента при использовании уравнения где =0,625 и Vp в матрице =18000 футов/с для песка,3) прогнозируют поровое давление на основании зарегистрированной скорости, полученной методом незатухающих гармонических колебаний, для чегоa) вычисляют суммарное вертикальное напряжение (напряжение в покрывающем слое) путем интегрирования плотностиSv=(z)gdz,где z - глубина; p(z) - объемная плотность породы пласта как функция глубины; d - плотность; g постоянная силы тяжести,b) получают эффективное вертикальное напряжение при использовании корреляции Холбрукаv=max(1-),где max ипредставляют собой постоянные параметры, зависящие от литологии, определенные Холбруком,с) определяют поровое давление при использовании закона TerzaghiPp=(Sv-v)/,где Рр - поровое давление; Sv - суммарное напряжение в покрывающем пласте, вычисленное на этапе 3 а; v - эффективное вертикальное напряжение, вычисленное на этапе 3b;- коэффициент Байота,равный 1-(К/К), где К - объемная сжимаемость, К - сжимаемость зерен,4) прогнозируют градиент давления гидравлического разрыва пласта с помощью модифицированного способа Холбрука при использовании скорости кусочка породы, полученной методом незатухающих гармонических колебаний, из уравненияSh min=[Pp+v(1-)],где Sh min - минимальное горизонтальное напряжение, которое представляет собой нижний предел градиента давления гидравлического разрыва пласта;- пористость, полученная на этапе 2; Рр - поровое давление, полученное на этапе 3(с); v - эффективное напряжение, вычисленное на этапе 3(b);показатель локальной прочности (для района разведочных работ используют 0,85 по умолчанию),5) или прогнозируют градиент давления гидравлического разрыва пласта с помощью модифицированного способа Cam Clay при использовании диэлектрических измерений кусочков породы (ДИКП),для чегоa) измеряют диэлектрическую проницаемость кусочка породы,b) вычисляют площадь поверхности (ПП) на основании диэлектрических измерений кусочков породы при использовании уравнения ПП=0,0274(ДИКП)1,9467,с) определяют безразмерный коэффициент М, характеризующий механическую прочность пласта,при использовании площади поверхности, полученной на основании диэлектрических измерений кусочков породы на этапе 5(b), из уравнения где A=1,54-1,36 Е-3(площадь поверхности),d) получают значения минимального горизонтального напряжения при использовании уравнения Краткое описание чертежей В дальнейшем изобретение поясняется описанием предпочтительного варианта выполнения со ссылками на сопровождающие чертежи, на которых фиг. 1 изображает схему экспериментальной установки, согласно изобретению; фиг. 2 - резонансный спектр частот для кусочка сланца толщиной 2,5 мм, согласно изобретению; фиг. 3 - диаграмму зависимости скоростей, полученных методом незатухающих гармонических колебаний, от частоты для контрольных материалов - люцита и соли, согласно изобретению; фиг. 4 - диаграмму зависимости скоростей, полученных в реальном времени методом незатухающих гармонических колебаний на скважине в Мексиканском заливе, согласно изобретению; фиг. 5 - спрогнозированные в реальном времени значения порового давления и градиента давления гидравлического разрыва пласта при использовании данных полученных кусочков породы методом незатухающих гармонических колебаний и диэлектрических измерений, в зависимости от глубины на выбранных интервалах скважины в Мексиканском заливе, согласно изобретению; фиг. 6 - диаграмму зависимости значений скоростей, полученных методом незатухающих гармонических колебаний и кабельным акустическим каротажем, от глубины на выбранных интервалах скважины в Мексиканском заливе, согласно изобретению; фиг. 7 - диаграмму зависимости спрогнозированных значений порового давления и градиента давления гидравлического разрыва пласта от глубины для скважины в Мексиканском заливе после бурения интервала и сбора данных кабельного акустического каротажа; прогнозы по методике незатухающих гармонических колебаний и на основе диэлектрических измерений кусочков породы сделаны в реальном времени во время бурения на том же самом интервале и приведены для сравнения, согласно изобретению; фиг. 8 - диаграмму зависимости спрогнозированных значений порового давления и градиента давления гидравлического разрыва пласта от глубины для другой скважины в Мексиканском заливе, полученных при использовании данных метода незатухающих гармонических колебаний и диэлектрических измерений кусочков породы. Подробное описание предпочтительных вариантов воплощения изобретения Использование заявленного способа согласно изобретению обеспечивает ряд преимуществ по сравнению с любым способом, известным в настоящее время. Прежде всего, обеспечивается возможность точного прогнозирования в реальном времени порового давления и градиента гидравлического разрыва пласта при использовании результатов измерений на образцах реальной породы на месте нахождения буровой установки. Способ также обеспечивает возможность прогнозирования прочности пласта в буро-3 005450 вой скважине без проведения акустического каротажа, что предоставляет операторам возможность получать все входные данные, необходимые для анализа устойчивости ствола скважины, и проводить анализ устойчивости ствола скважины, а для успешного бурения скважины давать рекомендации относительно оптимальных давлений в скважине. Потенциально эти особенности обеспечивают значительное снижение затрат на бурение, связанных с точной оптимизацией траектории скважины. Кроме того, результаты измерений могут быть использованы для сравнения с получаемыми позднее данными каротажа при бурении для анализа соответствия, при этом выявление различий может способствовать определению в реальном времени опасных зон, смягчению возникающих проблем, предотвращению неустойчивости ствола скважины и минимизации времени простоя при эксплуатации. Прибор для измерений методом незатухающих гармонических колебаний, используемый в изобретении для измерения характеристик небольших образцов кусочков породы на месте нахождения буровой установки, представляет собой портативное устройство, содержащее держатель образца с двумя пьезоэлектрическими ультразвуковыми преобразователями, расположенными навстречу друг другу. Один преобразователь используют для возбуждения акустических волн, авторой преобразователь используют для их обнаружения. Один из преобразователей закреплен на микрометрическом винте микрометра, что позволяет одновременно измерять длину образца. Образец располагают между двумя преобразователями. При подаче частоты возбуждения в пределах диапазона, соответствующего нескольким резонансам стоячих волн, в составном резонаторе создаются резонансы ультразвуковых стоячих волн. Огибающую принятого сигнала усиливают с помощью широкополосного приемника с регулируемым коэффициентом усиления, который подключают к переносному компактному персональному компьютеру. При поступлении нескольких входных параметров с помощью вспомогательной программы вычисляются скорости. Пример результатов измерений, выполненных для образца сланца из бурового шлама толщиной 2,5 мм(0,5 дюйма), показан на фиг. 2. Портативное устройство для измерений методом незатухающих гармонических колебаний, использованное в настоящем изобретении (фиг. 1), можно получить от ТЕМСО Company, Талса, штат Оклахома. Ниже приведены определения терминов, используемых в описании изобретения:- затухание звуковых колебаний, Нп/см; НГК - незатухающие гармонические колебания; МНГК метод незатухающих гармонических колебаний; f - частота, Гц, с-1; D - диаметр, м; L, z - толщина образца, м;- длина волны, м; ИГ СБР - измеренная глубина стола бурового ротора, м; ФВГфактическая вертикальная глубина; n - число половин длин волн, целое число; Vp - скорость звука, м/с, мкс/фут; Р первичная волна; S -поперечная/вторичная волна; v - вертикальное эффективное напряжение;- пористость; mах и- постоянные параметры, зависящие от литологии; Sv - суммарное напряжение в покрывающем пласте; (z) - объемная плотность породы как функция глубины; z - глубина; g - постоянная силы тяжести; d - плотность; Рр - прогнозируемое поровое давление; Рр - поровое давление, где- коэффициент Байота=1-(K/K), K - объемная сжимаемость, K -сжимаемость зерен. Нижние индексы:n - число половин длин волн; р - первичный. Устройство для измерения методом незатухающих гармонических колебаний (МНГК) основано на создании резонансов стоячих волн в составном резонаторе при подаче частоты f возбуждения в пределах диапазона, соответствующего нескольким резонансам стоячих волн. Посредством второго преобразователя принимается амплитудно-модулированный сигнал, включающий резонансы, показанные на фиг. 2. Каждый из n резонансов стоячих волн характеризуется резонансной частотой fn. На практике при наличии плоских волн предполагается, что отношение /D является достаточно малым, где- длина акустической волны, D - диаметр. Резонанс будет наблюдаться, по существу, на каждом интервале =v/f, чему соответствует L=n/2, где L - толщина образца или расстояние между поверхностями двух преобразователей, n - целое число. Фазовая скорость акустических волн для Р- или S-волн может быть выражена какVp = 2Lf, где f = fn+1-fn представляет собой разность частот двух последовательных резонансов. Например, для образца сланца с Vp=6550 футов/с и L=0,06 дюйма f будет составлять 0,67 МГц, а практическая частота будет находиться в пределах 1-10 МГц. При подготовке к измерениям оператор получает образец кусочка породы с двумя противоположными параллельными поверхностями, отстоящими друг от друга на расстоянии L (расстояние L соответствует длине образца кусочка породы), прикрепляет к каждой из двух поверхностей преобразователи Рили S-волн, в соответствующем частотном диапазоне непрерывно изменяет частоту синусоидального сигнала, выделяет и регистрирует огибающую принятого амплитудно-модулированного сигнала, измеряет f и подставляет значения в уравнение Vp=2L-f. Подготовка образца включает в себя формирование двух гладких и параллельных поверхностей и может быть выполнена шкуркой. Кусочки породы следует отбирать внимательно, чтобы уменьшить число ошибок, обусловленных неопределенностью глубины, неоднородностями и механическими повреждениями кусочков буровой породы и/или обрушениями. Можно подготавливать и измерять образцы с мелкими зернами толщиной в долю дюйма меньше, чем через 10 мин после помещения бурового шлама на вибрационное сито.-4 005450 Для проведения измерений образец зажимают между двумя преобразователями. Микрометр прикрепляют к металлической раме и используют для позиционирования преобразователя, тем самым обеспечивают возможность измерения длины образца под нагрузкой. Чтобы гарантировать достаточное качество сигнала без влияния на характеристики образца, на обе поверхности подготовленного образца наносят соответствующим образом выбранный гель для улучшения акустической связи. Требуемый диапазон подаваемой частоты обеспечивают посредством генератора частот. Огибающую принятого сигнала усиливают с помощью широкополосного приемника с регулируемым коэффициентом усиления, который подключают к переносному компактному персональному компьютеру. Скорости вычисляют с помощью вспомогательной программы компьютера при поступлении нескольких входных параметров. Использование способа согласно настоящему изобретению обеспечивает получение точных данных в реальном времени на образцах меньше дюйма. Однако было замечено, что, хотя значительные изменения точности отсутствуют, для образцов от 2,5 до 3 мм имеется возможность получать наиболее точные данные. Для образцов меньше 2,5 мм или образцов, которые толще 4-5 мм, иногда получают неверные данные. После определения значений скоростей для подготовленных образцов, имеющих предпочтительную длину, каждую точку данных скорости преобразуют в прогнозируемое поровое давление, используя модифицированный способ Холбрука. Модифицированный способ Холбрука представляет собой один из наиболее надежных способов для прогнозирования порового давления. Способ Холбрука выгодно отличается независимостью от выбора ориентации сланца, что требуется для большинства других известных способов. При использовании этого способа непосредственно учитываются литологические изменения при прогнозировании порового давления. Модифицированный способ Холбрука основан на трех основных уравнениях. Первое уравнение представляет собой определение вертикального эффективного напряжения, то есть закон эффективного напряжения Terzaghi где Sv - суммарное напряжение в покрывающем пласте; Рр -поровое давление;- коэффициент Байота,равный 1-(K/K), K - объемная сжимаемость, K - сжимаемость зерен. Sv обычно получают путем интегрирования плотности по данным каротажа, образуя второе уравнение Третье уравнение представляет собой эмпирическое уравнение, полученное на основании расширенного множества данных каротажа скважины и измерений пластового давления. Этим уравнением связь эффективного напряжения с пористостью пласта представлена в виде где max ипредставляют собой постоянные параметры, зависящие от литологии (Холбрук, 1996). Для уравнения (3) в качестве исходных данных требуется пористость. Пористость может быть измерена непосредственно в лаборатории при использовании образцов керна или может быть получена из разнообразных каротажных данных (сопротивления, проводимости или акустических), сейсмических данных или из данных лабораторных измерений кернов, или из скорости кусочков породы. Предпочтительно определять пористость по скорости, чтобы минимизировать влияние изменений температуры и солености на величины каротажных данных. Кроме того, небольшие изменения пористости, связанные с пластами, находящимися под аномальными давлениями, в значительной степени отражаются на акустических характеристиках. Применив метод незатухающих гармонических колебаний,скорости, полученные на основании измерений, можно использовать для определения пористости и объемной сжимаемости. Пористость может быть определена по скорости с помощью акустического пластового коэффициента согласно уравнению где =0,625 и Vp в матрице=18000 футов/с для песка. После определения пористости уравнение (3) используют для получения вертикального эффективного напряжения v. После нахождения вертикального эффективного напряжения уравнение (1) используют для прогнозирования порового давления. Во втором варианте осуществления изобретения можно спрогнозировать градиент давления гидравлического разрыва пласта путем выполнения двух различных вычислений. В первом случае градиент давления гидравлического разрыва пласта прогнозируют с помощью модифицированного способа Холбрука, используя скорость кусочков породы, определенную методом незатухающих гармонических колебанийSh min=[Pp+v(l-)] где Sh min - минимальное горизонтальное напряжение, которое характеризует нижний предел градиента гидравлического разрыва пласта;- пористость; Рp - поровое давление; v -эффективное напряже-5 005450 ние, вычисление которого рассматривалось выше;- коэффициент локального напряжения (0,85 по умолчанию для района разведочных работ). Альтернативно градиент давления гидравлического разрыва пласта может быть спрогнозирован с помощью модифицированного метода Cam Clay при использовании результатов диэлектрических измерений кусочков породы (ДИКП) следующим образом:a) измеряют диэлектрическую проницаемость кусочка породы,b) на основании результата диэлектрических измерений вычисляют площадь поверхности (ПП) ПП=0,0274- (ДИКП)1,9467,c) определяют безразмерный коэффициент М, характеризующий механическую прочность пласта,при использовании площади поверхности, полученной на основании диэлектрических измерений кусочка породы где А=1,54-1,36 Е-3(площадь поверхности),d) получают значения минимального горизонтального напряжения согласно уравнению После вычисления прогнозируемых значений порового давления и градиента давления гидравлического разрыва пласта значения могут быть нанесены на график в зависимости от глубины. Дополнительно это рассмотрено в примерах 2 и 3 и отображено на фиг. 5, 7 и 8. В настоящем изобретении метод незатухающих гармонических колебаний для измерения скоростей сочетается с новыми способами, разработанными заявителем, что обеспечивает возможность точного прогнозирования порового давления и градиента давления гидравлического разрыва пласта по скорости продольных волн. Способ согласно настоящему изобретению может быть использован для проверки прогноза порового давления на основании скоростей сейсмических волн в районе разведочных работ, а также для калибровки сейсмических данных до того, как станут доступны данные акустического каротажа. Настоящее изобретение будет также весьма полезным для обнаружения слабых зон, которые предрасположены к быстрому увеличению и утрате циркуляции, и для решения проблем устойчивости ствола скважины, например заклинивания труб и налипания породы на долото. Способ согласно настоящему изобретению является удобным, относительно недорогим и простым в реализации. Используя оборудование, которое может быть доставлено на буровую установку в портфеле,и используя образцы размером меньше дюйма, можно в течение нескольких минут определить скорости и преобразовать их в предсказанные значения порового давления и градиента давления гидравлического разрыва пласта. Опыт работы с прибором для измерений методом незатухающих гармонических колебаний на разведочно-эксплуатационной скважине во время работы буровой установки показал, что существует значительное отклонение получаемых результатов от интервала скоростей сейсмических волн, однако было обнаружено, что они намного лучше согласуются с данными кабельного акустического каротажа, когда указанные данные кабельного акустического каротажа собирались на более поздней стадии. При использовании изобретения на месторождении было замечено, что при интервалах глубин около 30 футов в большинстве случаев точность была достаточной, но для повышения точности кусочки породы можно собирать с более частыми интервалами, например в пределах глубин 10 футов. Кроме того, измерения выполняются при напряжении по одной оси, а не по трем осям, однако было обнаружено, особенно для сланцев, небольшое различие между результатами, полученными при напряжении по одной оси и по трем осям, не приводит к значительному отклонению от реального давления. Нижеследующие примеры служат иллюстрацией изобретения, раскрытого в настоящей заявке. Эксперимент Пример 1. Контрольные материалы исследовались в лаборатории и на месте нахождения буровой установки для проверки точности методики измерения скорости (методом незатухающих гармонических колебаний), которая является источником исходных данных для способа, рассмотренного в настоящем изобретении. Результаты контрольных исследований представлены на фиг. 3. Образцы соли, использованные при измерениях, имели примесь. Поэтому были получены два различных значения f, соответствующие чистой поваренной соли и примеси. На базе этих двух различных значений f зависимости скорости от частоты получены методом незатухающих гармонических колебаний для двух различных скоростей. Пример 2. Прибор проверяли в полевых условиях на месте нахождения буровой установки в Мексиканском заливе, а измерения в реальном времени методом незатухающих гармонических колебаний выполняли для отдельных интервалов глубин буровой скважины. Из-за больших размеров скважины акустический каротаж во время бурения не проводили. Поэтому скорости, определяемые в реальном времени методом незатухающих гармонических колебаний, были единственной основой во время бурения для калибровки порового давления, предсказанного до бурения на основании сейсмических данных. Полученные зави-6 005450 симости скорости от глубины для отдельных интервалов представлены на фиг. 4. Соответствующиедиаграммы порового давления и градиента давления гидравлического разрыва пласта, полученные на основании измерений скоростей методом незатухающих гармонических колебаний, представлены на фиг. 5, а полученные на основании измерений диэлектрической проницаемости при использовании кусочков породы, показаны на фиг. 6. По окончании бурения интервала собирали данные кабельного акустического каротажа. Сравнение скоростей, полученных с помощью кабельного зонда, и скоростей, измеренных методом незатухающих гармонических колебаний, представлено на фиг. 6. Сравнение итоговых значений порового давления и градиента давления гидравлического разрыва пласта, полученных на основании сейсмических данных, данных кабельного акустического каротажа,скоростей, измеренных методом незатухающих гармонических колебаний, и градиента давления гидравлического разрыва пласта, полученного на основании диэлектрических измерений кусочков породы,приведено на фиг. 7. На фиг. 7 видно, что спрогнозированные в реальном времени значения порового давления и градиента гидравлического разрыва пласта на основании измерений методом незатухающих гармонических колебаний и диэлектрических измерений кусочков породы являются точными и обеспечивают чрезвычайно полезную информацию для подтверждения глубины нахождения уступа в скважине для опоры башмака обсадной колонны, массы бурового раствора, для анализа устойчивости ствола скважины и другой буровой активности с целью безопасного бурения. Пример 3. Кусочки породы из другой скважины в Мексиканском заливе использовали для прогнозирования порового давления и градиента давления гидравлического разрыва пласта при использовании скорости,измеренной методом незатухающих гармонических колебаний, и данных диэлектрических измерений кусочков породы. Спрогнозированные значения порового давления и градиента давления гидравлического разрыва пласта, полученные на основании скоростей и данных диэлектрических измерений кусочков породы, представлены на фиг. 8, как и полученные из каротажа профили порового давления и градиента давления гидравлического разрыва пласта, а также реальные данные о давлении жидкости, измеренные при периодических испытаниях пласта, и данные испытаний на просачивание, уже проведенных при различных глубинах башмака обсадной колонны на соседних скважинах месторождения. Видно хорошее соответствие спрогнозированных значений порового давления, полученных на основании измерений методом незатухающих гармонических колебаний, данным, измеренным при периодических испытаниях пласта. ФОРМУЛА ИЗОБРЕТЕНИЯ 1. Способ прогнозирования порового давления во время бурения по кусочкам породы, заключающийся в том, чтоa) получают кусочки породы при бурении,b) вычисляют скорость по кусочкам породы, используя метод незатухающих гармонических колебаний,c) вычисляют пористостьпласта,d) используют пористостьдля вычисления вертикального эффективного напряжения v согласно уравнениюv=max(1-),где max и- постоянные параметры, зависящие от литологии,e) вычисляют суммарное напряжение Sv покрывающего пласта путем интегрирования плотности по данным каротажаSv=(z)gdz,где z - глубина; (z) - объемная плотность породы пласта как функция глубины; d - плотность; g - постоянная силы тяжести,f) прогнозируют поровое давление согласно уравнению Рр=(Sv-v)/,где Sv - суммарное напряжение в покрывающем пласте; Рр - поровое давление;- коэффициент Байота,равный 1-(K/K), K - объемная сжимаемость, K - сжимаемость зерен. 2. Способ по п.1, отличающийся тем, что осуществляют построение диаграммы зависимости спрогнозированного порового давления от глубины. 3. Способ по п.1, отличающийся тем, что используют метод незатухающих гармонических колебаний, который реализуют при использовании прибора, содержащего портативное устройство, состоящее из держателя образца с двумя пьезоэлектрическими ультразвуковыми преобразователями, обращенными поверхностями друг к другу, причем один преобразователь используют для возбуждения акустических волн, а второй преобразователь используют для обнаружения. 4. Способ по п.3, отличающийся тем, что дополнительноa) подготавливают кусочки породы для помещения в держатель образца с двумя противоположными параллельными поверхностями, разделенными расстоянием L, где L - длина образца кусочка породы,-7 005450b) размещают образец между преобразователями Р- и S-волн,c) непрерывно подают синусоидальный сигнал в пределах заданного частотного диапазона,d) извлекают и регистрируют огибающую принятого амплитудно-модулированного сигнала,e) измеряют ширину f полосы и осуществляют подстановку этого значения в уравнение v=2Lf для получения результата измерения скорости для каждого максимума амплитуды, где f = fn+1-fn представляет собой разность частот между двумя последовательными резонансными значениями. 5. Способ по п.2, отличающийся тем, что дополнительно осуществляют сравнение указанного спрогнозированного порового давления с прогнозами на основании скоростей сейсмических волн для калибровки данных до получения данных акустического каротажа. 6. Способ определения прогнозируемого порового давления при бурении по субдюймовым кусочкам породы, заключающийся в том, чтоa) получают результат измерения скорости кусочка породы путем подготовки образца с двумя противоположными параллельными поверхностями, отделенными расстоянием L (L - длина образца),b) помещают образец между преобразователями Р- и S-волн,c) непрерывно подают синусоидальный сигнал в пределах заданного частотного диапазона,d) определяют и регистрируют огибающую принятого амплитудно-модулированного сигнала, иe) измеряют f для каждого максимума амплитуды, где f = fn+1-fn представляет собой разность частот между двумя последовательными резонансными значениями, и путем подстановки f в уравнениеVp=2Lf вычисляют скорость Vp для каждого максимума амплитуды, где L - расстояние между параллельными поверхностями или длина кусочка породы,f) преобразуют значения Vp скорости в прогнозируемое поровое давление, для чего вычисляют пористость пласта,используют пористость и вычисляют вертикальное эффективное напряжение v согласно уравнениюv=max(1-),где max ипредставляют собой постоянные параметры, зависящие от литологии,g) вычисляют суммарное напряжение Sv покрывающего пласта путем интегрирования плотности по данным каротажаSv=(z)gdz,где z - глубина, (z) - объемная плотность породы пласта как функция глубины, d - плотность, g - постоянная силы тяжести,h) прогнозируют поровое давление согласно уравнениюPp=(Sv-v)/,где Sv - суммарное напряжение в покрывающем пласте, Рр - поровое давление,- коэффициент Байота,равный 1-(K/K), K - объемная сжимаемость, K - сжимаемость зерен,i) осуществляют построение кривой порового давления в зависимости от глубины. 7. Способ прогнозирования градиента давления гидравлического разрыва пласта по кусочкам породы при бурении, заключающийся в том, чтоa) получают результат измерения скорости Vp в кусочке породы, для чегоb) изготавливают образец с двумя противоположными параллельными поверхностями, разделенными расстоянием L (L - длина образца),c) помещают образец между преобразователями Р- и S-волн,d) непрерывно подают синусоидальный сигнал в пределах заданного частотного диапазона,e) определяют и регистрируют огибающую принятого амплитудно-модулированного сигнала,f) измеряют f - разность частот между двумя последовательными резонансными значениями для каждого максимума амплитуды и осуществляют подстановку f в уравнение Vp=2Lf,g) используют полученный результат измерения скорости для определения пористости по скорости с помощью акустического метода пластового коэффициента при использовании уравнения где =0,625 и Vp в матрице=18000 футов/с для песка,h) прогнозируют поровое давление на основании скорости, полученной методом незатухающих гармонических колебаний, для чегоi) вычисляют суммарное вертикальное напряжение (напряжение в покрывающем пласте) путем интегрирования плотностиSv=(z)gdz,где z - глубина; (z) - объемная плотность породы пласта как функция глубины; d - плотность; g - постоянная силы тяжести,j) получают эффективное вертикальное напряжение при использовании корреляции Холбрукаv=max(1-),ипредставляют собой постоянные параметры, зависящие от литологии, определенные Хол где max бруком,k) определяют поровое давление согласно закону TerzaghiPp=(Sv-v)/,где Рр - поровое давление; Sv - суммарное напряжение в покрывающем пласте, вычисленное на этапе i; v - эффективное вертикальное напряжение, вычисленное на этапе j;- коэффициент Байота, равный 1-(K/K), K объемная сжимаемость, K - сжимаемость зерен,1) прогнозируют градиент давления гидравлического разрыва пласта с помощью модифицированного способа Холбрука при использовании скорости кусочка породы, полученной методом незатухающих гармонических колебанийSh min=[Рр+v(1-)],где Sh min - минимальное горизонтальное напряжение, которое представляет собой нижний предел градиента давления гидравлического разрыва пласта;- пористость, полученная на этапе g; Рр - поровое давление, полученное на этапе k; v - эффективное напряжение, вычисленное на этапе j;- показатель локальной прочности (для района разведочных работ используют 0,85 по умолчанию). 8. Способ по п.7, отличающийся тем, что дополнительно осуществляют построение диаграммы зависимости спрогнозированного градиента давления гидравлического разрыва пласта от глубины. 9. Способ по п.7, отличающийся тем, что осуществляют частотные и амплитудные измерения кусочков породы с использованием прибора незатухающих гармонических колебаний, содержащего портативное устройство, состоящее из держателя образца с двумя пьезоэлектрическими ультразвуковыми преобразователями, расположенными поверхностями друг к другу, причем один преобразователь используют для возбуждения акустических волн, а второй преобразователь используют для обнаружения. 10. Способ прогнозирования в реальном времени градиента давления гидравлического разрыва пласта при бурении скважины при использовании субдюймовых кусочков породы, заключающийся в том,чтоa) получают результат измерения скорости Vp в кусочке породы, для чегоb) изготавливают образец с двумя противоположными параллельными поверхностями, разделенными расстоянием L,c) помещают образец между преобразователями Р- и S-волн,d) непрерывно подают синусоидальный сигнал в пределах заданного частотного диапазона,e) определяют и регистрируют огибающую принятого амплитудно-модулированного сигнала,f) измеряют f - разность частот между двумя последовательными резонансами для каждого максимума амплитуды и осуществляют подстановку f в уравнение Vp=2Lf,g) используют полученный результат измерения скорости для определения пористости по скорости с помощью акустического метода пластового коэффициента согласно уравнению где =0,625 и Vp в матрице=18000 футов/с для песка,h) прогнозируют поровое давление на основании скорости, полученной методом незатухающих гармонических колебаний, для чегоi) вычисляют суммарное вертикальное напряжение (напряжение в покрывающем пласте) путем интегрирования плотностиSv=(z)gdz,где z - глубина; (z) - объемная плотность породы пласта как функция глубины; d - плотность; g - постоянная силы тяжести,j) получают эффективное вертикальное напряжение при использовании корреляции Холбрукаv=max(1-),где max ипредставляют собой постоянные параметры, зависящие от литологии, определенные Холбруком,k) определяют поровое давление согласно закону Terzaghi Рр=(Sv-v)/,где Рp - поровое давление; Sv - суммарное напряжение в покрывающем пласте, вычисленное на этапе i; vэффективное вертикальное напряжение, вычисленное на этапе j;- коэффициент Байота, равный 1-(K/K), K объемная сжимаемость; K - сжимаемость зерен,l) прогнозируют градиент давления гидравлического разрыва пласта с помощью модифицированного способа Cam Clay при использовании диэлектрических измерений кусочков породы (ДИКП), для чегоm) измеряют диэлектрическую проницаемость кусочка породы,-9 005450n) вычисляют площадь поверхности (ПП) на основании диэлектрических измерений кусочков породы согласно уравнению ПП=0,0274(ДИКП)1,9467; о) определяют безразмерный коэффициент М, характеризующий механическую прочность пласта,при использовании площади поверхности, полученной на основании диэлектрического измерения кусочка породы на этапе n,где А=1,54-1,36 Е-3(площадь поверхности),р) получают значения для минимального горизонтального напряжения пласта согласно уравнению 11. Способ по любому из пп.1, 6, 7, 10, отличающийся тем, что осуществляют сбор кусочков породы на интервалах, при которых обеспечивается точность от 5 до 50 футов (от 1,524 до 15,24 м). 12. Способ по п.11, отличающийся тем, что кусочки породы собирают на интервалах, где обеспечивается точность 30 футов (9,144 м). 13. Способ по любому из пп.4, 6, 7, 10, отличающийся тем, что длина L образца находится в пределах от около 1,5 до 4 мм. 14. Способ по любому из пп.1, 7, 10, отличающийся тем, что кусочки породы насыщены. 15. Способ по любому из пп.1, 6, 7, 10, отличающийся тем, что дополнительно выполняют вычисления согласно компьютерной программы.