Прогнозирование показателей разработки месторождения

55,380

Таблица 6.2 - Расчет налогов, пошлин и платежей, млн. дол.

Прогнозируемый год	1	2	3	4
Акциз	Добыча газа, млн. м3	63,480	60,300	56,730	55,380
	Стоимость по ценам внутреннего рынка	0,646	0,614	0,578	0,564
	Итого	0,097	0,092	0,087	0,085
Роялти	Выручка (всего без НДС)	0,646	0,614	0,578	0,564
	Вычеты:
	Транспортные расходы	0,086	0,081	0,077	0,075
	Акциз	0,097	0,092	0,087	0,085
	Всего	0,183	0,173	0,163	0,159
	База для расчета	0,464	0,440	0,414	0,404
	Итого	0,074	0,070	0,066	0,065
Налог на МСБ	Выручка (всего без НДС)	0,646	0,614	0,578	0,564
	Вычеты:
	Транспортные расходы	0,086	0,081	0,077	0,075
	Акциз	0,097	0,092	0,087	0,085
	Всего	0,183	0,173	0,163	0,159
	База для расчета	0,464	0,440	0,414	0,404
	Итого	0,046	0,044	0,041	0,040
Налог на имущество	0,001	0,001	0,001	0,001
Налог на содержание жилого фонда	Выручка (всего без НДС)	0,646	0,614	0,578	0,564
	Вычеты:
	Акциз	0,097	0,092	0,087	0,085
	База для расчета	0,549	0,522	0,491	0,479
	Итого	0,008	0,008	0,007	0,007
Дорожный налог	Выручка (всего без НДС)	0,646	0,614	0,578	0,564
	Вычеты:
	Акциз	0,097	0,092	0,087	0,085
	База для расчета	0,549	0,522	0,491	0,479
	Итого	0,014	0,013	0,012	0,012
Отчисления с ФЗП	ФЗП	0,020	0,020	0,020	0,020
	Итого	0,008	0,008	0,008	0,008
ВСЕГО	0,249	0,237	0,223	0.218

Таблица 6.3 - Расчет экономических показателей, млн. дол.

Прогнозируемый год	1	2	3	4	ВСЕГО
Товарная добыча газа, млн. м3	63,480	60,300	56,730	55,380	235,890
Выручка	0,646	0,614	0,578	0,564	2,401
Транспортные расходы	0,086	0,081	0,077	0,075	0,318
Акциз	0,097	0,092	0,087	0,085	0,360
Эксплуатационные расходы
Текущие затраты	0,160	0,160	0,160	0,160	0,640
Роялти	0,074	0,070	0,066	0,065	0,276
Налог на МСБ	0,046	0,044	0,041	0,040	0,172
Дорожный налог	0,014	0,013	0,012	0,012	0,051
Отчисления с ФЗП	0,008	0,008	0,008	0,008	0,032
Амортизационные отчисления	0,020	0,020	0,020	0,020	0,080
Всего	0,322	0,316	0,308	0,305	1,251
Валовая прибыль	0,141	0,125	0,106	0,099	0,471
Налог на имущество	0,001	0,001	0,001	0,001	0,004
Налог на содержание жилого фонда	0,008	0,008	0,007	0,007	0,031
Балансовая прибыль	0,132	0,116	0,098	0,091	0,437
Налог на прибыль	0,040	0,035	0,029	0,027	0,131
Прибыль после выплаты налога	0,092	0,081	0,068	0,064	0,306
Чистый дисконти-рованный доход	0,092	0,071	0,052	0,042	0,257
Доходы бюджета	0,288	0,271	0,252	0,245	1,057

7. Безопасность и экологичность проекта

7.1 Основные направления обеспечения безопасности и экологичности эксплуатации газовых скважин в НГДУ “Оханефтегаз”.

Безопасность и экологичность инженерно-технических решений, предлагаемых в данном проекте, будут оцениваться по уровню безопасности, достигнутом в НГДУ “Оханефтегаз”.

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

1) природный газ в трубопроводах и в технологическом оборудовании находятся под высоким давлением;

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

3) необходимость обслуживания установок в любых метеорологических условиях на открытых площадках, а так же в ночное время;

4) проведение на промысле работ, связанных с электрическим током и электроприборами;

5) при исследовательских работах геофизической службой на скважине имеют место радиоактивные приборы;

6) аварийные ситуации в системе “скважина - шлейф - УПГ”.

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

Рассматриваемые условия труда относятся к вредным условиям. Если не предприняты меры для нейтрализации вредно действующих факторов, то они отрицательно влияют на здоровье персонала и приводят к профессиональным заболеваниям.

Реализация проекта связана со следующими видами воздействия на окружающую природную среду:

1) изъятие: полезных ископаемых (природного газа); водных ресурсов (на технологические и бытовые нужды);

2) привнос: химических веществ; твёрдых и жидких бытовых и производственных отходов; шума и вибрации.

Основными объектами явлются недра, атмосферный воздух, поверхностные и подземные воды, почва, флора и фауна, работающий персонал.

Для утилизации выбросов газа предусмотрены системы сброса на факел и на свечу.

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

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

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

Таким образом, основные направления направления обеспечения безопасности и экологичности эксплуатации газовых скважин в НГДУ “Оханефтегаз” -- это мероприятия по обеспечению производственной безопасности, мер по ликвидации последствий ЧС и охраны окружающей среды.

7.2 Мероприятия по обеспечению производственной безопасности

Наиболее опасными на промысле являются следующие объекты: площадки сепарации газа, котельная, замерной пункт газа, парк метанольных емкостей. Все эти объекты относятся к категории производств группы А, по классу взрыво- и пожароопасности (ПЭУ) - В-1г по категории и группе взрывоопасной смеси - 2 ТЗ по ГОСТу 12.1.011-78.

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

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

Для наружного пожаротушения предусмотрен кольцевой водовод. Для внутреннего пожаротушения в помещениях категорий производства А, В предусмотрены пожарные краны. Рабочие места оснащены ящиками с песком, совковыми лопатами, вёдрами, кошмами и огнетушителями (ОП-10, ОП-50). Использование пожарного оборудования и инвентаря для хозяйственных, производственных и других нужд, не связанных с пожаротушением, запрещается.

Взрыво- и пожароопасные вещества удаляются из помещений системой вытяжной вентиляции.

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

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

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

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

Мероприятия по защите от статического электричества:

1) технологическое оборудование и трубопроводы для предотвращения опасностей, связанных с искровыми разрядами статического электричества заземлены; максимальное сопротивление контура заземления от статического электричества не должно превышать 100 Ом;

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

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

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

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

На каждого работающего предусмотрен объем производственных помещений не менее 15 м³, а площадь не менее 4,5 м², высота производственных помещений не менее 3,2 м. Полы в помещениях ровные и нескользкие.

На рабочих местах установлена достаточная освещенность, равномерное распределение яркости, отсутствие резких теней, постоянство освещенности во времени.

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

Работники снабжаются спецодеждой и спецобувью общего назначения: костюм (комбинезон) хлопчатобумажный, плащ непромокаемый, сапоги кирзовые или резиновые, рукавицы брезентовые, костюм утепленный для метеоусловий.

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

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

7.3 Обеспечение мер по ликвидации чрезвычайных ситуаций

В условиях НГДУ “Оханефтегаз” возможно возникновение следующих видов ЧС:

- техногенного характера;

- природного характера;

- военно-политического характера.

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

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

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

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

ЧС природного характера. Месторождение Шхунное находится в Сахалинской области -- в северной части острова Сахалин. Климат района морской, довольно суровый. Зима продолжительная, холодная, с частыми буранами и сильными ветрами. Наиболее низкая температура (до минус 33 ^оС) отмечается в январе. Продолжительность морозного периода от 150 до 180 дней в году. Лето короткое, дождливое с максимальной температурой от 28 до 30 ^оС в июле-августе.

Северная часть Сахалина характеризуется повышенной сейсмической активностью.

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

Для предупреждения и предотвращения ЧС на предприятии действует отдел ГО и ЧС, который решает задачи выявления потенциальных источников ЧС на территории предприятия и риск их возникновения. На основе проведенного анализа с помощью специальных методик выявляются потенциально опасные производственные объекты и на основе этого прогнозируются последствия воздействия возможных ЧС на население и подведомственные территории. Отталкиваясь от полученных результатов, осуществляется выбор, обоснование и реализация направлений деятельности обеспечения защиты населения и территории предприятия. К ним относится:

1) осуществление комплекса профилактических мероприятий по предотвращению возникновения и снижению ущерба от ЧС;

2) организация защиты населения и его жизнеобеспечения в ЧС;

3) обеспечение устойчивости работы хозяйственных объектов в ЧС;

4) организация аварийно-спасательных и других неотложных работ в очагах поражения и зонах заражения.

7.4 Мероприятия по обеспечению охраны окружающей среды

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

Для аварийного или планового задавливания скважин перед капремонтом предусматривается установка на шлейфах в 100 м от устья специального кранового узла с задвижкой для подключения задавочного агрегата.

Минимизация загрязнения атмосферного воздуха, включая шумовое воздействие, при эксплуатации объектов обеспечивается путем:

1) использования в качестве топлива природного газа в котлоагрегатах;

2) размещения агрегатов и механизмов создающих повышенный шумовой фон в отдельных изолированных помещениях;

3) покрытия источников шума и подводящих трубопроводов противошумовой изоляцией;

4) соблюдения необходимых расстояний между жилой застройкой и промышленными объектами;

5) герметизации системы сбора, подготовки и транспорта продукции на всем пути движения;

6) отвода газа при срабатывании предохранительных клапанов и полного сжигание газа на свече факела.

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

Большинство рек, попадающих в зону влияния проекта, имеют протяженность до 200 км и квалифицируются как малые. Поэтому в соответствии с Постановлением Совета Министров РСФСР от 14.01.81 "Об усилении охраны малых рек РСФСР от загрязнения, засорения и истощения и о рациональном использовании их водных ресурсов" установлены специальные меры по охране малых рек с целью восстановления их водности и чистоты.

Так, в водоохранных зонах при производстве работ не допускается загрязнение поверхности земли, в частности, свалка мусора и отходов производства, а также площадки стоянок автотранспорта, ремонт строительных машин и механизмов и транспортных средств, размещение складов ГСМ и химреагентов.

Для предотвращения загрязнения и изменения режима водных объектов проводятся следующие работы:

1) очистка русла рек и пойм от загромождающих их предметов;

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

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

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

5) производится постоянный контроль за степенью загрязнения территории промысла и прилегающих водотоков.

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

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

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

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

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

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

2) исключение устройства временных опор и подмостей в русле рек;

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

Работы вблизи рек выполняются в строгом соответствии со сроками и рекомендациями специалистов СахНИРО при согласовании с рыбной инспекцией и областным Комитетом по охране природы.

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

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

При производстве работ в лесах проводятся противопожарных мероприятий.

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

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

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

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

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

1) профилактику сверхнормативного экологического ущерба;

2) формирование информационной базы состояния и изменений окружающей природной среды.

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

Комплекс работ по мониторингу включает в себя выполнение мероприятий по охране недр, систематический отбор и анализ проб поверхностных и подземных вод, донных грунтов и бентоса, почв, воздуха на объектах месторождения.

В данном разделе были рассмотрены мероприятия, принятые в НГДУ ”Оханефтегаз” по обеспечению пожаро- и взрывобезопасности, электробезопасности, комфортных условий труда работников, а также мероприятия по охране недр, атмосферного воздуха, поверхностных и подземных вод, животного мира, почвы и растительности.

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

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

Заключение

В заключение подытожим полученные результаты.

Итак, рассмотрено газовое месторождение Шхунное. По фактическим данным разработки уточнены режим работы залежи и запасы газа (по методу падения пластового давления). Режим залежи -- газовый; уточненные запасы оцениваются в 2788 млн. м³, что на 118 млн. м³ больше подсчитанных объемным методом.

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

На данной модели был воспроизведен процесс разработки за пять лет, причем задавались отборы по скважинам и “введение” их в разработку по фактическим данным. Относительная погрешность расхождения промысловых и расчетных данных не превышала 2,91 %, из чего был сделан вывод о пригодности данной модели для прогнозных расчетов. Следует отметить, что расчеты производились с учетом взаимодействия двух пластов залежи в тех скважинах, где они вскрыты перфорацией совместно.

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

Прогнозные расчеты выполнены на четыре года. По результатам расчетов за прогнозируемый период добыча газа составит 235,89 млн. м³, пластовое давление снизится с 5,86 до 5,19 МПа.

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

Технико-экономические показатели разработки были рассчитаны на прогнозируемый период. За этот период разработка месторождения будет характеризоваться следующими показателями:

- эксплуатационные расходы - 0,64 млн. дол.;

- поступления государству - 1,057 млн. дол.;

- средняя себестоимость добычи газа - 5,314 дол./тыс. м³;

- чистый доход - 0,306 млн. дол.;

- дисконтированный доход - 0,223 млн. дол.

Безопасность и экологичность инженерно-технических решений, предлагаемых в проекте, оценивались по уровню безопасности в НГДУ “Оханефтегаз”, где достигнут в данном направлении высокий уровень, что приводит к выводу о безопасности и экологичности данного проекта.

Список использованных источников

1. Азиз Х.,Сеттари Э. Математическое моделирование пластовых систем. - М.: Недра, 1982. - 407 с.

2. Закиров С.Н. Разработка газовых, газоконденсатных и нефтегазоконденсатных месторождений. - М.: Струна, 1998. - 628 с.

3. Многомерная и многокомпонентная фильтрация: Справочное пособие /С.Н. Закиров, Б.Е. Сомов, В.Я. Гордон и др. - М.: Недра, 1988. - 335 с.

4. Прогнозирование и регулирование разработки газовых месторождений /С.Н. Закиров, В.И. Васильев, А.И. Гутников и др. - М.: Недра, 1984. - 295 с.

5. Кричлоу Генри Б. Современная разработка нефтяных месторождений - проблемы моделирования. - М.: Недра, 1977. - 303 с.

6. Технологический режим работы газовых скважин /З.С. Алиев, С.А. Андреев, А.П. Власенко и др. - М.: Недра, 1978. - 279 с.

Приложение А

(обязательное)

Программа решения дифференциального уравнения неустановившейся фильтрации газа в неоднородном пласте

{$M 65520,0,655360}

program both;

uses Crt;

function z(Per:real):real;

begin

z:=0.000987437*Per*Per-0.02530907*Per+1.002522429;

end;

function power(x,st:real):real;

begin

if st=0 then power:=1

else power:=exp(st*ln(x));

end;

TYPE

tipR=real;

m=array[1..19,1..34]of real;

m1=array[1..16,1..34]of real;

procedure sip (i,j:integer;g,c,a,b,f,R,g_1,f_1,alfa,v2,v3:tipR;var ggg,fff,v1:tipR);

var

gg,cc,aa,bb,ff:tipR;

begin

bb:=f/(1+alfa*f_1);

cc:=c/(1+alfa*g_1);

aa:=a+alfa*(cc*g_1+bb*f_1)-cc*f_1-bb*g_1;

gg:=(g-alfa*cc*g_1)/aa;

ff:=(b-alfa*bb*f_1)/aa;

ggg:=gg;

fff:=ff;

v1:=(R-cc*v2-bb*v3)/aa;

end;

procedure sipp1(i,j:integer;g,c,a,b,f,R,b_1,f_1,alfa,v2,v3:tipR;var bbb,fff,v1:tipR);

var

gg,cc,aa,bb,ff:tipR;

begin

gg:=g/(1+alfa*b_1);

cc:=c/(1+alfa*f_1);

aa:=a+alfa*(gg*b_1+cc*f_1)-gg*f_1-cc*b_1;

bb:=(b-alfa*gg*b_1)/aa;

ff:=(f-alfa*cc*f_1)/aa;

bbb:=bb;

fff:=ff;

v1:=(R-cc*v2-gg*v3)/aa;

end;

procedure VII(giv:m;var gi:m;dt,t:real; var PcpVII,p5_2,p8_2,p10_2,p16_2,p75_2,p11,p15:real;

q5_1,q5_2,q8_1,q8_2,q10_1,q10_2,q16_1,q16_2,q75_1,q75_2,q11,q15:real)

const

N=19;{strok}

M=34;{stolbcov}

Ppl=7.6e6;{Pa}

k1=0.661e-12;{mkm2, IV block}

k2=0.44e-12;{mkm2, III block}

kp=0.32;

al1=0.77;

al2=0.76;

mu=0.011e-3;{Pa*sec}

dx=100;{m}

dy=100;{m}

pat=100000;

day=86400;

type

tipR=real;

mas1=array[1..20]of tipR;

mas2=array[1..N,1..M]of tipR;

dmas=^mas2;

VAR

dat,dat1:text;

al:mas1;

i,j,count,k,cell:integer;

qm,p:mas2;

g,c,a,b,f,d,R,sum,ZapN,ZapT,Qdob,Pcp,delta,Hcp,ggg,fff,bbb,v1:tipR;

b1,g1,f1,kh,V,dp,pp:dmas;

flag:boolean;

procedure alf(var al:mas1);

var

Продолжение приложения А

i,j:integer;

par,par1,par11,par2:tipR;

begin

par11:=100;

for i:=1 to N-1 do

for j:=1 to M-1 do

if h[i,j]<>0 then

if (h[i,j+1]<>0)and(h[i+1,j]<>0) then

begin

par:=dx*dx*(kh^[i,j+1]+kh^[i,j])/(kh^[i+1,j]+kh^[i,j])/dy/dy;

par1:=pi*pi/(2*N*N*(1+par));

par2:=pi*pi/(2*M*M*(1+1/par));

if par1<par2 then

begin

if par1<par11 then par11:=par1;

end

else

begin

if par2<par11 then par11:=par2;

end;

end;

al[1]:=1;

for i:=2 to 10 do al[i]:=1-power(par11*10,(i-1)/10);

end;

procedure skv;

begin

{5}

qm[10,13]:=q5_2/86.4;

{75}

qm[9,5]:=q75_2/86.4;

{10}

qm[6,15]:=q10_2/86.4;

{11}

qm[15,11]:=q11/86.4;

{8}

qm[10,26]:=q8_2/86.4;

{16}

qm[7,29]:=q16_2/86.4;

{15}

qm[14,23]:=q15/86.4;

end;

BEGIN

GetMem(kh,M*N*SizeOf(tipR));

GetMem(b1,M*N*SizeOf(tipR));

GetMem(pp,M*N*SizeOf(tipR));

GetMem(dp,M*N*SizeOf(tipR));

GetMem(V,M*N*SizeOf(tipR));

GetMem(g1,M*N*SizeOf(tipR));

GetMem(f1,M*N*SizeOf(tipR));

ZapN:=0;

for i:=1 to N do

begin

for j:=1 to M do

begin

if h[i,j]<>0 then

begin

if j<=21 then kh^[i,j]:=k2*h[i,j]

else kh^[i,j]:=k1*h[i,j];

p[i,j]:=giv[i,j];

ZapN:=ZapN+293/297*h[i,j]*dx*dy*kp*al1*Ppl/pat/z(Ppl/1e6)*z(pat/1e6);

end

else

begin

kh^[i,j]:=0;

p[i,j]:=0;

end;

pp^[i,j]:=p[i,j];

qm[i,j]:=0;

b1^[i,j]:=0;

g1^[i,j]:=0;

f1^[i,j]:=0;

V^[i,j]:=0;

dp^[i,j]:=0;

end;

end;

alf(al);

sum:=0;delta:=0;

count:=0;

k:=0;

skv;

repeat

count:=count+1;

if odd(count) then

begiN {obratnoe yporyadochivanie}

Продолжение приложения А

for i:=2 to N-1 do

for j:=M-1 downto 2 do

begin

g:=0;c:=0;b:=0;f:=0;a:=0;d:=0;R:=0;

if h[i,j]<>0 then

begin

if h[i,j-1]=0 then g:=0

else g:=(dx/dy)*pp^[i,j]*(kh^[i,j-1]+kh^[i,j])/mu/z(pp^[i,j]/1e6)/2;

if h[i-1,j]=0 then c:=0

else c:=(dy/dx)*pp^[i,j]*(kh^[i-1,j]+kh^[i,j])/mu/z(pp^[i,j]/1e6)/2;

if h[i+1,j]=0 then b:=0

else b:=(dy/dx)*pp^[i,j]*(kh^[i+1,j]+kh^[i,j])/mu/z(pp^[i,j]/1e6)/2;

if h[i,j+1]=0 then f:=0

else f:=(dx/dy)*pp^[i,j]*(kh^[i,j+1]+kh^[i,j])/mu/z(pp^[i,j]/1e6)/2;

a:=-(g+c+b+f)-dx*dy/dt*2*al1*kp*h[i,j]/z(pp^[i,j]/1e6);

d:=2*pat*qm[i,j]*297/293-dx*dy/dt*2*al1*kp*h[i,j]*p[i,j]/z(p[i,j]/1e6);

R:=d-(g*pp^[i,j-1]+c*pp^[i-1,j]+a*pp^[i,j]+b*pp^[i+1,j]+f*pp^[i,j+1]);

if odd(count) then k:=k+1;

if k>10 then k:=1;

sipp(i,j,g,c,a,b,f,R,g1^[i-1,j],f1^[i,j+1],al[k],V^[i-1,j],V^[i,j+1],ggg,fff,v1);

V^[i,j]:=v1;

g1^[i,j]:=ggg;

f1^[i,j]:=fff;

end;

end;

flag:=true;

for i:=N-1 downto 2 do

for j:=2 to M-1 do

if h[i,j]<>0 then

begin

dp^[i,j]:=V^[i,j]-g1^[i,j]*dp^[i,j-1]-f1^[i,j]*dp^[i+1,j];

if abs(dp^[i,j])>1e-4 then flag:=false;

pp^[i,j]:=pp^[i,j]+dp^[i,j];

end;

{} for i:=1 to N do

for j:=1 to M do begin g1^[i,j]:=0;f1^[i,j]:=0 end;

enD

eLSe

begiN

for i:=2 to N-1 do

for j:=2 to M-1 do

begin

g:=0;c:=0;b:=0;f:=0;a:=0;d:=0;R:=0;

if h[i,j]<>0 then

begin

if h[i,j-1]=0 then g:=0

else g:=(dx/dy)*pp^[i,j]*(kh^[i,j-1]+kh^[i,j])/mu/z(pp^[i,j]/1e6)/2;

if h[i-1,j]=0 then c:=0

else c:=(dy/dx)*pp^[i,j]*(kh^[i-1,j]+kh^[i,j])/mu/z(pp^[i,j]/1e6)/2;

if h[i+1,j]=0 then b:=0

else b:=(dy/dx)*pp^[i,j]*(kh^[i+1,j]+kh^[i,j])/mu/z(pp^[i,j]/1e6)/2;

if h[i,j+1]=0 then f:=0

else f:=(dx/dy)*pp^[i,j]*(kh^[i,j+1]+kh^[i,j])/mu/z(pp^[i,j]/1e6)/2;

a:=-(g+c+b+f)-dx*dy/dt*2*al1*kp*h[i,j]/z(pp^[i,j]/1e6);

d:=2*pat*qm[i,j]*297/293-dx*dy/dt*2*al1*kp*h[i,j]*p[i,j]/z(p[i,j]/1e6);

R:=d-(g*pp^[i,j-1]+c*pp^[i-1,j]+a*pp^[i,j]+b*pp^[i+1,j]+f*pp^[i,j+1]);

if odd(count) then k:=k+1;

if k>10 then k:=1;

sipp1(i,j,g,c,a,b,f,R,b1^[i,j-1],f1^[i-1,j],al[k],V^[i-1,j],V^[i,j-1],bbb,fff,v1);

V^[i,j]:=v1;

b1^[i,j]:=bbb;

f1^[i,j]:=fff;

end;

end;

flag:=true;

for i:=N-1 downto 2 do

for j:=M-1 downto 2 do

if h[i,j]<>0 then

begin

dp^[i,j]:=V^[i,j]-b1^[i,j]*dp^[i+1,j]-f1^[i,j]*dp^[i,j+1];

if abs(dp^[i,j])>1e-4 then flag:=false;

pp^[i,j]:=pp^[i,j]+dp^[i,j];

end;

{} for i:=1 to N do

for j:=1 to M do begin b1^[i,j]:=0;f1^[i,j]:=0 end;

enD;

until flag;

ZapT:=0;Pcp:=0;Hcp:=0;

for i:=1 to N do

begin

for j:=1 to M do

begin

gi[i,j]:=pp^[i,j];

ZapT:=ZapT+293/297*h[i,j]*dx*dy*kp*al1*p[i,j]/pat/z(p[i,j]/1e6)*z(pat/1e6);

Продолжение приложения А

Qdob:=Qdob+qm[i,j]*dt;

if h[i,j]<>0 then begin Pcp:=Pcp+pp^[i,j]*h[i,j];Hcp:=Hcp+h[i,j] end;

end;

end;

{ if (ZapN-ZapT-Qdob)/Qdob*100>delta then delta:=(ZapN-ZapT-Qdob)/Qdob*100;}

Pcp:=Pcp/Hcp;

PcpVII:=Pcp;

p5_2:=pp^[10,13];

p8_2:=pp^[10,26];

p10_2:=pp^[6,15];

p16_2:=pp^[7,29];

p75_2:=pp^[9,5];

p15:=pp^[14,23];

p11:=pp^[15,11];

FreeMem(kh,M*N*SizeOf(tipR));

FreeMem(b1,M*N*SizeOf(tipR));

FreeMem(pp,M*N*SizeOf(tipR));

FreeMem(dp,M*N*SizeOf(tipR));

FreeMem(V,M*N*SizeOf(tipR));

FreeMem(g1,M*N*SizeOf(tipR));

FreeMem(f1,M*N*SizeOf(tipR));

END;

procedure soprVII(gi:m;dt,t:real; var u5_2,u8_2,u10_2,u16_2,u75_2,u11,u15:real;

e5_1,e5_2,e8_1,e8_2,e10_1,e10_2,e16_1,e16_2,e75_1,e75_2,e11,e15:real);

const

N=19;{strok}

M=34;{stolbcov}

Ppl=7.6e6;{Pa} {!!!!!!!!!!!!!!!!!!}

k1=0.661e-12;{mkm2, IV block}

k2=0.44e-12;{mkm2, III block}

kp=0.32;

al1=0.77;

al2=0.76;

mu=0.011e-3;{Pa*sec}

dx=100;{m}

dy=100;{m}

pat=100000;

day=86400;

type

tipR=real;

mas1=array[1..20]of tipR;

mas2=array[1..N,1..M]of tipR;

dmas=^mas2;

VAR

dat,dat1:text;

al:mas1;

i,j,count,k,cell:integer;

qm,p,u:mas2;

g,c,a,b,f,d,R,sum,ZapN,ZapT,Qdob,Pcp,delta,Hcp,ggg,fff,bbb,v1:tipR;

b1,g1,f1,kh,V,dp,pp:dmas;

flag:boolean;

procedure alf(var al:mas1);

var

i,j:integer;

par,par1,par11,par2:tipR;

begin

par11:=100;

for i:=1 to N-1 do

for j:=1 to M-1 do

if h[i,j]<>0 then

if (h[i,j+1]<>0)and(h[i+1,j]<>0) then

begin

par:=dx*dx*(kh^[i,j+1]+kh^[i,j])/(kh^[i+1,j]+kh^[i,j])/dy/dy;

par1:=pi*pi/(2*N*N*(1+par));

par2:=pi*pi/(2*M*M*(1+1/par));

if par1<par2 then

begin

if par1<par11 then par11:=par1;

end

else

begin

if par2<par11 then par11:=par2;

end;

end;

al[1]:=1;

for i:=2 to 10 do al[i]:=1-power(par11*10,(i-1)/10);

end;

procedure skv;

begin

{5}

qm[10,13]:=e5_2{/86.4};

{75}

qm[9,5]:=e75_2{/86.4};

Продолжение приложения А

{10}

qm[6,15]:=e10_2{/86.4};

{11}

qm[15,11]:=e11{/86.4};

{8}

qm[10,26]:=e8_2{/86.4};

{16}

qm[7,29]:=e16_2{/86.4};

{15}

qm[14,23]:=e15{/86.4};

end;

BEGIN

GetMem(kh,M*N*SizeOf(tipR));

GetMem(b1,M*N*SizeOf(tipR));

GetMem(pp,M*N*SizeOf(tipR));

GetMem(dp,M*N*SizeOf(tipR));

GetMem(V,M*N*SizeOf(tipR));

GetMem(g1,M*N*SizeOf(tipR));

GetMem(f1,M*N*SizeOf(tipR));

for i:=1 to N do

begin

for j:=1 to M do

begin

if h[i,j]<>0 then

if j<=21 then kh^[i,j]:=k2*h[i,j]

else kh^[i,j]:=k1*h[i,j]

else kh^[i,j]:=0;

p[i,j]:=gi[i,j];

u[i,j]:=0;

pp^[i,j]:=u[i,j];

qm[i,j]:=0;

b1^[i,j]:=0;

g1^[i,j]:=0;

f1^[i,j]:=0;

V^[i,j]:=0;

dp^[i,j]:=0;

end;

end;

alf(al);

sum:=0;delta:=0;

count:=0;

k:=0;

skv;

repeat

count:=count+1;

if odd(count) then

begiN {obratnoe yporyadochivanie}

for i:=2 to N-1 do

for j:=M-1 downto 2 do

begin

g:=0;c:=0;b:=0;f:=0;a:=0;d:=0;R:=0;

if h[i,j]<>0 then

begin

if h[i,j-1]=0 then g:=0

else g:=(dx/dy)*p[i,j]*(kh^[i,j-1]+kh^[i,j])/mu/z(p[i,j]/1e6)/2;

if h[i-1,j]=0 then c:=0

else c:=(dy/dx)*p[i,j]*(kh^[i-1,j]+kh^[i,j])/mu/z(p[i,j]/1e6)/2;

if h[i+1,j]=0 then b:=0

else b:=(dy/dx)*p[i,j]*(kh^[i+1,j]+kh^[i,j])/mu/z(p[i,j]/1e6)/2;

if h[i,j+1]=0 then f:=0

else f:=(dx/dy)*p[i,j]*(kh^[i,j+1]+kh^[i,j])/mu/z(p[i,j]/1e6)/2;

a:=-(g+c+b+f)-dx*dy/dt*al1*kp*h[i,j];

d:=pat*qm[i,j]/dx/dy-dx*dy/dt*al1*kp*h[i,j]*p[i,j]/z(p[i,j]/1e6)*u[i,j];

R:=d-(g*pp^[i,j-1]+c*pp^[i-1,j]+a*pp^[i,j]+b*pp^[i+1,j]+f*pp^[i,j+1]);

if odd(count) then k:=k+1;

if k>10 then k:=1;

sipp(i,j,g,c,a,b,f,R,g1^[i-1,j],f1^[i,j+1],al[k],V^[i-1,j],V^[i,j+1],ggg,fff,v1);

V^[i,j]:=v1;

g1^[i,j]:=ggg;

f1^[i,j]:=fff;

end;

end;

flag:=true;

for i:=N-1 downto 2 do

for j:=2 to M-1 do

if h[i,j]<>0 then

begin

dp^[i,j]:=V^[i,j]-g1^[i,j]*dp^[i,j-1]-f1^[i,j]*dp^[i+1,j];

if abs(dp^[i,j])>1e-4 then flag:=false;

pp^[i,j]:=pp^[i,j]+dp^[i,j];

end;

{} for i:=1 to N do

for j:=1 to M do begin g1^[i,j]:=0;f1^[i,j]:=0 end;

Продолжение приложения А

enD

eLSe

begiN

for i:=2 to N-1 do

for j:=2 to M-1 do

begin

g:=0;c:=0;b:=0;f:=0;a:=0;d:=0;R:=0;

if h[i,j]<>0 then

begin

if h[i,j-1]=0 then g:=0

else g:=(dx/dy)*p[i,j]*(kh^[i,j-1]+kh^[i,j])/mu/z(p[i,j]/1e6)/2;

if h[i-1,j]=0 then c:=0

else c:=(dy/dx)*p[i,j]*(kh^[i-1,j]+kh^[i,j])/mu/z(p[i,j]/1e6)/2;

if h[i+1,j]=0 then b:=0

else b:=(dy/dx)*p[i,j]*(kh^[i+1,j]+kh^[i,j])/mu/z(p[i,j]/1e6)/2;

if h[i,j+1]=0 then f:=0

else f:=(dx/dy)*p[i,j]*(kh^[i,j+1]+kh^[i,j])/mu/z(p[i,j]/1e6)/2;

a:=-(g+c+b+f)-dx*dy/dt*al1*kp*h[i,j];

d:=pat*qm[i,j]/dx/dy-dx*dy/dt*al1*kp*h[i,j]*p[i,j]/z(p[i,j]/1e6)*u[i,j];

R:=d-(g*pp^[i,j-1]+c*pp^[i-1,j]+a*pp^[i,j]+b*pp^[i+1,j]+f*pp^[i,j+1]);

if odd(count) then k:=k+1;

if k>10 then k:=1;

sipp1(i,j,g,c,a,b,f,R,b1^[i,j-1],f1^[i-1,j],al[k],V^[i-1,j],V^[i,j-1],bbb,fff,v1);

V^[i,j]:=v1;

b1^[i,j]:=bbb;

f1^[i,j]:=fff;

end;

end;

flag:=true;

for i:=N-1 downto 2 do

for j:=M-1 downto 2 do

if h[i,j]<>0 then

begin

dp^[i,j]:=V^[i,j]-b1^[i,j]*dp^[i+1,j]-f1^[i,j]*dp^[i,j+1];

if abs(dp^[i,j])>1e-4 then flag:=false;

pp^[i,j]:=pp^[i,j]+dp^[i,j];

end;

{} for i:=1 to N do

for j:=1 to M do begin b1^[i,j]:=0;f1^[i,j]:=0 end;

enD;

until flag;

u5_2:=pp^[10,13];

u8_2:=pp^[10,26];

u10_2:=pp^[6,15];

u16_2:=pp^[7,29];

u75_2:=pp^[9,5];

u15:=pp^[14,23];

u11:=pp^[15,11];

FreeMem(kh,M*N*SizeOf(tipR));

FreeMem(b1,M*N*SizeOf(tipR));

FreeMem(pp,M*N*SizeOf(tipR));

FreeMem(dp,M*N*SizeOf(tipR));

FreeMem(V,M*N*SizeOf(tipR));

FreeMem(g1,M*N*SizeOf(tipR));

FreeMem(f1,M*N*SizeOf(tipR));

END;

procedure VIIa(giv1:m1;var gi1:m1;dt,t:real; var PcpVIIa,p5_1,p8_1,p10_1,p16_1,p75_1:real;

q5_1,q5_2,q8_1,q8_2,q10_1,q10_2,q16_1,q16_2,q75_1,q75_2:real);

const

N=16;{strok}

M=34;{stolbcov}

Ppl=7.6e6;{Pa}

k1=1.877e-12;{mkm2, IV block}

k2=2.062e-12;{mkm2, III block}

kp=0.3;

al1=0.78;

al2=0.76;

mu=0.011e-3;{Pa*sec}

dx=100;{m}

dy=100;{m}

pat=100000;

day=86400;

type

tipR=real;

mas1=array[1..10]of tipR;

mas2=array[1..N,1..M]of tipR;

dmas=^mas2;

VAR

dat,dat1:text;

al:mas1;

i,j,count,k,cell:integer;

qm,p,h:mas2;

g,c,a,b,f,d,R,sum,ZapN,ZapT,Qdob,Pcp,delta,Hcp,ggg,fff,bbb,v1:tipR;

Продолжение приложения А

b1,g1,f1,kh,V,dp,pp:dmas;

flag:boolean;

procedure alf(var al:mas1);

var

i,j:integer;

par,par1,par11,par2:tipR;

begin

par11:=100;

for i:=1 to N-1 do

for j:=1 to M-1 do

if h[i,j]<>0 then

if (h[i,j+1]<>0)and(h[i+1,j]<>0) then

begin

par:=dx*dx*(kh^[i,j+1]+kh^[i,j])/(kh^[i+1,j]+kh^[i,j])/dy/dy;

par1:=pi*pi/(2*N*N*(1+par));

par2:=pi*pi/(2*M*M*(1+1/par));

if par1<par2 then

begin

if par1<par11 then par11:=par1;

end

else

begin

if par2<par11 then par11:=par2;

end;

end;

al[1]:=1;

for i:=2 to 4 do al[i]:=1-power(par11,(i-1)/3);

end;

procedure skv;

begin

{5}

qm[9,13]:=q5_1/86.4;

{75}

qm[8,5]:=q75_1/86.4;

{10}

qm[5,15]:=q10_1/86.4;

{8}

qm[9,26]:=q8_1/86.4;

{16}

qm[6,29]:=q16_1/86.4;

end;

BEGIN

GetMem(kh,M*N*SizeOf(tipR));

GetMem(b1,M*N*SizeOf(tipR));

GetMem(pp,M*N*SizeOf(tipR));

GetMem(dp,M*N*SizeOf(tipR));

GetMem(V,M*N*SizeOf(tipR));

GetMem(g1,M*N*SizeOf(tipR));

GetMem(f1,M*N*SizeOf(tipR));

ZapN:=0;

for i:=1 to N do

begin

for j:=1 to M do

begin

if h[i,j]<>0 then

begin

if j<=21 then kh^[i,j]:=k2*h[i,j]

else kh^[i,j]:=k1*h[i,j];

p[i,j]:=giv1[i,j];

ZapN:=ZapN+h[i,j]*dx*dy*kp*al1*Ppl/pat/z(Ppl/1e6)*z(pat/1e6);

end

else

begin

kh^[i,j]:=0;

p[i,j]:=0;

end;

pp^[i,j]:=p[i,j];

qm[i,j]:=0;

b1^[i,j]:=0;

g1^[i,j]:=0;

f1^[i,j]:=0;

V^[i,j]:=0;

dp^[i,j]:=0;

end;

end;

alf(al);

sum:=0;Qdob:=0;delta:=0;

count:=0;

k:=0;

skv;

repeat

count:=count+1;

Продолжение приложения А

if odd(count) then

begiN {obratnoe yporyadochivanie}

for i:=2 to N-1 do

for j:=M-1 downto 2 do

begin

g:=0;c:=0;b:=0;f:=0;a:=0;d:=0;R:=0;

if h[i,j]<>0 then

begin

if h[i,j-1]=0 then g:=0

else g:=(dx/dy)*pp^[i,j]*(kh^[i,j-1]+kh^[i,j])/mu/z(pp^[i,j]/1e6)/2;

if h[i-1,j]=0 then c:=0

else c:=(dy/dx)*pp^[i,j]*(kh^[i-1,j]+kh^[i,j])/mu/z(pp^[i,j]/1e6)/2;

if h[i+1,j]=0 then b:=0

else b:=(dy/dx)*pp^[i,j]*(kh^[i+1,j]+kh^[i,j])/mu/z(pp^[i,j]/1e6)/2;

if h[i,j+1]=0 then f:=0

else f:=(dx/dy)*pp^[i,j]*(kh^[i,j+1]+kh^[i,j])/mu/z(pp^[i,j]/1e6)/2;

a:=-(g+c+b+f)-dx*dy/dt*2*al1*kp*h[i,j]/z(pp^[i,j]/1e6);

d:=2*pat*qm[i,j]-dx*dy/dt*2*al1*kp*h[i,j]*p[i,j]/z(p[i,j]/1e6);

R:=d-(g*pp^[i,j-1]+c*pp^[i-1,j]+a*pp^[i,j]+b*pp^[i+1,j]+f*pp^[i,j+1]);

if odd(count) then k:=k+1;

if k>4 then k:=1;

sipp(i,j,g,c,a,b,f,R,g1^[i-1,j],f1^[i,j+1],al[k],V^[i-1,j],V^[i,j+1],ggg,fff,v1);

V^[i,j]:=v1;

g1^[i,j]:=ggg;

f1^[i,j]:=fff;

end;

end;

flag:=true;

for i:=N-1 downto 2 do

for j:=2 to M-1 do

if h[i,j]<>0 then

begin

dp^[i,j]:=V^[i,j]-g1^[i,j]*dp^[i,j-1]-f1^[i,j]*dp^[i+1,j];

if abs(dp^[i,j])>1e-4 then flag:=false;

pp^[i,j]:=pp^[i,j]+dp^[i,j];

end;

for i:=1 to N do

for j:=1 to M do begin g1^[i,j]:=0;f1^[i,j]:=0 end;

enD

eLSe

begiN

for i:=2 to N-1 do

for j:=2 to M-1 do

begin

g:=0;c:=0;b:=0;f:=0;a:=0;d:=0;R:=0;

if h[i,j]<>0 then

begin

if h[i,j-1]=0 then g:=0

else g:=(dx/dy)*pp^[i,j]*(kh^[i,j-1]+kh^[i,j])/mu/z(pp^[i,j]/1e6)/2;

if h[i-1,j]=0 then c:=0

else c:=(dy/dx)*pp^[i,j]*(kh^[i-1,j]+kh^[i,j])/mu/z(pp^[i,j]/1e6)/2;

if h[i+1,j]=0 then b:=0

else b:=(dy/dx)*pp^[i,j]*(kh^[i+1,j]+kh^[i,j])/mu/z(pp^[i,j]/1e6)/2;

if h[i,j+1]=0 then f:=0

else f:=(dx/dy)*pp^[i,j]*(kh^[i,j+1]+kh^[i,j])/mu/z(pp^[i,j]/1e6)/2;

a:=-(g+c+b+f)-dx*dy/dt*2*al1*kp*h[i,j]/z(pp^[i,j]/1e6);

d:=2*pat*qm[i,j]-dx*dy/dt*2*al1*kp*h[i,j]*p[i,j]/z(p[i,j]/1e6);

R:=d-(g*pp^[i,j-1]+c*pp^[i-1,j]+a*pp^[i,j]+b*pp^[i+1,j]+f*pp^[i,j+1]);

if odd(count) then k:=k+1;

if k>4 then k:=1;

sipp1(i,j,g,c,a,b,f,R,b1^[i,j-1],f1^[i-1,j],al[k],V^[i-1,j],V^[i,j-1],bbb,fff,v1);

V^[i,j]:=v1;

b1^[i,j]:=bbb;

f1^[i,j]:=fff;

end;

end;

flag:=true;

for i:=N-1 downto 2 do

for j:=M-1 downto 2 do

if h[i,j]<>0 then

begin

dp^[i,j]:=V^[i,j]-b1^[i,j]*dp^[i+1,j]-f1^[i,j]*dp^[i,j+1];

if abs(dp^[i,j])>1e-4 then flag:=false;

pp^[i,j]:=pp^[i,j]+dp^[i,j];

end;

for i:=1 to N do

for j:=1 to M do begin b1^[i,j]:=0;f1^[i,j]:=0 end;

enD;

until flag;

ZapT:=0;Pcp:=0;Hcp:=0;

for i:=1 to N do

begin

for j:=1 to M do

begin

Продолжение приложения А

gi1[i,j]:=pp^[i,j];

ZapT:=ZapT+h[i,j]*dx*dy*kp*al1*p[i,j]/pat/z(p[i,j]/1e6)*z(pat/1e6);

Qdob:=Qdob+qm[i,j]*dt;

if h[i,j]<>0 then begin Pcp:=Pcp+pp^[i,j]*h[i,j];Hcp:=Hcp+h[i,j] end;

end;

end;

if (ZapN-ZapT-Qdob)/Qdob*100>delta then delta:=(ZapN-ZapT-Qdob)/Qdob*100;

Pcp:=Pcp/Hcp;

PcpVIIa:=Pcp;

p5_1:=pp^[9,13];

p8_1:=pp^[9,26];

p10_1:=pp^[5,15];

p16_1:=pp^[6,29];

p75_1:=pp^[8,5];

FreeMem(kh,M*N*SizeOf(tipR));

FreeMem(b1,M*N*SizeOf(tipR));

FreeMem(pp,M*N*SizeOf(tipR));

FreeMem(dp,M*N*SizeOf(tipR));

FreeMem(V,M*N*SizeOf(tipR));

FreeMem(g1,M*N*SizeOf(tipR));

FreeMem(f1,M*N*SizeOf(tipR));

END;

procedure soprVIIa(gi1:m1;dt,t:real; var u5_1,u8_1,u10_1,u16_1,u75_1:real;

e5_1,e5_2,e8_1,e8_2,e10_1,e10_2,e16_1,e16_2,e75_1,e75_2:real);

const

N=16;{strok}

M=34;{stolbcov}

Ppl=7.8e6;{Pa}

k1=1.877e-12;{mkm2, IV block}

k2=2.062e-12;{mkm2, III block}

kp=0.3;

al1=0.78;

al2=0.76;

mu=0.011e-3;{Pa*sec}

dx=100;{m}

dy=100;{m}

pat=100000;

day=86400;

type

tipR=real;

mas1=array[1..10]of tipR;

mas2=array[1..N,1..M]of tipR;

dmas=^mas2;

VAR

dat,dat1:text;

al:mas1;

i,j,count,k,cell:integer;

qm,p,h,u:mas2;

g,c,a,b,f,d,R,sum,ZapN,ZapT,Qdob,Pcp,delta,Hcp,ggg,fff,bbb,v1:tipR;

b1,g1,f1,kh,V,dp,pp:dmas;

flag:boolean;

procedure alf(var al:mas1);

var

i,j:integer;

par,par1,par11,par2:tipR;

begin

par11:=100;

for i:=1 to N-1 do

for j:=1 to M-1 do

if h[i,j]<>0 then

if (h[i,j+1]<>0)and(h[i+1,j]<>0) then

begin

par:=dx*dx*(kh^[i,j+1]+kh^[i,j])/(kh^[i+1,j]+kh^[i,j])/dy/dy;

par1:=pi*pi/(2*N*N*(1+par));

par2:=pi*pi/(2*M*M*(1+1/par));

if par1<par2 then

begin

if par1<par11 then par11:=par1;

end

else

begin

if par2<par11 then par11:=par2;

end;

end;

al[1]:=1;

for i:=2 to 4 do al[i]:=1-power(par11,(i-1)/3);

end;

procedure skv;

begin

{5}

qm[9,13]:=e5_1{/86.4};

{75}

qm[8,5]:=e75_1{/86.4};

{10}

Продолжение приложения А

qm[5,15]:=e10_1{/86.4};

{8}

qm[9,26]:=e8_1{/86.4};

{16}

qm[6,29]:=e16_1{/86.4};

end;

BEGIN

GetMem(kh,M*N*SizeOf(tipR));

GetMem(b1,M*N*SizeOf(tipR));

GetMem(pp,M*N*SizeOf(tipR));

GetMem(dp,M*N*SizeOf(tipR));

GetMem(V,M*N*SizeOf(tipR));

GetMem(g1,M*N*SizeOf(tipR));

GetMem(f1,M*N*SizeOf(tipR));

ZapN:=0;

for i:=1 to N do

begin

for j:=1 to M do

begin

if h[i,j]<>0 then

if j<=21 then kh^[i,j]:=k2*h[i,j]

else kh^[i,j]:=k1*h[i,j]

else kh^[i,j]:=0;

p[i,j]:=gi1[i,j];

u[i,j]:=0;

pp^[i,j]:=u[i,j];

qm[i,j]:=0;

b1^[i,j]:=0;

g1^[i,j]:=0;

f1^[i,j]:=0;

V^[i,j]:=0;

dp^[i,j]:=0;

end;

end;

alf(al);

sum:=0;Qdob:=0;delta:=0;

count:=0;

k:=0;

skv;

repeat

count:=count+1;

if odd(count) then

begiN {obratnoe yporyadochivanie}

for i:=2 to N-1 do

for j:=M-1 downto 2 do

begin

g:=0;c:=0;b:=0;f:=0;a:=0;d:=0;R:=0;

if h[i,j]<>0 then

begin

if h[i,j-1]=0 then g:=0

else g:=(dx/dy)*p[i,j]*(kh^[i,j-1]+kh^[i,j])/mu/z(p[i,j]/1e6)/2;

if h[i-1,j]=0 then c:=0

else c:=(dy/dx)*p[i,j]*(kh^[i-1,j]+kh^[i,j])/mu/z(p[i,j]/1e6)/2;

if h[i+1,j]=0 then b:=0

else b:=(dy/dx)*p[i,j]*(kh^[i+1,j]+kh^[i,j])/mu/z(p[i,j]/1e6)/2;

if h[i,j+1]=0 then f:=0

else f:=(dx/dy)*p[i,j]*(kh^[i,j+1]+kh^[i,j])/mu/z(p[i,j]/1e6)/2;

a:=-(g+c+b+f)-dx*dy/dt*al1*kp*h[i,j];

d:=pat*qm[i,j]/dx/dy-dx*dy/dt*al1*kp*h[i,j]*p[i,j]/z(p[i,j]/1e6)*u[i,j];

R:=d-(g*pp^[i,j-1]+c*pp^[i-1,j]+a*pp^[i,j]+b*pp^[i+1,j]+f*pp^[i,j+1]);

if odd(count) then k:=k+1;

if k>4 then k:=1;

sipp(i,j,g,c,a,b,f,R,g1^[i-1,j],f1^[i,j+1],al[k],V^[i-1,j],V^[i,j+1],ggg,fff,v1);

V^[i,j]:=v1;

g1^[i,j]:=ggg;

f1^[i,j]:=fff;

end;

end;

flag:=true;

for i:=N-1 downto 2 do

for j:=2 to M-1 do

if h[i,j]<>0 then

begin

dp^[i,j]:=V^[i,j]-g1^[i,j]*dp^[i,j-1]-f1^[i,j]*dp^[i+1,j];

if abs(dp^[i,j])>1e-4 then flag:=false;

pp^[i,j]:=pp^[i,j]+dp^[i,j];

end;

for i:=1 to N do

for j:=1 to M do begin g1^[i,j]:=0;f1^[i,j]:=0 end;

enD

eLSe

begiN

for i:=2 to N-1 do

Продолжение приложения А

for j:=2 to M-1 do

begin

g:=0;c:=0;b:=0;f:=0;a:=0;d:=0;R:=0;

if h[i,j]<>0 then

begin

if h[i,j-1]=0 then g:=0

else g:=(dx/dy)*p[i,j]*(kh^[i,j-1]+kh^[i,j])/mu/z(p[i,j]/1e6)/2;

if h[i-1,j]=0 then c:=0

else c:=(dy/dx)*p[i,j]*(kh^[i-1,j]+kh^[i,j])/mu/z(p[i,j]/1e6)/2;

if h[i+1,j]=0 then b:=0

else b:=(dy/dx)*p[i,j]*(kh^[i+1,j]+kh^[i,j])/mu/z(p[i,j]/1e6)/2;

if h[i,j+1]=0 then f:=0

else f:=(dx/dy)*p[i,j]*(kh^[i,j+1]+kh^[i,j])/mu/z(p[i,j]/1e6)/2;

a:=-(g+c+b+f)-dx*dy/dt*al1*kp*h[i,j];

d:=pat*qm[i,j]/dx/dy-dx*dy/dt*al1*kp*h[i,j]*p[i,j]/z(p[i,j]/1e6)*u[i,j];

R:=d-(g*pp^[i,j-1]+c*pp^[i-1,j]+a*pp^[i,j]+b*pp^[i+1,j]+f*pp^[i,j+1]);

if odd(count) then k:=k+1;

if k>4 then k:=1;

sipp1(i,j,g,c,a,b,f,R,b1^[i,j-1],f1^[i-1,j],al[k],V^[i-1,j],V^[i,j-1],bbb,fff,v1);

V^[i,j]:=v1;

b1^[i,j]:=bbb;

f1^[i,j]:=fff;

end;

end;

flag:=true;

for i:=N-1 downto 2 do

for j:=M-1 downto 2 do

if h[i,j]<>0 then

begin

dp^[i,j]:=V^[i,j]-b1^[i,j]*dp^[i+1,j]-f1^[i,j]*dp^[i,j+1];

if abs(dp^[i,j])>1e-4 then flag:=false;

pp^[i,j]:=pp^[i,j]+dp^[i,j];

end;

for i:=1 to N do

for j:=1 to M do begin b1^[i,j]:=0;f1^[i,j]:=0 end;

enD;

until flag;

u5_1:=pp^[9,13];

u8_1:=pp^[9,26];

u10_1:=pp^[5,15];

u16_1:=pp^[6,29];

u75_1:=pp^[8,5];

FreeMem(kh,M*N*SizeOf(tipR));

FreeMem(b1,M*N*SizeOf(tipR));

FreeMem(pp,M*N*SizeOf(tipR));

FreeMem(dp,M*N*SizeOf(tipR));

FreeMem(V,M*N*SizeOf(tipR));

FreeMem(g1,M*N*SizeOf(tipR));

FreeMem(f1,M*N*SizeOf(tipR));

END;

CONST

day=86400;

VAR

f,f1,shit:text;

i,j,iter:integer;

giv,gi:m;giv1,gi1:m1;

py5,py8,py10,py16,py75,py11,py15:real;

p5_1,p5_2,p8_1,p8_2,p10_1,p10_2,p16_1,p16_2,p75_1,p75_2,p11,p15:real;

q5_1,q5_2,q8_1,q8_2,q10_1,q10_2,q16_1,q16_2,q75_1,q75_2:real;

e5_1,e5_2,e8_1,e8_2,e10_1,e10_2,e16_1,e16_2,e75_1,e75_2,e11,e15:real;

u5_1,u5_2,u8_1,u8_2,u10_1,u10_2,u16_1,u16_2,u75_1,u75_2,u11,u15:real;

r5_1,r5_2,r8_1,r8_2,r10_1,r10_2,r16_1,r16_2,r75_1,r75_2,r11,r15:real;

q5,q8,q10,q16,q75,q11,q15,dob:real;

dt,t,PcpVII,PcpVIIa,Qdob,Qdob1,Sum,Sum1:real;

QQ,Q1,JJJ,JJ,J1,JJ1,J11,JJJ1,lam,gam,lam1,gam1:real;

pad:boolean;

function zab(py,q:real):real;

var

tet,cp,p3,pp3,tcp,s:real;

begin

q:=abs(q);

pp3:=py;

tcp:=(299-297)/ln(299/297);

repeat

p3:=pp3;

cp:=(p3+py)/2;

s:=0.03415*0.57*700/tcp/z(cp);

tet:=1.325e-12*0.025*z(cp)*z(cp)*tcp*tcp*(exp(2*s)-1)/

6.22e-2/6.22e-2/6.22e-2/6.22e-2/6.22e-2;

pp3:=sqrt(py*py*exp(-2*s)-tet*q*q*exp(-2*s));

until abs(pp3-p3)<0.05;

zab:=pp3;

end;

Продолжение приложения А

procedure Debit(p5_1,p5_2,p8_1,p8_2,p10_1,p10_2,p16_1,p16_2,p75_1,p75_2:real;

var q5,q8,q10,q16,q75:real;

var q5_1,q5_2,q8_1,q8_2,q10_1,q10_2,q16_1,q16_2,q75_1,q75_2:real);

const

ddp=0;

var

pc11,pc15,q,pc5_1,pc8_1,pc10_1,pc16_1,pc75_1,pc5_2,pc8_2,pc10_2,pc16_2,pc75_2:real;

begin

{5}

q:=q5/2;

repeat

q5_1:=q;

if q>0 then

pc5_1:=sqrt(p5_1*p5_1-0.0011*abs(q5_1)-0.000005*q5_1*q5_1)

else

pc5_1:=sqrt(p5_1*p5_1+0.0011*abs(q5_1)+0.000005*q5_1*q5_1);

pc5_2:=pc5_1-ddp{zab(pc5_1,q5_1)};

if pc5_2<p5_2 then

q5_2:=(-0.007+sqrt(0.007*0.007+4*0.000045*(p5_2*p5_2-pc5_2*pc5_2)))/(2*0.000045)

else

q5_2:=-(-0.007+sqrt(0.007*0.007+4*0.000045*(pc5_2*pc5_2-p5_2*p5_2)))/(2*0.000045);

q:=q5-q5_2;

until abs(q-q5_1)<0.0001;

{8}

q:=q8/2;

repeat

q8_1:=q;

if q>0 then

pc8_1:=sqrt(p8_1*p8_1-0.0005*abs(q8_1)-0.00001*q8_1*q8_1)

else

pc8_1:=sqrt(p8_1*p8_1+0.0005*abs(q8_1)+0.00001*q8_1*q8_1);

pc8_2:=pc8_1-ddp{zab(pc8_1,q8_1)};

if pc8_2<p8_2 then

q8_2:=(-0.0037+sqrt(0.0037*0.0037+4*0.000046*(p8_2*p8_2-pc8_2*pc8_2)))/(2*0.000046)

else

q8_2:=-(-0.0037+sqrt(0.0037*0.0037+4*0.000046*(pc8_2*pc8_2-p8_2*p8_2)))/(2*0.000046);

q:=q8-q8_2;

until abs(q-q8_1)<0.0001;

{10}

if q10=0 then begin q10_1:=0;q10_2:=0 end

else

begin

q:=q10/2;

repeat

q10_1:=q;

if q>0 then

pc10_1:=sqrt(p10_1*p10_1-0.0011*abs(q10_1)-0.000005*q10_1*q10_1)

else

pc10_1:=sqrt(p10_1*p10_1+0.0011*abs(q10_1)+0.000005*q10_1*q10_1);

pc10_2:=pc10_1-ddp{zab(pc10_1,q10_1)};

if pc10_2<p10_2 then

q10_2:=(-0.007+sqrt(0.007*0.007+4*0.000045*(p10_2*p10_2-pc10_2*pc10_2)))/(2*0.000045)

else

q10_2:=-(-0.007+sqrt(0.007*0.007+4*0.000045*(pc10_2*pc10_2-p10_2*p10_2)))/(2*0.000045);

q:=q10-q10_2;

until abs(q-q10_1)<0.0001;

end;

{16}

if q16=0 then begin q16_1:=0;q16_2:=0 end

else

begin

q:=q16/2;

repeat

q16_1:=q;

if q>0 then

pc16_1:=sqrt(p16_1*p16_1-0.0005*abs(q16_1)-0.00001*q16_1*q16_1)

else

pc16_1:=sqrt(p16_1*p16_1+0.0005*abs(q16_1)+0.00001*q16_1*q16_1);

pc16_2:=pc16_1-ddp{zab(pc16_1,q16_1)};

if pc16_2<p16_2 then

q16_2:=(-0.0037+sqrt(0.0037*0.0037+4*0.000046*(p16_2*p16_2-pc16_2*pc16_2)))/(2*0.000046)

else

q16_2:=-(-0.0037+sqrt(0.0037*0.0037+4*0.000046*(pc16_2*pc16_2-p16_2*p16_2)))/(2*0.00046);

q:=q16-q16_2;

until abs(q-q16_1)<0.0001;

end;

{75}

if q75=0 then begin q75_1:=0;q75_2:=0 end

else

begin

q:=q75/2;

repeat

q75_1:=q;

if q>0 then

Продолжение приложения А

pc75_1:=sqrt(p75_1*p75_1-0.0011*abs(q75_1)-0.000005*q75_1*q75_1)

else

pc75_1:=sqrt(p75_1*p75_1+0.0011*abs(q75_1)+0.000005*q75_1*q75_1);

pc75_2:=pc75_1-ddp{zab(pc75_1,q75_1)};

if pc75_2<p75_2 then

q75_2:=(-0.007+sqrt(0.007*0.007+4*0.000045*(p75_2*p75_2-pc75_2*pc75_2)))/(2*0.000045)

else

q75_2:=-(-0.007+sqrt(0.007*0.007+4*0.000045*(pc75_2*pc75_2-p75_2*p75_2)))/(2*0.000045);

q:=q75-q75_2;

until abs(q-q75_1)<0.0001;

end;

pc11:=sqrt(PcpVII*PcpVII/1e12-0.0011*abs(q11)-0.000005*q11*q11);

pc15:=sqrt(PcpVII*PcpVII/1e12-0.0011*abs(q15)-0.000005*q15*q15);

py5:=zab(pc5_1,q5);

py8:=zab(pc8_1,q8);

py10:=zab(pc10_1,q10);

py16:=zab(pc16_1,q16);

py75:=zab(pc75_1,q75);

py11:=zab(pc11,q11);

py15:=zab(pc15,q15);

if py5<4.5 then pad:=true;

if py8<4.5 then pad:=true;

if py10<4.5 then pad:=true;

if py16<4.5 then pad:=true;

if py75<4.5 then pad:=true;

if py11<4.5 then pad:=true;

if py15<4.5 then pad:=true;

end;

procedure Ytoch(var q5,q8,q10,q16,q75,q11,q15:real);

begin

WRiteln('Ydo=',(q5+q8+q10+q75+q16+q11+q15):1:13);

q5:=q5{q5_2}-gam*(r5_2*u5_2-lam){+q5_1}-gam1*(r5_1*u5_1-lam1);

if q5<10 then q5:=10

else if q5>100.9 then q5:=100.9;

q8:=q8{q8_2}-gam*(r8_2*u8_2-lam){+q8_1}-gam1*(r8_1*u8_1-lam1);

if q8<10 then q8:=10

else if q8>84.1 then q8:=84.1;

q10:=q10{q10_2}-gam*(r10_2*u10_2-lam){+q10_1}-gam1*(r10_1*u10_1-lam1);

if q10<10 then q10:=10

else if q10>69.1 then q10:=69.1;

q16:=q16{q16_2}-gam*(r16_2*u16_2-lam){+q16_1}-gam1*(r16_1*u16_1-lam1);

if q16<10 then q16:=10

else if q16>53.1 then q16:=53.1;

q75:=q75{q75_2}-gam*(r75_2*u75_2-lam){+q75_1}-gam1*(r75_1*u75_1-lam1);

if q75<10 then q75:=10

else if q75>100.7 then q75:=100.7;

q11:=q11-gam*(r11*u11-lam);

if q11<10 then q11:=10

else if q11>101.8 then q11:=101.8;

q15:=q15-gam*(r15*u15-lam);