УДК 556.536:681.3
АНАЛИЗ ПРОЦЕССА ИДЕНТИФИКАЦИИ ПАРАМЕТРОВ МОРФОМЕТРИИ РЕК И ВОДОХРАНИЛИЩ В КОМПЬЮТЕРНЫХ МОДЕЛЯХ ГИДРАВЛИЧЕСКИХ РАСЧЕТОВ (РЕЗУЛЬТАТИВНОСТЬ И ОДНОЗНАЧНОСТЬ КАЛИБРОВКИ)
© 2015 Л.К. Левит-Гуревич, Д.А. Никифоров
Институт водных проблем РАН, г. Москва
Поступила в редакцию 25.05.2015
Компьютерные модели гидравлических расчетов водохранилищ и участков рек предназначены для постоянной по запросам эксплуатации в инженерных и управленческих целях, оценки гидрологических явлений и пр. При подготовке компьютерных гидравлических моделей возникает необходимость их калибровки - корректировки первоначально заданных параметров морфометрии водного объекта с целью приближения результатов расчетов, полученных в определенных гидрологических условиях, к наблюденным фактически данным при этих условиях - уровням на створах наблюдений, а также к установленным характеристикам водохранилищ. Рассмотрен процесс калибровки по шагам. Доказано, что процесс калибровки однозначно приводит к искомой цифровой модели изучаемого участка реки.
Ключевые слова: река, гидравлический расчет, компьютерная программа, цифровая гидравлическая модель, калибровка
Компьютерной моделью гидравлических расчетов участка реки с водохранилищем называется структурированный набор цифровых данных в соответствии с используемой программой расчетов (HEC-RAS, MIKE и др.). Модель водного объекта является как бы эталонной формой заполнения гидрологических данных, что позволяет осуществлять многократные гидравлические расчеты с решением задач определенного вида, проводя замену данных в разных вариантах. Гидравлический расчет ведется по створам, нумерация створов i при подготовке данных идет от устья основного участка реки (водохранилища) против течения и отдельно от устья каждого притока, фиксируются расстояния между створами. Поперечные сечения каждого створа задаются в ориентации течения воды от условной топографической точки с высотной отметкой Zo - самой левой на сечении створа, каждая последующая топографическая точка сечения определяется отметкой и ее расстоянием от самой левой точки сечения. Отмечаются бровки русла, т.е. границы русла и пойм, задаются коэффициенты шероховатости nf n chanl, nright русла и пойм. Перечисленные данные являются базовой геометрической информацией цифровой гидравлической модели и используются во всех расчетах. Информация по гидрологии для изучаемого участка реки - это оперативные данные, предназначенные для конкретных расчетов. Таков состав данных цифровой гидравлической модели.
При расчете установившегося движения воды гидрологическая информация содержит
Левит-Гуревич Леонид Константинович, кандидат технических наук, старший научный сотрудник. E-mail: [email protected]
Никифоров Дмитрий Андреевич, младший научный сотрудник. E-mail: [email protected]
расходы Qi в створах и уровень воды в начале участка ; гидрологическая информация при расчете неустановившегося движения содержит гидрографы притоков Qi() , граничные условия в виде хода уровней ), - как правило, ) , - и
начальные условия Qi (0) (неустановившийся расчет требует также задания общего периода и интервала расчета ). В результате гидравлического расчета определяются уровни во всех створах, постоянные значения Zi при установившемся движении и ход уровней ) на каждом створе при неустановившемся движении. Но цифровая модель участка реки является лишь идеализацией гидравлических условий, не отражая всех гидравлических тонкостей реального морфологического строения реки в принципе.
Рис. 1. Схема расчетного гидравлического у частка с нумерацией створов
Программами гидравлических расчетов не учитываются многорукавность, косоструйность, изменения шероховатости в руслах и пойме между соседними створами и многое другое. И эта идеализация принципиальна. Исправление топографии поперечных сечений и коэффициентов шероховатости проводится корректировкой параметров
сечений с поверкой гидравлическими расчетами, т.е. путем калибровки базовой цифровой гидравлической модели. Принцип калибровки гидравлической модели следующий: это процесс, включающий этапы гидравлического расчета, сравнения полученных после расчета уровней с фактическими данными на постах наблюдений и оценка различий между ними, устранения этих различий путем изменения параметров сечений. Критерием калибровки является максимальное приближение полученных по расчетам уровней и объемов водохранилищ к фактическим уровням реки на створах - постах наблюдений в соответствующих гидрологических условиях и к выверенным объемам водохранилища [1,2]. Необходимость калибровки подчеркивали многие исследователи [3, 4], но ее формализация по сути началась с работ [5, 6], где даны общие принципы калибровки.
Калибровка является задачей, обратной гидравлическому расчету. Обратная задача по отношению к прямой задаче решения дифференциальных уравнений (для гидравлики это уравнения Сен-Венана и Шези [7, 8]), в принципе, некорректна по
Итак, калибровка цифровой модели участка реки с водохранилищем для гидравлических расчетов представляет собой процесс корректировок параметров модели (ширины русла и пойм, глубины, шероховатости, расстояний). При этом на каждом шаге калибровки выбираются перспективные
створы I и параметры х, их поперечных сечений, изменение которых целесообразно по сравнению с другими параметрами в смысле быстрого достижения минимума оценки близости рассчитанных и фактических уровней на створах наблюдений, как и близости к оценке объемов водохранилища, -функция F . Выбор осуществляется путем вычисления определенных оценок dF> 0 целесообразности корректировок. В этой статье процесс калибровки не рассматривается подробно, но обозначены шаги калибровки и доказывается результативность и конечность этого процесса. Кратко постановку задачи калибровки можно сформулировать следующим образом. Функция ^, оценивающая близость уровней
Fz )2/
Тихонову [9], и необходимо проводить регуляризацию задачи. Регуляризация задачи калибровки заключается во включении в ее формулировку (помимо требований максимального приближения расчетных уровней и объемов водохранилища к фактическим уровням и проектным объемам) условий обеспечения наименьшего отклонения исправленных, полученных в результате калибровки параметров описания поперечных сечений от первоначально заданных. Допустимо, однако, не включать эти условия непосредственно в формулировку задачи калибровки, а ограничиться требованием лишь малых изменений параметров модели при шагах калибровки, что равнозначно требованию небольших отклонений от начальных данных.
На профиле поперечного сечения по данным высотным отметкам земли можно проследить общее строение пойм и русла. Параметры схематизированных сечений: Н - глубина русла и пойм, Ь -ширина потока, 2 - отметки дна, В -ширина реки (рис. 2).
Fz (г) = (г,); УI(г,) = (1-1?*)2 /дг,
-а)
где - определенные величины, присущие каждому створу, позволяющие привести члены указанных сумм с разницами уровней на рассматриваемых постах наблюдений к сравнимым друг с другом величинам. Здесь I - створы постов наблюдений. При этом уровни на створах связаны друг с другом функцией гидравлического расчета
Н(г, г,.....г ), и каждый уровень зависит от
используемых параметров X(х,X,...х ), в принципе от всех параметров. Под параметрами X, ] = 1, т понимаем указанные выше ширину и
глубину русла, ширину пойм, коэффициенты шероховатости, соответствующие общему строению поперечных сечений. В зависимости от размеров
при большом, может быть, влиянии на уровень г, изменения некоторого параметра на -м створе влияние его на следующих створах уменьшается, зачастую уровень соседнего створа уже не
испытывает сильного влияния параметра -го створа.
Рис. 2. Схема поперечного сечений русла и пойм, левой и правой
Для водохранилища, входящего в состав рассматриваемого участка реки, задаются отметка нормативных уровней: нормальный подпорный 2погт (НПУ), уровень мертвого объема 2Леай (УМО), форсированный подпорный 2?огсе (ФПУ), и данные проектных или фактических выверенных
объемов водохранилища: полный упогт_и мертвый уйеай_ № объемы, полезный как разница полного и мертвого объемов. Обычно известны расстояния от створа плотины до выклинивания подпора водохранилища при разных уровнях у гидроузла, достаточно для НПУ и УМО.
пграсчет тгфакт \2 пграсчет тгфакт \2 Р _ У^УМО _ ' УМО ) VНПУ _ ' НПУ )
пграсчет Т/фа УФПУ _ ''ътг
ФПУ
)2
тгфакт УУМО
тгфакт * НПУ
тл факт * ФПУ
(2)
Т г расчет „ ^ - тгфакт -
где уумо - рассчитанный объем водохранилища при уровне мертвого объема; уумо - объем водохранилища при уровне мертвого объема фактический или принятый по проектным данным, аналогично для уровней НПУ и ФПУ.
Таким образом, калибровка модели гидравлических расчетов для участка реки формализуется в виде задачи минимизации функции, которая вычисляется в результате гидравлического расчета. Осуществляются вычислительные эксперименты, состоящие из гидравлических расчетов и вычислений Е - суммы квадратов относительных отклонений вычисленных уровней от фактических
на створах наблюдений Р2 плюс сумма относительных отклонений вычисленных объемов водохранилищ от заданных р . Выбираются параметры для корректировки модели, вносятся исправления, проводится новый гидравлический расчет.
Если в задачу непосредственно включаются
оценки изменения первоначально заданных ис-
*
ходных данных Х] в связи с регуляризацией, о чем говорилось выше, и ставится требование минимума отклонений откалиброванных топографических данных от заданных первоначально по расчетным створам параметров, то оценки близости этих параметров базовой гидравлической модели от исходных можно выписать как
2
/ _ тут х _ д.*)/х*] . Полная оценка качества
калибровки базовой цифровой модели реки и водохранилища для гидравлических расчетов в этом случае запишется:
Р _ Pz х Р + Ру х Ру + Р/ х /
где р - весовые коэффициенты, pz + ру + pf _ 1.
Будем рассматривать общий случай формализации калибровки базовой гидравлической модели для неустановившегося движения воды. Формализация калибровки модели при установившемся движении та же, но без временного фактора I. В частности, можно для определенности ориентироваться на гидравлические расчеты при оперативном управлении водохранилищами, тогда I - сутки. При установившимся движении воды калибровка ведется по уровням 2 в створах наблюдений, полученных в результате гидравлического расчета минимизируются разница (разбежка) рассчитанных и фактических уровней на створах
наблюдений. Калибровка при неустановившимся движении минимизирует среднеквадратические разбежки 2 (I) в расчетные интервалы времени или разницу максимума и минимума отклонений.
Итак, X(х,х2,х3,..,хт) - калибруемые параметры: коэффициент шероховатости, ширина и глубина русла на каждом створе и др. Оперируем полным списком параметров, не разделенных по створам, поскольку изменение каждого в отдельности или группы параметров влияют на уровни пропускаемых в реке расходов в принципе на всех выше и нижележащих створах. Назовем однонаправленным такое изменение параметров поперечного сечения, что изменение каждого из них ведет к одинаковому, одному и тому же изменению уровня на створе, повышению или уменьшению уровня.
По «физическим» соображениям всякое однонаправленное пропорциональное изменение параметров по некоторому векторному направлению в пространстве параметров X(х,х2,х3,..,х ) влечет за собой монотонное направленное только верх или вниз изменение уровней Zi в принципе на всех створах I, - не только на створах, где изменялись параметры. Возможно, пропорциональность изменения параметров и изменений при этом уровней как-то связаны, но лишь «физические» соображения не дают точных соотношений между ними. Следует заметить, что качественная картина изменения уровней и других гидравлических характеристик потока при оговоренных изменениях параметров одинакова при установившимся и неустановившимся движении. Заметим также, что локальный качественный анализ изменений уровня на одном створе можно оценочно проводить по формуле Шези. Таким образом, оперируем параметрами калибровки X (будем далее обозначать просто X), уровнями Z , которые меняются монотонно при однонаправленном изменении X и оценками качества калибровки Е -это строго выпуклая функция как сумма квадратов. При калибровке можно проводить разнонаправленные изменения параметров на створах, но будем рассматривать процесс как последовательные шаги однонаправленных изменений.
Исследование функции Р оценки калибровки гидравлической модели следует начать с обоснования квазивыпуклости функции Р как функции от параметров X . Понятие квазивыпуклости функции представлено в Википедии. Возьмем определение: минимизируемая функция квазивыпуклая тогда и только тогда, когда для любого допустимого её значения область аргумента с меньшими значениям функции выпуклая, т.е.: Р (X) квазивыпуклая, если ур * множество
{X} : Р(X) < Р* выпукло. Выпуклая функция -
квазивыпуклая; квазивыпуклая функция не обязательно выпуклая. Рассмотрим сначала функции одного переменного.
Утверждение: если функция Р(2) выпуклая, а 2 = 2 (X) монотонная, то Р(X) = Р(2 (X)) -квазивыпуклая функция.
Доказательство. Пусть Хр и Хе такие, что
F (X,) = F *
и
F(Хе) = F . Возьмем точку Хрд Х„ » X,
между двумя точками хр и :
Xp < XPQ < Xg . Поскольку 2 = 2(X) монотонная функция от X, то гр (Хр) < гре (Хре ) < ге (Хе ) , - не важно < или
>. Найдутся р > 0 и q > 0, р + q = 1, такие, что ¿рд (ХРв) = р х гр (Хр) + q х ^ (Хе) Из выпуклости Р(2) имеем:
F (Zpq) = F[ p х 2Р (XP) + q x 2q (XQ)] < p x F (2P (XP)) + q x F (ZQ (XQ)) =
= p x F (Xp) + q x F (Xg) = p x F * + q x F * = F *
(3)
т.е. F(Xpq)=F(2(Xpq)) <F* при F(X,)=F* f(Xq) = F, Xp < Xpq < Xq,
Рассмотрим теперь многомерную функцию
F=F г , г2, г3,.., гп) от набора уровней
г = {г, г, г,-, г } по створам , = 1, т, так что уровни зависят от множества параметров Х = {х, Х2, Х3,.., Хт }, и уровень на каждом створе
есть функция ^ = Zi (х, X, X,-, Хя). В этом случае
F = Fг ) = ЩХ1>Iг2(^х2,..хтI,г(^х2,..хт)] =
= F(Х2,"Хт )
Анализ поведения функции F(Х), исходя из свойств F=F (г) как квадратичной функции и монотонности г , = г ,(Х) по каждому параметру
Х] в отдельности показывает, что при бесконечном увеличении каждого параметра F стремится к определенному конечному значению, большему Fmin , поэтому ближе к Fmin , начиная от некоторых
значений параметров Х] F(Х) приобретает форму замкнутой фигуры (рис. 3).
Для справедливости утверждения о квазивыпуклости функции F(х1;х2,..хт) аналогично одномерному случаю недостаточно условие простой монотонности функций г I = (X1, X2,.., Хт), I = 1, т независимо друг от друга.
Введем понятие согласованной монотонности. При согласованной монотонности при линейном однонаправленном изменении параметров x{х1,х2,х3,..,хт} от Х к Т, в общем случае однонаправленного на нескольких створах, функции г, (Х), ,= 1,п меняются подобным образом:
, что и есть квазивыпуклость. ♦
Z(X) = Z(X) + r xR(X),i = \n; r = const ; R(X) монотонная, R( X') = 0.
Рис. 3. Форма функции F (Х) вблизи минимума
Утверждение: При согласованной монотонности для любых вариантов параметров Т и Т с соответствующими им наборами функций г иг, набору параметров X": Т < Т < Т
(т.е. х" < х) < х" , ] = 1, т) соответствует г" : г < Т" < г, и найдутся р>0, q>0, p+q=1, такие, что г (Т") = р х г (Т) + q х г (X) единственные для всех I = 1, п.
Доказательство. Если найдутся две функции тX гп(Х"'), их можно представить как Х")=рцх Х)+qцх гц(Т), ¿¡Г) = рчх г„(Т)+%х гч(Г),
р^>q^>0> р,+q,= 1; ргq¡>0, рп+%=1, то присо-
гласованнои
монотонности
ZM(X") = ZJyX') + Гмх R(X'), Iß'") = Z/X ) + rM x R(X"'), также как Zn(T) = Z„(X) + r„x R(X),
Z„(X") = Z7(X)+r^ x R(X). Тогда для первой функции^ rRT^p^m+q^xrRT) или R(X) = p^ x R(X)+^ x R(F), также
R(X') = ^ x R(X) + q x R(X*) для второй функции, -
одно и то же уравнение с разными р и q. Но по условиям к этим коэффициентам (в начале доказательства), необходимо p^ = p^ = p, q = ^ = q .♦ Заметим, если наоборот, определить согласованную монотонность, как одно для i = 1, n, условие
Z (Xй) = p x Z (X)+q x Zt (X) й
i i i , доказывается подобие
изменения Zi(X) на всех створах.
Согласованная монотонность интуитивно представляется из «физических» соображений: при однонаправленном изменении некоторых параметров меняются уровни потока, но выдерживается определенный порядок этих изменений. Невозможно, чтобы в условиях спокойного потока в одном створе уровень повышался, потом понижался. Превалирующее большинство рек имеют спокойное течение и в половодье. При бурном и только бурном потоке горной реки картина аналогична. Более сложно, если спокойный и бурный потоки перемежаются как на предгорных реках с порогами.
Переходим к основному Утверждению: Если F(Zl,Z2..,Zn) = 2^/ (Z), все fi - выпуклые, а Zi = Zi(x, x2,.., xm), i = 1, n , - согласованно монотонные функции, то при однонаправленном изменении параметров
F(X1,.., Xm ) = F[Z1(X1,.., X.. Zn ^^ Xm )] - ™а-
зивыпуклая функция.
Доказательство. Пусть Хр(xp,xp,xp,..,xpm) и
Q VQ VQ
Xq(xQ,xQ,x3Q,..,xm) такие, что F(Xp) = F
F (Xq ) = F *.
Q) = F . При
F (Xp) = 2 f ( Xp ) = 2 ff =F', F (Xq) = 2 f ( Xq ) = £ f =F *. между двумя точками
Возьмем xpq
XP
X,
Q :
XP < xPq < q .(по одним параметрам x может
быть
xp < xPQ < xQ
,-ß
по другим параметрам
функции Zi _ Zг (X) согласованно монотонные, найдутся р, Ч, Р + Ч _ 1, такие, что для всех V;' _ 1 , п функций можно записать ZPq(Xpq) _ р х ZP(X,) + ч х Zqi ). Из выпуклости / (Z) имеем тогда для всех г: / ^рЧ) _ /[ р х Zip (X,)+Ч х zq ^)] < р х / )+Ч х / (^) _ _ р х / (X,)+Ч х / (x0) _ р х /р*+Ч х /Ч*
Тогда
q\_
и
этом также
хр > хР6 > х6, указанный знак < условен, важно, что изменение параметров однонаправленное, т.е. изменение каждого параметра влияют на уровни одинаково). Поскольку каждая
^ _ zi (X) _ Zi (х, х ,.хи) монотонная функция по каждому из хг, то по однонаправленной совокупности хх2..,хт} Zp(X,) <ZIQ^) <Ze(Xe), -или Z/> (X]?) > ZPg (XPg) > ^ (Xg), не имеет значения, что рассматривать. Поскольку кроме того все
р Щ _ 2/ (zpq) < р х£ / (zp)+ч х 2г/ (zгq)
/ + Чх^/р _рхР + чхР _Р
Т.е. Р (X рч) _ Р (Z (Xpq)) < Р * при Р (X,) _ Р
Р) _ Р , ^ < < Xe, - это квазивыпуклость.
Таким образом, когда «текущие» перед следующим шагом процесса калибровки параметры X*{х*,х*,х*,..х*} делят все пространство параметров на октанты и в октанте, который соответствует нужному однонаправленному изменению параметров, оценочная функция Р(^) квазивыпуклая, при этом Р* _ Р(X*). Нужное изменение параметров определяется необходимостью поднять или опустить уровни ^ _ Zi (х, х, ., хт), ; _ 1, п , приближая их к фактическим уровням, заданным на постах наблюдений, а также увеличить или уменьшить расчетные значения характерных объемов водохранилища с целью приблизить их к фактическим (проектным) объемам.
Может рассматриваться один створ или некоторая совокупность створов, на которых необходимо поднять или опустить уровни, поэтому многомерная геометрическая картина проявления квазивыпуклости функции Р(X) в пространстве параметров может быть различной на разных шагах калибровки. Квазивыпуклость может проявляться в октанте в пространстве всех т параметров, и на многомерной плоскости (затрагиваются т-1 параметр), и на многомерном ребре с т-2 параметрами. Элемент многомерной геометрии, в пределах которого проявляется квазивыпуклость Р (X), зависит от числа параметров однонаправленного изменения, выбранных для уменьшения оценки Р (X).
Квазивыпуклая функция унимодальная, имеет единственный минимум. Рис. 4, на котором схематично даны изолинии равных значений Е(X) на множестве координат X(^, х2 ) , иллюстрирует шаги процесса калибровки гидравлической модели.
Выпуклое множество X: Р(X) < Р* (рис. 4) ограничено, характеризуется диаметром d и точка
и
X*{х', х', х',..х'} «текущих» параметров процесса калибровки по самой процедуре построения находится на границе множества. Шаг калибровки заключается в корректировке параметров на определенные величины Дг, ДВ, АН, Дп изменения высотных отметок, ширины и глубины русла и пойм, которые могут выражаться абсолютными значениями, при этом точность достижения решения при калибровке полагается на практике ДВ=3-10 м, Д, ДЯ=7-15 см. Изменение шероховатости Дп принимается не менее 0,0025.
F(Xе) = ¥ *. При этом ¥ (Хр) = 2 (Хр) = 2 ур =F,,
также F(X,) = 2, УЪ (Хе) = 2, УУ =F*. Возьмем Хрд
между двумя точками Хр и Х е : Хр < Хре < Хв .(по одним параметрам х может быть хр < хре < хе по другим параметрам
хр > хре > Xе ).
Поскольку каждая ^ = 2 (X) = 2 (X, X , ) монотонная функция по каждому из однонаправленно корректируемому Х], то по совокупности однонаправленно изменяемых {х1, х2.., хя} имеем гр (Хр) < г ре (Хре) < гв (Хв), или
(Хр) > ^(ХРе) > (Хд), не имеет значения. Поскольку каждая функция г 1 = (X) монотонная, найдутся р,, qi, р, + qi = 1, так, что для каждой из
них
можно
записать
Рис. 4. Изолинии равного значения Р (X) и схема действий по калибровке модели
Шаг калибровки перемещает текущую точку параметров X*{х*,х',х1*,..х1*} внутрь выпуклого
множества X: F(X) < F*, на следующем шаге
диаметр множества меньше предыдущего и так происходит до тех пор, пока с заданной точностью корректировок параметров мы не достигнем точки минимума. Процесс калибровки конечен, поскольку заданы конечные размеры изменений параметров.
Задача калибровки, таким образом, представляет собой поиск точки минимума в координатах калибровочных переменных (Xтп ) ; при этом совсем не обязательно этот минимум окажется равен нулю (Xтп ) = 0, - общий случай (Xтп ) > 0. Если (X) =0, достигается полная калибровка гидравлической модели, (X) >0 указывает на необходимость добавления новых расчетных створов на участке или/и добавления гидрологической информации.
Рассмотрим такой вопрос: в какой степени меняется свойство квазивыпуклости при отсутствии равномерной монотонности г = г (X) = г (X, X .хт), поскольку свойство равномерной монотонности основано на «физических» соображениях, что в природе может выполняться не абсолютно, а с определенной ограниченной неточностью.
г?(Хп) = р, х гр(Хр) + д1 х г*Т) и из выпуклости У (г) имеем для каждого г:
У (гГ) = У [ Л х гр (Хр) + q¡ х % (Хе)] < р, х У г) + q¡ х У ) =
= р, х У (Хр) + q¡ х У (Хе) = р, х У р + q¡ х Тогда
F (Хп) = 2 У (гл < 2[р, х У г) + г х У г )]=2[р, х + q¡ х у' ]=
= 2[( р + р) х Г,"' + + Ч) х ]=
= рх2 Ур + q х 2 У' +2[pУ"' + У ]=(р + q) х F' + 2[pУ"' + У ] = F' + 2[pУ"' + У ]
Здесь:
| п | п п п
р=-2р, q=12ъ , р = р - р ¡=Ъ - q, 2р =0' 2ч=0
п ,=1 п
1=1
1=1
,=1
из р + Ъ = 1, р, + ъ = 1 V,, Р+^ = 0 уи р<-2' ¡<-2. Тогда:
п п
¥(ХИ) < F * +2 А х (Ур - /■')■; Д = 2 Р. х (Ур - /У); 1=1 1=1
(9)
Итак, F(X ) < F" + А < F*, если А < 0.
Т.е. ¥ (X ) = ¥ (г (X )) < ¥" при ¥ (Хр) = ¥",
(Хв) = ¥ , Д < 0, - квазивыпуклость соблюдается.
При А > 0 необходим дополнительный анализ. Заметим сначала, что Д - небольшое число в сравнении с ¥ , так как содержит одни разности близких величин, а множители р ограничены: Р | < ^, и кроме
того
Швторяем, пусть Хр(х-р,Х2р,Хзр,..,х>я) и 2р = 0; 2(Ур ~У') = 0; При Д> 0 имеем Хе (хв, хв, хв,.., хЯ) такие, что ¥ (Хр) = ¥ * и '=- '=-
п
п
Р< Р + А, возможны два случая: Р^) < Р* < Р* + А и Р* < РX) < Р* + А. Первый случай означает, что выпуклое множество {X}: Р(X) < Р* существует, и это соответствует
квазивыпуклости Р(X). Во втором случае множества X}: Р^) < Р* нет, что может означать невозможность достичь результата калибровки, соответствующего заданной точности поиска, приведенной ранее: ДВ=3-10 м, Дг, ДН=7-15 см, Дп =0,0025. При такой ситуации необходимо рассмотреть большее число створов с дополнительными данными по гидрологии.
Выводы: практические расчеты по калибровке гидравлических моделей показывают, что приведенные выше ситуации имеют место, поскольку зачастую пользователь для облегчения калибровки экономит на числе створов и объеме гидрологических данных.
СПИСОК ЛИТЕРАТУРЫ:
1. Никифоров, Д.А. Калибровка гидравлических моделей рек и водохранилищ. Влияние корректировки параметров программной модели на точность гидравлических расчетов // Современные проблемы водохранилищ и их водосборов. Труды междунар. научно-практ. конф. Управление
водными ресурсами. Гидро- и геодинамические процессы. Том I. - Пермь, 2013. С. 253-260.
2. Никифоров, Д.А. Влияние корректировки параметров
программного комплекса HEC-RAS на точность гидравлических расчетов и формализация процесса калибровки гидравлических моделей рек и водохранилищ // Чистая вода России. Сб. материалов XII междун. научно-практ. симпозиума. -Екатеринбург, 2013. С.226-234.
3. Грушевский, М. С. Неустановившееся движение воды в реках и каналах. - Л.: Гидрометеоиздат,1982. 288 с.
4. Кюнж, Ж.А. Численные методы в задачах речной гидравлики: практ. применение / Ж.А. Кюнж, Ф.М. Холли, А. Вервей. - М.: Энергоиздат, 1984. 255 с.
5. Корень, В.И. Определение геометрических и гидравлических характеристик речного русла путем решения обратных задач для уравнений Сен-Венана / В.И. Корень, Л.С. Кучмент // Водные ресурсы. 1973. №4.С. 83-100.
6. Левит-Гуревич, Л.К. Математические основы калибровки компьютерных гидравлических моделей рек и водохранилищ / Л.К. Левит-Гуревич, Д.А. Никифоров // Чистая вода России. Сб. материалов XII междунар. научно-практ. симпозиума. -Екатеринбург, 2013. С. 190-198.
7. Агроскин, ИИ. Гидравлика. 4-ое изд. / ИИ. Агроскин, Г.Т. Дмитриев, Ф.И. Пикалов. - М.: Энергия, 964. 352 с.
8. Штеренлихт, Д.В. Гидравлика. Учебник для ВУЗов. 3-е изд. - М.: КолосС, 2004. 656 с.
9. Тихонов, А.Н. Методы решения некорректных задач. 2-ое изд. / А.Н. Тихонов, В.Я. Арсенин. - М.: Наука, Главная редакция физико-математической литературы, 1979. 564 с.
PROCESS IDENTIFICATION ANALYSIS OF MORPHOMETRY PARAMETERS OF RIVERS AND RESERVOIRS IN COMPUTER MODELS OF HYDRAULIC CALCULATIONS (EFFECTIVENESS AND UNIQUENESS OF CALIBRATION)
© 2015 L.K. Levit-Gurevich, D.A. Nikiforov
Institute of Water Problems RAS, Moscow
Computer models of hydraulic calculations of water reservoirs and river sections are designed for continuous operation at the request in engineering and management purposes, evaluation of hydrological phenomena and so forth. In the preparation of hydraulic computer models there is a need of calibration, adjustment to the original set parameters morphometry of the water body in order to bring the results of calculations obtained in certain hydrological conditions, according to the observations of fact, under these conditions, levels of cross-sections for observation and to established characteristics of reservoirs. The process of calibration steps. It is proved that the calibration process clearly leads to the desired digital model of the studied area of the river.
Key words: river, hydraulic calculation, computer software, digital hydraulic model calibration
Leonid Levit-Gurevich, Candidate of Technical Sciences, Senior Research Fellow. E-mail: [email protected]
Dmitriy Nikiforov, Minor Research Fellow. E-mail: [email protected]