Processing math: 100%

Lietard天然裂缝宽度预测模型求解新方法

彭浩, 李黔, 尹虎, 唐志强

彭浩, 李黔, 尹虎, 唐志强. Lietard天然裂缝宽度预测模型求解新方法[J]. 石油钻探技术, 2016, 44(3): 72-76. DOI: 10.11911/syztjs.201603013
引用本文: 彭浩, 李黔, 尹虎, 唐志强. Lietard天然裂缝宽度预测模型求解新方法[J]. 石油钻探技术, 2016, 44(3): 72-76. DOI: 10.11911/syztjs.201603013
PENG Hao, LI Qian, YIN Hu, TANG Zhiqiang. A New Solution Method for the Lietard Natural Fracture Width Prediction Model[J]. Petroleum Drilling Techniques, 2016, 44(3): 72-76. DOI: 10.11911/syztjs.201603013
Citation: PENG Hao, LI Qian, YIN Hu, TANG Zhiqiang. A New Solution Method for the Lietard Natural Fracture Width Prediction Model[J]. Petroleum Drilling Techniques, 2016, 44(3): 72-76. DOI: 10.11911/syztjs.201603013

Lietard天然裂缝宽度预测模型求解新方法

基金项目: 

国家科技重大专项课题"钻井工程设计和工艺软件"(编号:2011ZX05021-006)资助。

详细信息
    作者简介:

    彭浩(1988-),男,四川达州人,2011年毕业于西南石油大学石油工程专业,在读硕士研究生,主要从事井漏分析及钻井工程设计研究工作。E-mailpenghao323@qq.com。

  • 中图分类号: TE258

A New Solution Method for the Lietard Natural Fracture Width Prediction Model

  • 摘要: 采用图版法求解Lietard天然裂缝宽度预测模型时速度慢、易产生人为误差,为迅速准确地预测井漏时的天然裂缝宽度,给堵漏作业和屏蔽暂堵作业提供决策依据,开展了数值解法研究。研究Lietard模型发现,用Lietard模型计算得到的理论漏失特征曲线存在近似直线段,可作线性化处理,因此,基于井漏时实测得到的漏失数据,采用最小二乘法建立实钻漏失特征曲线近似直线段线性参数的计算模型,并运用自适应搜索法找出与实钻漏失特征曲线唯一对应的理论漏失特征曲线,根据该曲线的无因次有限侵入值反演出天然裂缝宽度。经过试验研究,建立了实钻漏失特征曲线近似直线段线性参数的计算模型,实现了无盲区自动匹配,将匹配误差控制在0.001%以内,克服了人为误差,并能快速完成Lietard模型求解。研究结果表明,Lietard模型的求解精度与实钻漏失特征曲线近似直线段的线性相关度呈正比关系,在同一相关度下,新方法比图版法求解速度更快、精度更高。
    Abstract: When the Lietard natural fracture width prediction model is solved by means of the chart method, the solving velocity is low and human errors tend to occur. In order to quickly and accurately predict the natural width with lost circulation and provide a decision-making basis for plugging operations and shield temporary plugging operations, the numerical solution method for Lietard model was studied in this paper. It was shown from the analysis on Lietard model that there was an approximate straight line section in the theoretical leakage characteristics curve obtained from Lietard model and it could be linearized. Based on the measured leakage data during the lost circulation, the calculation model for the linear parameters of approximate straight line section in the actual drilling leakage characteristics curve was established by means of the least square method. After the theoretical leakage characteristics curve which was the only one corresponding to the actual drilling leakage characteristics curve was identified by using the adaptive search method, natural fracture width could be inversed on the basis of the dimensionless finite invasion factor. Based on experimental studies, the calculation model for the linear parameters of approximate straight line section in the actual drilling leakage characteristics curve was built up and automatic matching without blind area was realized with matching error less than 0.001%. And furthermore, human errors were avoided and the Lietard mode could be solved quickly. It was shown that the solution accuracy of Lietard model was proportional to the linear correlation of approximate straight line section in the actual drilling leakage characteristics curve. And for the same correlation, new method was faster and more accurate than the chart method.
  • 目前,海上钻井几乎都采用随钻测井,陆上钻井采用随钻测井的比例也在不断提高[1]。其中,随钻电阻率成像测井仪器能提供高分辨率的井壁图像,从而可以较准确地判断出裂缝、孔洞等地质构造,实现对地层的准确评价和实时地质导向[26],在水平井钻井中发挥着重要作用。国外的随钻电阻率成像测井理论和测井仪器已经趋于成熟,Schlumberger、Halliburton和Baker Hughes等公司都推出了随钻电阻率成像测井仪器[713],这些成像测井仪器的钮扣电极在纵向上分布1~3排,而周向上钮扣电极布置的较少,造成仪器周向扫描时间较长。目前,国内在随钻电阻率成像测井理论研究和仪器研制方面尚处于起步阶段,有必要借助数值模拟手段对随钻电阻率成像测井进行研究。

    笔者设计了一种新的测井仪器钮扣电极系分布方案,并增加了测量侧向电阻率和钻头电阻率的功能,不仅缩短了测量时间,同时具有2种不同探测深度的电阻率成像、地层评价和地质导向功能,是一种高分辨率、多模式、多参数和近钻头的测量方案。为了分析该仪器方案的探测特征,借助有限元模拟平台,考察了其在复杂层状地层、周向异常体地层和水平井地层中的测井响应特征。

    随钻电阻率成像测井仪器有2种激励机制:一种是直接给电极加载电流;另一种是通过螺绕环激励在钻铤上产生等电位,以达到自动聚焦的作用。第二种方法在工艺上容易实现,因此被广泛应用[14],笔者的仪器方案也应用该原理。假定钻铤在井轴方向上无限长,井轴与柱面坐标系Z轴一致,发射螺绕环可以等效为长度磁矩的理想化磁环[15],如图1所示。

    图  1  发射螺绕环等效为理想化磁环示意
    Figure  1.  Schematic of a launching spiral ring that is equivalent to the idealized magnetic ring

    实际测量过程中由于测量频率低,可以忽略频率的影响,因此可以将螺绕环等效为延长的电压偶极子[1617]。此时,测量原理与传统侧向测井类似,采用欧姆定律对视电阻率进行标定。视电阻率的计算公式为:

    Ra=KUI (1)

    式中:Ra为视电阻率,Ω·m;K为仪器常数;U为螺绕环两端的电压,V;I为纽扣电极和接收螺绕环接收到的电流,A。

    根据电磁场原理,可得到特定仪器在空间均匀场内的响应,但是实际测井环境复杂,具有明显的非均质性,径向上由井眼、侵入带、过渡带和原状地层组成,而纵向上由目的层和围岩组成,很难利用解析方法求解如此复杂的地层模型,需要借助数值方法。因此,利用COMSOL Multiphysics有限元软件建立水平层状地层、异常体地层和水平井地层等3种地层模型,进行复杂地层的数值模拟。

    数值模拟验证的仪器由1个发射螺绕环、2排钮扣电极(R4、R5周向相隔90°,各分布4个钮扣电极)和2个接收螺绕环组成(见图2),可以测量不同深度的电阻率、侧向电阻率和钻头电阻率,对不同方位钮扣电极的测量结果进行加权平均可以获得浅侧向电阻率和中侧向电阻率。

    图  2  仪器结构示意
    Figure  2.  Structure of the instrument

    在确定源距和钮扣电极直径之前,需要考察二者对测量电流的影响,以确定最优的仪器结构参数。模拟时,发射螺绕环两端电压U取0.1 V,地层电阻率Rt的变化范围为0.1~1 000.0 Ω·m,钮扣电极与发射螺绕环之间的距离(源距)为0.10~1.50 m,钮扣电极直径为10.0 mm,不考虑井眼的影响,模拟结果如图3所示。从图3可以看出:随着源距增大,测量电流信号的变化幅度越来越小,最后基本趋于稳定;不同地层电阻率下的测量信号随源距变化趋势基本相同;源距相同时,测量电流与电阻率呈反比关系。

    图  3  源距对测量电流的影响
    Figure  3.  Effect of source distance on measured current

    同理,模拟了钮扣电极直径对测量电流信号的影响,源距为0.508 m,钮扣电极直径的变化范围为5.0~50.0 mm,其他模拟参数与图3相同,结果如图4所示。从图4可以看出:随着钮扣电极直径增大,测量电流信号在双对数坐标中呈线性增大趋势;不同地层电阻率下的测量信号随钮扣电极直径变化的趋势基本相同;纽扣电极直径相同时,测量电流与电阻率呈反比关系。

    图  4  钮扣电极直径对测量电流的影响
    Figure  4.  Effect of button electrode diameter on measured current

    对比图3图4可以发现,源距对测量电流信号的影响较小,因此可以灵活选取。钮扣电极直径对测量信号影响较大,可综合其他因素选取。国外测井仪器测量结果表明,钮扣电极直径较小时,其纵向分辨率较高,但只能探测电阻率为几百欧姆米的地层;适当增大钮扣电极直径,虽然降低了其纵向分辨率,但增大了其探测地层电阻率的范围,可以探测电阻率为几千欧姆米的地层。因此,综合考虑钮扣电极测量地层电阻率的范围、钮扣电极纵向分辨率和测量信号等3个因素,设计了2种不同直径的钮扣电极。

    综上,最终选取图2中的仪器结构和以下参数进行模拟:钮扣电极R4的直径为10.0 mm,为高分辨率钮扣电极,钮扣电极R5直径为25.4 mm,为标准钮扣电极。发射螺绕环与纽扣电极R4的距离LTR4为0.508 m,发射螺绕环与纽扣电极R5的距离LTR5为1.016 m,用于测量深侧向电阻率和钻头电阻率2个螺绕环间的距离Lr为0.381 m。由于侧向电阻率和钻头电阻率的测量原理和测井响应在文献[1617]中均有介绍,下面主要研究钮扣电极测量模式的测井响应特征。

    由于随钻电阻率成像测井仪钮扣电极的直径较小,因此可以分辨较薄的地层。为了研究上述结构仪器对地层的纵向分辨能力,建立了14层的水平层状地层,每层地层坐标、厚度和地层电阻率属性如表1所示。

    表  1  水平层状地层模型参数
    Table  1.  The model parameters of horizontally layered strata
    编号纵向坐标/m地层厚度/m地层电阻率/(Ω·m)
    1–100.000100.00010
    200.005100
    30.0050.00510
    40.0100.010100
    50.0200.01010
    60.0300.020100
    70.0500.02010
    80.0700.040100
    90.1100.04010
    100.1500.060100
    110.2100.06010
    120.2700.080100
    130.3500.08010
    140.43099.57010
    下载: 导出CSV 
    | 显示表格

    利用COMSOL Multiphysics 有限元软件模拟水平层状地层的结果如图5所示。模型第1层为巨厚层,因此没有显示,图5中只显示了从第2层到第13层及第14层的部分地层的测井响应。由于钮扣电极R5的直径为钮扣电极R4直径的2.54倍,明显地,钮扣电极R4对地层的分辨率高于钮扣电极R5。当地层厚度达到0.01 m时,钮扣电极R4的视电阻率接近模型值,而对于钮扣电极R5,当地层厚度达到0.02 m时,其视电阻率才开始接近模型值。当地层厚度大于0.02 m后,R4和R5均可以分辨地层,通过对比钮扣电极直径和其纵向分辨率可以发现,其对地层的分辨率大致为钮扣电极的直径尺寸。同时,由于模拟中只考虑了1个发射螺绕环的情况,没有对视电阻率进行补偿,因此模拟得到的测量曲线和地层模型不对称,在靠近上下地层界面处,电阻率出现“一高一低”的情况。

    图  5  测井仪器在水平层状地层的测井响应
    Figure  5.  Logging response of the logging instrument in horizontally layered strata

    为了考察仪器的周向探测特性,建立了含有方位地层的周向异常体地层模型(见图6),通过改变异常体张开角度来考察仪器的方位探测特性。计算模型由仪器结构、井眼、地层和异常体组成。异常体初始位置位于正北方向,张开角度θ的变化范围为0°~360°,异常体厚度为无限厚,分布于井眼之外。井眼直径Dh为215.9 mm,钻井液电阻率Rm为0.1 Ω·m,地层电阻率Rt为1.0 Ω·m,异常体电阻率Rb为100.0 Ω·m。

    图  6  含方向性异常体的地层模型
    Figure  6.  Stratigraphic model with directional anomalous bodies

    以钮扣电极R4为例,模拟异常体张开角度从0°变化到360°时不同方位的测井响应,结果如图7所示。图7中,RN4代表R4位于正北方向的钮扣电极,RE4、RS4、RW4分别代表R4位于正东、正南、正西方位上的钮扣电极。从图7可以看出:当位于正北方向的异常体张开角度从0°到90°增大时(从正北方向两侧对称增大),正北方位钮扣电极测量的视电阻率呈线性增大,从90°到135°缓慢接近异常体电阻率;考虑到方位钮扣电极分布的对称性,RE4和RS4视电阻率曲线重合,当异常体张开角度从0°增至135°时,该方位钮扣电极对异常体几乎没有识别能力;当异常体张开角度从135°增至270°时,正东方向钮扣电极的视电阻率基本呈线性增大,此后视电阻率随异常体张开角度增大保持不变。对比而言,位于正南方向的钮扣电极由于距离异常体较远,因此对异常体的识别度较低,当异常体角度大于315°时,其视电阻率才开始增大,并接近异常体电阻率。

    图  7  钮扣电极视电阻率与异常体张开角度的关系曲线
    Figure  7.  The relationship curve between the apparent resistivity of the button electrode and the anomalous body opening angle

    由于R2测量的侧向视电阻率和R3测量的钻头视电阻率没有方位探测特性,异常体张开角度为0°时(即不考虑异常体),二者的视电阻率接近地层真电阻率,为1.0 Ω·m(见图8),可以看出曲线略微受到井眼的影响,其中钻头视电阻率受井眼的影响较严重。当异常体张开角度从0°到360°变化,侧向和钻头的视电阻率均随异常体张开角度增大而增大,但侧向视电阻率略大于钻头视电阻率。当异常体张开角度增加到360°时,侧向和钻头的视电阻率接近异常体的电阻率(仍受到井眼的影响)。对比图7图8可以看出,钮扣电极与侧向电阻率测量电极、钻头电阻率测量电极对异常体的灵敏度不同,钮扣电极可以分辨较小张开角度的异常体,而侧向电阻率测量电极和钻头电阻率测量电极则无法检测较小张开角度的异常体,因此在测井解释方面,可以利用方位钮扣电极测量结果识别方位性高阻储层。

    图  8  深侧向与钻头视电阻率与异常体张开角度的关系曲线
    Figure  8.  The relationship curve between the apparent resistivity of deep laterolog/bit and the anomalous body opening angle

    随钻电阻率成像测井相比于常规电阻率成像测井的优势是其可以应用于大斜度井和水平井,为了考察仪器在水平井中的测井响应,建立了如图9所示的水平井地层模型。该模型由3层地层组成,上下层为围岩,电阻率Rs为1 Ω·m,中间层为目的层,电阻率Rt为10 Ω·m,仪器位于目的层中,目的层厚度H为2 m,仪器初始位置位于目的层中间,坐标Z为0,向上靠近地层界面Z值为正,向下靠近地层界面Z值为负。

    图  9  水平井数值模拟模型示意
    Figure  9.  The model of horizontal well numerical simulation

    钮扣电极R4测量的水平井中不同方位视电阻率与仪器距离地层界面距离的关系如图10所示,R4正北方向和正南方向的钮扣电极靠近地层界面,而正东和正西方向的钮扣电极与地层界面垂直。从图10可以看出:正北方向和正南方向钮扣电极的视电阻率曲线与仪器在直井中的测井响应曲线类似,当仪器靠近地层界面处时,由于电荷的累积,具有“犄角”现象,仪器离开地层界面时也是如此;仪器在地层上下界面处的测井响应不对称;相比而言,正东方向和正西方向钮扣电极的测量曲线几乎重合,具有良好的对称性。

    图  10  钮扣电极测量的视电阻率与仪器距地层界面距离的关系曲线
    Figure  10.  The relationship curve between the apparent resistivity measured by the buttonelectrode and the distance of the instrument to strata interface

    将R4和R5不同方位的钮扣电极测量的视电阻率进行加权平均,可以获得不同径向探测深度的浅侧向电阻率和深侧向电阻率,可以用于地层评价。计算结果表明,浅侧向R4、中侧向R5和深侧向R2视电阻率相差不大,且关于地层对称(见图11)。该结果与H. M. Wang等人[18]模拟的双侧向结果类似,测量的钻头视电阻率也关于地层模型对称,但是其测量值远远小于目的层的真实电阻率。

    图  11  侧向、钻头测量的视电阻率与仪器距地层界面距离的关系曲线
    Figure  11.  The relationship curve between the apparent resistivity measured by the laterolog/bit and the distance of the instrument to strata interface

    1)随钻电阻率成像测井纵向分辨率取决于钮扣电极的直径,并与钮扣电极的直径相当。测井仪器周向设计分布4个方位性钮扣电极,能够识别方位性高阻地层。

    2)不同方位的钮扣电极在水平井中的测井响应特征不同,靠近地层界面钮扣电极的测井曲线在地层界面处有明显的“犄角”现象,而与地层界面垂直的钮扣电极以及仪器侧向电阻率测量电极在地层界面处的测井响应与常规电缆侧向电阻率测井类似。

    3)作为仪器研发的先导,数值模拟可以有效缩短仪器的研发周期,但是其模拟环境大多为理想环境,与真实地层环境具有一定的差距,建议尽快研制出随钻电阻率测井仪器样机和建立地层模型,以验证该仪器方案的可行性。

  • [1] 李大奇,康毅力,曾义金,等.缝洞型储层缝宽动态变化及其对钻井液漏失的影响[J].中国石油大学学报(自然科学版),2011,35(5):76-81. LI Daqi,KANG Yili,ZENG Yijin,et al.Dynamic variation of fracture width and its effects on drilling fluid lost circulation in fractured vuggy reservoirs[J].Journal of China University of Petroleum(Edition of Natural Science),2011,35(5):76-81.
    [2] 王业众,康毅力,游利军,等.裂缝性储层漏失机理及控制技术进展[J].钻井液与完井液,2007,24(4):74-77. WANG Yezhong,KANG Yili,YOU Lijun,et al.Progresses in mechanism study and control:mud losses to fractured reservoirs[J].Drilling Fluid Completion Fluid,2007,24(4):74-77.
    [3] 蒋海军,鄢捷年.架桥粒子粒径与裂缝有效流动宽度匹配关系的试验研究[J].钻井液与完井液,2000,17(4):1-3,7. JIANG Haijun,YAN Jienian.Laboratory study on the compatibility between the diameter of bridging particles and the effective fracture width[J].Drilling Fluid Completion Fluid,2000,17(4):1-3,7.
    [4]

    LUTHI S M,SOUHAITE P.Fracture apertures from electrical borehole scans[J].Geophysics,1990,55(7):821-833.

    [5]

    HORNBY B E,JOHNSON D L,WINKLER K W,et al.Fracture evaluation using reflected Stoneley-wave arrivals[J].Geophysics,1989,54(10):1274-1288.

    [6]

    LAVROV A.Flow of truncated power-law fluid between parallel walls for hydraulic fracturing applications[J].Journal of Non-Newtonian Fluid Mechanics,2015,223:141-146.

    [7]

    OMOSEBI A O,ADENUGA K A.Pressure drop versus flow rate profiles for power-law and Herschel-Bulkley fluids[R].SPE 162999,2012.

    [8]

    MAJIDI R,MISKA S,ZHANG Jianguo.Fingerprint of mud losses into natural and induced fractures[R].SPE 143854,2011.

    [9]

    LAVROV A,TRONVOLL J.Numerical analysis of radial flow in a natural fracture:applications in drilling performance and reservoir characterization[R].SPE 103564,2006.

    [10]

    LIETARD O,SPIVEY J.Revisiting pressure transient testing of hydraulically fractured wells:a single,simple and exact analytical solution covering bilinear,linear,and transition in between flow regimes[R].SPE 139508,2011.

    [11] 金业权,胡满,吴谦,等.Macondo深水井漏油事故防喷器系统失效原因分析[J].石油钻探技术,2014,42(4):53-58. JIN Yequan,HU Man,WU Qian,et al.Analysis of deepwater BOP failure in the Macondo Well accident[J].Petroleum Drilling Techniques,2014,42(4):53-58.
    [12] 臧艳彬,王瑞和,张锐.川东北地区钻井漏失及堵漏措施现状分析[J].石油钻探技术,2011,39(2):60-64. ZANG Yanbin,WANG Ruihe,ZHANG Rui.Current situation analysis of circulation lost and measures in Northeast Sichuan Basin[J].Petroleum Drilling Techniques,2011,39(2):60-64.
    [13] 林英松,蒋金宝,秦涛.井漏处理技术的研究及发展[J].断块油气田,2005,12(2):4-7. LIN Yingsong,JIANG Jinbao,QIN Tao.The development of well loss processing technology[J].Fault-Block Oil Gas Field,2005,12(2):4-7.
    [14] 邹德永,赵建,郭玉龙,等.渗透性砂岩地层漏失压力预测模型[J].石油钻探技术,2014,42(1):33-36. ZOU Deyong,ZHAO Jian,GUO Yulong,et al.A model for predicting leak-off pressure in permeable-sandstone formations[J].Petroleum Drilling Techniques,2014,42(1):33-36.
    [15] 刘加杰,钟颖,张浩.利用钻井液漏失资料预测裂缝宽度[J].中国井矿盐,2014,45(3):20-22. LIU Jiajie,ZHONG Ying,ZHANG Hao.Forecast crack width with drilling fluid leakage data[J].China Well and Rock Salt,2014,45(3):20-22.
    [16] 李大奇,康毅力,刘修善,等.裂缝性地层钻井液漏失动力学模型研究进展[J].石油钻探技术,2013,41(4):42-47. LI Daqi,KANG Yili,LIU Xiushan,et al.Progress in drilling fluid loss dynamics model for fractured formations[J].Petroleum Drilling Techniques,2013,41(4):42-47.
    [17]

    LIETARD O,UNWIN T,GUILLOT D,et al.Fracture width LWD and drilling mud/LCM selection guidelines[R].SPE 36832,1996.

    [18]

    VERGA F,CARUGO C,CHELINI V,et al.Detection and characterization of fractures in naturally fractured reservoirs[R].SPE 63266,2000.

    [19]

    GUILLOT D.A digest of rheological equations[J].Developments in Petroleum Science,1990,28:A1-A8.

    [20] 李秀敏,江卫华.相关系数与相关性度量[J].数学的实践与认识,2006,36(12):188-192. LI Xiuming,JIANG Weihua.Research on linear correlation and dependence measure[J].Mathematics in Practice and Theory,2006,36(12):188-192.
    [21] 严丽坤.相关系数与偏相关系数在相关分析中的应用[J].云南财贸学院学报,2003,19(3):78-80. YAN Likun.Application of correlation coefficient and based correlation coefficient in related analysis[J].Journal of Yunnan University of Finance and Economics,2003,19(3):78-80.
  • 期刊类型引用(7)

    1. 柳杰,李立伟,田小强,许林,沈雄伟,苏方波,谢培. 随钻聚焦电阻率测井仪响应特性分析与应用. 石油管材与仪器. 2023(03): 46-51 . 百度学术
    2. 尤嘉祺,陈刚,陈思嘉,阳质量,许月晨,唐章宏. 随钻侧向电阻率仪器钮扣电极扫描成像影响因素分析. 测井技术. 2022(01): 11-14+28 . 百度学术
    3. 刘冀昆,王成虎,乾增珍,周昊,谷佳诚. 超声波钻孔成像技术在岩体质量定量评价中的应用. 工程勘察. 2022(12): 6-11+46 . 百度学术
    4. 康正明,柯式镇,李新,倪卫宁,李飞. 基于随钻电阻率成像测井仪的洞穴评价模型理论研究. 中南大学学报(自然科学版). 2021(05): 1542-1551 . 百度学术
    5. 林昕,苑仁国,秦磊,刘素周,苏朝博,卢中原,于忠涛,谭伟雄. 地质导向钻井前探技术现状及进展. 特种油气藏. 2021(02): 1-10 . 百度学术
    6. 康正明,柯式镇,李新,倪卫宁,李飞. 随钻电阻率成像测井仪定量评价地层界面探究. 石油钻探技术. 2020(04): 124-130 . 本站查看
    7. 李新,倪卫宁,米金泰,康正明,闫立鹏,宋朝晖. 一种基于非接触耦合原理的新型随钻微电阻率成像仪器. 中国石油大学学报(自然科学版). 2020(06): 46-52 . 百度学术

    其他类型引用(3)

计量
  • 文章访问数:  2895
  • HTML全文浏览量:  112
  • PDF下载量:  2650
  • 被引次数: 10
出版历程
  • 收稿日期:  2015-07-28
  • 修回日期:  2016-03-03
  • 刊出日期:  1899-12-31

目录

/

返回文章
返回