Abstract and keywords
Abstract (English):
The article describes suggested new method for finding influence zone of sample, which is based on the mathematical modeling of geological parameter variability anisotropy. This influence zone is area bounded by pedal curve, and it is limited by correlation radius that depends on direction.

Keywords:
mineralization heterogeneity, geological parameter, sample, variability, anisotropy, correlation, vector, matrix, tensor, pedal curve
Text
Publication text (PDF): Read Download

Оценки геологических показателей и их погрешностей зависят от  геологических и технологических факторов, а достоверность оценок от адекватного отражения используемым математическим аппаратом геолого-технологических факторов.

Геологические факторы представляют собой особенности внутреннего строения месторождений полезных ископаемых, связанные с периодами формирования залежей и воздействием диагенеза, эпигенеза, метаморфизма, тектонических деформаций, процессов трещинообразования. Внутреннее строение залежей характеризуется варьированием от непрерывного оруденения до крайне неравномерной дискретной минерализации.

Непрерывность оруденения участков недр может быть изотропно-однородной. В этом случае зона влияния пробы простирается на некоторое одинаковое расстояние во всех направлениях (окружность, сфера – частный случай подеры, подероида). Но чаще всего оруденение анизотропно-неоднородное и зона влияния пробы зависит от направления в пространстве. При  крайне неравномерном оруденении зона влияния пробы практически ограничивается её объёмом.

Технологические факторы представляются формой, размерами проб и расстояниями между пробами. Большим размерам проб отвечают большие зоны их влияния.

Количественное описание свойств оруденения недр выражается понятием изменчивости геологических показателей. Под изменчивостью понимают принятие геологическим показателем в разных точках пространства различных значений. В зависимости от характеристик оруденения изменчивость проявляется как случайная, либо коррелированная, либо координированная (функциональная). Причём коррелированная изменчивость представляется коррелированной изотропной и коррелированной анизотропной изменчивостями.

Крайне неравномерному оруденению, либо большим расстояниям между пробами отвечает случайная изменчивость геологических показателей. Для описания случайной изменчивости применяют аппарат математической статистики следующим образом. Данные опробования отождествляют со взаимно-независимыми значениями случайной величины, которая полностью определяется законом распределения вероятности и основными характеристиками- математическим ожиданием и дисперсией.

Непрерывному оруденению в пределах зон влияния проб соответствуют коррелированная, либо координированная изменчивости геологических показателей. Математическое моделирование коррелированной изотропной изменчивости осуществляется на основе стационарных эргодических случайных функций, основной характеристикой которых является корреляционная функция. Корреляционная функция отражает связи между значениями переменной, устанавливает предельные расстояния распространения этих связей и на её основе определяют оценки и погрешности оценок геологических показателей.

Коррелированную анизотропную изменчивость описывают математическими нестационарными анизотропными моделями путём  отождествления значений геологического показателя с реализацией нестационарной случайной функцией, которая в каждом фиксированном направлении представляет собой стационарную эргодическую функцию [1]. При этом корреляционная функция выражается тензором второго ранга, который описывает анизотропную структуру связей поля значений геологического показателя.

В случае координированной изменчивости значения геологического показателя интерполируют полиномами различных степеней. Погрешности интерполирования называют ошибками аналогий и выражают их разностями на порядок выше степени интерполяционного многочлена [1].

Исходя  из изложенного, актуальна проблема определения зон влияния проб. Ввиду того, что изотропная изменчивость является частным случаем анизотропной изменчивости показателей, будем решать эту проблему для коррелированной анизотропной изменчивости.

Не ограничивая общности, рассмотрим решение данной проблемы в плоскости, что отвечает погоризонтному изучению геологических показателей и отработке полезных ископаемых. Согласно [1] анизотропная изменчивость геологического показателя выражается в векторно-тензорной форме

 

 

 

Kh;α=nKhnT =cosα;sinαKxx(h)Kxy(h)Kyx(h)Kyy(h)cosαsinα,                                  (1)

 

где  Kh;α  – корреляционная функция, описывающая силу связей между значениями геологического показателя в зависимости от шага h  и направления в плоскости α;  Kh  − тензор второго ранга, выраженный в выбранной на месторождении системе координат квадратной неособенной симметричной матрицей; n, nT − соответственно единичный и транспонированный единичный векторы направляющих косинусов;  Kxxh, Kyy(h)  – корреляционные функции, отражающие связи между значениями геологического показателя в зависимости от шага h  соответственно по направлениям осей Ox  и Oy ; Kxyh=Kyx(h)  – зависящие от шага h  ковариационные функции, отражающие связь значений геологического показателя между направлениями в плоскости.

Главные (собственные) оси тензора Kh  являются главными осями анизотропии коррелированной изменчивости. Приведём тензор Kh  к главным осям и получим [1] направления наибольших и наименьших связей

tgφ1,2h=λ1,2h-Kxx(h)Kxy(h)=Kxy(h)λ1,2h-Kyy(h),  (2)

и корреляционные функции, отражающие экстремальные значения связей по этим направлениям

 

λ1,2h=Kxxh+Kyy(h)2±Kxxh-Kyy(h)22+Kxy2h .                                 (3)

 

Тогда корреляционная функция (1) в системе главных осей Ox'y'  примет наиболее простую каноническую форму

 

Kh;α'=nα'K'hnTα'=cosα';sinα'λ1h00λ2hcosα'sinα' ,                  (4)

 

где  α' – определяющий направление угол в системе координат главных осей; λ1hλ2h  − корреляционные функции в направлениях главных осей, отражающие наибольшие и наименьшие связи по этим направлениям.

Обозначим через  rα'  значение радиуса корреляции в направлении α'.  Радиус корреляции должен удовлетворять двум условиям:

  1. Корреляционные функции λ1h и λ2h  соответственно определяют радиусы корреляций r1  и r2 , которые принимают наибольшее и наименьшее значения и устанавливают предельные расстояния распространения связей между значениями геологического показателя по направлениям осей (2). Поэтому

minr1;r2rα'maxr1;r2.          (5)

  1. Из симметричности связей (4) относительно главных осей следует симметрия значений радиусов корреляции, т.е.

Kh;α'=Kh;-α'rα'=r-α'.    (6)

Получим выражение для радиуса корреляции rα' . Поставим в соответствие экстремальным значениям r1  и r2  вектор

r=r1; r2.                               (7)

Найдём проекцию вектора (7) на направление α' , которая выражает значение радиуса корреляции в направлении α'  и равна скалярному произведению векторов r  и nα'=cosα';sinα'

rα'=npα'r=rnα'=r1cosα'+ r2sinα'.  (8)

Проекция (8) точно отражает значения радиусов корреляции в направлениях главных осей изменчивости r0=r1  и rπ2=r2 . Однако по другим направлениям 0<α'<π2   приближённо и не удовлетворяет условию (5). Действительно, при  α'=arctgr2r1  получим rα'=r12+r22 , что превышает и r1<r12+r22 , и  r2<r12+r22 .

Тем не менее, на базе векторного представления радиуса корреляции (7), (8) перейдём к более сложному выражению радиуса корреляции, основанному на тензоре второго ранга, который может, удовлетворяя свойствам (5), (6), адекватно выразить радиус корреляции в зависимости от направления α' . Возведём в квадрат проекцию (8)

r2α'=r12cos2α'+r1r2sin2α'+ r22sin2α',  (9)

который также содержит указанную в проекции (8) погрешность. Исключим её следующим образом. Потребуем, чтобы равенство (9) удовлетворяло свойству (6). Для этого вычислим его в симметричном относительно главной оси Ox'  направлении -α'

r2-α'=r12cos2α'-r1r2sin2α'+ r22sin2α' , (10)

вычтем его из равенства (9) и в силу требования свойства (6) получим

2∙r1r2sin2α'=r2α'-r2-α'≡0.     (11)

Тогда из равенства (9) с учётом (11) определим выражение квадрата радиуса корреляции в зависимости от направления α'

r2α'=r12cos2α'+ r22sin2α'.         (12)

Выражение (12) целесообразно представить в тензорной (матричной) форме. Обозначим через R'  тензор корреляции, представленный в главных осях Ox'y'  тензора Kh  неособенной диагональной матрицей, и единичный вектор направляющих  косинусов  nTα' , полученный транспонированием орта nα' :

R'=r1200r22  и  nTα'=cosα'sinα'.       (13)

Тогда из равенства (12) выразим радиус корреляции в зависимости от направления α'  через скалярное произведение тензора корреляции R'  и единичных векторов направляющих косинусов nα', nTα' :

rα'=cosα';sinα'r1200r22cosα'sinα'12 . (14)

Легко увидеть, что радиус корреляции (12), (14) удовлетворяет свойству (6) rα'=r-α' , а из теории матриц [2] следует, что стоящие на главной диагонали элементы r12  и r22  тензора корреляции  R'  (13), (14) являются экстремальными значениями выражений (12), (14), т.е. выполняется условие (5).

Формулы (12), (14) выражают зависимость радиуса корреляции rα'  от направления α'  в системе главных осей Ox'y'  тензора Kh  и, следовательно, тензора R'.  При этом график радиуса корреляции образует замкнутую кривую, ограничивающую зону влияния пробы. Установим вид этой кривой и изобразим её форму.

Предварительно заметим, что система главных осей Ox'y'  тензоров Kh  и R'  является декартовой ортогональной системой координат и образует угол  φ  с выбранной на месторождении системой координат Oxy .

Действительно, произведение левых и правых  частей равенств (2) с учётом (3) даёт соотношение

tgφ1htgφ2h=λ1h-Kxx(h)λ2h-Kxx(h)Kxy2=-1,  (15)

свидетельствующее о взаимной ортогональности главных осей Ox'  и Oy' . Тогда обозначим через φ1h=φ  – угол между осями Ox  и Ox'  и получим φ2h=φ+π2  – угол между осями Ox  и Oy' . Откуда следует, что φ – угол между принятой на месторождении системой координат Oxy  и системой Ox'y'  главных осей тензора R' .

Формулы радиуса корреляции (12), (14) являются в системе Ox'y'  главных осей тензора R'  уравнениями подеры (от греч. – нога) эллипса с полуосями r1, r2  относительно его центра [1]. Кривая-подера представляет собой геометрическое место оснований перпендикуляра, опущенных из начала координат на касательные к эллипсу. На рис. 1 приведены подеры, ограничивающие зоны влияния проб, для коррелированной изотропной и коррелированной анизотропной изменчивостей геологических показателей, отвечающие разным значениям коэффициента анизотропии  k=r1r2 .

Рис. 1. Кривые подер, ограничивающие зоны влияния проб, отвечающие изменчивостям:

1 – коррелированной изотропной r1=r2k=1 ;

2 – коррелированной анизотропной

 r1=2r2k=2 ;  3 – коррелированной

 анизотропной r2=0,  k=∞

 

Для удобства практического  применения радиуса корреляции (14), выраженного в системе Ox'y'  главных осей тензора R' , целесообразно его представить в принятой на месторождении системе координат Oxy .

Переход от системы координат Ox'y'  к принятой на месторождении системе Oxy  осуществляется с помощью матрицы поворота [2]

Q=cosφ-sinφsinφcosφ ,                   (16)

где что φ – угол между системами Ox'y'  и Oxy . При этом матрица Q  – ортогональная и удовлетворяет равенствам

 

QQT=QTQ=EQT=Q-1=cosφsinφ-sinφcosφ ,                                  (17)

 

где  E  − единичная матрица, а QT, Q-1  – соответственно транспонированная и обратная матрицы.

 

 

 

Теперь тензор R'  (13) и радиус корреляции rα'  (14) в принятой на месторождении системе координат Oxy  с учётом (16), (17) примут соответственно вид

 

R=QR'Q-1=r12+r22tg2φr12-r22tgφr12-r22tgφr12tg2φ+r22cos2φ ,                          (18)

r(α)=cosα;sinαr12+r22tg2φr12-r22tgφr12-r22tgφr12tg2φ+r22cosαsinα12cosφ,               (19)

 

где  α  – определяющий направление угол в системе координат Oxy.

Полученные выражения радиусов корреляции (19), (14) в зависимости от направления в пространстве и их геометрическая интерпретация (см. рис.) является важным результатом, так как только в пределах зон корреляции проб правомерна интерполяция, аппроксимация и экстраполяция разведочных данных [3].

References

1. Red'kin G.M. Nestacionarnoe anizotropnoe matematicheskoe modelirovanie neodnorodnostey sistem mineral'nogo syr'ya. M.: Izd-vo Associacii stroitel'nyh vuzov, 2007. 500 s.

2. Gantmaher F.R. Teoriya matric. M.: Nauka, 1967. 575 s.

3. Kazhdan A.B. Metodologicheskie osnovy razvedki poleznyh iskopaemyh. M.: Nedra, 1974. 270 s.


Login or Create
* Forgot password?