A MATLAB program for 1D strain rate inversion

  • Authors:

  • Hai-Bin Song;Lin Chen;Jiong Zhang;Chang-Yu Zhao;Chong-Zhi Dong

  • Affiliations:

  • Key Laboratory of Petroleum Resources Research, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China;Key Laboratory of Petroleum Resources Research, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China and Graduate University of Chinese Academy of Sciences, Beij ...;School of Geophysics & Information Technology, China University of Geosciences, Beijing 100083, China;School of Geophysics & Information Technology, China University of Geosciences, Beijing 100083, China;Key Laboratory of Petroleum Resources Research, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China and Graduate University of Chinese Academy of Sciences, Beij ...

  • Venue:
  • Computers & Geosciences
  • Year:
  • 2010

Quantified Score

Hi-index 0.00

Visualization

Abstract

This paper presents a MATLAB program designed to invert 1D strain rate from subsidence data. In forward modeling, we use an implicit finite difference scheme to solve the heat conduction equation with an advective term. In the inversion, we adopt the Powell algorithm to continually search for the optimal values of strain rate until the fit, defined by the difference between the calculated subsidence and the observed subsidence, is satisfactory. Synthetic datasets are generated, and one of them is used to test the inversion algorithm. The results show that the calculated subsidence fit the theoretical subsidence quite well, and the inverted strain rate oscillates around the true value and is a good approximation to the original strain rate variation. The program is applied to the northern continental margin of the South China Sea, and the inverted strain rate from WC1411 well reveals the multiple rifting events that occurred in this region. The inverted strain rate can be used to evaluate the stretching factor and provides constraints for dynamic modeling of lithospheric deformation.