Дослідження густини крейдяної товщі методом Монте-Карло на проммайданчику Рівненської АЕС із застосуванням 3D моделей

Abstract

Для дослідження густини товщі крейдяних відкладів на території проммайданчика Рівненської АЕС на основі спектрального розкладу розроблено нову ефективну методику статистичного моделювання випадкових 3D полів у застосуванні до задач геофізичного моніторингу довкілля. Розглянуто задачу чисельного моделювання реалізацій випадкового поля “шуму” для даних зміни густини крейдяної товщі у тривимірному просторі. Наведено статистичні моделі для гаусcових однорідних ізотропних випадкових полів у тривимірному прос¬торі, заданих статистичними характеристиками. На основі таких моделей, сформульованих алгоритмів і побудованих програм згенеровано реалізації випадкових 3D полів із кореляційними функціями різних типів. Проведено ефективне порівняння точності середньоквадратичної апроксимації моделювання запропонованими методами та методом ТВМ (turning band method).

Для исследования плотности толщи меловых отложений на территории промплощадки Ривненской АЭС на основе спектрального разложения разработана новая эффективная методика статистического моделирования случайных 3D полей в применении к задачам геофизического мониторинга окружающей среды. Рассмотрена задача численного моделирования реализаций случайного поля “шума” для данных изменения плотности меловой толщи в трёхмерном пространстве. Приведены статистические модели для гауcсовых однородных изотропных случайных полей в трёхмерном пространстве, которые заданы статистическими характеристиками. На базе таких моделей, сформулированных алгоритмов и построенных программ сгенерированы реализации случайных 3D полей с корреляционными функциями разных типов. Проведено эффективное сравнение точности среднеквадратической аппроксимации моделирования предложенными методами и методом ТВМ (turning band method).

The article is devoted to the application of the theory and methods of 3D random fields statistical simulation (Monte Carlo methods) to environmental geophysical monitoring problems. To investigate chalk layer density on the Rivne NPP industrial site a new effective statistical technique has been devised to simulate random fields in 3D space, based on spectral decomposition. The 2D data were selected from 3D density data of chalk rock strata at three depth levels (28, 29, 30 m from the surface). At each level, the data were presented as the sum of deterministic and random components. The deterministic 2D trend surface was constructed using spline interpolation. The random component (the so called “noise”) is a 2D homogeneous isotropic random field. The authors considered the problem of statistical simulation of “noise” for chalk layer density realizations as random fields in 3D space. Statistical models have been constructed for the gauss random fields in three-dimensional space given by their statistical characteristics. Using these models, formulated algorithms and created programs, the authors have obtained 3D random fields realization with difference Bessel types, Cauchy types, and hole effect with certain parameters values. 300 additional values were simulated in the intervals between observation points for each level by constructing original programs Spectr 3 and Spectr 3_1 based on the chosen statistical models. The authors compared mean-square errors of simulation made by the proposed methods and the ТВМ (turning band method) method. Statistical simulation method of random processes and fields in 3D space was I ntroduced based on spectral decompositions, in order to enhance map accuracy with chalk layer density data. The paper suggests a universal method of statistical simulation of geophysical data to generate random 3D fields’ realizations on grids with required accuracy and regularity.