All issues

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

In this study we discuss methods to solve optimal control problems based on neural network techniques. We study hierarchical dynamical two-level system for surface water quality control. The system consists of a supervisor (government) and a few agents (enterprises). We consider this problem from the point of agents. In this case we solve optimal control problem with constraints. To solve this problem, we use Pontryagin’s maximum principle, with which we obtain optimality conditions. To solve emerging ODEs, we use feedforward neural network. We provide a review of existing techniques to study such problems and a review of neural network’s training methods. To estimate the error of numerical solution, we propose to use defect analysis method, adapted for neural networks. This allows one to get quantitative error estimations of numerical solution. We provide examples of our method’s usage for solving synthetic problem and a surface water quality control model. We compare the results of this examples with known solution (when provided) and the results of shooting method. In all cases the errors, estimated by our method are of the same order as the errors compared with known solution. Moreover, we study surface water quality control problem when no solutions is provided by other methods. This happens because of relatively large time interval and/or the case of several agents. In the latter case we seek Nash equilibrium between agents. Thus, in this study we show the ability of neural networks to solve various problems including optimal control problems and differential games and we show the ability of quantitative estimation of an error. From the numerical results we conclude that the presence of the supervisor is necessary for achieving the sustainable development.

Объяснимый искусственный интеллект (Explainable AI, XAI) представляет собой область искусственного интеллекта, направленную на создание методов и инструментов для генерации интерпретируемых и понятных для человека объяснений решений ИИ. Актуальность объяснимости моделей возрастает по мере внедрения искусственного интеллекта в критически важные сферы (медицина, финансы, юриспруденция), где непрозрачность алгоритмов может приводить к серьезным последствиям для пользователей и общества. В работе представлен аналитический обзор современного состояния области XAI, охватывающий теоретические основы, методологию и практические применения.

Рассматриваемые методы объяснимого ИИ были отобраны и систематизированы на основе многоуровневой классификации методов XAI по постановке задачи (цель, целевая аудитория, тип данных), методологии (стадия применения, модель-специфичность, методы, масштаб) и форме результата (представление, презентация, метрики оценки).

Проведен сравнительный анализ методов объяснимого ИИ для различных областей применения. Для классического машинного обучения детально рассмотрены SHAP и LIME с выявлением их теоретических оснований, вычислительных характеристик и ограничений. Для компьютерного зрения систематизированы градиентные методы (SmoothGrad, Integrated Gradients), методы визуализации активаций (Grad-CAM, Grad-CAM++), методы на основе возмущений (RISE, Occlusion) и концептуальные объяснения (TCAV, Network Dissection). Особое внимание уделено специфике применения XAI к обработке естественного языка и большим языковым моделям, включая анализ достоверности цепочек размышлений (Chain-of-Thought), естественно-языковых объяснений и методов на основе графов атрибуции. Выделены фундаментальные ограничения существующих подходов к объяснимости LLM и определены направления дальнейших исследований.

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

Explainable Artificial Intelligence (XAI) is a field of artificial intelligence aimed at creating methods and tools for generating interpretable and human-understandable explanations of AI decisions. The relevance of model explainability increases with the deployment of artificial intelligence in critical domains (healthcare, finance, law), where algorithmic opacity can lead to serious consequences for users and society. This work presents an analytical review of the current state of the XAI field, covering theoretical foundations, methodology, and practical applications.

The examined explainable AI methods were selected and systematized based on a multi-level classification of XAI methods by problem formulation (goal, target audience, data type), methodology (application stage, model-specificity, methods, scale), and result form (representation, presentation, evaluation metrics).

A comparative analysis of explainable AI methods for various application domains is conducted. For classical machine learning, SHAP and LIME are examined in detail, revealing their theoretical foundations, computational characteristics, and limitations. For computer vision, gradient-based methods (SmoothGrad, Integrated Gradients), activation visualization methods (Grad-CAM, Grad-CAM++), perturbation-based methods (RISE, Occlusion), and conceptual explanations (TCAV, Network Dissection) are systematized. Special attention is paid to the specifics of applying XAI to natural language processing and large language models, including analysis of the faithfulness of Chain-of-Thought reasoning, natural language explanations, and attribution graph methods. Fundamental limitations of existing approaches to LLM explainability are identified and directions for future research are defined.

The review results demonstrate that XAI methods have reached significant maturity in classical machine learning and computer vision, however, their application to large language models remains an open research problem requiring the development of new explanation paradigms.

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

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

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

The article presents a method for linearization the effective function of material instantaneous deformation in order to generalize the torsional vibration equation to the case of nonlinearly deformable rheologically active shafts. It is considered layered and structurally heterogeneous, on average isotropic shafts made of nonlinearly viscoelastic components. The technique consists in determining the approximate shear modulus by minimizing the root-mean-square deviation in approximation of the effective diagram of instantaneous deformation.

The method allows to estimate analytically values of natural frequencies of layered and structurally heterogeneous nonlinearly viscoelastic shaft. This makes it possible to significantly reduce resources in vibration analysis, as well as to track changes in values of natural frequencies with changing geometric, physico-mechanical and structural parameters of shafts, which is especially important at the initial stages of modeling and design. In addition, the paper shows that only a pronounced nonlinearity of the effective state equation has an effect on the natural frequencies, and in some cases the nonlinearity in determining the natural frequencies can be neglected.

As equations of state of the composite material components, the article considers the equations of nonlinear heredity with instantaneous deformation functions in the form of the Prandtl’s bilinear diagrams. To homogenize the state equations of layered shafts, it is applied the Voigt’s hypothesis on the homogeneity of deformations and the Reuss’ hypothesis on the homogeneity of stresses in the volume of a composite body. Using these assumptions, effective secant and tangential shear moduli, proportionality limits, as well as creep and relaxation kernels of longitudinal, axial and transversely layered shafts are obtained. In addition, it is obtained the indicated effective characteristics of a structurally heterogeneous, on average isotropic shaft using the homogenization method previously proposed by the authors, based on the determination of the material deformation parameters by the rule of a mixture for the Voigt’s and the Reuss’ state equations.

В работе рассматриваются развитие и применение метода учета заполненности прямоугольных ячеек материальной средой, в частности жидкостью для повышения гладкости и точности конечно-разностного решения задач гидродинамики со сложной формой граничной поверхности. Для исследования возможностей предлагаемых разностных схем рассмотрены две задачи вычислительной гидродинамики — пространственно-двумерного течения вязкой жидкости между двумя соосными полуцилиндрами и переноса веществ между соосными полуцилиндрами. Аппроксимация задач по времени выполнена на основе схем расщепления по физическим процессам. Дискретизация операторов диффузии и конвекции выполнена на основе интегроинтерполяционного метода с учетом заполненности ячеек и без ее учета. Для решения задачи диффузии – конвекции при больших сеточных числах Пекле предложено использовать разностную схему, учитывающую функцию заполненности ячеек, и схему, построенную на основе линейной комбинации разностных схем «кабаре» и «крест» с весовыми коэффициентами, полученными в результате минимизации погрешности аппроксимации при малых числах Куранта. Для оценки точности численного решения в качестве эталона используется аналитическое решение, описывающее течение Куэтта – Тейлора. В случае непосредственного использования прямоугольных сеток (ступенчатой аппроксимации границ) относительная погрешность расчетов достигает 70 %, при тех же условиях использование предлагаемого метода позволяет уменьшить погрешность до 6%. Показано, что дробление прямоугольной сетки в 2–8 раз по каждому из пространственных направлений не приводит к такому же повышению точности, которой обладают численные решения, полученные с учетом заполненности ячеек. Предложенные разностные схемы, построенные на основе линейной комбинации разностных схем «кабаре» и «крест» с весовыми коэффициентами 2/3 и 1/3 соответственно, полученные в результате минимизации порядка погрешности аппроксимации, для задачи диффузии – конвекции обладают меньшей сеточной вязкостью и, как следствие, точнее описывают поведение решения в случае больших сеточных чисел Пекле.

The paper discusses the development and application of the accounting rectangular cell fullness method with material substance, in particular, a liquid, to increase the smoothness and accuracy of a finite-difference solution of hydrodynamic problems with a complex shape of the boundary surface. Two problems of computational hydrodynamics are considered to study the possibilities of the proposed difference schemes: the spatial-twodimensional flow of a viscous fluid between two coaxial semi-cylinders and the transfer of substances between coaxial semi-cylinders. Discretization of diffusion and convection operators was performed on the basis of the integro-interpolation method, taking into account taking into account the fullness of cells and without it. It is proposed to use a difference scheme, for solving the problem of diffusion – convection at large grid Peclet numbers, that takes into account the cell population function, and a scheme on the basis of linear combination of the Upwind and Standard Leapfrog difference schemes with weight coefficients obtained by minimizing the approximation error at small Courant numbers. As a reference, an analytical solution describing the Couette – Taylor flow is used to estimate the accuracy of the numerical solution. The relative error of calculations reaches 70% in the case of the direct use of rectangular grids (stepwise approximation of the boundaries), under the same conditions using the proposed method allows to reduce the error to 6%. It is shown that the fragmentation of a rectangular grid by 2–8 times in each of the spatial directions does not lead to the same increase in the accuracy that numerical solutions have, obtained taking into account the fullness of the cells. The proposed difference schemes on the basis of linear combination of the Upwind and Standard Leapfrog difference schemes with weighting factors of 2/3 and 1/3, respectively, obtained by minimizing the order of approximation error, for the diffusion – convection problem have a lower grid viscosity and, as a corollary, more precisely, describe the behavior of the solution in the case of large grid Peclet numbers.

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

For modeling and statistical analysis of data characterized by cyclicity (periodicity) in various areas of science are used circular or wrapped distribution models. The phase distribution function of a harmonic and phase-shift-keying signal in case additive white Gaussian noise is considered. Algorithms for modeling random phases sample of harmonic and modulated signals with specified parameters and correlation function are presented. Expressions for the phase distribution density of the phase-shift-keying signal are given. It is shown that the phase probability density function of the phase-shift-keying signal becomes multimodal. In addition, the probability density function under consideration is a periodic function, which means that the trigonometric Fourier basis can be used to decompose it into a series. In paper for the first time, analytical expressions for the coefficients of the Fourier series when decomposing the density under consideration into a harmonic basis are obtained, and the derivation of the corresponding expressions are presented. Examples of computer modeling and corresponding graphical materials of calculating Fourier coefficients of the phase probability density function for harmonic and phase-shift-keying signals are presented. A formula for the cumulative distribution function and its decomposition into a Fourier series are also obtained. Based on the representation of the phase probability density function in the form of a Fourier series, a comparison is made with other circular distributions often used in practical problems, the Mises distribution and the wrapped normal distribution. The results obtained in this work are of theoretical and practical interest for modeling and statistical analysis of signal phases in various applied problems in area radio engineering, digital communication, radar, etc. In particular, in the problems of estimating the signal-to-noise ratio, the bit error rate, as well as the reliability of demodulator solutions, i. e. soft demodulation of phase-shift-keying signals. Analytical expressions for the Fourier series coefficients can be used to estimate the empirical probability density function.

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

The increased potential for jamming satellite navigation signals, which provide the highest positioning accuracy for moving objects, makes the development of alternative navigation systems comparable in accuracy but autonomous extremely important. One of the most effective such approaches is the integration of inertial and optical navigation systems (NS), as they are more resistant to artificial interference. Among NS data, one can distinguish systems that use a velocity field calculation method when processing optical flow, allowing the determination of the linear and angular velocity vectors of an object in the absence of terrain maps and reference points. However, a serious drawback of this method is the high computational cost of determining the velocity field (optical flow parameters), which is difficult to implement onboard an object. In this regard, the article considers an approach that allows for constructing a tightly coupled autonomous inertial-optical navigation scheme for objects moving along known (programmed) trajectories. This scheme utilizes a navigation algorithm that is easily implemented in onboard computers due to the discovered possibility of estimating the navigation vector without preliminary calculation of optical flow parameters, as well as by utilizing functional dependencies of navigation variables arising on analytical (orthodromic) trajectories. To illustrate the generality of the solution, cases of a rigidly mounted video camera on the object and its two-degree stabilization are studied. The navigation algorithm, providing a stochastic estimate of the full vector of linear and angular motion parameters based on measurements of an integrated tightly coupled inertial-optical NS, is built on the basis of an extended Kalman filter for correlated noise of the object and the observer. A numerical experiment illustrating the effectiveness of the proposed approach is conducted.

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

The semiclassical approximation method is applied for solution construction of the Fokker–Planck equation with quadratic nonlocal nonlinearity and various coefficients in models of asset returns estimation. Analitical expressions determining nonlinear evolution operator are obtained in semiclasical approximation.

Работа посвящена проблеме анализа близости популяционной системы к опасным границам, при пересечении которых в системе разрушается устойчивое сосуществование взаимодействующих популяций. В качестве причины такого разрушения рассматриваются случайные возмущения, неизбежно присутствующие в любой живой системе. Это исследование проводится на примере известной модели взаимодействия популяций хищника и жертвы, учитывающей как стабилизирующий фактор конкуренции хищника за отличные от жертвы ресурсы, так и дестабилизирующий фактор насыщения хищника. Для описания насыщения хищника используется трофическая функция Холлинга второго типа. Динамика системы исследуется в зависимости от коэффициента, характеризующего насыщение хищника, и коэффициента конкуренции хищника за отличные от жертвы ресурсы. В работе дается параметрическое описание возможных режимов динамики детерминированной модели, исследуются локальные и глобальные бифуркации и выделяются зоны устойчивого сосуществования популяций в равновесном и осцилляционном режимах. Интересной математической особенностью данной модели, впервые рассмотренной Базыкиным, является глобальная бифуркация рождения цикла из петли сепаратрисы. В работе исследуется воздействие шума на равновесный и осцилляционный режимы сосуществования популяций хищника и жертвы. Показано, что увеличение интенсивности случайных возмущений может привести к значительным деформациям этих режимов вплоть до их разрушения. Целью данной работы является разработка конструктивного вероятностного критерия близости этой стохастической системы к опасным границам. Основой предлагаемого математического подхода является техника функций стохастической чувствительности и метод доверительных областей — доверительных эллипсов, окружающих устойчивое равновесие, и доверительных полос вокруг устойчивого цикла. Размеры доверительных областей пропорциональны интенсивности шума и стохастической чувствительности исходных детерминированных аттракторов. Геометрическим критерием выхода популяционной системы из режима устойчивого сосуществования является пересечение доверительных областей и соответствующих сепаратрис детерминированной модели. Эффективность данного аналитического подхода подтверждается хорошим соответствием теоретических оценок и результатов прямого численного моделирования.

The paper is devoted to the analysis of the proximity of the population system to dangerous boundaries. An intersection of these boundaries results in the collapse of the stable coexistence of interacting populations. As a reason of such destruction one can consider random perturbations inevitably presented in any living system. This study is carried out on the example of the well-known model of interaction between predator and prey populations, taking into account both a stabilizing factor of the competition of predators for another than prey resources, and also a destabilizing saturation factor for predators. To describe the saturation of predators, we use the second type Holling trophic function. The dynamics of the system is studied as a function of the predator saturation, and the coefficient of predator competition for resources other than prey. The paper presents a parametric description of the possible dynamic regimes of the deterministic model. Here, local and global bifurcations are studied, and areas of sustainable coexistence of populations in equilibrium and the oscillation modes are described. An interesting feature of this mathematical model, firstly considered by Bazykin, is a global bifurcation of the birth of limit cycle from the separatrix loop. We study the effects of noise on the equilibrium and oscillatory regimes of coexistence of predator and prey populations. It is shown that an increase of the intensity of random disturbances can lead to significant deformations of these regimes right up to their destruction. The aim of this work is to develop a constructive probabilistic criterion for the proximity of the population stochastic system to the dangerous boundaries. The proposed approach is based on the mathematical technique of stochastic sensitivity functions, and the method of confidence domains. In the case of a stable equilibrium, this confidence domain is an ellipse. For the stable cycle, this domain is a confidence band. The size of the confidence domain is proportional to the intensity of the noise and stochastic sensitivity of the initial deterministic attractor. A geometric criterion of the exit of the population system from sustainable coexistence mode is the intersection of the confidence domain and the corresponding separatrix of the unforced deterministic model. An effectiveness of this analytical approach is confirmed by the good agreement of theoretical estimates and results of direct numerical simulations.

Безопасное проведение аварийно-восстановительных работ на аварийных морских газоконденсатных скважинах возможно при учете опасных факторов, препятствующих проведению противофонтанных мероприятий. Одним из таких факторов является загазованность района работ вследствие выхода из водной толщи большого количества легкого, по сравнению с воздухом, природного газа, а также паров более тяжелых компонентов газового конденсата (ГК). Для оценки распределения взрывоопасных концентраций паров нефтепродукта в приводном слое атмосферы необходимо определить характеристики источника загазованности. На основании анализа теоретических работ, посвященных формированию поля скорости в верхнем слое моря вследствие выхода на поверхность большого количества газа, предложена аналитическая модель для расчета размеров области, в которой происходит испарение значительного количества поступающего на поверхность ГК при авариях на мелководных скважинах. Рассматривается стационарный режим истечения пластового продукта при открытом фонтанировании газонефтяных скважин морского базирования при подводном расположении их устья. Построена малопараметрическая модель испарения нефтепродуктов из пленок различной толщины. Показано, что размер зоны интенсивного испарения ГК при подводном выбросе на мелководных скважинах определяется объемным потоком жидкой фракции ГК, его фракционным составом и выбранным порогом для оценки потока паров нефтепродукта в атмосферу. В контексте данной работы мелководными называются скважины при дебите газа от 1 до 20 млн м³ на глубинах порядка 50–200 метров. В этом случае струя пластового флюида из устья скважины на морском дне трансформируется в пузырьковый шлейф, типичная для летне-осеннего периода стратификация водной толщи не ограничивает выход шлейфа на поверхность моря, а скорость подъема пузырьков позволяет не принимать во внимание процесс растворения газа. Проведенный анализ был ограничен условиями близкими к штилевым. Такие условия благоприятны для проведения морских операций, однако неблагоприятны с точки зрения рассеяния высоких концентраций паров нефтепродуктов в приводном слое атмосферы над морем. В результате проведенной работы предложено аналитическое соотношение для приближенной оценки зоны интенсивного испарения ГК.

Safe carrying out of emergency recovery operations at emergency offshore gas condensate wells is possible when taking into account the hazardous factors that prevent anti-fontanning measures. One of such factors is the gassiness of the operation zone due to the release from the water column of a large amount of light, as compared to air, natural gas, as well as vapours of heavier components of gas condensate. To estimate the distribution of explosive concentration of petroleum product vapours in the near surface layer of the atmosphere, it is necessary to determine the characteristics of the source of the contamination. Based on the analysis of theoretical works concerning to the formation of the velocity field in the upper layer of the sea as a result of large amounts of gas coming to the surface, an analytical model is proposed to calculate the size of the area in which a significant amount of gas condensate coming to the surface is vaporised during accidents at shallow-water wells. The stationary regime of reservoir fluid flow during fountaining of offshore gas and oil wells with an underwater location of their mouths is considered. A low-parametric model of oil product evaporation from films of different thickness is constructed. It is shown that the size of the zone of intensive evaporation at shallow-water wells is determined by the volume flow of liquid fraction, its fractional composition and selected threshold for estimation of oil product vapour flow into the atmosphere. In the context of this work shallow water wells are wells with gas flow rate from 1 to 20 million cubic meters at sea depths of about 50–200 metres. In this case, the formation fluid jet from the wellhead on the seabed is transformed into a bubble plume, the stratification of the water column, typical for the summer-autumn period, does not limit the plume’s exit to the sea surface, and the velocity of bubble rise allows the gas dissolution process to be disregardded. The analysis was limited to almost calm hydrometeorological conditions. Such conditions are favourable for offshore operations, but unfavourable from the point of view of dispersion of high concentrations of oil product vapours in the near surface layer of the atmosphere. As a result of this work, an analytical dependence for an approximate assessment of the zone of intensive evaporation of gas condensate is proposed.

Рассматриваются вопросы моделирования течения в дальних областях осесимметричного безымпульсного турбулентного следа, возникающего при обтекании тела. Основными особенностями исследуемого течения являются его практически бессдвиговый характер и близость к автомодельному режиму, что подтверждается имеющимися экспериментальными данными. Для описания течения в дальнем следе привлекается полуэмпирическая $e-\varepsilon$-модель турбулентности, основанная на алгебраической модели рейнольдсовых напряжений Роди. С помощью теоретико-группового анализа получена автомодельная редукция привлекаемой модели турбулентности к системе обыкновенных дифференциальных уравнений. Проведено построение приближенного автомодельного решения второго рода соответствующей краевой задачи для редуцированной системы уравнений с помощью метода, основанного на сшивании усечений асимптотических разложений ее решения, найденных в окрестностях граничных точек. Привлекаемый подход также позволил осуществить выделение приближенного значения показателя автомодельности и других неизвестных параметров задачи в процессе выполнения процедуры сшивания асимптотических разложений. Сопоставление сконструированного в виде рядов Пюизё приближенного решения поставленной краевой задачи, включающего 16 аппроксимирующих слагаемых, с численным решением, найденным с помощью метода стрельбы, показало, что максимальная относительная погрешность между решениями не превосходит 2%. Установленное в данной работе приближенное значение показателя автомодельности задачи хорошо согласуется с точным значением, полученным ранее другими авторами в результате аналитического решения соответствующей нелинейной задачи на собственные значения. Проведено сопоставление приближенного решения с доступными экспериментальными данными по осесимметричному безымпульсному турбулентному следу, полученными в ИГиЛ СО РАН, имеет место удовлетворительное соответствие. Построенные в работе приближенные формулы имеют практическое значение и могут быть использованы для быстрой оценки значений гидродинамических величин в дальних областях безымпульсного осесимметричного турбулентного следа. Развиваемый метод построения приближенных решений краевых задач свободной турбулентности может найти применение и в других прикладных областях, которые не связаны с теорией турбулентности.

In this paper we consider the problem of modeling of flow in the far regions of an axisymmetric momentumless turbulent wake that occurs when a body is flowed. The main features of the studied flow are its almost shear-free nature and its proximity to the self-similar regime, which is confirmed by the available experimental data. A semi-empirical $e-\varepsilon$ model of turbulence based on the Rodi algebraic model of Reynolds stresses is invoked to describe a flow in the far wake. Using group-theoretic analysis, similarity reduction of the invoked model of turbulence to a system of ordinary differential equations is obtained. An approximate self-similar solution of the second kind of the corresponding boundary value problem for the reduced system of equations was constructed using a method based on matching the truncations of asymptotic expansions of its solution found in neighborhoods of the boundary points. The involved approach also allowed for the extraction of an approximate value of the self-similarity exponent and other unknown parameters of the problem during the procedure of matching of the asymptotic expansion. Comparison of the approximate solution of the given boundary value problem constructed in the form of Puiseux series, which includes 16 approximating terms, with the numerical solution found using the shooting method, showed that the maximum relative error between the solutions does not exceed 2%. The approximate value of the self-similarity exponent of the problem established in this work is in a good agreement with the exact value obtained earlier by other authors as a result of the analytical solution of the corresponding nonlinear eigenvalue problem. The approximate solution was compared with the available experimental data on the axisymmetric momentumless turbulent wake obtained at the Lavrentyev Institute of Hydrodynamics of the Siberian Branch of the Russian Academy of Sciences, and a satisfactory agreement was found. The approximate formulas developed in this work have practical significance and can be used for quick estimation of the values of hydrodynamic quantities in the far regions of the axisymmetric momentumless turbulent wake. The developed method for constructing approximate solutions to the boundary value problems of free turbulence can also be used in other applied fields that are not related to turbulence theory.