Lithologic Reservoirs ›› 2026, Vol. 38 ›› Issue (4): 63-76.doi: 10.12108/yxyqc.20260406

• PETROLEUM EXPLORATION • Previous Articles     Next Articles

Rock physics modeling method for overpressure reservoirs of Miocene Huangliu Formation in Ledong slope, Yinggehai Basin

CHEN Hao1,2(), HUANG Handong1,2(), ZHANG Tieming3, CUI Gang4, PENG Jiahui1,2   

  1. 1 State Key Laboratory of Petroleum Resources and ProspectingChina University of Petroleum (Beijing), Beijing 102249, China
    2 College of GeophysicsChina University of Petroleum (Beijing)Beijing 102249, China
    3 Science & Technology Information DepartmentCNPC Greatwall Drilling Engineering Co., Ltd., Beijing 100101, China
    4 No. 1 Production PlantPetroChina Huabei Oilfield Company, Renqiu 062550Hebei, China
  • Received:2025-10-17 Revised:2025-12-26 Online:2026-07-01 Published:2026-07-06
  • Contact: HUANG Handong E-mail:ch19up@163.com;webhhd@163.com

Abstract:

Existing rock physics models cannot accurately characterize the variation of elastic parameters in overpressure reservoir, which limits the precision of overpressure reservoir characterization. Taking overpressure re-servoirs of Miocene Huangliu Formation in Ledong slope of Yinggehai Basin as an example, a multi-porosity rock physics model for overpressure reservoirs was proposed, and its application performance was analyzed. Based on the proposed model, the influence factors of elastic parameters in overpressure reservoirs were analyzed, and a rock physics template incorporating formation effective stress and porosity parameters were established. The results show that: (1) The multi-porosity overpressure rock physics modeling method is proposed based on the multi-porosity theory by comprehensively considering effects of soft pores, stiff pores, and bound-water. By introducing the formation pressure coefficient and effective stress into the pore-space stiffness theory, and employing a gradient-based optimization algorithm to adaptively optimize pore aspect ratios, the model captures the variation of elastic parameters with effective stress under complex pore structures. (2) The predicted P-wave and S-wave velocities of Huangliu Formation in Ledong slope area by this model are in high agreement with well logging interpretations, with prediction errors below 8%. Compared with conventional models, the proposed model achieves the best fitting performance, with the highest R2 values of 0.920 and 0.937, respectively. (3) In overpressure reservoirs, rock elastic parameters are jointly controlled by five factors, including effective stress, soft-pore proportion, clay content, bound-water porosity, and movable-fluid porosity. In overpressure formations, with increasing pore pressure and decreasing effective stress, the P-wave to S-wave velo-city ratio (vp/vs) increases. Compared with sandstone, mudstone exhibits lower velocity and higher vp/vs values. The vp/vs ratio is the sensitive parameter for distinguishing sandstone from mudstone in overpressure clastic reservoirs. In the study area, the discriminant threshold of vp/vs for sandstone and mudstone is approximately 1.64 in the normal-pressure section, and 1.70-1.75 in the overpressure section. (4) The sandstone reservoir distribution predicted by the vp/vs section inverted on the basis of the proposed multi-porosity overpressure rock physics model is basically identical with the gas logging interpretation conclusion.

Key words: overpressure formation, effective stress, elastic parameter, pore-space stiffness, soft-pore proportion, rock physics modeling, rock physics template, Ledong slope, Yinggehai Basin

CLC Number: 

  • P631

Fig. 1

Structural units of Ledong slope (a) and comprehensive stratigraphic column of Cenozoic (b), Yinggehai Basin"

Table 1

Correspondence between soft/stiff pores and conventional pore types"

对照项 传统孔隙类型划分 本文孔隙类型划分
分类依据 孔隙成因与几何形态(粒间孔、粒内孔、溶蚀孔、微裂缝等) 孔隙纵横比及其在有效
应力作用下的力学响应
核心判据 孔隙尺度、形态与成因 压力敏感性与压缩行为
与粒间孔
关系
粒间孔为碎屑岩中
主要储集空间
孔隙纵横比较低、弱胶结的
粒间孔通常表现为柔性孔;
孔隙纵横比较高、强胶结的
残余粒间孔多表现为刚性孔
与粒内孔
关系
粒内孔尺度较小、
分布离散
多数情况下孔隙纵横比较高,表现为相对刚性孔;在弱胶
结或连通条件下亦可呈现
一定柔性
与黏土孔/
微裂缝关系
黏土片间孔、微裂缝
尺度小、形态扁平
通常孔隙纵横比较低,对
有效应力高度敏感,
表现为典型柔性孔
分类目的 描述孔隙几何特征
及成因类型
描述孔隙对有效应力变化
的力学响应及其对弹性
参数的影响

Fig. 2

Simulation the variation of stiff pores and soft pores with the depth"

Fig. 3

Thin section photographs of Miocene Huangliu Formation at different burial depths of well A in Ledong slope, Yinggehai Basin"

Fig. 4

Variation of sample bulk modulus (a) and shear modulus (b) with normal effective stress under different mudstone contents of Miocene Huangliu Formation in Ledong slope, Yinggehai Basin"

Fig. 5

Simulation of variation of soft-pore proportion with depth in normal-pressure and overpressure formations"

Fig. 6

Experimental simulation of coupled confining pressure-pore pressure effects in sandstone: Variations of P-wave and S-wave velocities (a) and vp/vs ratio (b) with pore pressure under different confining pressures"

Table 2

Formulas of four pore-space stiffness theoretical models"

作者 Kd模型 $\frac{{K}_{\varphi }}{{K}_{m}}$ μd模型
Russell[13] ${K}_{d}=\frac{{K}_{m}}{1+\frac{\varphi }{k}}$ $k=\frac{{K}_{\varphi }}{{K}_{m}}=f\left({P}_{e}\right)=M+N\cdot ln\left({P}_{e}\right)$ ${\mu }_{d}=\frac{{\mu }_{m}{K}_{d}}{{K}_{m}}$
Dinh等[14] ${K}_{d}=\frac{{K}_{m}}{1+\frac{\varphi }{k}}$ $k=\frac{{K}_{\varphi }}{{K}_{m}}=f\left({P}_{e}\right)=L+H\cdot {P}_{e}$ μd=P+Q⋅Pe
刘仕友等[16] ${K}_{d}=\frac{{K}_{m}}{1+\frac{\varphi }{k}}$ $k=\frac{{K}_{\varphi }}{{K}_{m}}=f\left({P}_{e}\right)=R+S\cdot ln\left({P}_{e}\right)$ ${\mu }_{d}=\frac{{\mu }_{m}{K}_{d}}{T{K}_{m}}$
本文模型 ${K}_{d}=\frac{{K}_{m}}{1+\frac{\varphi }{k}}$ $k=\frac{{K}_{\varphi }}{{K}_{m}}=f\left({P}_{e}, {P}_{c}\right)=\left(D+E\cdot {e}^{F\cdot {P}_{e}}\right)\cdot {P}_{c}$ ${\mu }_{d}=\frac{{\mu }_{m}\cdot {K}_{d}\cdot ln(G\cdot {P}_{c})}{{K}_{m}}$

Fig. 7

Workflow of the GBO-based multi-porosity overpressure rock physics model"

Fig. 8

Prediction of pore pressure and formation pressure coefficient of Miocene Huangliu Formation in well A, Ledong slope"

Table 3

Input parameters of framework minerals of Miocene Huangliu Formation in well A, Ledong slope"

矿物成分 体积模量/GPa 剪切模量/GPa 密度/(g·cm-3
石英 39.000 30 2.65
黏土 32.000 19 2.82
方解石 27.000 20 2.85
2.650 0 1.03
0.001 0 0.38

Fig. 9

Comparison of prediction results from different models of Miocene Huangliu Formation in well A, Ledong slope"

Fig. 10

Prediction errors comparison of different rock physics models of Miocene Huangliu Formation in well A, Ledong slope"

Fig. 11

Prediction performance and error distribution histograms of the GBO-based multi-porosity rock physics model of Miocene Huangliu Formation in well A, Ledong slope"

Table 4

Statistics of velocity prediction errors of different rock physics models on the same test dataset of Miocene Huangliu Formation in well A, Ledong slope"

模型 RMSE/(m·s-1) R2
纵波
速度
横波
速度
纵波
速度
横波
速度
Xu-White模型
(未考虑地层压力影响)
334.69 192.01 -0.824 0.147
传统孔隙空间刚度理论模型
(未考虑束缚水影响)
163.81 93.094 0.563 0.799
改进孔隙空间刚度理论模型
(未考虑束缚水影响)
146.49 74.366 0.651 0.872
多重孔隙结构超压岩石物理
模型(考虑束缚水影响)
91.938 58.919 0.862 0.919
GBO优化多重孔隙结构超压
岩石物理模型(最终优化模型)
70.220 52.176 0.920 0.937

Table 5

Parameter settings for sensitivity analysis of the multi-porosity overpressure rock physics model"

输入参数 取值
泥质含量占骨架比例/% 10~90
压力系数 0.4~2.3
有效应力/MPa 10~80
可动流体孔隙度/% 0.1~10.0
束缚水孔隙度/% 0.1~10.0
柔性孔隙占比/% 10~90
束缚水孔隙纵横比 0.05
柔性孔隙纵横比 0.08
刚性孔隙纵横比 0.50

Fig. 12

Variation of elastic parameters with soft-pore proportion and effective stress simulated by the multi-porosity overpressure rock physics model"

Fig. 13

Variation of elastic parameters with mudstone content and effective stress simulated by the multi-porosity overpressure rock physics model"

Fig. 14

Variation of elastic parameters with bound-water porosity and movable-fluid porosity simulated by the multi-porosity overpressure rock physics model"

Table 6

Statistics of effective stress and formation pressure coefficient of Miocene Huangliu Formation, Ledong slope"

有效应力/
MPa
70 65 60 55 50 45 40 35 30 25 20
地层压力
系数
0.70 0.85 0.98 1.12 1.25 1.38 1.51 1.64 1.77 1.90 2.04

Fig. 15

Rock physics template of overpressure sandstone-mudstone reservoirs in Miocene Huangliu Formation of Ledong slope"

Fig. 16

vp/vs inversion section of Miocene Huangliu Formation across well A in Ledong slope"

Fig. 17

Porosity inversion section of Miocene Huangliu Formation across well A in Ledong slope"

[1] 陈鸣, 余立文, 秦瑞, 等. 泥页岩地层孔隙压力预测新方法[J]. 复杂油气藏, 2021, 14(3):27-30.
CHEN Ming, YU Liwen, QIN Rui, et al. A new method for prediction of pore pressure in shale formation[J]. Complex Hydrocarbon, 2021, 14(3):27-30.
[2] 陈红果, 张凤奇, 江青春, 等. 松辽盆地徐家围子断陷白垩系沙河子组超压形成机制及其演化特征[J]. 岩性油气藏, 2025, 37(1):102-114.
CHEN Hongguo, ZHANG Fengqi, JIANG Qingchun, et al. Overpressure-generating mechanism and its evolution characteristics of Cretaceous Shahezi Formation in Xujiaweizi Fault Depression,Songliao Basin[J]. Lithologic Reservoirs, 2025, 37(1):102-114.
[3] 赵军, 李勇, 文晓峰, 等. 基于斑马算法优化支持向量回归机模型预测页岩地层压力[J]. 岩性油气藏, 2024, 36(6):12-22.
ZHAO Jun, LI Yong, WEN Xiaofeng, et al. Prediction of shale formation pore pressure based on Zebra Optimization Algorithm-optimized support vector regression[J]. Lithologic Reservoirs, 2024, 36(6):12-22.
[4] RUIZ F, CHENG A. A rock physics model for tight gas sand[J]. The Leading Edge, 2010, 29(12):1484-1489.
[5] 武陈月, 印兴耀, 印林杰, 等. 考虑压力的深部页岩储层地震岩石物理建模方法及应用[J]. 石油地球物理勘探, 2023, 58(4):893-901.
WU Chenyue, YIN Xingyao, YIN Linjie, et al. Rock physical modelling for deep shale reservoirs involving the influence of pressure[J]. Oil Geophysical Prospecting, 2023, 58(4):893-901.
[6] XU Shiyu, WHITE R E. A new velocity model for clay‐sand mixtures[J]. Geophysical Prospecting, 1995, 43(1):91-118.
[7] 原宏壮. 各向异性介质岩石物理模型及应用研究[D]. 青岛: 中国石油大学(华东), 2007.
YUAN Hongzhuang. Study of anisotropic rock physics model and application[D]. Qingdao: China University of Petroleum(East China), 2007.
[8] 王璞, 吴国忱. 基于自相容近似的致密储层岩石物理建模[J]. 地球物理学进展, 2015, 30(5):2233-2238.
WANG Pu, WU Guochen. The rock physics modeling for tight reservoir based on the self-consistent approximation[J]. Progress in Geophysics, 2015, 30(5):2233-2238.
[9] 朱伟. 气饱和超压页岩的岩石物理建模研究[J]. 地球物理学进展, 2020, 35(3):1116-1121.
ZHU Wei. Rock physics modeling of gas-saturated over-pressured shale[J]. Progress in Geophysics, 2020, 35(3):1116-1121.
[10] 仲保温. 碎屑岩储层压力相关岩石物理建模方法研究[D]. 青岛: 中国石油大学(华东), 2022.
ZHONG Baowen. Research on pressure related rock physical modeling method of clastic rock reservoir[D]. Qingdao: China University of Petroleum (East China), 2022.
[11] RIAHI M A, FAKHARI M G. Pore pressure prediction using seismic acoustic impedance in an overpressure carbonate reservoir[J]. Journal of Petroleum Exploration and Production Technology, 2022, 12(12):3311-3323.
[12] MAVKO G, MUKERJI T. Seismic pore space compressibility and Gassmann’s relation[J]. Geophysics, 1995, 60(6):1743-1749.
[13] RUSSELL B. A Gassmann-consistent rock physics template[J]. CSEG Recorder, 2013, 38(6):22-30.
[14] DINH H, MIRKO V D B, RUSSELL B. Pore space stiffness approach for a pressure-dependent rock-physics model[R]. Dallas: Society of Exploration Geophysicists, 2016.
[15] ZONG Zhaoyun, JIANG Man, XU Miaomiao. Influencing factor analysis of the elastic properties of shale with rock-physical model including pressure effects[J]. Interpretation, 2020, 8(3):T515-T524.
[16] 刘仕友, 仲保温, 周凡, 等. 高压背景砂岩储层岩石物理建模方法与应用[J]. 地球物理学进展, 2022, 37(4):1567-1575.
LIU Shiyou, ZHONG Baowen, ZHOU Fan, et al. Rock physical modeling method and application of sandstone reservoir in high pressure background[J]. Progress in Geophysics, 2022, 37(4):1567-1575.
[17] 姜曼. 考虑压力作用的地震岩石物理建模方法研究[D]. 青岛: 中国石油大学(华东), 2021.
JIANG Man. Rock physical models involving the influence of the pressure[D]. Qingdao: China University of Petroleum (East China), 2021.
[18] 杨骐羽, 李景叶, 吴凡, 等. 考虑软矿物纵横比的页岩岩石物理建模及其应用[J]. 物探与化探, 2024, 48(1):98-104.
YANG Qiyu, LI Jingye, WU Fan, et al. A petrophysical model of shales considering soft-mineral aspect ratios and its application[J]. Geophysical and Geochemical Exploration, 2024, 48(1):98-104.
[19] CHEN Fubin, ZONG Zhaoyun, YIN Xingyao. Rock elasticity and crack property descriptions using an improved acoustoelasticity model with nonlinear stress-strain relationship[J]. Geoenergy Science and Engineering, 2024, 241:213211.
[20] 刘为, 杨希冰, 张秀苹, 等. 莺歌海盆地东部黄流组重力流沉积特征及其控制因素[J]. 岩性油气藏, 2019, 31(2):75-82.
LIU Wei, YANG Xibing, ZHANG Xiuping, et al. Characteristics and controlling factors of gravity flow deposits of Huangliu Formation in eastern Yinggehai Basin[J]. Lithologic Reservoirs, 2019, 31(2):75-82.
[21] 郝芳, 李思田, 龚再升, 等. 莺歌海盆地底辟发育机理与流体幕式充注[J]. 中国科学D辑:地球科学, 2001, 31(6):471-476.
HAO Fang, LI Sitian, GONG Zaisheng, et al. Mechanisms of diapir development and episodic fluid charging in the Yinggehai Basin[J]. Science in China Series D:Earth Science, 2001, 31(6):471-476.
[22] 刘平, 童传新, 曾小明, 等. 莺歌海盆地莺东斜坡乐东区中深层地层水特征及其与油气关系[J]. 海洋地质前沿, 2019, 35(8):44-50.
LIU Ping, TONG Chuanxin, ZENG Xiaoming, et al. Characte-ristics of stratigraphic water in the middle and deep layers of Ledong District,Yingdong Slope,Yinggehai Basin and its relationship with oil and gas[J]. Marine Geology Frontiers, 2019, 35(8):44-50.
[23] 李伟, 刘平, 艾能平, 等. 莺歌海盆地乐东地区中深层储层发育特征及成因机理[J]. 岩性油气藏, 2020, 32(1):19-26.
LI Wei, LIU Ping, AI Nengping, et al. Development characteristics and genetic mechanism of mid-deep reservoirs in Ledong area,Yinggehai Basin[J]. Lithologic Reservoirs, 2020, 32(1):19-26.
[24] 范彩伟, 贾茹, 柳波, 等. 莺歌海盆地中央坳陷带成藏体系的盖层评价及控藏作用[J]. 岩性油气藏, 2023, 35(1):36-48.
FAN Caiwei, JIA Ru, LIU Bo, et al. Caprock evaluation and its reservoir control of different accumulation systems in central depression zone of Yinggehai Basin[J]. Lithologic Reservoirs, 2023, 35(1):36-48.
[25] 杨计海, 黄保家. 莺歌海凹陷东斜坡L气田天然气成因及运移模式[J]. 石油勘探与开发, 2019, 46(3):450-460.
YANG Jihai, HUANG Baojia. Origin and migration model of natural gas in L gas field,eastern slope of Yinggehai Sag,China[J]. Petroleum Exploration and Development, 2019, 46(3):450-460.
[26] 李绪深, 杨计海, 范彩伟, 等. 南海北部海域高温超压天然气勘探新进展与关键技术:以莺歌海盆地乐东斜坡带为例[J]. 中国海上油气, 2020, 32(1):23-31.
LI Xushen, YANG Jihai, FAN Caiwei, et al. New progress and key technologies for high temperature and overpressure natural gas exploration in the northern part of South China Sea:Taking the Ledong slope belt of Yinggehai Basin as an example[J]. China Offshore Oil and Gas, 2020, 32(1):23-31.
[27] 王立锋, 宋瑞有, 陈殿远, 等. 莺歌海盆地D13区新近系黄流组大型海底扇地震识别及含气性预测[J]. 岩性油气藏, 2022, 34(4):42-52.
WANG Lifeng, SONG Ruiyou, CHEN Dianyuan, et al. Seismic identification and gas-bearing prediction of large-scale submarine fans of Neogene Huangliu Formation in D13 area of Yinggehai Basin[J]. Lithologic Reservoirs, 2022, 34(4):42-52.
[28] SHAPIRO S A, TROYAN V N. Stress dependences of velocities in porous and fractured rocks[R]. Horten,64th EAGE Conference & Exhibition, 2002.
[29] ATHY L F. Density,porosity,and compaction of sedimentary rocks[J]. AAPG Bulletin, 1930, 14(1):1-24.
[30] SHAPIRO S A. Elastic piezosensitivity of porous and fractured rocks[J]. Geophysics, 2003, 68(2):482-486.
[31] CHRISTENSEN N I, WANG H F. The influence of pore pressure and confining pressure on dynamic elastic properties of Berea sandstone[J]. Geophysics, 1985, 50(2):207-213.
[32] AGHLI G, MOGHADDAS N H, SABERI M R, et al. The effect of varying pore pressure on carbonate reservoirs elastic response available to purchase[J]. Geophysics, 2025, 90(3):M45-M58.
[33] 吴亚宁, 黄捍东, 徐海, 等. 基于Wolfe搜索机制的随机梯度地震反演方法[J]. 地球物理学报, 2024, 67(11):4309-4324.
WU Yaning, HUANG Handong, XU Hai, et al. Seismic inversion method based on SGO-WS[J]. Chinese Journal of Geophysics, 2024, 67(11):4309-4324.
[34] 刘海钰, 胡林, 刘兵, 等. 琼东南盆地深水区中新统梅山组超压成因与天然气成藏特征[J]. 岩性油气藏, 2025, 37(5):70-82.
LIU Haiyu, HU Lin, LIU Bing, et al. Origin of overpressure and natural gas accumulation characteristics of Miocene Mei-shan Formation in the deepwater area,Qiongdongnan Basin[J]. Lithologic Reservoirs, 2025, 37(5):70-82.
[35] 李中, 刘和兴, 李磊, 等. 基于改进的Bowers法预测南海陵水深水高温地层异常压力[J]. 西安石油大学学报(自然科学版), 2019, 34(6):60-66.
LI Zhong, LIU Hexing, LI Lei, et al. Prediction of abnormal pressure in Lingshui deep-water high temperature formation based on improved bowers method[J]. Journal of Xi’an Shiyou University (Nature Science Edition), 2019, 34(6):60-66.
[1] REN Yufei, YAN Jianping, YAN Ke, HUANG Lisha, WANG Min, GENG Bin. Quantitative characterization method for clastic rock components based on improved electrical imaging logging mineral probability spectrum [J]. Lithologic Reservoirs, 2026, 38(4): 91-100.
[2] CHEN Zhihong, LYU Zhengxiang, HU Gaowei, GAO Mengtian, CHEN Yabing, JIN Feng, MENG Hongyu. Hydrocarbon accumulation condition and main controlling factors of Miocene Sanya Formation in Yingdong slope area, Yinggehai Basin [J]. Lithologic Reservoirs, 2025, 37(6): 1-12.
[3] FAN Caiwei, JIA Ru, LIU Bo, FU Xiaofei, HOU Jingxian, JIN Yejun. Caprock evaluation and its reservoir control of different accumulation systems in central depression zone of Yinggehai Basin [J]. Lithologic Reservoirs, 2023, 35(1): 36-48.
[4] YANG Kaile, HE Shenglin, YANG Zhaoqiang, WANG Meng, ZHANG Ruixue, REN Shuangpo, ZHAO Xiaobo, YAO Guangqing. Diagenesis characteristics of tight sandstone reservoirs with high temperature,overpressure and high CO2 content: A case study of Neogene Meishan-Huangliu Formation in LD10 area,Yinggehai Basin [J]. Lithologic Reservoirs, 2023, 35(1): 83-95.
[5] WANG Lifeng, SONG Ruiyou, CHEN Dianyuan, XU Tao, PAN Guangchao, HAN Guangming. Seismic identification and gas-bearing prediction of large-scale submarine fans of Neogene Huangliu Formation in D13 area of Yinggehai Basin [J]. Lithologic Reservoirs, 2022, 34(4): 42-52.
[6] LI Wei, LIU Ping, AI Nengping, SHAO Yuan, HOU Jingxian. Development characteristics and genetic mechanism of mid-deep reservoirs in Ledong area,Yinggehai Basin [J]. Lithologic Reservoirs, 2020, 32(1): 19-26.
[7] LIU Wei, YANG Xibing, ZHANG Xiuping, DUAN Liang, SHAO Yuan, HAO Defeng. Characteristics and controlling factors of gravity flow deposits of Huangliu Formation in eastern Yinggehai Basin [J]. Lithologic Reservoirs, 2019, 31(2): 75-82.
[8] Han Guangming,Pan Guangchao,Fu Chen,Luo Qi,Shao Yuan,Wang Rui. Influence of velocity changing of gas reservoir and seal on seismic response and AVO type [J]. Lithologic Reservoirs, 2016, 28(2): 107-113.
[9] QU Jianhua, WANG Zesheng, REN Benbing, BAI Yu, WANG Bin. Genetic mechanism analysis and prediction method of abnormal high pressure in Mahu slope area, Junggar Basin [J]. Lithologic Reservoirs, 2014, 26(5): 36-39.
[10] HAN Guangming, ZHOU Jiaxiong, PEI Jianxiang, WANG Lifeng. Essence of diapir and its relationship with natural gas accumulation in Yinggehai Basin [J]. Lithologic Reservoirs, 2012, 24(5): 27-31.
[11] ZHENG Hongjun, CAO Zhenglin, YAN Cunfeng, XU Ziyuan, SUN Songling. Using seismic rock physics modeling to predict accumulation unit of biogenetic gas in Sanhu area [J]. Lithologic Reservoirs, 2010, 22(4): 20-24.
[12] LIU Wenlin. Geophysical response characteristics of coal bed methane [J]. Lithologic Reservoirs, 2009, 21(2): 113-115.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!
TRENDMD: