The Ordovician carbonate rocks of the Yingshan formation in the Tarim Basin have a complex pore structure owing to diagenetic and secondary structures. Seismic elastic parameters（e.g., wave velocity） depend on porosity and pore structure. We estimated the average specific surface, average pore-throat radius, pore roundness, and average aspect ratio of carbonate rocks from the Tazhong area. High P-wave velocity samples have small average specific surface, small average pore-throat radius, and large average aspect ratio. Differences in the pore structure of dense carbonate samples lead to fluid-related velocity variability. However, the relation between velocity dispersion and average specific surface, or the average aspect ratio, is not linear. For large or small average specific surface, the pore structure of the rock samples becomes uniform, which weakens squirt fl ow and minimizes the residuals of ultrasonic data and predictions with the Gassmann equation. When rigid dissolved（casting mold） pores coexist with less rigid microcracks, there are significant P-wave velocity differences between measurements and predictions.
The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.
