WaveRock軟件是研究3D數字巖石寬（全）頻帶地震響應，建立儲層巖石巖性、物性參數及含油氣性與地震屬性之間定量關系的一種新的有效工具，可有效提高地震巖性識別、儲層預測和流體檢測的精度和可靠性。
　WaveRock軟件具有下列優越的性能：
1、采用獨特的變密度交錯網格優化高階有限差分算法和多核并行算法以及PML吸收邊界條件，模擬精度高，邊界吸收效果理想，計算速度快；
2、采用2D/3D全彈性波方程、2D/3D孔隙介質Biot/BISQ方程；
3、采用多種3D數字巖石建模方法（多組分規則形狀均勻分布、隨機分布、非均勻分形分布、巖心薄片X-CT掃描圖）以及礦物和巖石模型庫來快速建模；
4、快速、高精度測定碎屑巖模型、多孔隙（含裂隙）介質模型、多相非均勻飽和孔隙介質模型的寬（全）頻帶速度和Q值?？梢耘c巖石物理測試結果、測井數據以及各種理論公式對比校驗，并建立數據庫；
5、能夠高精度模擬分析孔隙尺度下不同孔隙結構和流體分布的微波場特征，提取巖性、物性和流體敏感地震屬性；
6、能夠高精度模擬不同尺度數字巖石跨頻帶彈性波響應，優選敏感地震屬性，建立儲層巖性、物性及含油氣性與地震屬性之間的定量預測模型。

北方林泰公司與中科院地質-地球物理研究所合作取得重要成果，該成果發表在

2014年地震學報（英文版）?！禘arthquake Science》，June 2014, Volume 27, Issue 3, pp 285-299

Poroelastic finite-difference modeling for ultrasonic waves in digital porous cores
Li-Yun Fu  Yan Zhang Zhenglin Pei  Wei Wei  Luxin Zhang

Abstract

Scattering attenuation in short wavelengths has long been interesting to geophysicists. Ultrasonic coda waves, observed as the tail portion of ultrasonic wavetrains in laboratory ultrasonic measurements, are important for such studies where ultrasonic waves interact with small-scale random heterogeneities on a scale of micrometers, but often ignored as noises because of the contamination of boundary reflections from the side ends of a sample core. Numerical simulations with accurate absorbing boundary can provide insight into the effect of boundary reflections on coda waves in laboratory experiments. The simulation of wave propagation in digital and heterogeneous porous cores really challenges numerical techniques by digital image of poroelastic properties, numerical dispersion at high frequency and strong heterogeneity, and accurate absorbing boundary schemes at grazing incidence. To overcome these difficulties, we present a staggered-grid high-order finite-difference (FD) method of Biot’s poroelastic equations, with an arbitrary even-order (2L) accuracy to simulate ultrasonic wave propagation in digital porous cores with strong heterogeneity. An unsplit convolutional perfectly matched layer (CPML) absorbing boundary, which improves conventional PML methods at grazing incidence with less memory and better computational efficiency, is employed in the simulation to investigate the influence of boundary reflections on ultrasonic coda waves. Numerical experiments with saturated poroelastic media demonstrate that the 2L FD scheme with the CPML for ultrasonic wave propagation significantly improves stability conditions at strong heterogeneity and absorbing performance at grazing incidence. The boundary reflections from the artificial boundary surrounding the digital core decay fast with the increase of CPML thicknesses, almost disappearing at the CPML thickness of 15 grids. Comparisons of the resulting ultrasonic coda Qsc values between the numerical and experimental ultrasonic S waveforms for a cylindrical rock sample demonstrate that the boundary reflection may contribute around one-third of the ultrasonic coda attenuation observed in laboratory experiments.

Keywords
Digital porous cores – Ultrasonic coda – Poroelastic finite-difference modeling – Unsplit convolutional PML absorbing boundary