Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,222

ЧИСЛЕННОЕ РЕШЕНИЕ УРАВНЕНИЙ РАСТВОРЕНИЯ И ВЫНОСА СОЛЕЙ ПРИ ФИЛЬТРАЦИИ В ТРЕЩИНОВАТЫХ ПОРОДАХ

Магомедова А.В. 1 Иванов В.В. 1
1 ФГБОУ ВПО «Дагестанский государственный технический университет»
В статье рассматривается численное решение уравнений растворения и выноса солей при фильтрации в трещиноватых породах. Решена краевая задача двухмерной конвективной диффузии разностным методом. Решение позволяет прогнозировать растворение и вынос солей в основаниях гидротехнических сооружений. Приведен алгоритм решения краевой задачи двухмерной конвективной диффузии в прямоугольной области и приведены численные результаты расчета распределений концентраций раствора в заполненной трещине. По приведенному алгоритму для вычисления распределения концентраций по трещине составлена и от­лажена программа на языке TPascal 7.1. Для оценки интенсивности раскрытия трещины при раство­рении ее стенок рассмотрены средние по длине тре­щины приращения раскрытия и получены соответствующие выражения. Приведенные приемы позволяют прибли­женно учесть изменение размеров трещины в процессе раство­рения.
фильтрация
трещина
растворение солей
конвективная диффузия
порода
гидротехнические сооружения
1. Баламирзоев А.Г. Развитие теории и методов прогнозирования суффозионных деформаций при фильтрации в трещиноватых основаниях гидротехнических сооружений: дис. ... д-ра техн.наук. – Махачкала, 2006. – 397 с.
2. Баламирзоев А.Г., Косиченко Ю.М., Магомедова А.В. Расчет суффозионной осадки при неравномерном распределении растворимых солей в грунтах оснований гидротехнических сооружений // Изв. вузов. Сев.-кав. регион. Техн. науки. – 2002. – № 3. – С. 82–84.
3. Beригин Н.Н. Некоторые вопросы химической гидродинамики, представ-ляющие интерес для мелиорации и гидротехники.// Известия АН СССР, ОТН. – 1953. – № 10.
4. Орадовская А.Е. Фильтрационное выщелачивание дисперсно- распределенного гипса из песчано-глинистых пород // Растворение и выщелачивание горных пород. – М.: Госстройиздат, 1957.
5. Патpашев А.Н., Аpyтюнян Н.Х. Диффyзия солей пpи одномеpной фильтpации // Изв. НИИГ. – 1941. – т. 30. – С. 64–77.

При движении воды по трещинам пород, содержащих растворимые включения, происходит увеличение раскрытия трещин за счет раст­ворения их стенок. Процесс выноса солей включает как диффузион­ный транспорт их от поверхности растворения, так и конвективный перенос фильтрующимся потоком. Рассмотрим процесс для наиболее простой схемы одиночной заполненной дисперсным материалом трещи­ны длиной l и шириной 2h; скорость фильтрации считается постоян­ной и равной v (рисунок). Такого рода задача была исследована А.Н. Патрашевым [5] для условий стационарной диффузии и при пре­небрежении диффузионным переносом в направлении конвекции. Счи­талось также, что коэффициент диффузии одинаков по всем направ­лениям. Следующее аналитическое решение, свободное от указанных допущений, при­ведено в работе [2].

pic_34.tif

Расчетная схема растворения заполненной трещины

В связи с тем, что v = const, коэффициенты продольной и попе­речной конвективной диффузии D1 и D2 также будут постоянными, но различными по величине. Процесс растворения и выноса солей из трещины в таких условиях будет описываться уравнением нестацио­нарной конвективной диффузии в прямоугольной области 0 ≤ x ≤ l, 0 ≤ y ≤ 2h:

Eqn49.wmf (1)

с краевыми условиями:

C(x,y,0) = C(x,0,t) = C0; (2)

Eqn50.wmf (3)

. Eqn51.wmf (4)

Применяя к (1) интегральное преобразование Лапласа по переменной t с учетом начального условия (2), получим

Eqn52.wmf (5)

причем

Eqn53.wmf,

где p – параметр преобразования.

Граничные условия (2)–(4) также могут быть представ­лены в форме изображений:

Eqn54.wmf (6)

Eqn55.wmf (7)

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

Eqn56.wmf.

При этом параметр преобразования может принимать ряд значений, определяемых из уравнения cosλnh = 0, т. е.

Eqn57.wmf (8)

Определив изображения постоянной Eqn58.wmf и второй производной Eqn59.wmf, включающие гранич­ные условия (6), получим следующее одномерное уравнение для u(x, λn, p):

Eqn60.wmf (9)

Это уравнение должно решаться при граничных условиях

Eqn61.wmf (10)

Представим общее решение (9) в виде

Eqn62.wmf (11)

где

Eqn63.wmf

Подстановка (11) в условия (10) позволяет определить постоянные

Eqn64.wmf

С учетом этих значений постоянных частное решение (9) записывается в виде

Eqn65.wmf (12)

где Eqn66.wmf

Обратный переход к функции с осуществляется по формуле обращения

Eqn67.wmf

Имея в виду, что

Eqn68.wmf

получим

Eqn69.wmf (13)

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

Eqn70.wmf (14)

Здесь

Eqn71.wmf Eqn72.wmf

Eqn73.wmf Eqn74.wmf

Eqn75.wmf Eqn76.wmf

Eqn77.wmf Eqn78.wmf

– критерий Пекле; Eqn79.wmf – критерий Струхаля или кратность обмена жидкости в трещине.

В этих обозначениях перепишем уравнение (1) в следующем виде [1]

Eqn80.wmf (15)

Умножая обе части уравнения (15) на l2/D1, получим

Eqn81.wmf

или

Eqn82.wmf (16)

где Eqn83.wmfC = C/C0.

Для получения стационарного распределения концентраций в трещине уравнение (16) решаем разностным методом, называемым продольно-­поперечной разностной схемой. В этой схеме переход от слоя n к слою n + 1 осуществляется в два этапа. На первом этапе опреде­ляют промежуточные значения Eqn84.wmf из системы уравнений

Eqn85.wmf (17)

а на втором этапе, пользуясь найденными значениями Yij, находят Eqn86.wmf из системы уравнений

Eqn87.wmf (18)

Алгоритм решения уравнений (17) и (18). Перепишем уравнение (17) в виде

,Eqn88.wmf (19)

где

Eqn89.wmf

Eqn90.wmf

Eqn91.wmf

Eqn92.wmf

Уравнение (19) решается при каждом фиксированном j = 1, 2,..., N2 – 1 методом прогонки по переменному i. Чтобы применить прогонку, надо знать граничные значения Eqn93.wmf Eqn94.wmf В данном случаеEqn95.wmf, Eqn96.wmf. После того как все Eqn97.wmf найдены, решается уравнение (18). Переписывая это уравнение подробнее:

Eqn98.wmf (20)

где

Eqn99.wmf

Eqn100.wmf

видим, что при каждом фиксированном i = 1, 2,...N1 – 1 его можно ре­шить с помощью одномерной прогонки. Граничные условия в данной задаче Eqn101.wmf, Eqn102.wmf По приведенному алгоритму для определения распределения концентраций по трещине составлена и от­лажена программа на языке ТPascal 7.1.

В таблице приведены численные результаты расчета распределений концентраций раствора в заполненной трещине по указанной программе в сопоставлении с численными результатами, полученными по методу А.С. Малышева.

Распределение концентраций раствора в заполненной трещине при h0 = 0,1, D0 = 0,1, x = 1

y0\x0

По методу А.С. Малышева

По предлагаемому методу

0,0

0,2

0,4

0,6

0,8

 

0,0

0,2

0,4

0,6

0,8

 

0,0

0,2

0,4

0,6

0,0

0,0

0,0

0,0

1,0

0,896

0,658

0,591

1,0

0,935

0,857

0,805

1,0

0,927

0,869

0,810

1,0

0,831

0,695

0,603

 

0,0

0,0

0,0

0,0

1,0

0,795

0,646

0,614

1,0

0,943

0,866

0,794

1,0

0,891

0,786

0,715

1,0

0,827

0,712

0,660

 

Сравнительный анализ численных результатов показывает удов­летворительное соответствие между полученными численными резуль­татами по приведенной программе и результатами, полученными по методу А.С. Малышева.

Стационарное распределение концентраций в трещине (при ϑ → ∞) определяется зависимостью (14) без последнего чле­на [2].

В тех случаях, когда скорость фильтрации v значительна, коэффициент конвективной диффузии может определяться из соотношения D1,2 = λ1,2v. Критерий Пекле, в свою очередь, может быть представлен в виде ξ = l/2λ1. Так как параметр λ1 по величине близок к среднему размеру зерен заполнителя трещины, представляет интерес получение зависимостей, при­годных при больших значениях критерия Пекле (x ≥ 10). Оценивая для этого аргумент гиперболических функций в (13), найдем, что в области 0,1 < х0 < 0,9 с незначительной погреш­ностью можно принимать shz = ez/2. Тогда (13) упрощается и имеет вид

Eqn103.wmf (21)

Переход к оригиналу и введение принятых безразмерных величин дает следующую зависимость для распределения кон­центраций при больших значениях критерия Пекле (x ≥ 10).

Eqn104.wmf (22)

Здесь

Eqn105.wmf

– диффузионный критерий Фурье.

Из выражений (14) и (22) можно определить диффузи­онный поток солей от растворяющейся стенки трещины и, сле­довательно, приращение раскрытия трещины 2d в любом ее сечении на некоторый момент времени t. Очевидно, что

Eqn106.wmf (23)

Подставляя (14) в (23), получим после некоторых пре­образований следующую зависимость для d:

Eqn107.wmf (24)

Для оценки интенсивности раскрытия трещины при раство­рении ее стенок удобно рассматривать средние по длине тре­щины приращения раскрытия. Из (24) для этой цели может быть получено выражение

Eqn108.wmf (25)

Имея в виду, что при больших значениях ξ можно поло­жить Eqn109.wmf, а также заменяя гиперболические функции показательными Eqn110.wmf, получим для стационар­ного режима растворения

Eqn111.wmf (26)

Величина раскрытия трещины в процессе растворения

Eqn112.wmf, (27)

определяемая по зависимости (24), не остается постоянной, а меняется как по длине трещины, так и во времени. Следова­тельно, задача, рассмотренная выше при h = const, дает доста­точно точные результаты только при δ < 0,2h. Если это условие не соблюдается, расчет должен проводиться в две стадии. Сна­чала для заданного времени t по (25) или (26) определяется среднее раскрытие трещины <2h>, которое затем подставляет­ся во все итоговые зависимости. Этот прием позволяет прибли­женно учесть изменение размеров трещины в процессе раство­рения.

Рецензенты:

Агаханов Э.К., д.т.н., профессор, зав. кафедрой «Автомобильные дороги, основания и фундаменты», ФГБОУ ВПО «Дагестанский государственный технический университет», г. Махачкала;

Баламирзоев А.Г., д.т.н., профессор кафедры прикладной математики Махачкалинского филиала МАДИ, г. Махачкала.

Работа поступила в редакцию 22.04.2013.


Библиографическая ссылка

Магомедова А.В., Иванов В.В. ЧИСЛЕННОЕ РЕШЕНИЕ УРАВНЕНИЙ РАСТВОРЕНИЯ И ВЫНОСА СОЛЕЙ ПРИ ФИЛЬТРАЦИИ В ТРЕЩИНОВАТЫХ ПОРОДАХ // Фундаментальные исследования. – 2013. – № 6-3. – С. 546-550;
URL: http://fundamental-research.ru/ru/article/view?id=31544 (дата обращения: 25.04.2019).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.252