Scientific journal
Advances in current natural sciences
ISSN 1681-7494
"Перечень" ВАК
ИФ РИНЦ = 0,775

STATISTICAL PROCESSING TECHNIQUE AND INTERPRETATION OF THE RESULTS OF EXPERIMENTAL RESEARCHES OF FILTRATION CONSOLIDATION

Ageykina O.V. 1 Vorontsov V.V. 2 Sufyanov R.R. 1, 3
1 Industrial University of Tyumen
2 Department of Education Administration of the City of Tyumen
3 Research and Design Institute «Neftegazproekt» Tyumen Industrial University
A characteristic feature of the construction and operation of buildings and structures on water-saturated soils at the base is the unstable stress-strain state of the soil massif, including as a result of filtration consolidation, the formation and dispersion of excess pore pressure. Knowledge of the patterns of pore pressure development in water-saturated clay arrays allows you to calculate the parameters of the foundations of construction sites, organize monitoring of their stability, monitor the course of filtration consolidation, and improve the efficiency of geomechanical calculations. At the same time, the diversity of soils led to the emergence of a wide range of mathematical models. Observing and measuring the characteristics of an object as part of a mathematical description, it is important to correctly process and interpret the experimental data, which will allow evaluating and comparing the results for testing hypotheses to identify the essential properties and patterns of the filtration consolidation process. The numerical values ​​of the parameters obtained as a result of experimental studies can be further used only if they are reliable, that is, they correspond to certain levels and criteria of mathematical statistics. The article describes an iterative approach to determining the parameters of a mathematical model of filtration consolidation of water-saturated clay soils when comparing the excess pore pressure in the soil base under load with a hard stamp with a function describing its development over time. Based on the processing of experimental data, the numerical values ​​of the model parameters were obtained. Experimentally confirmed the conclusion about the loss of bearing capacity in the underlying layers of water-saturated bases due to a jump in excess pore pressure at maximum values ​​of the filtration path length. The necessity of developing a criterion for assessing the quality and completeness of the iterative process in processing the results of the experiment, which will allow to transfer the entire data processing process to automatic mode, is shown.
water saturated clay soils
filtration consolidation
model parameters
statistical processing
experimental data
1. Kosterin A.V., Skvortsov E.V. Seepage Consolidation under Plane Deformation of Elastic Half-Space // Fluid Dynamics. 2018. Т. 53. № 2. Р. 270–276. DOI: 10.1134/S0015462818020106.
2. Podgornova N.N. Filtration consolidation of clay soil in the base of road construction // Oil and Gas Studies. 2010. № 1.
P. 96–100 (in Russian).
3. Aleksandrov A.S., Kalinin A.L., Tsyguleva M.V. Distributing ability of sandy soils reinforced with geosynthetics // Magazine of Civil Engineering. 2016. № 6 (66). Р. 35–48 (in Russian).
4. Ageikina O.V., Ageikin V.N. Determination of parameters of traffic flows on commercial roads of Western Siberia for environmental calculations // Materialy dokladov na XIX mezhdunarodnoy nauchno-prakticheskoy konferentsii (g. Tyumen’, 17 marta 2017 g.). Tyumen’: Izd-vo TIU, 2017. P. 13–17 (in Russian).
5. Bazykin I.V., Litvinenko G.I. The use of multivariate statistical analysis in processing the results of measuring the strength of concrete of hydraulic structures // Gidrotekhnicheskoye stroitel’stvo. 2011. № 6. Р. 34–35 (in Russian).
6. Ageikin V.N., Vorontsov V.V., Podgornova N.N. Forecast of excess pore pressures in waterlogged soils // Geotekhnika: aktual’nyye teoreticheskiye i prakticheskiye problemy: Mezhvuzovskiy tematicheskiy sbornik trudov. SPb., 2006. Р. 79–83 (in Russian).
7. Bulychev A.A., Verkhoturov V.N., Gulyaev B.A. Modern methods of biophysical research. Praktikum po biofizike. Uchebnoye posobiye dlya biol. spets. Vuzov / Pod red. A.B. Rubina. M.: Vysshaya shkola, 1988. 359 р. (in Russian).

Характерной особенностью строительства и эксплуатации зданий и сооружений на водонасыщенных грунтах в основании является нестабильное напряженно-деформированное состояние грунтового массива, в том числе вследствие фильтрационной консолидации, образования и рассеивания избыточного порового давления. Знание закономерностей развития порового давления в глинистых водонасыщенных массивах позволяет рассчитать параметры фундаментов объектов строительства, организовывать мониторинг их устойчивости, производить наблюдения за ходом фильтрационной консолидации и повысить эффективность геомеханических расчетов. При этом многообразие грунтов обусловило появление широкого спектра математических моделей [1, 2].

Наблюдая и измеряя характеристики объекта в рамках того или иного математического описания, важно грамотно обработать и интерпретировать экспериментальные данные, что позволит оценить и сопоставить результаты для проверки гипотез, для выявления существенных свойств и закономерностей процесса фильтрационной консолидации. Полученные в результате экспериментальных исследований численные значения параметров могут быть далее использованы лишь в том случае, если они достоверны, т.е. соответствуют определенным уровням и критериям математической статистики [3–5].

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

Материалы и методы исследования

Ранее предложена математическая модель, описывающая фильтрационную консолидацию водонасыщенного грунта, а именно [2, 6]:

ageyk01.wmf (1)

ageyk02.wmf ageyk03.wmf

где N(t) – избыточное поровое давление, a – напряжения, необходимые для преодоления начального градиента фильтрации, C(t) – пористость грунта, k – коэффициент пропорциональности; kw – модуль объемного сжатия поровой воды, kf – коэффициент фильтрации, γw – удельный вес воды, L(y,z) – длина пути фильтрации воды

В представленном описании величины и характер избыточного порового давления зависит от знака дискриминанта характеристического уравнения:

ageyk04.wmf (2)

Конкретный случай расчета можно представить в следующем виде:

ageyk05.wmf (3)

Для параметра A принято обозначение

ageyk07.wmf

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

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

При подготовке статьи были использованы результаты исследований, проведенные на кафедре строительной механики ФГБОУ ВПО ТюмГАСУ в 2004–2006 гг. и в отделе криогенных ресурсов ООО «НИПИ Нефтегазпроект» при выполнении инженерных изысканий грунтовых карьеров на Западно-Чатылькинском лицензионном участке в 2018–2019 гг.

В ходе экспериментальных исследований напряженно-деформированного состояния водонасыщенного суглинка в условиях плоской деформации на грунтовых моделях варьировались следующие факторы: ширина загружающей пластины 50 и 20 см; величина полосовой распределенной нагрузки 4,6 кПа, 7,25кПа и 19,2 кПа; глубина заложения датчиков относительно полосы нагружения от 0,30 до 0,8 м. Грунтовое основание было представлено водонасыщенным суглинком нарушенной структуры, который послойно укладывался в лоток. Одновременно с укладкой грунта были установлены датчики для определения полных (общих) и избыточных поровых давлений в грунтовом массиве. Для этого использовали тензорезисторные мембранные месдозы и пьезометры. Сигналы датчиков регистрировали, независимо друг от друга, с использованием автоматического измерителя деформаций и измерительного регистрирующего комплекса [6].

Результаты исследования
и их обсуждение

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

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

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

D = 1,2734 – 0,00862*P – 0,6106*L/B, (4)

где P – равномерно-распределенная нагрузка, приложенная к жесткому штампу, Па; L – длина пути фильтрации, B – ширина жесткого штампа.

Достоинством рассматриваемой модели фильтрационной консолидации водонасыщенных грунтов является четкое отражение физической сущности реального процесса при относительной простоте математической формулировки. Кроме того, структура математического описания удобна для расчетов в конкретной пространственной точке модели при использовании подхода, реализованного в Delphi 6. При этом обеспечена возможность определения неизвестных параметров в автоматическом режиме, в интерактивном, а также прямым перебором в заданном интервале. Необходимо заметить, что определение показателей параметров модели базируется на подходах восстановления показателей экспонент и реализовано на основе известных идей [7].

Статистическая значимость уравнения проверена с помощью коэффициента детерминации и критерия Фишера, Fкp = 3,37 < F = 30,81. Установлено также, что параметры модели статистически значимы. Кроме того расчет коэффициента ранговой корреляции Спирмена наряду с тестом Голдфелда – Квандта позволили подтвердить гипотезу об отсутствии гетероскедастичности, расчет критерия Дарбина – Уотсона подтвердил отсутствие автокорреляции, а расчетное значение RS-критерия, 2,7 < RS = 3,656 < 3,7, показало адекватность модели по нормальности распределения остаточной компоненты.

Однако метод экспериментального определения параметров модели в целом для объекта потребовал более внимательного рассмотрения и детальной статистической обработки результатов. Дело в том, что в модели, как отмечалось ранее, заложены три расчетных случая в зависимости от знака дискриминанта характеристического уравнения (1–3). Поэтому при аппроксимации экспериментальных данных достаточно велик риск ошибки. В этом отношении показателен пример обработки показаний датчика № 4 при распределенной нагрузке 4,6 кПа.

В рассматриваемом случае, при первичной обработке результатов эксперимента с использованием программы в автоматическом режиме, определен случай D = 0, рис. 1, а. Однако интерпретация расчетов по всей совокупности показаний датчиков свидетельствует о том, что речь идет о случае – случай D < 0, рис. 1, б.

ag1.tif

а)

ag2.tif

б)

Рис. 1. Аппроксимация экспериментальных данных для датчика № 4, при значении дискриминанта D = 0 (а) и D < 0 (б)

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

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

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

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

Заключение

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

ag3a.tif ag3b.tif

ag3c.tif ag3d.tif

Рис. 2. Теоретические и экспериментальные значения параметров в модели фильтрационной консолидации: а) параметр λ1; б) параметр λ2 в) параметр В; г) параметр С

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

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

Показана необходимость исследований в части разработки и обоснования критерия оценки качества и завершенности итерационного процесса при обработке экспериментальных данных.