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

СПОСОБ И СИСТЕМА ОПРЕДЕЛЕНИЯ ХАРАКТЕРИСТИК ГИДРОРАЗРЫВА С ИСПОЛЬЗОВАНИЕМ ДАННЫХ О МИКРОСЕЙСМИЧЕСКИХ СОБЫТИЯХ В изобретении представлены системы и способы определения характеристик гидроразрыва с использованием данных о микросейсмических событиях для выявления ориентаций подземных трещин, расстояний между подземными трещинами и углов падения плоскостей подземных трещин.(71)(73) Заявитель и патентовладелец: ЛЭНДМАРК ГРАФИКС КОРПОРЕЙШН (US) Перекрестные ссылки на родственные заявки Не применимо. Указания на исследования, спонсируемые из федерального бюджета Не применимо. Область техники, к которой относится изобретение Настоящее изобретение относится к системам и способам определения характеристик гидроразрыва с использованием данных о микросейсмических событиях. В частности, настоящее изобретение относится к определению характеристик гидроразрыва с использованием данных о микросейсмических событиях для выявления ориентаций подземных трещин, расстояний между подземными трещинами и углов падения плоскостей подземных трещин. Уровень техники Данные о микросейсмических событиях иногда получают при проведении мероприятий по осуществлению гидроразрыва пласта через скважину, пробуренную для освоения углеводородных месторождений. Данные получают при помощи низкомагнитудных (микро)сейсмических событий, которые связаны с трещинами, образующимися при гидроразрыве. Суть гидроразрыва подземного пласта состоит в создании в этом пласте трещин путем нагнетания в структуру горной породы флюидов и проппантов под высоким давлением с целью интенсификации добычи из этого пласта углеводородов. За местоположением точек, в которых происходят микросейсмические события, обычно ведется наблюдение в реальном времени, причем местоположения этих точек, отображаемые в трехмерном (3D) виде, могут проверяться по мере появления этих точек. Местоположения точек, в которых происходят микросейсмические события, также доступны для анализа по завершении операций гидроразрыва и, таким образом, доступны для сравнения с результатами, полученными в соседних скважинах. Микросейсмические события обычно происходят вдоль подземных трещин или рядом с подземными трещинами,которые могут представлять собой искусственно созданные разрывы или разрывы, образовавшиеся ранее естественным путем и повторно раскрытые при осуществлении гидроразрывных операций. На ориентацию трещин большое влияние оказывает режим напряжений, существующий на текущий момент, а также наличие систем трещин, появившихся в различные периоды времени в прошлом, когда ориентация напряжений отличалась от текущей. Каждое регистрируемое и анализируемое отдельно взятое микросейсмическое событие, отличное от других микросейсмических событий, является результатом разрыва, осуществленного через скважину и имеющего ориентацию, магнитуду, местоположение и другие характеристики, которые могут быть определены при помощи наклономера или на основе сейсмических данных, полученных от датчиков. Разрыв может характеризоваться и другими параметрами, например длиной, шириной, высотой и давлением. В отношении местоположения каждой точки, в которой происходит микросейсмическое событие, существует неопределенность. Неопределенность в направлении плоскости координатных осей XY отличается от неопределенности в вертикальном направлении глубинной области (координатной оси Z). Неопределенность местоположения каждой точки, в которой происходит событие, может быть представлена в виде удлиненного эллипсоида вращения. В некоторых случаях ориентация микросейсмических событий и расстояния между микросейсмическими событиями известны, что соответствует классическим моделям представления разрыва в виде трещины с двумя крыльями, часто используемым при математических описаниях анализа разрывов. В других случаях плотное облако данных, представляющее объемную область, охватывающую все точки микросейсмических данных, указывает на сложный характер распределения искусственно созданных трещин или повторно раскрытых естественных трещин. В этих случаях анализ микросейсмических данных становится в большой степени субъективным и допускает разные интерпретации, но даже в этих случаях в облаке данных выделяют закономерности, которые могут указывать на форму трещин, имеющихся в подземном пласте. Область напряжений, существующая в текущий момент времени, может отличаться от области напряжений, имевшейся на исходный момент образования разрывов. На ориентацию трещин, искусственно создаваемых в текущий момент времени при осуществлении гидроразрывов, большое влияние оказывает состояние напряжений в подземном пласте. При рассмотрении напряжений в вертикальном и двух горизонтальных направлениях в подземном пласте выделяют наличие некоторой анизотропии напряжений,которая имеется всегда. Чем выше эта анизотропия, тем более плоской структурой характеризуются трещины, вызванные гидроразрывом, и тем больше они соответствуют традиционной модели трещины с двумя крыльями. Чем больше проницаемость породы, тем более плоскую структуру будут иметь трещины. Чем более изотропный характер носит режим напряжений, тем интенсивнее фронты распространения трещин могут отражаться от границ нарушения сплошности горной породы и создавать сложные системы трещин. На сегодняшний день существует несколько методов определения характеристик трещин, используемых для апробирования и выявления ориентаций искусственно созданных или образованных естественным путем подземных трещин, углов падения плоскостей таких трещин и расстояний между такими трещинами. Один из этих методов заключается в определении облака данных, охватывающего все точки микросейсмических данных, для построения моделируемой объемной области коллектора (от англ. simulatedreservoir volume - SRV) или расчетной моделируемой объемной области коллектора (от англ. estimatedreservoir volume - ERV). На основе данной информации оценивается объем горной породы, в который нагнетались флюиды и проппанты. Следует отметить, что лишь малое количество энергии, сообщенной грунту, доходит до поверхности и регистрируется как микросейсмические события. Фактически, энергетические характеристики множества микросейсмических событий лежат за пределами чувствительности контрольно-измерительной аппаратуры. Таким образом, при применении этих методов определения характеристик трещин использование геометрических параметров облака данных дает лишь приблизительные результаты, которые могут оказаться недействительными. Интерпретация характера распределения трещин на основе системы данных о микросейсмических событиях относится к другим методам, описанным в публикации "Optimizing Horizontal Completion Technologies in the Barnett Shale Using MicroseismicFracture Mapping" (Оптимизация технологии заканчивания горизонтальных скважин в формации глинистого сланца "Барнетт" с использованием карты распределения микросейсмических данных о трещиноватости) авторами М.K. Fisher и др. Еще один метод определения характеристик трещин предполагает получение и анализ данных наземной сейсморазведки (активных и пассивных). В некоторых случаях может иметь место поляризация поперечных волн, связанная с наличием в подземном пласте трещин. Скорость распространения поперечных волн в одних направлениях может быть выше, нежели в других направлениях. В методах определения характеристик трещин также используются обменные волны (например, преобразующиеся из продольных в поперечные). Кроме того, существует несколько различных методов определения характеристик трещин, согласно которым данные о микросейсмических событиях математически связываются с моделью подземного пласта, и создается дискретная система трещин (от англ. discrete fracture network - DFN) или совокупность предполагаемых характеристик трещин. К таким методам относятся, например, методы, описанные в американских патентных заявках 2010/0307755 и 2011/0029291. Еще один метод автоматического определения характеристик трещин, позволяющий вести численный анализ облака данных и выявлять главные плоскости по точкам микросейсмических данных для моделирования коллектора углеводородов, был предложен в публикации "Quantitative Interpretation of MajorPlanes From Microseismic Event Locations with Application in Production Prediction" (Количественная интерпретация главных плоскостей ориентаций на основе местоположения точек, в которых происходят микросейсмические события, применительно к прогнозированию дебита) авторами М.J. Williams и др. Данный метод представляет подход к решению задачи определения направлений трещин по вероятностным структурам, но позволяет определить ориентацию только одной трещины. Краткое описание изобретения Предложенные настоящим изобретением системы и способы определения характеристик гидроразрыва с использованием данных о микросейсмических событиях для выявления ориентаций подземных трещин, расстояний между подземными трещинами и углов падения плоскостей подземных трещин удовлетворяют вышеописанные потребности и позволяют устранить один или несколько недостатков известного уровня техники. В одном варианте осуществления настоящего изобретения предложен способ определения характеристик гидроразрыва с использованием данных о микросейсмических событиях, содержащий: i) выявление при помощи процессора компьютера одного или более углов простирания по усредненным точкам микросейсмических данных, относящихся к одной или нескольким совокупностям плоскостей, и неопределенности их ориентаций, причем указанные один или более углов простирания имеют большее количество усредненных точек микросейсмических данных при соответствующем угле простирания, чем количество усредненных точек микросейсмических данных при других соответствующих углах простирания; ii) определение угла падения каждой совокупности плоскостей с учетом весового коэффициента контакта между каждой точкой микросейсмических данных и каждой плоскостью в каждой совокупности плоскостей; iii) определение местоположения каждой плоскости в каждой совокупности плоскостей; и iv) выполнение статистического анализа для каждой совокупности плоскостей с общим выявленным углом простирания для выявления неопределенности, связанной с ориентацией соответствующей совокупности подземных трещин, расстоянием между подземными трещинами соответствующей совокупности подземных трещин и углом падения плоскостей соответствующей совокупности подземных трещин. В другом варианте осуществления настоящего изобретения предложен носитель программного обеспечения для длительного хранения данных, физически содержащий выполняемые компьютером команды для определения характеристик гидроразрыва по точкам микросейсмических данных, причем данные команды при их выполнении обеспечивают следующее:i) выявление одного или более углов простирания по усредненным точкам микросейсмических данных, относящихся к одной или нескольким совокупностям плоскостей, и неопределенности их ориентаций, причем указанные один или более углов простирания имеют большее количество усредненных точек микросейсмических данных при соответствующем угле простирания, чем количество усредненных точек микросейсмических данных при других соответствующих углах простирания;ii) определение угла падения каждой совокупности плоскостей с учетом весового коэффициента контакта между каждой точкой микросейсмических данных и каждой плоскостью в каждой совокупности плоскостей;iii) определение местоположения каждой плоскости в каждой совокупности плоскостей; иiv) выполнение статистического анализа для каждой совокупности плоскостей с общим выявленным углом простирания для выявления неопределенности, связанной с ориентацией соответствующей совокупности подземных трещин, расстоянием между подземными трещинами соответствующей совокупности подземных трещин и углом падения плоскостей соответствующей совокупности подземных трещин. Дополнительные аспекты, преимущества и варианты осуществления настоящего изобретения будут понятны специалисту после ознакомления с нижеприведенным описанием различных вариантов осуществления настоящего изобретения и соответствующими чертежами. Краткое описание чертежей Ниже приведено описание настоящего изобретения со ссылками на прилагаемые чертежи, на которых сходные элементы имеют сходные позиционные обозначения. На фиг. 1 А показана блок-схема алгоритма, иллюстрирующая вариант осуществления способа,предложенного настоящим изобретением. На фиг. 1 В показано продолжение блок-схемы алгоритма способа, иллюстрируемого на фиг. 1 А. На фиг. 2 А проиллюстрирован один формат отображения усредненных точек микросейсмических данных при каждой ориентаций, причем указанное отображение осуществляется на этапе 108 блоксхемы алгоритма, приведенной на фиг. 1 А. На фиг. 2 В проиллюстрирован другой формат отображения усредненных точек микросейсмических данных при каждой ориентаций, причем указанное отображение осуществляется на этапе 108 блоксхемы алгоритма, приведенной на фиг. 1 А. На фиг. 2 С проиллюстрирован еще один формат отображения усредненных точек микросейсмических данных при каждой ориентаций, причем указанное отображение осуществляется на этапе 108 блоксхемы алгоритма, приведенной на фиг. 1 А. На фиг. 3 А проиллюстрирован формат отображения данных, приведенный на фиг. 2 А, после применения стандартного статистического анализа по методу Монте-Карло. На фиг. 3 В проиллюстрирован формат отображения данных, приведенный на фиг. 2 В, после применения стандартного статистического анализа по методу Монте-Карло. На фиг. 3 С проиллюстрирован формат отображения данных, приведенный на фиг. 2 С, после применения стандартного статистического анализа по методу Монте-Карло. На фиг. 4 показан пример двухмерной карты оптимальных ориентаций плоскостей и расстояний между некоторыми вертикальными плоскостями падения, представляющими совокупность подземных трещин, обрабатываемую на этапе 113 блок-схемы алгоритма, приведенной на фиг. 1 В. На фиг. 5 показана таблица, в которой приведены статистические данные для трех семейств плоскостей, показанных на фиг. 4. На фиг. 6 показана структурная схема, иллюстрирующая вариант осуществления системы, предложенной настоящим изобретением. Подробное описание предпочтительных вариантов осуществления изобретения. Объект настоящего изобретения описан на основе конкретных вариантов осуществления изобретения, однако данное описание никоим образом не ограничивает объем настоящего изобретения. Таким образом, объект настоящего изобретения может быть реализован по-иному и может включать в себя различные этапы или различные сочетания этапов, схожих с описанными в данном документе, вкупе с другими технологиями. Кроме того, слово "этап" может использоваться в данном документе для описания моментов и периодов времени, относящихся к выполнению различных процессов, происходящих в рамках предложенных способов, однако данное слово не подразумевает какой-либо определенной последовательности действий, если это прямым текстом не указано в отношении какого-то конкретного порядка действий. Нижеприведенное описание изобретения относится к нефтегазовой отрасли промышленности,однако системы и способы, предложенные настоящим изобретением, не ограничиваются применением только в этой отрасли промышленности, и могут в равной степени применяться в других отраслях промышленности. Отдельные микросейсмические события, вызываемые гидроразрывом, характеризуются данными о микросейсмических событиях, которые могут использоваться для определения ориентаций созданных при гидроразрыве и образованных естественным путем трещин, имеющихся в подземном пласте, углов падения плоскостей таких трещин и расстояний между такими трещинами. В объемной области задают и вводят плоскости. Для задания плоскости требуется по меньшей мере три связанные точки микросейсмических событий. Может быть рассчитано пересечение каждой плоскости и соответствующего удлиненного эллипсоида вращения, выражающего неопределенность местоположения в направлениях координатных осей X, Y и Z, выходящих из каждой найденной центральной точки данных о микросейсмическом событии. Плоскости могут систематически поворачиваться на основе входных данных так, чтобы по ориентации данных о микросейсмических событиях, падению плоскости таких данных и расстоянию между такими данными можно было определить наиболее вероятную совокупность плоскостей, соответствующих входным данным. С любой из ориентаций может быть связана любая трещина, которая может располагаться на пересечении двух и более семейств трещин. Вероятность того, что при помощи данных о микросейсмических событиях была определена истинная плоскость трещины в подземном слое, возрастает с увеличением количества событий, лежащих на плоскости. Количество изначально имеющихся систем трещин не оценивается, так как множество геологических событий могло привести к появлению нескольких групп перекрывающихся систем трещин в любой взятой области. Ниже раскрыт способ, позволяющий определить наиболее вероятную ориентацию каждой отдельной трещины и семейства трещин, расстояния между трещинами и семействами трещин и падение плоскости отдельной трещины и семейства трещин, причем данный способ также может использоваться для определения любых корреляций с данными о микросейсмических событиях. Описание способа. На фиг. 1 показана блок-схема, иллюстрирующая один вариант осуществления способа 100, предложенного настоящим изобретением. На этапе 101 осуществляют ввод данных о микросейсмических событиях при помощи пользовательского интерфейса/видеоинтерфейса, описанного со ссылкой на фиг. 6. Данные о микросейсмических событиях, вызванных естественными геологическими факторами или гидроразрывом, поступают от датчиков. Данные о микросейсмических событиях, таким образом, могут содержать несколько массивов данных, которые могут представлять отдельные события гидроразрыва, относящиеся к одной или нескольким скважинам. Данные о микросейсмических событиях включают в себя точки микросейсмических данных и информацию, относящуюся к местоположению каждой точки сейсмических данных в координатах X, Y и Z и к неопределенности местоположения, связанной с каждой из координат X, Y и Z. Данные о микросейсмических событиях также могут содержать моментную магнитуду и тензор момента,связанные с каждой точкой микросейсмических данных. Касательно неопределенности местоположения для каждой точки микросейсмических данных, неопределенность местоположения может содержать значения неопределенности для каждой точки микросейсмических данных по каждой из координат X, Y и Z в отдельности или оценочное значение этой неопределенности, которое может применяться ко всем точкам микросейсмических данных. Дополнительно данные о микросейсмических событиях могут содержать другую информацию, относящуюся к местоположению датчиков, их свойствам и скоростным характеристикам. Данные о микросейсмических событиях могут поступать от любых доступных источников данных, в том числе, но не ограничиваясь таковым, от наклономеров, наземных датчиков или внутрискважинных датчиков, установленных в одной или нескольких наблюдательных скважинах. Кроме того, данные о микросейсмических событиях могут дополняться информацией, полученной из скважины,в том числе, но не ограничиваясь таковым, из кернов, дипольных зондов, приборов акустического каротажа или приборов фото- и видеокаротажа. По данным ориентированного дипольного акустического каротажа, например, можно определить разницу между напряжением в вертикальном направлении и максимальным и минимальным напряжениями в горизонтальных направлениях. Для осуществления этапов способа 100 могут использоваться и другие данные. На этапе 102 при помощи известных методов и на основе данных о микросейсмических событиях,введенных на этапе 101, выявляют местоположение всех точек микросейсмических данных и их соответствующую неопределенность местоположения. На этапе 103 задают ориентацию (например, угол простирания и угол падения) плоскости. Например, угол простирания и угол падения могут быть установлены в (исходные) значения 0 и 0, или 0 и 360, или 0 и 90, или 360 и 90 соответственно. На этапе 104 создают плоскость, имеющую вектор нормали (sincos , sinsin , cos ) и проходящую через среднее местоположение всех точек микросейсмических данных для ориентации плоскости(угла простирания и угла падения). Кроме того, известными методами рассчитывается кратчайшее расстояние от каждой точки микросейсмических данных до плоскости. Угол псивыражает угол простирания, угол тетавыражает угол падения. Кратчайшее расстояние между каждой точкой микросейсмических данных и плоскостью может быть рассчитано по формуле Таким образом, кратчайшее расстояние между каждой точкой микросейсмических данных и плоскостью имеет положительный знак, если точка находится с одной стороны плоскости, и отрицательный знак, если точка находится с другой стороны плоскости. Каждая точка микросейсмических данных характеризуется кратчайшим расстоянием, имеющим знак . Изменяющееся количество точек микросейсмических данных, таким образом, будет связано с вероятными плоскостями, каждая из которых отстоит от этих точек на расстояния, имеющие знак. Такие плоскости будут объединены в совокупности и могут быть связаны с возможными плоскостями трещин. На этапе 105 при помощи известных методов создают совокупности плоскостей, параллельных плоскости, созданной на этапе 104, причем каждая новая плоскость пересекает объемную область, пред-4 025073 ставленную неопределенностью местоположения каждой точки микросейсмических данных. На этапе 106 подсчитывают количество плоскостей, созданных на этапе 105, и количество точек микросейсмических данных, контактирующих с этими плоскостями. Для ориентации плоскостей также при помощи известных методов рассчитываются усредненные точки микросейсмических данных, контактирующие с каждой плоскостью. На этапе 107 угол простирания или угол падения для плоскости, созданной на этапе 104, увеличивают на 1, после чего повторяют этапы 104-106, пока угол простирания не станет равным 0 или 360, а угол падения не станет равным 0 или 90. После повторения этапов 104-106, пока угол простирания не станет равным 0 или 360, а угол падения не станет равным 0 или 90, переходят к этапу 108 способа 100. На этапе 108 отображают усредненные точки микросейсмических данных, рассчитанные при каждой ориентации на этапе 106, для определения известными методами главных ориентаций простирания усредненных точек микросейсмических данных, относящихся к совокупности плоскостей и их ориентаций. Кроме того, отображают неопределенность ориентаций для анализа. Таким образом, тщательно обрабатывают, анализируют и отображают результаты, полученные на этапах 104-107 для различных сочетаний углов простирания и углов падения. На данном этапе данные могут отображаться в любом из множества предпочтительных форматов, однако способ 100 не ограничивается каким-либо определенным форматом отображения этих данных, и на примерах, приведенных на фиг. 2 А-2 С и 3 А-3 С, проиллюстрированы различные форматы, в которых могут отображаться указанные данные. На фиг. 2 А проиллюстрирован один формат осуществляемого на этапе 108 отображения усредненных точек микросейсмических данных при каждой ориентации. На фиг. 2 А показан пример графика простирания и падения, соответствующих усредненным точкам микросейсмических данных при каждой ориентаций. Точки микросейсмических данных в координатах XY показаны в виде затемнений серого цвета. Найдены три семейства плоскостей с главными углами простирания 112, 132 и 178. Кроме того,семейство плоскостей трещин с углом простирания 112 имеет предпочтительный угол падения приблизительно 84, в то время как семейства плоскостей трещин с углами простирания 132 и 178 в большей степени отклонены от вертикали и имеют предпочтительные углы падения, в среднем 74. На фиг. 2 В проиллюстрирован другой формат отображения усредненных точек микросейсмических данных при каждой ориентаций. На фиг. 2 В показана гистограмма усредненных точек микросейсмических данных при каждой ориентаций в зависимости от угла простирания. На фиг. 2 С проиллюстрирован еще один формат отображения усредненных точек микросейсмических данных при каждой ориентаций. На фиг. 2 В показан график радиального распределения усредненных точек микросейсмических данных при каждой ориентаций. На фиг. 3 А показаны результаты анализа данных, отображаемых в приведенном на фиг. 2 А формате, на примере применения стандартного статистического анализа по методу Монте-Карло. Другими словами, отображается неопределенность местоположения или главный угол простирания по усредненным точкам микросейсмических данных, относящимся к совокупности плоскостей и их неопределенности ориентаций. Таким образом, при помощи статистического анализа по методу Монте-Карло может быть оценена достоверность определения главных углов простирания в усредненных точках микросейсмических данных, относящихся к совокупности плоскостей и их неопределенности ориентации. Неопределенность каждого микросейсмического события в направлениях координатных осей X, Y, Z подчиняется распределению Гаусса. Случайный процесс накладывается на координату каждого микросейсмического события по координатной оси X в виде функции х+ xN(0;1), где х - радиус неопределенности в направлении координатной оси X, a N(0;1) - нормальное распределение. Аналогичный процесс накладывается на координаты микросейсмического события в направлениях координатных осей Y и Z. После наложения этих процессов на все события этапы 104-107 повторяются для осуществления сбора усредненных точек микросейсмических данных и отображения их в пространстве значений углов простирания и падения. После семплирования множества создаваемых величин на карте остаются синфазные компоненты, а противофазные компоненты взаимно исключаются. Таким образом, главные углы простирания могут быть определены с большей четкостью. На фиг. 3 В показаны результаты анализа данных, отображаемых в приведенном на фиг. 2 В формате, на примере применения стандартного статистического анализа по методу Монте-Карло. На фиг. 3 С показаны результаты анализа данных, отображаемых в приведенном на фиг. 2 С формате, на примере применения стандартного статистического анализа по методу Монте-Карло. На этапе 109 повторяют этапы 104-106 при каждом главном угле простирания для каждого соответствующего угла падения на основе микросейсмических данных, введенных на этапе 101. При этом с учетом неопределенности местоположения и магнитуды каждого микросейсмического события может быть выполнен подробный анализ каждой плоскости при главном угле простирания. На этапе 110 в каждой совокупности созданных плоскостей находят плоскость, имеющую максимальное количество пересечений с точками микросейсмических данных, и определяют ее местоположение при помощи известных методов. Эта плоскость представляет собой плоскость с наилучшим соответствием относительно точек микросейсмических данных, введенных на этапе 101. На этапе 111 рассчитывают весовой коэффициент контакта между плоскостью, найденной на этапе 110, и пересекающимися точками микросейсмических событий, причем данный расчет основывается на следующих параметрах: i) неопределенность местоположения каждой точки микросейсмических данных в направлениях координатных осей X, Y, Z; ii) кратчайшее расстояние от центра каждой точки микросейсмических данных до плоскости; iii) числовое взвешенное значение, выражающее моментную магнитуду каждого отдельного микросейсмического события; iv) корреляция между точками микросейсмических данных и плоскостью. Каждая точка микросейсмических данных выражает объемную область с наибольшей степенью определенности местоположения в центре и наименьшей степенью определенности на периферии. Точки микросейсмических данных с большей моментной магнитудой имеют большую степень определенности местоположения, нежели точки микросейсмических данных с меньшей моментной магнитудой, регистрация которых может быть затруднительна. Другими словами, низкая моментная магнитуда свидетельствует о том, что данные точки микросейсмических данных характеризуются наибольшей неопределенностью их местоположения. Таким образом, высокая моментная магнитуда для конкретной точки микросейсмических данных, измеренная в ее центре, может лежать в диапазоне 0,51,0. Для всех точек, расположенных на периферической части сжатого у полюсов сфероида, ограниченного неопределенностями в направлениях координатных осей X, Y, Z, присваивается весовой коэффициент контакта со значением 0,05. Значение весового коэффициента линейно убывает от центра местоположения к периферии и выражает неопределенность местоположения. Численное значение весового коэффициента контакта, применяемого к каждой точке сейсмических данных, является произведением следующих элементов: i) значения, соответствующего степени достоверности определения каждого микросейсмического события, согласно представляемому масштабированию моментной магнитуды; ii) неопределенности расположения каждого микросейсмического события в направлении координатных осей X, Y, Z; iii) кратчайшего расстояния от центра точки микросейсмических данных до пересекающей плоскости. Статистическая корреляция между точками микросейсмических данных и плоскостью далее может быть проанализирована на основе обратной величины расстояния между плоскостью и точками микросейсмических данных. Если точки микросейсмических данных лежат в одной плоскости трещины, сумма этих расстояний должна быть минимальной. На этапе 112 точки микросейсмических данных, контактирующие с плоскостью, найденной на этапе 110, удаляют, и повторяют этапы 109-112, пока не останется ни одной определенной на этапе 110 плоскости, контактирующей с тремя или большим количеством точек микросейсмических данных. После определения каждой из плоскостей, контактирующих с тремя или большим количеством точек микросейсмических данных, переходят к этапу 113 способа 100. На этапе 113 выполняют статистический анализ для каждой совокупности плоскостей (выражающей совокупность подземных трещин) с общим главным углом простирания с учетом весовых коэффициентов контакта, рассчитанных на этапе 111, для нахождения неопределенности, относящейся к ориентаций каждой совокупности подземных трещин, расстоянию между подземными трещинами в каждой совокупности подземных трещин и падению плоскости каждой совокупности подземных трещин. Анализ может содержать, например, основные статистические данные (среднее значение, срединное значение, среднеквадратическое отклонение простирания, падения и длины), рассчитываемые для каждой совокупности плоскостей, имеющих общий главный угол простирания. На этапе 114 отображают результаты, полученные на этапе 113. Неопределенность, относящаяся к ориентации всех потенциальных подземных трещин, расстоянию между потенциальными подземными трещинами и падению плоскостей всех потенциальных подземных трещин и связанная с микросейсмическими событиями, которыми сопровождается гидроразрыв, может быть отображена и представлена для анализа специалистом, принимающим решение по принятию или отклонению результатов анализа на основе геологической истории области, в которой могут образоваться трещины. Могут иметь место трещины, образованные в ранние геологические периоды, когда режим напряжений отличался от текущего режима напряжений. Таким образом, специалист-аналитик может учитывать сведения о геологической истории рассматриваемой области при определении всех возможных сочетаний простирания, падения и расстояний между трещинами, которые могут не быть признаны действительными ввиду того, что данные о микросейсмических событиях не дают полного описания трещин гидроразрыва, в частности потому, что маломощные микросейсмические события могут не улавливаться приборами с недостаточно высокой чувствительностью для регистрации таких микросейсмических событий. Кроме того, совокупности плоскостей, определенные способом 100 и имеющие наилучшее математическое соответствие, могут в геологическом смысле быть неприемлемыми, а совокупности плоскостей с другой степенью математического соответствия могут быть более уместными для общего анализа и оценки трещин. На фиг. 4 приведена для примера двухмерная карта оптимальных углов простирания и расстояний для некоторых вертикальных плоскостей падения, выражающих совокупность подземных трещин, анализируемых на этапе 113. Массив данных, представленных на карте, охватывает 284 микросейсмических события. Для иллюстрации результатов анализа, проводимого на этапах 109-115, были выбраны три главных угла простирания (290, 312, 358). Для большей наглядности дополнительные ориентации не были выбраны. Всего найдена 31 плоскость. Семь плоскостей составляют семейство с углом простира-6 025073 ния 290, 10 плоскостей составляют семейство с углом простирания 312 и 14 плоскостей составляют семейство с углом простирания 358. Эти плоскости содержат 221 микросейсмическое событие, что составляет 78% от общего количества микросейсмических событий, зарегистрированных в этом массиве данных. Для визуализации этих плоскостей все микросейсмические события (обозначенные точками серого цвета) спроецированы на плоскость с осями "север-восток" (плоскость с координатными осями XY), данное изображение является видом сверху. На рисунке показаны только 12 найденных плоскостей(обозначены отрезками серого цвета) с углом падения более 85: семейство из трех содержащих микросейсмические события плоскостей с углом простирания 290, семейство из четырех содержащих микросейсмические события плоскостей с углом простирания 312 и семейство из пяти содержащих микросейсмические события плоскостей с углом простирания 358. В таблице, приведенной на фиг. 5, показаны статистические данные для этих трех семейств плоскостей (показанных на фиг. 4). В данной таблице для каждого семейства плоскостей указаны среднее значение и среднеквадратическое отклонение ориентаций (угла простирания), падения и расстояния между трещинами. Таким образом, способ 100 обеспечивает более точное выявление различных семейств направлений трещин. В области, в отношении которой предпринимаются действия по интенсификации добычи углеводородов, могут содержаться множества семейств трещин, каждое из которых может иметь свою ориентацию, падение и расстояние между трещинами, что в большей степени характерно для областей, в которых выявлены сложные системы трещин. Таким образом, способ 100 может использоваться для следующих целей: 1) предоставление входных данных для модели сложной системы трещин, при помощи которой моделируются параметры семейств трещин и характер их изменения в одной скважине и на отрезке между двумя скважинами в некоторой области; 2) предоставление входных данных для механической модели геологической среды, содержащей напряжения, параметры горной породы, имеющей слоистую структуру и характеризующейся прочей информацией о структуре, стратиграфических характеристиках и свойствах подземного пласта; 3) предоставление входных данных для моделей эксплуатационных характеристик скважины и коллектора для анализа одной или нескольких скважин на предмет производительности, анализа дренирования и анализа общей величины дебита; 4) анализ в реальном времени i) обнаруженных систем трещин по результатам первой операции из ряда операций по интенсификации добычи углеводородов для выявления возможного закрытия обнаруженных и повторно открытых трещин и возможного пересечения новых направлений трещин с трещинами гидроразрыва в скважине, а также ii) данных о семействе трещин, обнаруженных в других соседних скважинах, для выявления возможного образования трещин с тем же направлением при текущей операции по интенсификации добычи углеводородов; и 5) анализ данных о микросейсмических событиях для определения связи между моментной магнитудой и тензорами момента отдельных микросейсмических событий и выявленными ориентациями семейств плоскостей трещин, расстояний между трещинами в семействах плоскостей трещин и падением семейств плоскостей трещин. Описание системы. Настоящее изобретение может быть реализовано посредством выполняемой компьютером программы, содержащей, например, программные модули, обычно называемые программными приложениями или прикладными программами, выполняемыми компьютером. Программные средства могут содержать, например, подпрограммы, программы, объекты, компоненты и блоки данных, выполняющие определенные задания или реализующие конкретные абстрактные типы данных. Программными средствами формируется интерфейс для взаимодействия компьютера с источниками входных сигналов. Программные средства также могут взаимодействовать с другими кодами для выполнения различных заданий в ответ на получение данных от соответствующих источников. Программные средства могут храниться и/или переноситься на различных запоминающих устройствах, например на компактных дисках,магнитных дисках, запоминающих устройствах на цилиндрических магнитных доменах и запоминающих устройствах на полупроводниковых приборах (например, на постоянном запоминающем устройстве или оперативном запоминающем устройстве). Кроме того, программные средства и результаты их действия могут передаваться при помощи различных средств передачи данных, например по оптическому волокну, металлическим проводам и/или посредством любой компьютерной сети, например через Интернет. Кроме того, специалисту должно быть понятно, что настоящее изобретение может быть реализовано различными компьютерными системами, в том числе переносными устройствами, многопроцессорными системами, средствами микропроцессорной или программируемой пользователем техники, миникомпьютерами, большими универсальными компьютерами и т. п. При использовании настоящего изобретения может применяться любое количество компьютерных систем и компьютерных сетей. Настоящее изобретение может быть реализовано в средах распределенных вычислений, в которых задания выполняются удаленными обрабатывающими устройствами, соединенными между собой в сеть передачи данных. В среде распределенных вычислений программные модули могут находиться как на локальном,-7 025073 так и на удаленном компьютерном запоминающем устройстве. Настоящее изобретение, таким образом,может осуществляться при использовании различных аппаратных средств, программных средств в любом их сочетании посредством как компьютерной системы, так и другой вычислительной системы. На фиг. 6 показана структурная схема, иллюстрирующая один вариант осуществления предложенной системы, предназначенной для реализации настоящего изобретения на компьютере. В систему входит вычислительный блок, также называемый вычислительной системой, который содержит запоминающее устройство, прикладные программы, интерфейс пользователя, видеоинтерфейс и процессор. Вычислительный блок представляет собой лишь один пример подходящей компьютерной среды, при этом область применения или функциональность настоящего изобретения не ограничивается использованием такого вычислительного блока. В запоминающем устройстве хранятся прикладные программы, которые также могут быть охарактеризованы как программные модули, содержащие выполняемые на компьютере команды, которые могут выполняться вычислительным устройством для осуществления настоящего изобретения, описанного в данном документе и проиллюстрированного на фиг. 1-5. Таким образом, запоминающее устройство содержит модуль определения характеристик трещин, посредством которого реализуются способы, описанные и проиллюстрированные со ссылками на фиг. 1 А, 1 В. Результаты, полученные модулем определения характеристик трещин, могут использоваться в модели сложной системы трещин, однако вместо модуля определения характеристик трещин могут применяться другие прикладные интерфейсные средства, работающие автономно. На фиг. 6 показано, что вычислительный блок имеет универсальное запоминающее устройство, однако вычислительный блок обычно содержит несколько носителей считываемой компьютером информации. Например, носители считываемой компьютером информации могут представлять собой (но не ограничиваются таковыми) устройства хранения считываемой компьютером информации. Запоминающее устройство вычислительной системы может представлять собой устройство хранения считываемой компьютером информации в виде энергозависимого и/или энергонезависимого запоминающего устройства,например постоянного запоминающего устройства (ПЗУ) или оперативного запоминающего устройства(ОЗУ). В ПЗУ обычно хранится базовая система ввода-вывода (от англ. BIOS- basic input/output system),содержащая основные команды, обеспечивающие передачу информации между элементами вычислительного блока, например при запуске. В ОЗУ обычно содержатся данные и/или программные модули, к которым имеет произвольный доступ вычислительный блок и/или которые выполняются вычислительным блоком в текущий момент времени. Например, вычислительный блок может содержать (но не ограничивается таковым) операционную систему, прикладные программы, прочие программные модули и программные данные. Компоненты, показанные в запоминающем устройстве, также могут содержаться в других стационарных/сменных, энергозависимых/энергонезависимых устройствах хранения считываемой компьютером информации или могут быть задействованы в вычислительном блоке при помощи интерфейса прикладного программирования (от англ. API -application program interface), который может находиться в отдельном вычислительном блоке, соединенном через компьютерную систему или сеть. Например, считывание данных из стационарного энергонезависимого магнитного носителя информации и запись данных в стационарный энергонезависимый магнитный носитель информации может осуществляться через накопитель на жестких дисках; считывание данных из сменного энергонезависимого магнитного диска и запись данных в сменный энергонезависимый магнитный диск может осуществляться через накопитель на магнитных дисках; считывание данных из сменного энергонезависимого оптического диска, например компакт-диска или другого оптического носителя информации, и запись данных в такой сменный энергонезависимый оптический диск может осуществляться через накопитель на оптических дисках. К другим сменным/стационарным энергозависимым/энергонезависимым устройствам хранения считываемой компьютером информации, которые могут использоваться в приведенной для примера операционной среде, относятся (но не ограничиваются таковыми) кассеты с магнитной лентой, карты флэш-памяти,универсальные цифровые диски, видеокассеты, твердотельные ОЗУ, твердотельные ПЗУ и т. п. Устройства ввода носителей и соответствующие устройства хранения считываемой компьютером информации,указанные выше, обеспечивают хранение считываемых компьютером команд, структур данных, программных модулей и других данных, используемых в вычислительном блоке. Пользователь может вводить команды и информацию в вычислительный блок через интерфейс пользователя, который может представлять собой устройство ввода, например клавиатуру и указательное устройство в виде компьютерной мыши, шарового манипулятора или сенсорной панели. Устройство ввода может представлять собой микрофон, джойстик, спутниковую тарелку, сканер и т. п. Эти и другие устройства ввода обычно соединяются с процессором через системную шину, но также могут соединяться через другие интерфейсные и шинные устройства, например через параллельный порт или универсальную последовательную шину (от англ. USB - universal serial bus). Через специальный интерфейс, например видеоинтерфейс, к системной шине может подключаться монитор или иное устройство отображения информации. Совместно с видеоинтерфейсом также может использоваться графический интерфейс пользователя (от англ. GUI- graphical user interface) для получе-8 025073 ния команд от интерфейса пользователя и для передачи команд в процессор. Помимо монитора в состав компьютера могут входить различные периферийные устройства вывода, например акустические колонки и принтер, соединяющиеся через интерфейс для периферийных устройств вывода. На прилагаемых чертежах не приведены все возможные внутренние компоненты вычислительного блока, однако такого рода компоненты и соединения между ними хорошо известны. Настоящее изобретение описано на примере предпочтительных вариантов его осуществления, однако специалисту булет понятно, что данные варианты осуществления изобретения не ограничивают объем настоящего изобретения. Таким образом, подразумевается, что при реализации настоящего изобретения могут использоваться другие варианты его осуществления, а также в описанные варианты осуществления изобретения могут вноситься изменения, не противоречащие сути и объему настоящего изобретения, что определено формулой настоящего изобретения и эквивалентными признаками настоящего изобретения. ФОРМУЛА ИЗОБРЕТЕНИЯ 1. Способ определения характеристик гидроразрыва с использованием данных о микросейсмических событиях и их соответствующих местоположениях, содержащий этапы, на которых выявляют при помощи процессора компьютера один или более углов простирания по усредненным точкам микросейсмических данных, относящихся к одной или нескольким совокупностям плоскостей, и неопределенности их ориентаций, причем указанные один или более углов простирания имеют большее количество усредненных точек микросейсмических данных при соответствующем угле простирания, чем количество усредненных точек микросейсмических данных при других соответствующих углах простирания; определяют местоположение каждой плоскости в каждой совокупности плоскостей; и выполняют статистический анализ для каждой совокупности плоскостей с общим выявленным углом простирания для выявления неопределенности, связанной с ориентацией соответствующей совокупности подземных трещин, расстоянием между подземными трещинами соответствующей совокупности подземных трещин и углом падения плоскостей соответствующей совокупности подземных трещин. 2. Способ по п.1, отличающийся тем, что дополнительно включает в себя этап, на котором задают исходную ориентацию плоскости, причем эта ориентация плоскости содержит угол простирания и угол падения. 3. Способ по п.2, отличающийся тем, что включает в себя этап, на котором выявляют выявленные углы простирания по усредненным точкам микросейсмических данных, относящихся к одной или нескольким совокупностям плоскостей, и неопределенности их ориентации, причем для этого выполняют следующие этапы: создают плоскость, имеющую вектор нормали и проходящую через среднее местоположение всех точек микросейсмических данных для ориентаций плоскости; рассчитывают кратчайшее расстояние от каждой точки микросейсмических данных до плоскости; создают совокупность плоскостей, параллельных указанной плоскости, причем каждая плоскость в указанной совокупности плоскостей пересекает объемную область, представленную неопределенностью местоположения каждой точки микросейсмических данных; подсчитывают количество плоскостей в указанной совокупности плоскостей и точки микросейсмических данных, контактирующие с плоскостями в указанной совокупности плоскостей; рассчитывают усредненные точки микросейсмических данных, контактирующие с каждой плоскостью в указанной совокупности плоскостей, для каждой ориентаций плоскостей в указанной совокупности плоскостей; увеличивают угол простирания или угол падения указанной плоскости на 1; и повторяют все предыдущие этапы, пока угол простирания не станет равным 0 или 360, а угол падения не станет равным 0 или 90. 4. Способ по п.3, отличающийся тем, что определяют угол падения для каждой совокупности плоскостей и местоположение каждой плоскости в каждой совокупности плоскостей, причем для этого выполняют следующие этапы: повторяют все этапы по п.3 при каждом выявленном угле простирания для каждого соответствующего угла падения; находят в каждой совокупности созданных плоскостей плоскость, имеющую максимальное количество пересечений с точками микросейсмических данных, и определяют ее местоположение; рассчитывают весовой коэффициент контакта между точками микросейсмических событий и плоскостью, имеющей максимальное количество пересечений с точками микросейсмических данных; удаляют точки микросейсмических данных, контактирующие с плоскостью, имеющей максимальное количество пересечений с точками микросейсмических данных; и повторяют все предыдущие этапы до тех пор, пока не останется ни одной плоскости, имеющей максимальное количество пересечений с точками микросейсмических данных и контактирующей с тремя или большим количеством точек микросейсмических данных. 5. Способ по п.1, отличающийся тем, что неопределенность ориентации одной или нескольких совокупностей плоскостей определяют при помощи статистического анализа по методу Монте-Карло. 6. Способ по п.4, отличающийся тем, что неопределенность, связанную с ориентацией каждой совокупности подземных трещин, расстоянием между подземными трещинами в каждой совокупности подземных трещин и углом падения плоскости каждой совокупности подземных трещин, рассчитывают с учетом весового коэффициента контакта для каждой совокупности плоскостей и точек микросейсмических данных, связанных с каждой совокупностью плоскостей. 7. Способ по п.4, отличающийся тем, что расчет весового коэффициента контакта основывается на следующих параметрах: неопределенность местоположения каждой точки микросейсмических данных,кратчайшее расстояние от каждой точки микросейсмических данных до плоскости, числовое взвешенное значение, выражающее моментную магнитуду, связанную с каждой точкой микросейсмических данных 8. Способ по п.4, отличающийся тем, что весовой коэффициент контакта назначают каждой соответствующей точке микросейсмических данных. 9. Способ по п.1, отличающийся тем, что дополнительно включает в себя этап, на котором отображают статистический анализ. 10. Носитель программного обеспечения для длительного хранения данных, физически содержащий выполняемые компьютером команды для определения характеристик гидроразрыва с использованием данных о микросейсмических событиях и их соответствующих местоположениях, причем данные команды при их выполнении выявляют один или более углов простирания по усредненным точкам микросейсмических данных,относящихся к одной или нескольким совокупностям плоскостей и неопределенности их ориентации,причем указанные один или более углов простирания имеют большее количество усредненных точек микросейсмических данных при соответствующем угле простирания, чем количество усредненных точек микросейсмических данных при других соответствующих углах простирания; определяют местоположение каждой плоскости в каждой совокупности плоскостей; и выполняют статистический анализ для каждой совокупности плоскостей с общим выявленным углом простирания для выявления неопределенности, связанной с ориентацией соответствующей совокупности подземных трещин, расстоянием между подземными трещинами соответствующей совокупности подземных трещин и углом падения плоскостей соответствующей совокупности подземных трещин. 11. Носитель программного обеспечения по п.10, согласно командам которого задают исходную ориентацию плоскости, причем эта ориентация плоскости содержит угол простирания и угол падения. 12. Носитель программного обеспечения по п.11, обеспечивающий выявление выявленных углов простирания по усредненным точкам микросейсмических данных, относящихся к одной или нескольким совокупностям плоскостей, и неопределенности их ориентации, путем выполнения этапов, на которых создают плоскость, имеющую вектор нормали и проходящую через среднее местоположение всех точек микросейсмических данных для ориентации плоскости; рассчитывают кратчайшее расстояние от каждой точки микросейсмических данных до плоскости; создают совокупность плоскостей, параллельных указанной плоскости, причем каждая плоскость в указанной совокупности плоскостей пересекает объемную область, представленную неопределенностью местоположения каждой точки микросейсмических данных; подсчитывают количество плоскостей в указанной совокупности плоскостей и точки микросейсмических данных, контактирующие с плоскостями в указанной совокупности плоскостей; рассчитывают усредненные точки микросейсмических данных, контактирующие с каждой плоскостью в указанной совокупности плоскостей, для каждой ориентаций плоскостей в указанной совокупности плоскостей; увеличивают угол простирания или угол падения указанной плоскости на 1; и повторяют все предыдущие этапы, пока угол простирания не станет равным 0 или 360, а угол падения не станет равным 0 или 90. 13. Носитель программного обеспечения по п.12, обеспечивающий определение угла падения для каждой совокупности плоскостей и местоположения каждой плоскости в каждой совокупности плоскостей, путем выполнения этапов, на которых повторяют все этапы по п.12 при каждом выявленном угле простирания для каждого соответствующего угла падения; находят в каждой совокупности созданных плоскостей плоскость, имеющую максимальное количество пересечений с точками микросейсмических данных, и определяют ее местоположение; рассчитывают весовой коэффициент контакта между точками микросейсмических событий и плоскостью, имеющей максимальное количество пересечений с точками микросейсмических данных; удаляют точки микросейсмических данных, контактирующие с плоскостью, имеющей максимальное количество пересечений с точками микросейсмических данных; и повторяют все предыдущие этапы, пока не останется ни одной плоскости, имеющей максимальное количество пересечений с точками микросейсмических данных и контактирующей с тремя или большим количеством точек микросейсмических данных. 14. Носитель программного обеспечения по п.10, в котором неопределенность ориентации одной или нескольких совокупностей плоскостей определяют согласно содержащимся на нем программам при помощи статистического анализа по методу Монте-Карло. 15. Носитель программного обеспечения по п.13, в котором неопределенность, связанную с ориентацией каждой совокупности подземных трещин, расстоянием между подземными трещинами в каждой совокупности подземных трещин и углом падения плоскости каждой совокупности подземных трещин,рассчитывают согласно содержащимся на нем программам с учетом весового коэффициента для каждой совокупности плоскостей и точек микросейсмических данных, связанных с каждой совокупностью плоскостей. 16. Носитель программного обеспечения по п.13, в котором расчет весового коэффициента контакта основывается на следующих параметрах: неопределенность местоположения каждой точки микросейсмических данных, кратчайшее расстояние от каждой точки микросейсмических данных до плоскости,числовое взвешенное значение, выражающее моментную магнитуду, связанную с каждой точкой микросейсмических данных. 17. Носитель программного обеспечения по п.13, в котором весовой коэффициент контакта назначают каждой соответствующей точке микросейсмических данных. 18. Носитель программного обеспечения по п.10 дополнительно включает в себя этап, на котором отображают статистический анализ.