留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

南海深水超浅层含水合物气层测井识别与饱和度估算方法

汤翟 裴健翔 赵军 骆玉虎 郭璃 贾将

汤翟,裴健翔,赵军,等. 南海深水超浅层含水合物气层测井识别与饱和度估算方法[J]. 地质科技通报,2025,44(6):317-329 doi: 10.19509/j.cnki.dzkq.tb20240082
引用本文: 汤翟,裴健翔,赵军,等. 南海深水超浅层含水合物气层测井识别与饱和度估算方法[J]. 地质科技通报,2025,44(6):317-329 doi: 10.19509/j.cnki.dzkq.tb20240082
TANG Di,PEI Jianxiang,ZHAO Jun,et al. Logging identification and saturation estimation method for hydrate-bearing gas layers in the deep water and ultra-shallow strata of the South China Sea[J]. Bulletin of Geological Science and Technology,2025,44(6):317-329 doi: 10.19509/j.cnki.dzkq.tb20240082
Citation: TANG Di,PEI Jianxiang,ZHAO Jun,et al. Logging identification and saturation estimation method for hydrate-bearing gas layers in the deep water and ultra-shallow strata of the South China Sea[J]. Bulletin of Geological Science and Technology,2025,44(6):317-329 doi: 10.19509/j.cnki.dzkq.tb20240082

南海深水超浅层含水合物气层测井识别与饱和度估算方法

doi: 10.19509/j.cnki.dzkq.tb20240082
基金项目: 海南省科技重点研发项目“莺−琼盆地浅层天然气富集机理与有利勘探方向研究”(ZDYF2024GXJS031)
详细信息
    作者简介:

    汤翟:E-mail:tang_di@126.com

    通讯作者:

    E-mail:zhaojun_70@126.com

  • 中图分类号: P631.8;P618.13

Logging identification and saturation estimation method for hydrate-bearing gas layers in the deep water and ultra-shallow strata of the South China Sea

Funds: Supported by the Science and Technology Special Fund of Hainan Province and the Hainan Provincial Key R&D Project "Study on the Enrichment Mechanism and Favorable Exploration Directions of Shallow Gas in the Ying-Qiong Basin" (No. ZDYF2024GXJS031).
More Information
  • 摘要:

    在海域天然气水合物的勘探过程中,含水合物气层因同时赋存天然气水合物与浅层气,电性测井响应特征异常复杂,导致其定性识别与定量评价面临较大挑战。为解决这一问题,本研究充分利用测井多物理场在水合物和浅层气中的差异响应,创新性地提出了基于纵波测井与电阻率测井多信息源的联合反演方法,以提高含水合物气层饱和度的计算精度。基于南海深水超浅层疏松砂岩含水合物气层“同增同减”的测井响应特征,综合筛选低自然伽马、低泥质含量、高孔隙度、稳定厚度的砂质层段作为储层,并结合孔隙度差值法、中子−密度曲线重叠法、剪切模量法等定性识别方法进行层位判定。随后,采用中子−密度联合法计算孔隙度,在三相Biot方程与阿尔奇公式均考虑含气饱和度的基础上,运用循环迭代法同步反演纵波速度与电阻率,优化联合误差,最终求解含水合物气层的饱和度。研究结果表明,在储层段内,通过综合运用含烃或水合物的测井指示方法,并结合流体性质判别表,可有效识别含水合物气层。含水合物气层的典型特征包括:电阻率绝对值法和正演电阻率曲线重叠法指示含烃或水合物,中子−密度曲线重叠法指示含浅层气,剪切模量高于背景值等。采用纵波与电阻率测井的联合反演方法计算含水合物气层饱和度,验证其可行性与可靠性。应用于L区块Z井的联合反演计算结果与岩心饱和度的吻合度达81.25%,L区块Y井的计算结果与单独使用水合物或浅层气计算模型的吻合度接近85%。研究成果可为现场含水合物气层的识别及饱和度计算提供重要参考,为深水区域水合物资源的精细化评价奠定技术基础。

     

  • 图 1  含水合物气层测井参数交会图

    Figure 1.  Crossplot of log parameters for hydrate-bearing gas layer

    图 2  L区块X井三相混合层声电响应特征图(1 in=25.4 mm,下同)

    Figure 2.  Acoustic-electric response characteristics of three-phase mixed layer in Well X, Block L

    图 3  L区块含水合物气层测井识别图(a. X井;b. Y井. 电阻率单位为Ω·m;中子单位为%;密度单位为g/cm3;剪切模量单位为GPa;下同)

    Figure 3.  Logging-based identification diagram of hydrate-bearing gas layer in Block L

    图 4  纵波和电阻率测井联合反演计算饱和度流程图

    Viv. 反演的纵波速度;Vp. 测井纵波速度;Riv. 反演的电阻率;Rt. 测井电阻率

    Figure 4.  Flowchart of joint inversion calculation of saturation using P-wave logging and resistivity logging

    图 5  L区块取样的糊状(a)和块状(b)泥砂样品

    Figure 5.  Viscous (a) and granular (b) mud and sand samples collected in Block L

    图 6  L区块阿尔奇公式系数交会图

    Figure 6.  Crossplot of Archie coefficients in Block L

    图 7  L区综合解释成果图(a. Z井;b. X井;c. Y井)

    Figure 7.  Comprehensive interpretation map of Block L

    表  1  琼东南盆地L区块含水合物气层测井参数统计结果

    Table  1.   Statistical results of logging parameters for gas hydrate reservoir in Block L, Southeast Qiongdong Basin

    流体性质 自然伽马/
    GAPI
    电阻率/
    (Ω·m)
    中子孔隙度/
    %
    纵波时差/
    (μs·ft−1)
    密度/
    (g·cm−3)

    83 1.5 57 175 2.02
    85 1.6 56 178 2.02
    81 1.7 49 176 2.04
    84 1.6 56 175 2.09
    85 1.6 52 178 2.07
    83 1.6 61 177 2.12
    83 1.8 61 172 2.05



    81 2.6 55 136 1.99
    82 2.9 69 168 1.98
    82 2.8 66 168 2.02
    83 3.5 71 124 1.97
    82 3.8 71 114 1.98
    83 4.0 64 132 2.01
    83 2.9 78 137 2.02





    78 1.8 55 164 1.91
    90 1.9 52 187 1.91
    85 2.1 49 192 1.92
    78 4.1 48 202 1.93
    78 4.4 49 206 1.97
    74 2.6 50 196 1.87
    73 2.3 50 179 1.9



    81 2.2 46 179 1.83
    78 2.0 42 180 1.76
    81 1.9 40 172 1.77
    79 1.9 40 178 1.69
    80 1.9 39 182 1.69
    81 2.0 39 181 1.67
    85 1.9 38 181 1.69
    注:1 ft=12 in=0.3048 m,下同
    下载: 导出CSV

    表  2  基于电阻率特征的流体性质判别

    Table  2.   Discrimination of fluid properties based on resistivity characteristics

    电阻率 $ {\mathbf{\varphi }}_{\rm{s}}-{\mathbf{\varphi }}_{\rm{w}} $ $ {\mathbf{\varphi }}_{\rm{D}}-{\mathbf{\varphi }}_{\rm{w}} $ TZ−TJ Gz−Gb 判别结果
    不变=0=0≤0=0水层
    增大<0≈0≤0>0水合物层
    增大<0>0>0>0含水合物气层水合物>浅层气
    ≈0>0水合物≈浅层气
    >0>0水合物<浅层气
    增大>>0>>0>0=0浅层气层
      注:$\varphi_{\mathrm{s}} $. 纵波计算孔隙度;φw. 含水孔隙度;φD. 密度计算孔隙度;TZ. 中子孔隙度;TJ. 中子孔隙度基值;Gz. 剪切模量;Gb. 剪切模量背景值
    下载: 导出CSV

    表  3  井下泥砂样品信息

    Table  3.   Subsurface sand and mud samples information

    序号 井深/m 黏土 石盐 石英 钾长石 斜长石 方解石 白云石
    φB/%
    1 1651.5 42.1 2.0 31.4 4.7 5.6 14.3 0
    2 1653.6 30.2 0.9 23.5 29.2 4.7 11.5 0
    3 1724.4 33.1 0 38.9 4.3 13.3 7.0 3.4
    4 1728.4 41.0 0 33.6 3.2 11.4 7.6 3.1
    下载: 导出CSV

    表  4  琼东南盆地L区Z井含水合物气层饱和度预测结果对比

    Table  4.   Comparison of predicted saturation in the hydrate-bearing gas layer for Well Z, Block L, Qiongdongnan Basin

    序号深度/m水合物饱和度预测值/%岩心实测饱和度/%误差值/%
    11630.5135.235.30.1
    21632.3433.433.60.2
    31646.5319.831.311.5
    41647.3521.222.10.1
    51649.5933.233.20
    61650.8242.232.69.6
    71653.8230.130.00.1
    81654.6222.422.30.1
    91655.1216.822.15.3
    101668.9231.330.90.4
    111670.2128.929.10.2
    121671.5925.925.80.1
    131674.2131.331.50.2
    141675.1129.930.20.3
    151677.5327.928.40.5
    161707.2618.218.50.3
    下载: 导出CSV
  • [1] SLOAN E D. Fundamental principles and applications of natural gas hydrates[J]. Nature, 2003, 426: 353-359. doi: 10.1038/nature02135
    [2] 李宁, 孙文杰, 李心童, 等. 天然气水合物饱和度测井解释模型及方程[J]. 石油勘探与开发, 2022, 49(6): 1073-1079.

    LI N, SUN W J, LI X T, et al. Gas hydrate saturation model and equation for logging interpretation[J]. Petroleum Exploration and Development, 2022, 49(6): 1073-1079. (in Chinese with English abstract
    [3] 马玲, 褚梦凡, 包锐. 浅谈南海水合物分解、甲烷转化及衍生碳的埋藏[J/OL]. 地质科技通报: 1-13[2024-04-03]. https://link.cnki.net/doi/10.19509/j.cnki.dzkq.tb20230591.

    MA L, CHU M F, BAO R. Hydrate decomposition, methane conversion and burial of methane-derived carbon in the South China Sea[J/OL]. Bulletin of Geological Science and Technology: 1-13[2024-04-03]. https://link.cnki.net/doi/10.19509/j.cnki.dzkq.tb20230591.(in Chinese with English abstract
    [4] 单文昊, 王琳, 吴祥恩, 等. 功能化多壁碳纳米管与L-亮氨酸复配体系中甲烷水合物动力学特征[J]. 地质科技通报, 2024, 43(5): 161-169.

    SHAN W H, WANG L, WU X E, et al. Kinetic characteristics of methane hydrate in functionalized multi-walled carbon nanotubes and L-leucine compounding system[J]. Bulletin of Geological Science and Technology, 2024, 43(5): 161-169. (in Chinese with English abstract
    [5] 郁桂刚, 欧文佳, 吴翔, 等. 天然气水合物分解动力学研究进展[J]. 地质科技通报, 2023, 42(3): 175-188.

    YU G G, OU W J, WU X, et al. Research advances on the dissociation dynamics of natural gas hydrates[J]. Bulletin of Geological Science and Technology, 2023, 42(3): 175-188. (in Chinese with English abstract
    [6] 邓勇. 琼东南盆地天然气水合物地震识别与饱和度预测[J]. 地球科学, 2024, 49(5): 1865-1875.

    DENG Y. Seismic identification and saturation prediction of natural gas hydrate in the Qiongdongnan Basin[J]. Earth Science, 2024, 49(5): 1865-1875. (in Chinese with English abstract
    [7] 李丽松, 苗琦. 天然气水合物勘探开发技术发展综述[J]. 天然气与石油, 2014, 32(1): 67-71.

    LI L S, MIAO Q. Review on natural gas hydrate exploration and development technology[J]. Natural Gas and Oil, 2014, 32(1): 67-71. (in Chinese with English abstract
    [8] GUERIN G, GOLDBERG D, MELTSER A. Characterization of in situ elastic properties of gas hydrate-bearing sediments on the Blake Ridge[J]. Journal of Geophysical Research (Solid Earth), 1999, 104(B8): 17781-17795. doi: 10.1029/1999JB900127
    [9] LEE M W, COLLETT T S. Gas hydrate and free gas saturations estimated from velocity logs on hydrate ridge, offshore oregon, U. S. A[M]//TREHU A M, BOHRMANN G, TORRES M E, et al. Proceedings of the ocean drilling program, scientific results. [S. l.]: [s. n.], 2006: 1-25.
    [10] MILKOV A V, DICKENS G R, CLAYPOOL G E, et al. Co-existence of gas hydrate, free gas, and brine within the regional gas hydrate stability zone at Hydrate Ridge (oregon margin): Evidence from prolonged degassing of a pressurized core[J]. Earth and Planetary Science Letters, 2004, 222(3/4): 829-843.
    [11] KANG D J, LU J A, ZHANG Z J, et al. Fine-grained gas hydrate reservoir properties estimated from well logs and lab measurements at the Shenhu gas hydrate production test site, the northern slope of the South China Sea[J]. Marine and Petroleum Geology, 2020, 122: 104676. doi: 10.1016/j.marpetgeo.2020.104676
    [12] LEE M W. Elastic velocities of partially gas-saturated unconsolidated sediments[J]. Marine and Petroleum Geology, 2004, 21(6): 641-650. doi: 10.1016/j.marpetgeo.2003.12.004
    [13] HUI G G, LI S Z, GUO L L, et al. Source and accumulation of gas hydrate in the northern margin of the South China Sea[J]. Marine and Petroleum Geology, 2016, 69: 127-145. doi: 10.1016/j.marpetgeo.2015.10.009
    [14] QIAN J, WANG X J, COLLETT T S, et al. Downhole log evidence for the coexistence of structure Ⅱ gas hydrate and free gas below the bottom simulating reflector in the South China Sea[J]. Marine and Petroleum Geology, 2018, 98: 662-674. doi: 10.1016/j.marpetgeo.2018.09.024
    [15] QIN X W, LU J G, LU H L, et al. Coexistence of natural gas hydrate, free gas and water in the gas hydrate system in the Shenhu area, South China Sea[J]. China Geology, 2020, 3(2): 210-220. doi: 10.31035/cg2020038
    [16] 谢莹峰, 陆敬安, 匡增桂, 等. 南海神狐海域水合物三相混合层测井评价方法研究[J]. 现代地质, 2022, 36(1): 182-192.

    XIE Y F, LU J A, KUANG Z G, et al. Well logging evaluation for three-phase zone with gas hydrate in the Shenhu area, South China Sea[J]. Geoscience, 2022, 36(1): 182-192. (in Chinese with English abstract
    [17] ZHAN L S, LU H L. Introducing realistic distributions of gas bubble size and inclusion aspect ratio to model the coexistence of gas hydrate and free gas[J]. Journal of Applied Geophysics, 2023, 213: 105023. doi: 10.1016/j.jappgeo.2023.105023
    [18] 高兴军, 于兴河, 李胜利, 等. 地球物理测井在天然气水合物勘探中的应用[J]. 地球科学进展, 2003, 18(2): 305-311.

    GAO X J, YU X H, LI S L, et al. Applicaton of geophysical well logging technology in exploration of gas hydrate[J]. Advance in Earth Sciences, 2003, 18(2): 305-311. (in Chinese with English abstract
    [19] 赵军, 武延亮, 周灿灿, 等. 天然气水合物的测井评价方法综述[J]. 测井技术, 2016, 40(4): 392-398.

    ZHAO J, WU Y L, ZHOU C C, et al. Review on logging technology evaluation methods of nature gas hydrate[J]. Well Logging Technology, 2016, 40(4): 392-398. (in Chinese with English abstract
    [20] 王秀娟, 王吉亮, 李伟, 等. 基于声波阻抗估算印度安达曼海天然气水合物饱和度[J]. 海洋地质与第四纪地质, 2013, 33(2): 163-168.

    WANG X J, WANG J L, LI W, et al. Gas hydrate saturation estimated from acoustic impedance in the Andaman Sea, India[J]. Marine Geology & Quaternary Geology, 2013, 33(2): 163-168. (in Chinese with English abstract
    [21] 康冬菊, 刘俊东, 李海燕, 等. 南海神狐海域天然气水合物储层参数定量评价方法[J]. 现代地质, 2024, 38(2): 385-397.

    KANG D J, LIU J D, LI H Y, et al. A quantitative evaluation method for the natural gas hydrate reservoir parameters in the Shenhu area, South China Sea[J]. Geoscience, 2024, 38(2): 385-397.
    [22] 张磊. 改进的三相Biot模型在水合物地球物理测井评价中的研究[D]. 北京: 中国地质大学(北京), 2015.

    ZHANG L. Research on improved three-phase Biot model in gas hydrate geophysical logging evaluation [D]. Beijing: China University of Geosciences (Beijing), 2015. (in Chinese with English abstract
    [23] 何彬, 匡邓晖, 王大魁. 海域天然气水合物测井响应的研究进展[J]. 科技信息, 2013(14): 3-4.

    HE B, KUANG D H, WANG D K. Research progress on logging response of natural gas hydrate in sea area[J]. Science & Technology Information, 2013(14): 3-4. (in Chinese)
    [24] 史謌, 杨东全. 岩石波速和孔隙度、泥质含量之间的关系研究[J]. 北京大学学报(自然科学版), 2001, 37(3): 379-384.

    SHI G, YANG D Q. The regression analysis study on velocity and porisity, and clay content of rocks[J]. Acta Scicentiarum Naturalum Universitis Pekinesis, 2001, 37(3): 379-384. (in Chinese with English abstract
    [25] 谢莹峰, 任金锋, 邓炜, 等. 琼东南盆地渗漏型天然气水合物识别与差异成藏特征[J]. 海洋地质与第四纪地质, 2024, 44(6): 1-11.

    XIE Y F, REN J F, DENG W, et al. Identification and the differential accumulation of leakage-typed gas hydrate in the Qiongdongnan Basin, northern South China Sea[J]. Marine Geology & Quaternary Geology, 2024, 44(6): 1-11. (in Chinese with English abstract
    [26] 贾培文, 吕瑶瑶, 苏丕波, 等. 神狐海域天然气水合物测井-地震响应特征[J]. 地质学报, 2024, 98(9): 2654-2665.

    JIA P W, LYU Y Y, SU P B, et al. Logging-seismic response characteristics of natural gas hydrates in the Shenhu area, South China Sea[J]. Acta Geologica Sinica, 2024, 98(9): 2654-2665. (in Chinese with English abstract
    [27] 赵军, 郑超, 裴健翔, 等. 深水浅层松散沉积物水合物赋存状态识别与饱和度计算[J]. 天然气地球科学, 2024, 35(10): 1713-1723.

    ZHAO J, ZHENG C, PEI J X, et al. Hydrate occurrence identification of shallow loose sediments in deep water and its saturation calculation[J]. Natural Gas Geoscience, 2024, 35(10): 1713-1723. (in Chinese with English abstract
    [28] 周吉林. 天然气水合物和游离气储层的测井与地震响应特征研究[D]. 山东青岛: 中国科学院大学(中国科学院海洋研究所), 2022.

    ZHOU J L. The characteristics of gas hydrate-bearing and free gas-bearing sediments identified from well log and seismic data[D]. Qingdao Shandong: Institute of Oceanology, Chinese Academy of Sciences, 2022. (in Chinese with English abstract
    [29] 赵军, 刘仕鑫, 李元平, 等. 基于测井资料的南海天然气水合物储层识别与参数计算[J]. 地球物理学进展, 2022, 37(6): 2508-2517.

    ZHAO J, LIU S X, LI Y P, et al. Logging identification and evaluation of gas hydrate reservoirs in the South China Sea[J]. Progress in Geophysics, 2022, 37(6): 2508-2517. (in Chinese with English abstract
    [30] 王圣宜, 邹长春, 彭诚, 等. 海域孔隙型储层天然气水合物赋存模式定量化表征: 声波和电阻率测井的约束[J]. 现代地质, 2023, 37(1): 127-137.

    WANG S Y, ZOU C C, PENG C, et al. Quantitative characterization of hydrate occurrence mode in marine pore-filling gas hydrate reservoirs: Constraints from acoustic and resistivity log data[J]. Geoscience, 2023, 37(1): 127-137. (in Chinese with English abstract
    [31] 肖昆, 邹长春, 卢振权, 等. 基于声波测井的冻土区孔隙型水合物储层水合物饱和度估算方法[J]. 地质学报, 2020, 94(5): 1664-1674.

    XIAO K, ZOU C C, LU Z Q, et al. The acoustic log method of estimating gas hydrate saturation in gas hydrate reservoirs[J]. Acta Geologica Sinica, 2020, 94(5): 1664-1674. (in Chinese with English abstract
    [32] 孙建孟, 罗红, 焦滔, 等. 天然气水合物储层参数测井评价综述[J]. 地球物理学进展, 2018, 33(2): 715-723.

    SUN J M, LUO H, JIAO T, et al. Review of the well logging assessment of natural gas hydrate reservoir parameters[J]. Progress in Geophysics, 2018, 33(2): 715-723. (in Chinese with English abstract
    [33] 刘瑞文, 李春山, 管加强. 孔隙流体类型对地层声波速度的影响[J]. 石油钻探技术, 2007, 35(4): 32-34.

    LIU R W, LI C S, GUAN J Q. The impact of pore fluid types on formations acoustic velocity[J]. Petroleum Drilling Techniques, 2007, 35(4): 32-34. (in Chinese with English abstract
    [34] 高红艳, 钟广法, 梁金强, 等. 应用改进的Biot-Gassmann模型估算天然气水合物的饱和度[J]. 海洋地质与第四纪地质, 2012, 32(4): 83-89.

    GAO H Y, ZHONG G F, LIANG J Q, et al. Estimation of gas hydrate saturation with modified biotgassmann theory: A case from northern South China Sea[J]. Marine Geology & Quaternary Geology, 2012, 32(4): 83-89. (in Chinese with English abstract
    [35] 周建, 宋延杰, 姜艳娇, 等. 海洋天然气水合物测井评价研究进展[J]. 西南石油大学学报(自然科学版), 2020, 42(2): 85-93.

    ZHOU J, SONG Y J, JIANG Y J, et al. The research progress of well logging evaluation of marine natural gas hydrate[J]. Journal of Southwest Petroleum University (Science & Technology Edition), 2020, 42(2): 85-93. (in Chinese with English abstract
    [36] 杨克兵, 王竞飞, 马凤芹, 等. 阿尔奇公式的适用条件分析及对策[J]. 天然气与石油, 2018, 36(2): 58-63.

    YANG K B, WANG J F, MA F Q, et al. Analysis and countermeasures about applicable conditions of Archie's formula[J]. Natural Gas and Oil, 2018, 36(2): 58-63. (in Chinese with English abstract
    [37] 赵军, 廖文海, 汤翟, 等. 深水浅层松散细粒水合物储层新型含水饱和度模型建立及应用[J]. 天然气地球科学, 2025, 36(2): 197-208.

    ZHAO J, LIAO W H, TANG Z, et al. Establishment and application of a new water saturation model for shallow and loose fine-grained hydrate reservoirs[J]. Natural Gas Geoscience, 2025, 36(2): 197-208. (in Chinese with English abstract
    [38] 胡高伟, 业渝光, 张剑, 等. 南海沉积物中天然气水合物饱和度与声学特性的关系[J]. 石油学报, 2013, 34(6): 1112-1118.

    HU G W, YE Y G, ZHANG J, et al. Relationship between gas hydrate saturation and acoustic properties of sediments in the South China Sea[J]. Acta Petrolei Sinica, 2013, 34(6): 1112-1118. (in Chinese with English abstract
    [39] 梁劲, 王明君, 王宏斌, 等. 南海神狐海域天然气水合物声波测井速度与饱和度关系分析[J]. 现代地质, 2009, 23(2): 217-223.

    LIANG J, WANG M J, WANG H B, et al. Relationship between the sonic logging velocity and saturation of gas hydrate in Shenhu area, northern slope of South China Sea[J]. Geoscience, 2009, 23(2): 217-223. (in Chinese with English abstract
    [40] 孔秀, 李琦, 姜晓虹, 等. 印度克里希纳−戈达瓦里盆地天然气水合物饱和度特征[J]. 应用海洋学学报, 2019, 38(3): 299-306.

    KONG X, LI Q, JIANG X H, et al. Characteristics of gas hydrate saturations in Krishna-Godavari Basin, India[J]. Journal of Applied Oceanography, 2019, 38(3): 299-306. (in Chinese with English abstract
  • 加载中
图(7) / 表(4)
计量
  • 文章访问数:  395
  • PDF下载量:  25
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-03-07
  • 录用日期:  2024-05-11
  • 修回日期:  2024-04-22
  • 网络出版日期:  2024-05-11

目录

    /

    返回文章
    返回