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

Abstract

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

The relevance. Sinking vertical mine shafts using artificial ground freezing requires controlling the state of the frozen wall formed around the shaft under construction. The most widespread and informative method of experimental control of the frozen wall state is borehole thermometry. The borehole thermometry data is used to adjust the parameters of the model (thermophysical properties of rocks) to recreate the temperature field in the entire volume of the frozen rock mass by means of mathematical modeling. Adjustment of the model parameters under certain conditions can lead to ambiguity in determining the model parameters, the reason for which is currently not studied. The main aim is to study the conditions for ambiguity occurrence in the calculated thermophysical properties of the rock mass as a result of the adjusting the parameters of the model according to the measurement data in the control boreholes. The object of the research is a frozen water-saturated rock mass around two mine shafts of the potash mine of the Nezhinsky mining and processing complex under construction. Research methods included the collection and processing of borehole thermometry data, the formulation of a mathematical model of heat transfer in the frozen rock mass, the numerical solution of the direct and inverse Stefan problems and the theoretical interpretation of the mismatch between the measured and calculated temperatures of the rocks. Results. We studied the mismatch between the measured and numerically calculated temperatures in the control boreholes on the phase plane of the adjusted parameters of the model - the thermal conductivities of frozen and unfrozen rocks. It is shown that the minimum of the mismatch functional under certain conditions moves over time or is not unique. The analysis of the arising ambiguity is carried out on the basis of a simplified mathematical model of the heat balance at the boundary of the phase transition front. It was found that at the ice growing stage, an unreliable determination of thermal conductivity in the frozen rock is possible due to the high temperature gradient in the rock mass near the control boreholes and the error of borehole inclinometry. At the ice holding stage, an unreliable determination of both thermal conductivities is possible, since the solution of the problem depends mainly on their ratio.