Формулировка математической модели глубоководного нефтеподъемника в условиях действия внутреннего потока газожидкостной смеси Текст научной статьи по специальности «Физика»

ОА Торопова. Формулировка математической модели глубоководного нефтеподъемника

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

Библиографический список

1. Мирзаджанзаде А.Х., Аметов И.М., Ковалев А.Г. Физика нефтяного и газового пласта. М.; Ижевск: Ин-т комп. исслед., 2005.

2. Голубева О.В. Курс механики сплошных сред. М.: Высш. шк., 1972.

3. Толпаев В.А., Ледовской В.И. Уравнения линейной двумерной фильтрации в искривленных пластах конечной толщины // ОПиПМ. 2004. Т.11, вып.1. С.143-146.

4. Толпаев В.А., Ледовской В.И. Математическое моделирование фильтрационных течений несжимаемой жидкости в искривленных пластах конечной толщи-

ны // ОПиПМ. 2005. Т. 12, вып. 2. С. 524-527.

5. Толпаев В.А., Ледовской В.И. Двумерные математические модели течений жидкости в круговом коническом слое постоянной толщины // Современные проблемы прикладной математики и математического моделирования: Материалы конф. Воронеж: Воронежск. гос. техн. акад., 2005. С. 222.

6. Толпаев В.А., Ледовской В.И. Оценки точности расчета дебитов скважин в искривленных пластах // Нефтепромысловое дело. 2004. № 12. С. 9-13.

УДК 539.3

ФОРМУЛИРОВКА МАТЕМАТИЧЕСКОЙ МОДЕЛИ ГЛУБОКОВОДНОГО НЕФТЕПОДЪЕМНИКА В УСЛОВИЯХ ДЕЙСТВИЯ ВНУТРЕННЕГО ПОТОКА ГАЗОЖИДКОСТНОЙ СМЕСИ

О.А. Торопова

Саратовский государственный технический университет, Институт социально-производственного менеджмента, кафедра информационных систем E-mail: toropova@sstu.ru

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

The Formulation of the Mathematical Model Deep-Water of the Oil Raiser, Interacting with Internal Flow of Gas-Liquid Mixture

O.A. Toropova

A new mathematical model for calculating the characteristics of a stress-deformated state of a deepwater raiser that interacts with an inner two-phase stream of the gas-liquid mixture, is formulated. The singular-disturbance type problem of nonlinear hydroelasticity is being solved.

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

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

Сформулируем модель установившегося движения газожидкостной смеси в вертикальном трубопроводе. Характер движения (эрлифта) газожидкостной смеси определяется, как правило, значением фактора газонасыщенности и (отношением объема свободного газа в смеси к ее объему, 0 < и < 1).

При постепенном увеличении этого фактора происходит последовательный переход от эмульсионного состояния смеси к диспергированному:

1) движение смеси в состоянии эмульсии, когда пузырьки газа более или менее равномерно распределены в жидкости;

2) слоистое движение, когда в трубопроводе перемещаются, заполняя почти все сечение, цилиндрики газа значительной длины, разделенные цилиндриками жидкости меньшей длины;

© О.А. Торопова, 2007

Известия Саратовского университета. 2007. Т. 7. Сер. Математика. Механика. Информатика, вып. 2

3) газ с частицами жидкости движется в осевой части трубы, а жидкость — тонким слоем по ее стенкам;

4) диспергированное движение жидкости (ш ~ 1) — это движение газа, заполняющего все сечение трубопровода и увлекающего за собой частицы жидкости.

Последний случай в реальных условиях эксплуатации нефтяных скважин практически не встречается, так как он связан с недостаточной сепарацией газа от жидкости при низких давлениях. Далее, предельное количество газа, которое можно растворить в единице объема жидкости с увеличением давления, возрастает. Поэтому существует значение давления P*, ниже которого начинается выделение газа (оно называется давлением насыщения). Движение газированных жидкостей при значениях давлений P < P* происходит в смеси со свободным газом, а по мере снижения давления (т.е. с ростом ш) происходит постепенное увеличение его массы в смеси. Рассмотрим случай установившегося движения смеси, принимая для газа изотермический закон состояния, а его растворимость в жидкости, подчиняющейся закону Генри [1], Vg = а0р 0Vf, где Vg — объем стандартного (т.е. при нормальном атмосферном давлении и при некоторой нормальной температуре) газа, растворенного в жидкости, Vf — объем жидкости, a0 — объемный коэффициент растворимости газа (эмпирическая безразмерная функция давления, в расчетах обычно полагают a0 = const), p 0=р/атм. Тогда

где рс — плотность газожидкостной смеси, pg, pf — плотности стандартного газа и жидкости. Если пренебречь эффектом усадки (т.е. увеличением удельного веса жидкости по мере снижения количества растворенного в смеси газа), то в формуле (1) pg = const, pf = const, ш = ш(р0): 0 < ш < 1 при

0 < р 0 < р*, ш = 0 при р 0 > р*.

Определим вектор распределенных гидродинамических нагрузок со стороны внутреннего потока гидросмеси и внешнего потока окружающей жидкости. Будем считать, что райзер (гибкий стержень кольцевого поперечного сечения с внутренним и внешним диаметрами d0, D0) в недеформированном состоянии имеет прямолинейную форму. Тогда формулами перехода от ортов локальной подвижной системы координат т, n к единичным базисным векторам ii, i2 осей декартовой правой системы отсчета х\Ох2 (начало которой выбрано на дне океана глубиной H) являются [2]

при этом йт/йв = кп, йп/йв = — кт. Здесь 5 — эйлерова дуговая координата, начало которой (в = 0) выбрано в нижнем сечении райзера, р — угол отклонения касательной к осевой линии деформированного райзера по отношению к вертикальной координате х2, к — кривизна осевой линии. В осях связанной системы координат вектор распределенных гидродинамических нагрузок (т), действующих на элемент райзера, можно представить в виде

где <Г/, Ти, — векторы сил взаимодействия элемента райзера с внутренним и внешним потоками жидкости, т — погонная масса райзера в воздухе. Вектор сил взаимодействия между частицами жидкости и газа (/) не входит в выражение (3). Действительно, пусть У/ = У/т — вектор скорости движения центра масс жидкости в элементарном объеме смеси, а Уд = Удт — вектор скорости движения центра масс свободного газа относительно скорости движения центра масс жидкости в том же объеме. Динамические уравнения Эйлера движения элементов жидкости и газа [1] запишутся для рассматриваемого случая в виде

т/(йУ//б£) = — (рР/т)/в — т/ д ею — / + /; тд(й(Уд + У/)/&) = — (рРдт)/в — тд д ею — /.

Складывая эти уравнения, находим / = —т/(йУ//б£) — тд(й(Уд + У/)/б£) — (Рст)/в — тс д ею. Здесь т/ = р/Р/; тд = рдРд; тс = рсР0; Рс = рР0; ею = соб рт — бій рп; Р/; Рд — части площади поперечного сечения райзера Р0 = 0.25пй2, занимаемые жидкостью и газом: Р0 = Р/ + Ри. Так как (Г/ ± т и / = (/ ,т)т + т х (/ х т), (/ х т) = 0, то

рс = ир 0pg + (І - u)pj-,

(І)

т = — sin ф i1 + cos ф i2,

if = — cos ф i1 — sin ф i2,

(2)

q = qf + qw — mg cos ф т + mg sin ф П,

(З)

Pc/s + ((mf (dvf /dt) + mg (d(Vg + Vf )/dt)),r) + mc g(ew ,т) = О,

qf = qhc + qac, qhc = —kPcf — mcgr x (ew x т) = —(kPc — mcg sin ф)П, qac = —т x ((mf (dvf /dt) + mg(d(Vg + Vf )/dt)) x т).

(4)

(5)

(6)

ОА Торопова. Формулировка математической модели глубоководного нефтеподъемника

Формулы (5)-(6) представляют собой обобщенные выражения для векторов нормальной силы гидростатического давления на стенки райзера (qhc) и нормальной силы инерции присоединенной массы газожидкостной гидросмеси (qac) (при отсутствии свободного газа в формулах (4)-(6) следует положить Pc = Pf, mg = ш = vg = 0 и они переходят в аналогичные выражения для нагрузок, полученные в [2]).

Преобразуем скалярные компоненты формул (4), (6). Имеем:

mf (dvf /dt) = mf ((vf T)/s)vf = mf Vf (vf /s)t + kmf v'jn, mg (d(vg + vf)/dt) = mg (vf + vg )((vg + vf)/s)q + kmg (vg + vf Для недиспергированного движения жидкости обычно vg ~ const, кроме того, нетрудно показать, что Fg ~ ш^0, Ff ~ (1 — шЩ, поэтому соотношения (4)-(6) преобразуются окончательно к виду

р/s + vf /s(pcvf + Pg(шгиg) + pcgcos у = 0, (7)

qhc = —(кр — Pcg sin p)F0ri, qac = —k(pf (1 — ш)v2 + шpg (vg + vf )2 )F0n, (8)

Pc = шр0 pg + (1 — ш)pf.

Пусть на райзер действует стационарный поток подводных течений со скоростью vc = vciT+vc2n = = —(vc sin yr + vc cos yn), vc = vc(x2), 0 < x2 < H. Вектор сил взаимодействия элемента райзера с внешним потоком окружающей жидкости Тш состоит из векторов нормальных сил гидростатического давления на стенки (qhw), инерции присоединенной массы жидкости (qaw) и гидродинамического сопротивления (qn)[2]:

qhw = (kPw — mw g sin y)n = (kmw g(H — X2) — mw g sin y)n,

Qaw = —т x ((mw(dvaw/dt)) x т) = —kmwv2c sin2 yn, qn = —0.5cnpwD0vc cos y|vc cos y| n (mw = pwFw,Fw = 0.25nD2).

Таким образом,

q = —mgcosy т + [(m + mc — mw)gsinф — k(pF0 — mwg(H — x2) + (pf (1 — u)v‘2 + +upg (vg + Vf )2 )F0 + mw v2 sin2 ф) — 0.5cn pw Dovc cos ф | vc cos ф | ]f.

(9)

Для определения зависимости входящего в (9) фактора газонасыщенности ш от давления р0 воспользуемся уравнениями неразрывности для потока гидросмеси. Они представляют собой условия постоянства весовых расходов стандартной жидкости (Qf)(pf — а0р0pg) x Qf = const и газа (Qg)pg(Qg + а0р0pfQf) = const, что позволяет сформулировать выражение для весового газового фактора (при эсплуатации нефтяных скважин его определяют обычно из замеров дебитов газа и нефти после сепарации).

Г = pg («0 + а0 р °)/(pf — а0р 0 pg) = const, «0 = Qg/Qf,

откуда

a0 = (pf /pg )Г — а0р 0 (Г + 1). (10)

С другой стороны,

«0 = р 0ш(vf + vg)/vf (1 — ш) = (р 0ш/(1 — ш))(1 + vg/vf). (11)

Из (10)—(11), учитывая, что vf = vf0/(1 — ш) (vf0 > 0 — известное значение начальной скорости жидкости на входе s = 0), получаем уравнение

в0 ш + ш/(1 — ш) = 01, 00 = vg/vf 0, 01 = а0/р 0 • (12)

Его корень, удовлетворяющий условию 0 < ш < 1,

201

ш =-------------------, = (13)

1 + 00 + 0i + \/(1 + e0 — 0i)2 + 401

Механика 55

Известия Саратовского университета. 2007. Т. 7. Сер. Математика. Механика. Информатика, вып. 2

дает искомую, нелинейную зависимость ш = ш(р0), 0 < р0 < р,0. Давление насыщения р,0 определяется из (10), принимая там а0 = 0:

р,0 = (pf /(а0 pg))(Г/(Г+1)). (14)

В этом случае в (13) в1 =0 и ш = 0. Видно также, что при р0 ^ +0 (в1 ^ +го) ш(р 0) ^ 1.

Сформулируем математическую модель. В задачах расчета морских нефтеподъемников удобно использовать уравнения их равновесия [2] в проекциях на оси декартовой системы координат Ox1 х2. Учитывая соотношения (2) и связь между независимыми переменными ds = dx2/ cos у, имеем

T' = kN/ cos у — (q, т)/ cos у, N' = —kT/ cos у — (q, n)/ cos у,

EIk' = —N/ cos у, у' = k/ cos у, u' = — tg у.

Здесь T, N — продольные и перерезывающие усилия, EI = const — изгибная жесткость, и — смещение осевой линии в горизонтальном направлении, (■)' = d(-)/dx (x = x2) .

Подстановка вектора внешних нагрузок (g) и использование уравнения (7) дает следующую систему:

T' = ш + kN/ cos ш,

N' = — (m + mc — mw)g tg у — k/ cos у[Г — pFo — (pf (1 — ш^2 + шpg(vg + vf )2)F0+

(15)

+mwg(H — x) — mwv^ sin2 у] + 0.5cnpwD0vc | vc cos у |,

EIk' = —N/ cos у, у' = k/ cos у, и' = — tg у, р' + v'f (pcvf + шpgvg) + pcg = 0.

Выберем в качестве характерных масштабов длины, погонного веса, площади и скорости соответственно величины H [м], w = pt(Fw — F0)g [kH/м], F0 [м2], v [м/c], (v = 1 м/с). После перехода к безразмерным переменным и параметрам

x0 = x/H, T0 = T/(wH), N0 = N/(wH/j,), k0 = kH, u0 = u/H, v0{) = v()/v,

V = л/EI/wH3, Y1 =0.01pggH, Y2 = pfd0/(pt x (D0 — d2)^ Y3 = pwD2/(pt(D2 — d0)),

74 = 0.01pt gH (d2 — d2 )/d0, Y5 = Y2 (v‘2/(gH)), Y6 = (pg/pt)lb, Y7 = Y3 (v‘2/(gH)),

7s = 0.5cnpwD0v2/w, 79 = 0.01 pgv2, jw = (pf /pg)Г, Y11 = 74а0 x (Г + 1),

р, = 710 /711, .р0 = 74р0, в0 = vg/vf 0,

(опуская верхний индекс ноль в обозначениях безразмерных величин), система уравнений (15) с учетом (13) запишется окончательно в виде

T' = 1 + jkN/ cos у, jN' = —(1 + 7цшр + 72(1 — ш) — 7з) tg у — k/ cos у x [T + 73(1 — x) — р—

— v20(75(1 — ш)-1 + 7бш((1 — ш)-1 + в0)2) — 77v2 sin2 у] + Ysvc | vc cos у |, jk' = —N/ cos у, у' = k/ cos у, и' = — tg у,

' Y1 шр + 72(1 — ш) «0 (16)

р' = -----^2------------------------------------------, в1 = ------------,

p2(1+0о(1°_ш)2) [ш(79р(1 — ш)-1 + Y6e0) + y5(1 — ш)] — 1 74р

ao = [710 — 711р, 0 <р < р-, ш Л 1+e°+ft+V(21+e°_ei)2+4e.! 0 <р < p*,

I0, P>P,, [0, р>р,.

При известном давлении гидросмеси на входе р0 и при шарнирном закреплении граничных сечений райзера краевыми условиями для системы (16) являются:

k(0) = u(0) = 0, р(0) = р0, k(1) = u(1) = 0, T = T1. (17)

А.Е. Федоров и др. О моделировании кожи человека

Здесь Т > 1 — заданное значение тягового усилия на платформе.

Таким образом, предложенная здесь математическая модель глубоководного райзера, взаимодействующего с внутренним двухфазным потоком газожидкостной смеси (16)-(17), относится к классу связанных задач нелинейной гидроупругости сингулярно возмущенного типа (так как при Н > 1000 м значение параметра ^ в (16) ^ < 10-3). Ее асимптотическое решение можно построить на основе использования алгоритма, предложенного в [3].

Библиографический список

1. Лойцанский Л.Г. Механика жидкости и газа. М.: Наука, 1970. 904 с.

2. Кузнецов В.В. Моделирование статических характеристик глубоководных трубопроводов морских гидротехнических комплексов. Ч.1. Формулировка уравнений установившегося пространственного движения //

Изв. вузов. Строительство. 1994. № 3. С. 69-72.

3. Кузнецов В.В. Моделирование статических характеристик глубоководных трубопроводов морских гидротехнических комплексов. Ч.2. Асимптотический анализ модельных уравнений // Изв. вузов. Строительство. 1994. № 4. С. 57-60.

УДК 501.1

О МОДЕЛИРОВАНИИ КОЖИ ЧЕЛОВЕКА

А.Е. Федоров*, В.А. Самарцев** , А.А. Адамов***

* Пермский государственный технический университет, кафедра теоретической механики,

** Пермская государственная медицинская академия им. акад. Е.А. Вагнера, кафедра общей хирургии,

***Институт механики сплошных сред Уральского отделения Российской академии наук,

E-mail: * spider@perm.raid.ru, *** adamov@icmm.ru

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

On Modelling of the Human Skin

A.E. Fedorov, V.A. Samartsev, A.A. Adamov

The mechanical properties of the skin at small and finite deformations are investigated. The models are suggested approximating viscoelastic properties of the skin. The constitutive relations forthree-dimensional stress-strain state are developed. Models parameters are determined by experiments in uniaxial tension of the skin specimens.

ВВЕДЕНИЕ

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

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

Первые исследования кожи с точки зрения ее механических свойств были сделаны Дюпюитреном (Эириу^еп) [1] в 1831 году во время расследования самоубийства в Париже. Он провел эксперименты на трупе и доказал, что раны на коже, сделанные конусообразным предметом, становятся похожими на линии и направления этих линий в различных частях тела отличаются.

Первая статья, посвященная анатомическому обобщению линий рассечений на коже, была написана Лангером (Langer) [2] в 1861 году. Им были обнаружены линии натяжения кожи (линии Лангера), а также был установлен факт, что кожа на теле человека имеет некоторое начальное натяжение.

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

© А.Е. Федоров, В.А. Самарцев, А.А. Адамов, 2007