Lg波衰减模型中建模误差的定量评估

Quantitative evaluation of modeling error in Lg-wave attenuation model

  • 摘要: 对Lg波衰减模型中建模误差的统计特征进行了系统研究,并建立了地壳二维Lg波衰减模型。由于Lg波振幅可能受到几何扩散函数的强烈影响,合理评估反演过程的误差对于能否使用最小二乘意义下的反演非常重要。通过在川滇及其邻近地区收集的建模误差样本,使用K-S数值检验方法、Q-Q图和正态分布图形检验方法对Lg波衰减层析成像反演的输入数据中建模误差的分布特征进行了统计分析。采用奇异值分解(SVD)和反投影方法,分别获得了川滇地区的QLg模型,定量计算模型的协方差矩阵和分辨率矩阵,定量评估了QLg模型中每个格点的分辨率和误差。结果表明:在一阶近似条件下建模误差服从正态分布;通过开发的数据筛选程序,可以产生一个接近完美正态分布的数据集;与反投影方法相比,利用SVD方法获得的地壳Q值的分辨率更高;在射线覆盖较好的区域,QLg模型的分辨率达到100 km,相对误差小于3%。

     

    Abstract: In this study, we performed comprehensive statistical studies for the characteristics of the modeling error in Lg-wave attenuation model and obtained a two-dimensional crustal Lg-wave attenuation model. The amplitudes of Lg-wave can be strongly influenced by the geometrical spreading function and it is essential to reasonably evaluate the error of the inversion process if we use the inverted method under the meaning of least squares. Based on the modeling error samples collected in the Sichuan-Yunnan region and its adjacent areas, we use three statistical test methods, namely the K-S test, the Q-Q plot, and the normal distribution plot, to analyze the distribution characteristics of modeling errors in the input data of Lg- wave attenuation tomography inversion. We used both the SVD-based tomographic method and the back-projection method to invert for the Lg Q model of the Sichuan-Yunnan region and its adjacent area respectively. Then, we calculated the covariance matrix and the resolution matrix of the model and evaluated the resolution and error of each grid point in the Lg Q model quantitatively. Our results indicate that the modeling errors in the input data obeyed a normal distribution under first-order approximation. By using the data screening technique, we generated a new dataset with nearly perfect normally distributed for Q tomography. Compared with the result of back-projection, the SVD-based method could allow for obtaining a finer resolution of crustal Q value in the areas with dense ray path coverages, where the model resolution can be resolved within a range of 100 km with a relative error of less than 3%.

     

/

返回文章
返回