СПОСОБ ОПРЕДЕЛЕНИЯ КАРСТОВОЙ ОБЛАСТИ

1. Способ определения карстовой области, измененной явлением карстообразования, содержащий этапы, на которых
a) принимают геологическую модель (100, 300) карстовой области, содержащую ячейки, причем упомянутая модель содержит по меньшей мере одну исходную зону (101), содержащую ячейки модели, и по меньшей мере одну целевую зону (102), содержащую ячейки модели, и каждая из ячеек модели имеет координаты;
b) определяют исходную ячейку (103) из указанных ячеек указанной по меньшей мере одной исходной зоны;
c) определяют целевую ячейку (104) из указанных ячеек указанной по меньшей мере одной целевой зоны;
d) если расстояние между исходной ячейкой и целевой ячейкой превышает заданное расстояние (304):
d1) определяют (305а) зону (106) бифуркации на основе по меньшей мере координат исходной ячейки, координат целевой ячейки, распределения расстояний, углового распределения,
d2) определяют (305b) ячейку бифуркации (105а) в упомянутой зоне бифуркации;
d3) повторяют этап d) для ячейки (105b) исходя из того, что упомянутая ячейка (105а) бифуркации определена в качестве исходной ячейки;
e) определяют в геологической модели траекторные ячейки на основе исходной ячейки, целевой ячейки и ячеек бифуркации, определенных на этапе d);
f) для каждой текущей ячейки в траекторных ячейках
f1) определяют (307) канальные характеристики на основе заданного распределения плотности каналов, заданного распределения проходов каналов и траекторных ячеек, связанных с упомянутой текущей ячейкой;
f2) определяют характеристики ячеек на основе канальных характеристик упомянутой текущей ячейки.
2. Способ по п. 1, в котором этапы b)-f) многократно повторяют с одной и той же исходной зоной и целевой зоной.
3. Способ по п. 1, в котором упомянутая модель содержит множество исходных зон и множество целевых зон, а этапы b)-f) повторяют многократно с различными комбинациями исходной зоны из указанного множества исходных зон и целевой зоны из указанного множества целевых зон.
4. Способ по п. 1, в котором исходную ячейку определяют на этапе b) в соответствии со стохастическим процессом и согласно заданному распределению в указанной по меньшей мере одной исходной зоне.
5. Способ по п. 1, в котором целевую ячейку определяют на этапе с) в соответствии со стохастическим процессом и согласно заданному распределению в указанной по меньшей мере одной целевой зоне.
6. Способ по п. 1, в котором ячейка (105b), определенная на основе упомянутой ячейки (105а) бифуркации, является упомянутой ячейкой бифуркации.
7. Способ по п. 1, в котором ячейка (105b), определенная на основе упомянутой ячейки (105а) бифуркации, является ячейкой, соответствующей проекции ячейки бифуркации на линию (110а), определяемую указанными исходной ячейкой и целевой ячейкой.
8. Способ по п. 1, в котором определяют траекторные ячейки посредством выполнения этапов, на которых
упорядочивают указанные ячейки бифуркации в списке ячеек в соответствии с порядком определения на этапе d);
добавляют указанную исходную ячейку в первую позицию в упомянутом списке ячеек;
добавляют указанную целевую ячейку в последнюю позицию в упомянутом списке ячеек;
определяют сегменты в модели между двумя следующими друг за другом ячейками в указанном списке ячеек;
определяют ячейки модели на основе алгоритма линии аппроксимации.
9. Способ по п. 8, в котором алгоритм линии аппроксимации является алгоритмом, выбранным из группы, содержащей алгоритм Брезенхэма, алгоритм DDA и вычислительный алгоритм Run-Slice.
10. Способ по п. 1, в котором определяют траекторные ячейки посредством выполнения этапов, на которых
упорядочивают ячейки бифуркации в списке ячеек в соответствии с порядком определения на этапе d);
добавляют указанную исходную ячейку в первую позицию в упомянутом списке ячеек;
добавляют указанную целевую ячейку в последнюю позицию в упомянутом списке ячеек;
определяют кривую интерполяции, аппроксимирующую список ячеек в указанной модели;
определяют ячейки модели на основе алгоритма кривой аппроксимации.
11. Способ по п. 1, в котором при определении канальных характеристик
c
для ячейки, являющейся двумерной ячейкой, для каждого направления пространства, перпендикулярного по меньшей мере одному краю текущей ячейки, если край упомянутой текущей ячейки, перпендикулярный упомянутому направлению, является общим с траекторной ячейкой, связанной с упомянутой текущей ячейкой, определяют по меньшей мере один канал в упомянутой текущей ячейке в упомянутом направлении, при этом радиус r упомянутого канала определен в упомянутом заданном распределении проходов каналов, или
c
для ячейки, являющейся трехмерной ячейкой, для каждого направления пространства, перпендикулярного по меньшей мере одной стороне текущей ячейки, если сторона упомянутой текущей ячейки, перпендикулярная упомянутому направлению, является общей с траекторной ячейкой, связанной с упомянутой текущей ячейкой, определяют по меньшей мере один канал в упомянутой текущей ячейке в упомянутом направлении, при этом радиус r упомянутого канала определен в упомянутом заданном распределении проходов каналов.
12. Способ по п. 1, в котором при определении канальных характеристик
для ячейки, являющейся двумерной ячейкой, для каждого направления пространства, перпендикулярного по меньшей мере одному краю текущей ячейки, если край упомянутой текущей ячейки, перпендикулярный упомянутому направлению, является общим с траекторной ячейкой, связанной с упомянутой текущей ячейкой,
c
определяют количество N каналов, которые должны быть определены на основе упомянутого заданного распределения плотности каналов;
c
определяют N каналы в упомянутой текущей ячейке в упомянутом направлении; или
для ячейки, являющейся трехмерной ячейкой, для каждого направления пространства, перпендикулярного по меньшей мере одной стороне текущей ячейки, если сторона упомянутой текущей ячейки, перпендикулярная упомянутому направлению является общей с траекторной ячейкой, связанной с упомянутой текущей ячейкой,
c
определяют количество N каналов, которые должны быть определены на основе упомянутого заданного распределения плотности каналов;
c
определяют N каналы в упомянутой текущей ячейке в упомянутом направлении.
13. Способ по п. 11 или 12, в котором при определении характеристик ячейки выполняют этапы, на которых
m
получают начальную характеристику пористости ϕ текущей ячейки;
eq
m
устанавливают эквивалентную пористость ϕ равной ϕ;
по меньшей мере для одного текущего канала в текущей ячейке в упомянутом направлении
temp
c
вычисляют пористость ϕ текущей ячейки на основе выражения , где S - сечение текущего канала и S - сечение упомянутого края или стороны упомянутой текущей ячейки, которые перпендикулярны упомянутому направлению и являются общими с траекторной ячейкой;
eq
temp
устанавливают эквивалентную пористость ϕ равной ϕ.
14. Способ по п. 11 или 12, в котором при определении характеристик ячейки дополнительно выполняют этапы, на которых
m
получают характеристику K начальной проницаемости текущей ячейки;
eq
m
устанавливают эквивалентную проницаемость K равной K;
по меньшей мере для одного текущего канала в текущей ячейке в упомянутом направлении
temp
c
c
вычисляют проницаемость K текущей ячейки на основе выражения , где S - сечение текущего канала, S - сечение упомянутого края или стороны упомянутой текущей ячейки, которые перпендикулярны упомянутому направлению и являются общими с траекторной ячейкой, и K - проницаемость текущего канала;
eq
temp
устанавливают эквивалентную проницаемость K равной K.
15. Способ по п. 14, в котором проницаемость текущего канала вычисляют на основе выражения .
16. Считываемый компьютером носитель для хранения данных, имеющий хранящуюся на нем компьютерную программу, содержащую программные команды, причем компьютерная программа пригодна для загрузки в блок обработки данных и выполнена с возможностью вызывать выполнение блоком обработки данных этапов любого из п.п. 1-15 формулы изобретения, когда компьютерная программа выполняется на устройстве обработки данных.