3-D paleotemperature modeling of the geothermal regime of sedimentary basins: Example of the Lunskaya depression,Sakhalin Island

The task of 3-D modeling of the thermal field of a sedimentary basin during sedimentation is considered. The aim of the modeling is to determine the temperature at any point of the basin at a given moment of geological time. The mathematical model is based on a system of equations of thermal conductivity for a heterogeneous layered medium with dynamic boundaries. The conditions of the continuous temperature and thermal flow are given at the boundaries of the adjacent layers. The temperature values, which are determined by the values of the secular course of the earth temperature, are given at the upper boundary coinciding with the sedimentation surface. The thermal flow value is considered to be given at the lower boundary. The medium is approximated using a vertical triangle prism, which is accepted in algorithms of interpretation of the gravitation field and characterized by random upper and lower basements and given values of the thermal physical parameters. The equations of thermal conductivity are solved on the basis of potential theory. The precision of this algorithm is demonstrated by calculation of a test example. The thermal evolution of the sedimentary complexes and dynamics of the major zone of oil formation are reconstructed and possible errors of paleotemperature interpretations caused by ignored 3-D modeling medium are determined on the example of the sedimentary basin of the Lunskaya depression of Sakhalin.