Search

Article

x

留言板

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

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

Density functional theory calculation of structure and electronic properties in N-methane hydrate

Luo Qiang Yang Heng Guo Ping Zhao Jian-Fei

Jiang Guang-Yu, Sun Chao, Xie Lei, Liu Xiong-Hou. Influence of surface duct on the vertical spatial characteristics of wind-generated noise in deep ocean. Acta Phys. Sin., 2019, 68(2): 024302. doi: 10.7498/aps.68.20181794
Citation: Jiang Guang-Yu, Sun Chao, Xie Lei, Liu Xiong-Hou. Influence of surface duct on the vertical spatial characteristics of wind-generated noise in deep ocean. Acta Phys. Sin., 2019, 68(2): 024302. doi: 10.7498/aps.68.20181794

Density functional theory calculation of structure and electronic properties in N-methane hydrate

Luo Qiang, Yang Heng, Guo Ping, Zhao Jian-Fei
Article Text (iFLYTEK Translation)
PDF
HTML
Get Citation
  • As a clean and efficient unconventional energy source, natural gas hydrate has been highly valued and vigorously developed by many countries in recent years. In order to solve the problem that the existing hydrate structure symmetry is not high, which leads the theoretical research to be restricted, it is imperative to explore a new type of methane hydrate structure with high symmetry. Using the first-principles method which is based on the density functional theory (DFT), the structure and electronic properties of N-methane hydrate are calculated in the generalized gradient approximation (GGA) for Grimme dispersion correction. The obtained results are shown below. 1) The water cage structure of N-methane hydrate is a truncated octahedron (4668), which is composed of 8 regular hexagons and 6 squares, and the average length of the hexagons and the average length of the squares are both 2.723 Å. The average bond length of water molecules is optimized to be 1.056 Å, and the average bond angle of water molecules is 107.738°. The average bond length of methane molecules is 1.0973 Å. The average distance from methane molecules to water molecules is 4.2831 Å that is longer than the distance in the I- methane hydrate. So N-methane hydrate can accommodate larger volumes of gas molecules. The symmetric group is IMˉ3M for N-methane hydrate, which has a simple and strict periodic stable structure. 2) The lattice parameter of N-methane hydrate is 7.70 Å, and the density is 0.903 g/cm3, which is greater than I-, II- and H-type hydrate density. 3) The x-ray diffraction(XRD) pattern of N-methane hydrateis calculated and is close to that of of I-methane hydrate, while the water cage of N-methane hydrate is larger. 4) The interaction between methane molecules and the water cage is van der Waals force, and the formation energy of N- methane hydrate is –0.247 eV, which indicates that the N-methane hydrate is easy to form. Both the density of states and partial density of states indicate that the interaction between methane and water cage is weak, and it relies on molecular force. 5) In addition, N-methane hydrate is an insulator material with the energy gap greater than 5 eV.
      Corresponding author: Luo Qiang, luoqiang@swpu.edu.cn ; Guo Ping, guopingswpi@vip.sina.com
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos.51875091,61701455) and the Science and Technology Project of Sichuan Province, China (Grant No.2018JY0174).

    水下风成噪声是海洋环境噪声的主要组成成分之一[1]. 对声呐设备而言, 风成噪声是限制其性能的重要干扰[2]. 另一方面, 风成噪声中包含许多环境和声场信息, 如海底底质[3,4]、海面粗糙度和海面风速[5,6]等. 深入了解海洋环境对风成噪声特性(包括噪声级、噪声方向性、噪声相关性等)的影响规律, 有利于提高声呐性能, 扩大其探测范围, 并且有利于提高利用风成噪声反演海洋环境参数的准确性.

    风成噪声建模是研究海洋环境对其影响规律的基础. 20世纪60年代初, Cron和Sherman[7]建立了最早的噪声模型, 即C/S模型, 该模型假设海深无限深, 噪声源均匀分布在海洋表面并具有cosmα(通常m=1m=2, α是以垂直向下为0°记的俯仰角)的指向性, 仅适用于海深较深或海底作用不明显的情况; Kuperman和Ingenito[8]假设具有一定强度和相关性的单极子源随机均匀地分布在海面以下某一深度的无穷大平面上, 结合波动理论, 考虑海底的作用, 给出了一般分层海洋环境下的噪声模型, 简称K/I 模型; Carey等[9]将抛物方程声传播模型应用于噪声建模当中, 建立了距离有关环境下的噪声模型; Perkins 和Kuperman等[10]利用二维声传播模型模拟三维环境下的噪声传播, 将噪声模型由二维拓展到了三维.

    在以上噪声模型的基础上, 业内学者开展了大量关于海洋环境对噪声特性影响的研究. 在浅海环境下, Hamson[11]研究了声速剖面和海底底质对噪声方向性和噪声级的影响; Yang和Yoo[12]给出了噪声垂直方向性的简正波表达式, 分析了海底底质和声速剖面对噪声垂直方向性的影响; Rouseff和Tang[13]以及江鹏飞等[14]结合实验数据分析了内波引起的噪声垂直方向性变化, 周建波等[15]研究了海面随机起伏对噪声级、噪声垂直方向性和垂直相关性的影响. 在深海环境下, Carey等[9]研究了斜坡地形对噪声垂直方向性的影响; Buckingham[16]研究了声吸收损失对噪声方向性和相关性的影响; 刘珊琪和李风华[17]以及王璟琰和李风华等[18]分析了不同海底参数下的噪声垂直相关性.

    本文主要关注表面声道对深海风成噪声垂直空间特性(包括噪声垂直方向性和噪声垂直相关性)的影响. 一般来说, 深海的声速分布呈现典型的“三层结构”, 即表面混合层、主跃变层和深海等温层[19]. 其中, 表面混合层的厚度与纬度和季节之间存在明显的依赖关系[20]. 在表面混合层内, 海水温度几乎不随深度变化, 声速分布呈现微弱的正梯度. 已有研究表明, 当存在表面混合层且声源位于混合层内部时, 一部分声能量将被限制在混合层内传播, 该声道称为表面声道[21], 同时在表面声道以下更大的深度范围内, 声传播特性也将会发生一定的变化. 由于噪声源靠近海面, 表面声道的存在将导致噪声传播特性发生改变, 进而影响噪声场特性. Urick[22]指出, 存在表面声道时其内部的噪声级上升, 水平方向上的噪声能量增强. 但在表面声道以下更大的深度范围内, 表面声道引起的噪声特性变化仍待研究.

    针对上述问题, 本文选用K/I模型描述噪声源的几何分布, 以Pekeris分支割线[23]下的简正波理论描述噪声的传播过程, 研究了存在表面声道时, 表面声道以下深度上噪声垂直方向性和噪声垂直相关性的变化, 并分析了表面声道参数变化的影响. 最后, 利用简正波声传播理论解释了表面声道引起噪声垂直空间特性变化的原因.

    为准确描述近场和远场噪声的传播过程, 本文选用Pekeris分支割线下的简正波声传播模型. 在分层介质中, 点源声场的波函数ψ(r,z)满足亥姆霍兹方程

    [2+k2(z)]ψ(r,z)=fs(z,ω)δ(r)2πr,
    (1)

    其中表示拉普拉斯算子, k(z)为波数, c(z)为海水声速, ω为信号角频率, r,z分别为距离和深度, fs(z,ω)为源函数, δ(r)为Dirac 函数.

    则波函数的积分表达式为[24]

    ψ(r,z)=12ψ(kr,z)H(1)0(krr)krdkr,
    (2)

    其中 ψ(r,z)为Hankel变换后的波函数, H(1)0为零阶第一类Bessel函数.

    场的这个积分表示可以用回路定理来解算. 根据复变函数中的Cauchy定理, 回路积分等于回路中所围的极点留数之和及绕分支割线左右岸的积分. 在分层介质传播理论中, 一般有两种分支割线取法, 即EJP(Ewing-Jardetsky-Press)割线[25]和Pekeris割线. 本文选用Pekeris割线, 其示意图见图1.

    图 1 Pekeris分支割线示意图\r\nFig. 1. Pekeris branch cut map.
    图 1  Pekeris分支割线示意图
    Fig. 1.  Pekeris branch cut map.

    图1中, 虚线表示Pekeris分支割线, 实心圆极点和空心圆极点分别对应于“正规”简正波和“非正规”简正波. Bartberger[26]指出, 当信号频率高于100 Hz时, Pekeris割线积分的贡献可以忽略不计, 声场格林函数的简正波表达式(即仅计入“正规”和“非正规”简正波, 忽略Pekeris割线积分)在近场和远场均可给出准确的计算结果.

    G(r,z,z0)=Mm=12πkrmrΨm(z0)Ψm(z)ei(krm+iαm)r=Mm=1am(r)Ψm(z),am(r)2πkrmrΨm(z0)ei(krm+iαm)r,
    (3)

    其中 z0为声源深度, ρ为海水密度, i=1为虚数单位, Ψm(z)为模态函数, krm为第m阶简正波的水平波数, αm为第m阶简正波的吸收系数, am为模态系数, M为简正波总阶数, 其中包含“正规”和“非正规”简正波. 理论上, 存在无穷多阶“非正规”简正波, 因此M应为. 但在数值计算中, 计入无穷多阶“非正规”简正波是不可能的. 在本文的数值仿真中, 计入尽可能多的“非正规”简正波(设置简正波搜寻的相速度上限为100000 m/s), 即M取足够大的值, 以保证计算结果的准确性.

    本文采用K/I模型来描述噪声源的几何分布, 即假设无限多单极子噪声源分布于海面以下某一无限大平面上, 噪声源深度小于λ/4(λ为波长), 且各噪声源之间互不相关. 图2为噪声源几何分布示意图.

    图 2 噪声源几何分布示意图\r\nFig. 2. Geometric distribution of the noise sources.
    图 2  噪声源几何分布示意图
    Fig. 2.  Geometric distribution of the noise sources.

    考虑半径为r, 宽度为dr的圆环, 假设该圆环内所有噪声源到深度为z的接收水听器的声场传递函数均为G(r,z,z0), 则圆环内所有噪声源在深度z处产生的声压为

    p(r,z)=G(r,z,z0)N(r)n=1qneiφn=G(r,z,z0)Qeiφ(r),
    (4)

    其中, N(r)为圆环内噪声源个数, qn为单个噪声源激发的随机声压振幅, Q为圆环内所有噪声源叠加后的声压振幅, φn为单个噪声源的随机相位, φ(r)为圆环内所有噪声源叠加产生的随机相位. 假定单位面积内的噪声源强度为S(r)(单位为dB/μPa2/Hz/m2), 则(4)式可改写为

    p(r,z)=G(r,z,z0)2πrdrS(r)eiφ(r).
    (5)

    噪声的垂直方向性描述了垂直面内各方向上的噪声能量分布, 可由垂直阵通过常规波束形成测得. 在某一俯角θ上, 垂直阵的常规波数形成响应即代表该方向上的噪声能量. 假设垂直阵中各阵元深度为zj,j=1,2,,J, 则对于某一圆环内的噪声源所产生的噪声场, 垂直阵的常规波束形成响应为

    B(r,θ)=|Jj=1p(r,zj)eikzjsinθ|2,
    (6)

    其中表示集合平均. 将(5)式代入(6)式, 有

    B(r,θ)=2πrS(r)|Jj=1G(r,zj,z0)eikzjsinθ|2dr.
    (7)

    由于各圆环的贡献是互不相关的, 垂直阵对噪声场的常规波束形成响应为

    B(θ)=02πrS(r)|Jj=1G(r,zj,z0)eikzjsinθ|2dr.
    (8)

    将(3)式代入(8)式中, 有

    B(θ)=0|Mm=1am(r)2πrS(r)Jj=1Ψm(z0,zj)eikzjsinθ|2dr.
    (9)

    bm(θ)Jj=1Ψm(z0,zj)eikzjsinθ,
    (10)
    ρm1m202πrS(r)am1(r)am2(r)dr,
    (11)

    则(9)式可改写为

    B(θ)=Mm1=1Mm2=1ρm1m2bm1(θ)bm2(θ),
    (12)

    式中当m1=m2时, bm(θ)bm(θ)为第m阶简正波的垂直方向性; ρm1m2为右端各项的叠加系数. 在浅海环境下, 各阶简正波之间的相干程度较弱, (11) 式中交叉项的叠加系数ρm1m2,m1m2可近似为0[8,12], 因此交叉项Mm1=1Mm2=1ρm1m2bm1(θ)bm2(θ),m1m2可忽略不计. 在深海环境下, 相邻阶简正波的相干程度增强, ρm1m2,m1m2不可近似为0[27], 则交叉项Mm1=1Mm2=1ρm1m2bm1(θ)bm2(θ),m1m2也不能忽略.

    将(3)式 代入(11)式可得

    ρm1m2=4π2krm1krm2Ψm1(z0)Ψm2(z0)×0S(r)ei(krm1krm2)r(αm1αm2)rdr.
    (13)

    由(13)式可知, Ψm(z0)直接影响与第m阶简正波相关的叠加系数ρmmρmm的大小, 进而影响(12) 式中右端与第m阶简正波相关的叠加项的强度, 并最终决定第m阶简正波在其来波方向上对噪声场能量贡献的大小. 因此, 本文中定义Ψm(z0)为第m阶简正波的简正波强度.

    噪声的垂直相关性描述了垂直两接收点上噪声的相关程度. 对于深度为z1z2的两个接收点, 噪声的垂直相关系数定义为归一化后的噪声互谱密度. 噪声互谱密度的表达式为

    S12=p(z1)p(z2).
    (14)

    仅考虑半径为r, 宽度为dr的圆环内的噪声源的贡献, 有

    dS12=p(r,z1)p(r,z2).
    (15)

    将(5)式代入(15)式, 有

    dS12=2πrS(r)drG(r,z1,z0)G(r,z2,z0).
    (16)

    由于各圆环的贡献是互不相关的, (14)式可改写为

    S12=02πrS(r)G(r,z1,z0)G(r,z2,z0)dr.
    (17)

    将(3)式代入(17)式, 有

    S12=02πrS(r)Mm1=1Mm2=1am1(r)am2(r)×Ψm1(z1)Ψm2(z2)dr.
    (18)

    将(11)式代入(18)式, 可得

    S12=Mm1=1Mm2=1ρm1m2Ψm1(z1)Ψm2(z2).
    (19)

    在浅海环境下, 由于ρm1m2,m1m2可近似为0, (19)式退化为右端非交叉项的叠加之和, 此时S12的表达式与Kuperman和Ingenito[8]给出的表达式一致. 但是, 在深海环境下, ρm1m2,m1m2不可近似为0, 则(19)式右端的交叉项Mm1=1Mm2=1ρm1m2Ψm1(z1)Ψm2(z2),m1m2也不能忽略.

    由(13)式可知, 简正波强度Ψm(z0)决定了(19)式右端与第m阶简正波相关的叠加项的强度, 表征第m阶简正波对噪声互谱密度S12的影响大小.

    S12进行归一化处理, 得到噪声的垂直相关系数为

    Γ12 = S12S11S22.
    (20)

    d=|z1z2|, 并做变量代换z=z1, 则Γ12可写作dz的函数, 即Γ12=Γ(z,d). 下文中, 称Γ(z,d)为深度z上的噪声垂直相关函数Γ(d).

    在第2节的理论基础上, 本节将对有无表面声道及表面声道参数变化时的噪声垂直方向性和垂直相关性进行仿真分析.

    首先, 图3给出了仿真中无表面声道时的海洋环境. 其中, 声速剖面为标准Munk曲线[28], 声道轴深度1300 m, 临界深度4800 m. 存在表面声道时, 保持其他环境参数不变, 仅改变表面声道内的声速分布. 仿真中设定表面声道厚度100 m, 声道内海水呈等温分布, 受静水压力的影响声道内的声速梯度为0.0167s1[29]. 定义存在表面声道时, 临界深度为进入表面声道以下传播的声能量的最小翻转深度. 则存在表面声道时, 临界深度上的声速等于表面声道下界深度上的声速. 仿真条件下, 存在表面声道时的临界深度为4210 m. 有无表面声道情况下, 声速剖面对比图如图4所示. 仿真中, 假设噪声频率为500 Hz, 设定噪声源深度为z0=0.5 m, 单位面积上的噪声源强度为S(r)=0dB/μPa2/Hz/m2. 如无特殊说明, 下文中保持各项仿真参数不变.

    图 3 无表面声道时的海洋环境\r\nFig. 3. Underwater environment in absence of surface duct.
    图 3  无表面声道时的海洋环境
    Fig. 3.  Underwater environment in absence of surface duct.
    图 4 有无表面声道时的声速剖面\r\nFig. 4. Sound speed profiles in absence and in presence of surface duct.
    图 4  有无表面声道时的声速剖面
    Fig. 4.  Sound speed profiles in absence and in presence of surface duct.

    仿真中噪声垂直方向性由一阵元个数J=20, 阵元间距为1 m的垂直阵通过常规波束形成计算得到. 使用(12)式计算得到有无表面声道时, 不同深度上的噪声垂直方向性B(θ), 如图5所示. 图中0表示海面方向, 180表示海底方向. 由于本文仅关注表面声道以下的噪声垂直空间特性, 下文中默认讨论的深度范围为表面声道以下.

    图 5 有无表面声道时不同深度上的噪声垂直方向性 (a) 无表面声道; (b) 有表面声道\r\nFig. 5. Vertical directionality of the noise at different depths in absence and in presence of surface duct: (a) In absence of surface duct; (b) in presence surface duct.
    图 5  有无表面声道时不同深度上的噪声垂直方向性 (a) 无表面声道; (b) 有表面声道
    Fig. 5.  Vertical directionality of the noise at different depths in absence and in presence of surface duct: (a) In absence of surface duct; (b) in presence surface duct.

    图5(a)图5(b)可以看出, B(θ)在有无表面声道时具有某些共同特征: 1)在临界深度以上, B(θ)在水平方向上始终存在一个水平凹槽, 且凹槽宽度随着深度远离声道轴而逐渐减小; 2)在临界深度以下, 水平凹槽消失. 此外, B(θ)在有无表面声道时存在明显的差异: 1)在临界深度以上, 有无表面声道时, B(θ)在水平凹槽的边缘靠近海底的一侧始终存在一个峰值, 如图中黑色虚线所示, 存在表面声道时, 这一峰值的高度明显升高, 即该方向上的噪声能量明显增强; 2)在临界深度以下, 存在表面声道时, 水平方向上的噪声能量明显增强. 由B(θ)的变化可知, 存在表面声道时, 各深度上来自某一方向的噪声能量增强, 而其他方向的噪声能量几乎不变. 因此, 表面声道的存在同时引起各深度上的噪声级升高.

    图6(a)图6(b)分别给出了1300 m (声道轴深度, 临界深度以上)和5000 m (临界深度以下)深度上有无表面声道时的B(θ), 图中黑色虚线标明了90方向, 即水平方向. 由图6(a)可以看出, 在1300 m深度上, B(θ)90方向上存在一个凹槽, 且在105方向上存在一个峰值, 存在表面声道时, 这一峰值的高度上升了大约3.2 dB. 由图6(b)可以看出, 在5000 m深度上, B(θ)90方向上的凹槽消失, 存在表面声道时90方向上的噪声能量上升了大约2.4 dB.

    图 6 1300 和5000 m深度上有无表面声道时的噪声垂直方向性 (a) 1300 m; (b) 5000 m\r\nFig. 6. Vertical directionality of the noise at 1300 and 5000 m in absence and in presence of surface duct: (a) 1300 m; (b) 5000 m.
    图 6  1300 和5000 m深度上有无表面声道时的噪声垂直方向性 (a) 1300 m; (b) 5000 m
    Fig. 6.  Vertical directionality of the noise at 1300 and 5000 m in absence and in presence of surface duct: (a) 1300 m; (b) 5000 m.

    利用(19)式和(20) 式计算得到有无表面声道时1300 和5000 m深度上的噪声垂直相关函数Γ(d), 分别如图7(a)图7(b)所示. 由于Γ(d)随着垂直距离d增加振荡衰减, 当d较大时, 可近似认为Γ(d)=0, 即噪声是不相关的. 因此, 本文只关注0d/λ5时的Γ(d).

    图 7 1300 和5000 m深度上有无表面声道时的噪声垂直相关函数 (a) 1300 m; (b) 5000 m\r\nFig. 7. Vertical coherence function of the noise at 1300 and 5000 m in absence and in presence of surface duct: (a) 1300 m; (b) 5000 m.
    图 7  1300 和5000 m深度上有无表面声道时的噪声垂直相关函数 (a) 1300 m; (b) 5000 m
    Fig. 7.  Vertical coherence function of the noise at 1300 and 5000 m in absence and in presence of surface duct: (a) 1300 m; (b) 5000 m.

    可以看出, 在临界深度以上和临界深度以下, 表面声道对噪声垂直相关性的影响规律也不相同. 在1300 m深度上, 存在表面声道时的Γ(d)相较于无表面声道时的Γ(d)随着垂直距离d的增加先后周期地向+1和–1方向偏移, 即噪声分别更趋近于正相干和负相干. 图7(a)中, 有无表面声道时Γ(d)取值相等的点将0d/λ5划分为三个区间, 如图中黑色虚线所示. 可以看出, 存在表面声道时, Γ(d)在区间I和区间III内向正相干方向偏移, 在区间II内向负相干方向偏移. 但是, 在5000 m深度上, 存在表面声道时的Γ(d)相较于无表面声道时的Γ(d)整体向正相干方向偏移. 总体来说, 表面声道对噪声垂直相关性的影响较小, 但表面声道的厚度增加将导致其对噪声垂直相关性的影响增强. 该现象将在下文中进一步讨论.

    由以上分析可知, 表面声道的存在将对深海噪声的垂直空间特性产生一定的影响. 本节进一步分析该影响的大小随表面声道参数变化的规律. 由于表面声道对B(θ)的影响更加明显, 本小节中选择B(θ)的变化大小来表征表面声道对噪声垂直空间特性的影响程度. 由Cox[30]的研究可知, 噪声垂直相关性可表示为噪声垂直方向性的函数, 因此表面声道对噪声垂直方向性的影响越大, 其对噪声垂直相关性的影响也越大.

    表面声道的参数主要有表面声道的厚度D和表面声道内的声速梯度G. 其中D与纬度和季节之间存在一定的相关关系, 其变化范围大约为0—300 m. 仿真中考虑较常见的情况, 取D=0 (无表面声道), D=100 m, D=150 m. 保持其他参数不变, 计算得到1300 和 5000 m深度上的噪声垂直方向性B(θ)分别如图8(a)图8(b)所示.

    图 8 表面声道厚度为0 (无表面声道), 50, 100 和150 m 时, 1300和5000 m 深度上的噪声垂直方向性 (a) 1300 m; (b) 5000 m\r\nFig. 8. Vertical directionality of the noise at 1300 and 5000 m with the thickness of the surface duct varing from 0 to 50, 100, and 150 m: (a) 1300 m; (b) 5000 m.
    图 8  表面声道厚度为0 (无表面声道), 50, 100 和150 m 时, 1300和5000 m 深度上的噪声垂直方向性 (a) 1300 m; (b) 5000 m
    Fig. 8.  Vertical directionality of the noise at 1300 and 5000 m with the thickness of the surface duct varing from 0 to 50, 100, and 150 m: (a) 1300 m; (b) 5000 m.

    图8(a)图8(b)可以看出, 在1300 m深度上, B(θ)105左右的峰值高度随着D增大明显升高; 在5000 m 深度上, 水平方向上的噪声能量随着D增大明显增强. 因此, D增大不改变表面声道对B(θ)的影响规律, 但其对B(θ)的影响程度明显增强. 相同地, D增大也不改变表面声道对Γ(d)的影响规律, 但其对Γ(d)的影响程度明显增强.

    当海洋表面的风浪搅拌作用较强时, 表面声道内的海水温度几乎不随深度变化, 此时表面声道内的声速梯度主要由静水压力引起, 大小约为0.0167s1. 当风浪搅拌作用较弱时, 表面声道内的声速梯度也随之减小. 因此, 实际海洋环境中, 表面声道内的声速梯度常小于0.0167s1[31]. 仿真中选取较常见的声速梯度0.002,0.0090.0167s1作为仿真参数, 保持其他参数不变, 计算得到1300和5000 m深度上的B(θ)分别如图9(a)图9(b)所示.

    图 9 无表面声道及表面声道内的声速梯度为0.002, 0.009和0.0167 s−1时, 1300和5000 m 深度上的噪声垂直方向性 (a) 1300 m; (b) 5000 m\r\nFig. 9. Vertical directionality of the noise at 1300 and 5000 m in absence of surface duct and with the sound speed gradient in the 100 m-thick surface duct varying from 0.002 to 0.009 and 0.0167 s−1: (a) 1300 m; (b) 5000 m.
    图 9  无表面声道及表面声道内的声速梯度为0.002, 0.009和0.0167 s−1时, 1300和5000 m 深度上的噪声垂直方向性 (a) 1300 m; (b) 5000 m
    Fig. 9.  Vertical directionality of the noise at 1300 and 5000 m in absence of surface duct and with the sound speed gradient in the 100 m-thick surface duct varying from 0.002 to 0.009 and 0.0167 s−1: (a) 1300 m; (b) 5000 m.

    图9(a)图9(b)可以看出, G取不同值时B(θ)的变化十分有限, 几乎可以忽略不计. 相同地, G的变化对Γ(d)的影响也很小, 几乎可以忽略不计.

    根据以上结果可以认为, 表面声道对噪声垂直空间特性的影响程度仅取决于D, 而几乎与G无关.

    无表面声道时, 噪声源在深海声道中激发的简正波可分为三类, 即无效简正波(invalid mode, IM), 折射简正波(refracted mode, Refr-M)和反射简正波(reflected mode, Refl-M). 其中, IM指对声场没有能量贡献的简正波; Refr-M指到达海底之前已经发生翻转, 不与海底作用的简正波; Refl-M指与海底发生作用的简正波. 存在表面声道时, 一部分Refr-M在到达表面声道下边界之前已经翻转, 进入表面声道内传播, 称为陷获简正波(trapped mode, TM). 本文对TM进行单独讨论, 因此,存在表面声道时, 噪声源激发的简正波分为IM, TM, Refr-M和Refl-M四类.

    定义无表面声道时, 海面声速为cN-SD-s, 海底声速为cN-SD-b; 有表面声道时, 海面声速为cSD-s, 海底声速为cSD-b, 表面声道下边界处的声速为cSD-i. 由图4可知, 各声速始终满足cSD-s<cSD-i<cN-SD-s<cN-SD-b=cSD-b. 则根据模态函数的WKB 近似, 有无表面声道时, 由噪声源激发的各类简正波的水平波数krm满足表1所列的条件. 根据各声速的相对大小, 表1同时给出了存在表面声道时, 各类简正波的水平波数krm的区间大小变化情况.

    表 1  有无表面声道时各类简正波的水平波数krm所处区间
    Table 1.  The intervals of krm for different kinds of normal modes in absence and in presence of surface duct.
    简正波类型Refl-MRefr-MTMIM
    无表面声道krm<ωcN-SD-bωcN-SD-b<krm<ωcN-SD-sωcN-SD-s<krm
    有表面声道krm<ωcSD-bωcSD-b<krm<ωcSD-iωcSD-i<krm<ωcSD-sωcSD-s<krm
    krm区间大小变化不变增大减小
    下载: 导出CSV 
    | 显示表格

    根据表1内各类简正波的水平波数krm所处的区间, 图10(a)图10(b)分别给出了无表面声道和有表面声道时, 各类简正波的本征值在复平面上的分布示意图. 由于存在表面声道时不同类简正波的水平波数krm的区间增大、减小或保持不变, 其阶数也将相应地增多、减少或保持不变. 如图10(b)所示, 存在表面声道时, IM阶数减少, Refr-M阶数增加, Refl-M 阶数基本不变, 同时一部分简正波进入表面声道传播, 即出现TM.

    图 10 有无表面声道时, 各类简正波的本征值在复平面上的分布 (a) 无表面声道; (b) 有表面声道\r\nFig. 10. Distribution of different kinds of normal modes on the complex plane in absence and in presence of surface duct: (a) In absence of surface duct; (b) in presence of surface duct.
    图 10  有无表面声道时, 各类简正波的本征值在复平面上的分布 (a) 无表面声道; (b) 有表面声道
    Fig. 10.  Distribution of different kinds of normal modes on the complex plane in absence and in presence of surface duct: (a) In absence of surface duct; (b) in presence of surface duct.

    表2给出了仿真条件下有无表面声道时各类简正波的阶数. 可以看出, 与无表面声道的情况相比, 存在表面声道时出现24阶TM, 同时IM减少154阶, Refr-M增加131阶, 而Refl-M的阶数基本不变.

    表 2  仿真条件下有无表面声道时各类简正波的阶数
    Table 2.  The indexes of different kinds of normal modes in absence and in presence of surface duct under the simulation environment.
    简正波类型Refl-MRefr-MTMIM
    无表面声道751—3611阶(共2861)590—750阶(共161阶)1—589 阶(共589阶)
    有表面声道752—3611阶(共2860阶)460—751阶(共292阶)436—459 阶(共24阶)1—435 阶(共435阶)
    下载: 导出CSV 
    | 显示表格

    图11给出了仿真条件下, 有无表面声道时, 各阶简正波的简正波强度Ψ(z0). 图中实线和虚线分别表示无表面声道和有表面声道; 黑色、绿色和蓝色分别表示IM, TM, Refr-M 和Refl-M. 由模态函数的WKB近似可知, IM的模态函数在到达噪声源所在深度z0前已经发生指数衰减, 因此IM的Ψ(z0)在有无表面声道时均很小, 如图11中黑色曲线所示. 对于其他类简正波, 其模态函数的第一个极大值点随着简正波阶数增加逐渐靠近海面, 因此其Ψ(z0)随着简正波阶数增加逐渐增强, 并在简正波阶数较高时趋于某一较大的值. 对比图11中红色实线和红色虚线可以看出, 存在表面声道时Refr-M的阶数增加, 因此Refr-M的Ψ(z0)在更低的阶数开始增强, 导致存在表面声道时Refr-M的Ψ(z0)明显强于无表面声道的情况. 由图11中蓝色实线和蓝色虚线可以看出, Refl-M的阶数较高, 其Ψ(z0)已增强到一个较大值, 且Ψ(z0)随阶数增强的速度减缓, 因此有无表面声道对Refr-M的Ψ(z0)影响较小, 几乎可忽略不计.

    图 11 有无表面声道时各阶简正波的简正波强度$\scriptstyle \varPsi(z_0) $ (a) 整体图; (b) 局部放大图\r\nFig. 11. Modal intensity $\scriptstyle \varPsi(z_0) $ of noise-sources-generated normal modes in absence and in presence of surface duct: (a) The overall plot ; (b) zoom in Fig. (a).
    图 11  有无表面声道时各阶简正波的简正波强度Ψ(z0) (a) 整体图; (b) 局部放大图
    Fig. 11.  Modal intensity Ψ(z0) of noise-sources-generated normal modes in absence and in presence of surface duct: (a) The overall plot ; (b) zoom in Fig. (a).

    综上, 仅考虑对声场有能量贡献的简正波, 表面声道对各类简正波的影响主要体现在两个方面, 一是出现TM; 二是使Refr-M阶数增加, 简正波强度增强.

    当表面声道的厚度D增加时cSD-scSD-i发生同等程度的减小, 其他声速不变. 结合以上分析, 此时TM和Refr-M的水平波数krm的区间[ωcSD-i,ωcSD-s][ωcSD-b,ωcSD-i]均增大, 相应地TM 和Refr-M的阶数增加, 简正波强度增强. 当表面声道内的声速梯度G增加时, cSD-s减小, 其他声速不变. 此时仅TM的水平波数krm的区间[ωcSD-i,ωcSD-s]增大, 则TM的阶数增加, 简正波强度增强, 而Refr-M的阶数和简正波强度几乎不变.

    图12(a)图12(b)分别给出了保持其他参数不变, DG取不同值时, 各阶简正波的简正波强度Ψ(z0), 图中黑色、绿色、红色和蓝色分别表示IM, TM, Refr-M和Refl-M, 图12(a)中实线、点划线和虚线分别表示D = 50, 100 和 150 m, 图12(b)中实线、点划线和虚线分别表示G=0.002,0.0090.0167s1. 为方便观察, 图中仅给出关注区域的局部放大图. 可以看出, D增大时, TM和Refr-M的阶数增多, 简正波阶数增强; G增大时, 仅TM的阶数增多, 简正波强度增强, Refr-M的阶数和简正波强度均无明显变化.

    图 12 $\scriptstyle D $和$\scriptstyle G $取不同值时各阶简正波的简正波强度$\scriptstyle \varPsi(z_0) $ (a) $\scriptstyle D $取不同值; (b) $\scriptstyle G $取不同值\r\nFig. 12. Modal intensity $\scriptstyle \varPsi(z_0) $ of noise-sources-excited normal modes with different $\scriptstyle D $ and different $\scriptstyle G $: (a) Different $\scriptstyle D $; (b) different $\scriptstyle G $.
    图 12  DG取不同值时各阶简正波的简正波强度Ψ(z0) (a) D取不同值; (b) G取不同值
    Fig. 12.  Modal intensity Ψ(z0) of noise-sources-excited normal modes with different D and different G: (a) Different D; (b) different G.

    由以上分析可知, 存在表面声道时, 噪声源激发的各类简正波中对声场有能量贡献的有TM, Refr-M和Refl-M. 本节依次分析这三类简正波对噪声垂直方向性的影响.

    4.2.1   TM对噪声垂直方向性的影响

    图13给出了仿真条件下, 存在表面声道时, 仅计入TM计算得到的传播损失. 可以看出, 这类简正波在传播过程中会发生能量“泄漏”, 导致一部分声能量进入表面声道以下传播. 由于“泄漏”声能量的出射角度较小, 所有“泄漏”能量都在临界深度附近翻转. 由其传播特征可知, TM将对表面声道以下, 临界深度以上的噪声垂直方向性产生影响.

    图 13 存在表面声道时, 仅考虑TM计算得到的传播损失\r\nFig. 13. Transmission loss calculated involving only TM in presence of surface duct.
    图 13  存在表面声道时, 仅考虑TM计算得到的传播损失
    Fig. 13.  Transmission loss calculated involving only TM in presence of surface duct.

    图14给出了仿真条件下, 存在表面声道时, 仅考虑TM由(12)式计算得到的1300 m深度上的噪声垂直方向性B(θ)TM. 由图14可以看出, 在1300 m上, B(θ)TM75105左右存在两个峰值, 峰值的高度分别为1.73和0 dB左右. 对比图6(a)图14可以看出, 在75105左右, B(θ)TM的峰值高度比B(θ)低20 dB以上. 因此, TM对噪声场的能量贡献很小, 不是表面声道引起噪声垂直方向性变化的主要因素.

    图 14 存在表面声道时, 仅考虑TM计算得到1300 m深度上的噪声垂直方向性\r\nFig. 14. Vertical directionality of the noise calculated involving only TM at 1300 m depth in presence of surface duct.
    图 14  存在表面声道时, 仅考虑TM计算得到1300 m深度上的噪声垂直方向性
    Fig. 14.  Vertical directionality of the noise calculated involving only TM at 1300 m depth in presence of surface duct.
    4.2.2   Refr-M对噪声垂直方向性的影响

    由模态函数的WKB近似可知, 第m阶简正波对应于两列俯仰角分别为π/2±arccos(krm/k(z))的声波, 相应地, 其对噪声场的能量贡献也集中在这两个方向上. 由表1可知, 在临界深度以上深度z处, 无表面声道和有表面声道时Refr-M的来波俯仰角分别覆盖[π2arccosc(z)cNoSD-b,π2arccosc(z)cNoSD-s][π2+arccosc(z)cNoSD-s,π2+arccosc(z)cNoSD-b][π2arccosc(z)cSD-b,π2arccosc(z)cNoSD-i][π2+arccosc(z)cSD-i,π2+arccosc(z)cSD-b]. 由于cN-SD-b=cSD-b, cN-SD-s<cSD-i, 存在表面声道时, Refr-M的俯仰角度覆盖范围将向90方向(即水平方向)展宽. 在临界深度以下深度z处, krm>k(z)的Refr-M已经发生翻转, 则无表面声道和有表面声道时, Refr-M的来波俯仰角分别覆盖[π2arccosc(z)cNoSD-b,π2+arccosc(z)cNoSD-b][π2arccosc(z)cSD-b,π2+arccosc(z)cSD-b]. 由cN-SD-b=cSD-b可知, 在临界深度以下, 有无表面声道时Refr-M的来波俯仰角范围不变. 另外, 由于表面声道的存在导致Refr-M的简正波强度增强, Refr-M在其来波方向上对噪声场的能量贡献也随之增强.

    图15(a)图15(b)分别给出了仿真环境下, 有无表面声道时, 仅考虑Refr-M计算得到的1300 和5000 m深度上的噪声垂直方向性B(θ)Refr-M.

    在仿真条件下, 由Refr-M的水平波数计算可知, 无表面声道和有表面声道时, Refr-M在1300 m深度上的来波俯仰角分别覆盖[74,75.6][104.4,106][74,77.2][102.8,106]. 因此, 图15(a)中, 有无表面声道时B(θ)Refr-M均在75105左右存在两个峰值, 但峰的宽度在存在表面声道时向水平方向发生一定的展宽. 同时, 由于Refr-M的简正波强度在存在表面声道时增强, 图15(a)B(θ)Refr-M75105左右的峰值高度在存在表面声道时均上升了4.5 dB 左右.

    类似地, 在仿真条件下, 由Refr-M的水平波数计算可知, 有无表面声道时, Refr-M在5000 m 深度上的来波俯仰角覆盖范围均为[84.05,95.95]. 相应地, 图15(b)中, B(θ)Refr-M在有无表面声道时均在90左右存在一个峰值, 且峰的宽度在有无表面声道时无明显变化. 但由于存在表面声道时Refr-M的简正波强度增强, 图15(b)B(θ)Refr-M的峰值高度在存在表面声道时上升了大约2.6 dB.

    图 15 有无表面声道时, 仅考虑Refr-M计算得到的1300和5000 m深度上的噪声垂直方向性 (a) 1300 m; (b) 5000 m\r\nFig. 15. Vertical directionality of the noise calculated involving only Refr-M in absence and in presence of surface duct at 1300 and 5000 m: (a) 1300 m; (b) 5000 m.
    图 15  有无表面声道时, 仅考虑Refr-M计算得到的1300和5000 m深度上的噪声垂直方向性 (a) 1300 m; (b) 5000 m
    Fig. 15.  Vertical directionality of the noise calculated involving only Refr-M in absence and in presence of surface duct at 1300 and 5000 m: (a) 1300 m; (b) 5000 m.
    4.2.3   Refl-M对噪声垂直方向性的影响

    4.1小节中已经指出, 表面声道的存在对Refl-M的水平波数krm的区间范围和简正波强度的影响很小, 几乎可以忽略不计. 因此, Refl-M的来波方向及其在来波方向上贡献的噪声能量在有无表面声道时无明显变化.

    图16(a)图16(b)分别给出了仿真条件下, 有无表面声道时, 仅考虑Refl-M计算得到的1300和5000 m深度上的噪声垂直方向性B(θ)Refl-M. 可以看出, 在1300和5000 m 深度上, 存在表面声道时B(θ)Refl-M均未发生明显的变化.

    图 16 有无表面声道时, 仅考虑Refl-M计算得到的1300和5000 m深度上的噪声垂直方向性 (a) 1300 m; (b) 5000 m\r\nFig. 16. Vertical directionality of the noise calculated involving only Refl-M in absence and in presence of surface duct at depth 1300 and 5000 m: (a) 1300 m; (b) 5000 m.
    图 16  有无表面声道时, 仅考虑Refl-M计算得到的1300和5000 m深度上的噪声垂直方向性 (a) 1300 m; (b) 5000 m
    Fig. 16.  Vertical directionality of the noise calculated involving only Refl-M in absence and in presence of surface duct at depth 1300 and 5000 m: (a) 1300 m; (b) 5000 m.

    对比图6(a)图15(a)图16(a)可以看出, 在1300 m深度上, Refl-M对噪声垂直方向性B(θ)的影响很小, 存在表面声道时Refr-M在105左右对噪声场的能量贡献增强导致B(θ)在该方向上的峰值高度上升, 但Refr-M在75左右的噪声能量贡献增强被来自海面的直达噪声能量所掩盖, B(θ)在该方向上变化并不明显. 相同地, 由图6(b)图15(b)图16(b)可以看出, 在5000 m深度上, Refl-M对B(θ)的影响很小, 存在表面声道时B(θ)在水平方向上的高度上升是由Refr-M在水平方向上的噪声能量贡献增强引起的.

    分析表明, IM和TM对表面声道以下噪声场的能量贡献很小, 几乎可以忽略不计. 因此, 噪声垂直相关性受IM和TM的影响可近似忽略, 主要由Refr-M和Refl-M决定.

    图17(a)图17(b)分别给出了仿真条件下, 有无表面声道时, 1300和5000 m深度上考虑全波场得到的噪声垂直相关函数Γ(d), 仅考虑Rrefr-M得到的噪声垂直相关函数Γ(d)Refr-M和仅考虑Refl-M得到的噪声垂直相关函数Γ(d)Refl-M. 图中Γ(d)Γ(d)Refr-MΓ(d)Refl-M分别以黑色、蓝色和红色曲线给出; 且实线和虚线分别表示无表面声道和有表面声道.

    图 17 有无表面声道时, 1300和5000 m深度上的$\scriptstyle \varGamma(d) $,$\scriptstyle \varGamma(d)_{\rm Refr{\text{-}} M} $和$\scriptstyle \varGamma(d)_{\rm Refl{\text{-}} M} $ (a) 1300 m; (b) 5000 m\r\nFig. 17. $\scriptstyle \varGamma(d) $, $\scriptstyle \varGamma(d)_{\rm Refr{\text{-}} M} $ and $\scriptstyle \varGamma(d)_{\rm Refl{\text{-}} M} $ in absence and in presence of surface duct at 1300 and 5000 m: (a) 1300 m; (b) 5000 m.
    图 17  有无表面声道时, 1300和5000 m深度上的Γ(d),Γ(d)Refr-MΓ(d)Refl-M (a) 1300 m; (b) 5000 m
    Fig. 17.  Γ(d), Γ(d)Refr-M and Γ(d)Refl-M in absence and in presence of surface duct at 1300 and 5000 m: (a) 1300 m; (b) 5000 m.

    图17(a)中, 蓝色实线和蓝色虚线基本重合, 表明表面声道的存在对Γ(d)Refl-M几乎没有影响, 原因是有无表面声道时Refl-M的阶数和简正波强度几乎没有发生变化. 对比图17(a)中红色实线和红色虚线可以看出, 表面声道的存在导致Γ(d)Refr-M的振荡周期变大, 幅度衰减变快. 结合(18)式可知, 该现象是由Refr-M的阶数增加和简正波强度增强引起的. 需要指出的是, 存在表面声道时Refr-M的阶数增加, 简正波强度增强, 导致其对Γ(d)的影响变大, 使Γ(d)Γ(d)Refr-M方向发生更大幅度的偏移, 是存在表面声道时Γ(d)发生变化的主要原因. 对比图17(a)中所有实线可以看出, 黑色实线和蓝色实线基本重合, 即无表面声道时Γ(d)Γ(d)Refl-M基本一致, Refr-M对Γ(d)几乎没有影响. 对比图17(a)中所有虚线可以看出, 当红色虚线位于0以上时, 黑色虚线偏离蓝色虚线向+1(即正相干)方向靠近, 而当红色虚线位于0以下时, 黑色虚线偏离蓝色虚线向−1(即负相干)方向靠近. 因此, 存在表面声道时, Refr-M对Γ(d)的影响增强, 导致存在表面声道时的Γ(d)相较于无表面声道时的Γ(d)随着d的增加先后周期地向正相干和负相干方向偏移.

    图17(b)中的结果与图17(a)是类似的, 不同的是, 在临界深度以下, 有能量贡献的Refr-M已经临近其翻转深度, 其模态函数随深度振荡的周期增大, 导致Γ(d)Refr-M的振荡周期也随之变大. 如图17(a)中红色实线和红色虚线所示, 有无表面声道时Γ(d)Refr-M0d/λ5范围内几乎均大于0. 因此, 作为Refr-M 和Refl-M共同作用的结果, 有无表面声道时Γ(d)均偏离Γ(d)Refl-M向正相干方向靠近, 但在有表面声道时, Refr-M的影响更大, 导致Γ(d)偏离Γ(d)Refl-M的幅度更大. 因此, 存在表面声道时的Γ(d)相较于无表面声道时的Γ(d)整体向正相干方向偏移.

    由4.2和4.3节的分析可知, 表面声道对噪声垂直空间特性的影响主要是由存在表面声道时噪声源激发的Refr-M阶数增加, 简正波强度增强引起的. 因此, 当表面声道的厚度D增加时, Refr-M的阶数增加, 简正波强度增强, 表面声道对噪声垂直空间特性的影响程度也随之增强;当表面声道内的声速梯度G增加时, Refr-M的阶数和简正波强度无明显变化, 则表面声道对噪声垂直空间特性的影响程度也无明显变化.

    本文以Pekeris割线下的简正波模型描述噪声信号的传播过程, 研究了深海环境下, 表面声道对表面声道以下风成噪声垂直空间特性(包括噪声垂直方向性和垂直相关性)的影响. 研究表明, 在临界深度以上和临界深度以下, 表面声道对噪声垂直空间特性的影响规律不同: 在临界深度以上, 表面声道的存在导致噪声垂直方向性在水平凹槽边缘靠近海底方向上的峰值升高, 噪声垂直相关性随垂直距离增加先后周期地向正相干和负相干方向偏移; 在临界深度以下, 表面声道的存在导致水平方向上的噪声能量增强, 噪声垂直相关性整体向正相干方向发生偏移. 当表面声道的厚度增加时, 表面声道对噪声垂直空间特性的影响明显增强; 当表面声道内的声速梯度增加时, 表面声道对噪声垂直空间特性的影响程度无明显变化. 进一步分析表明, 表面声道的存在引起噪声源激发的折射简正波的阶数增加、简正波强度增强, 导致折射简正波在其来波方向上对噪声场的能量贡献增强,对噪声垂直相关性的影响增大, 是表面声道引起噪声垂直空间特性变化的主要原因.

    感谢西北工业大学航海学院周建波博士后的讨论.

    [1]

    Holder G, Kamath V, Godbole S 1984 Annu. Rev. Energy 9 427

    [2]

    Max M, Lowrie A 1996 J. Petrol. Geol. 19 41Google Scholar

    [3]

    Collett T S 2002 AAPG Bull. 86 1971

    [4]

    Kvamme B, Kuznetsova T, Sapate A, Qorbani K 2016 J. Nat. Gas Sci. Eng. 35 1594Google Scholar

    [5]

    Zhao J, Chen X, Song Y, Zhu Z, Yang L, Tian Y, Wang J, Yang M, Zhang Y 2014 Energy Procedia 61 75Google Scholar

    [6]

    Khlebnikov V, Antonov S, Mishin A, Bakulin D, Khamidullina I, 梁萌, Vinokurov V, Gushchin P A 2016 天然气工业 36 40Google Scholar

    Khlebnikov V, Antonov S, Mishin A, Bakulin D, Khamidullina I, Liang M, Vinokurov V, Gushchin P A 2016 Nat. Gas Ind. 36 40Google Scholar

    [7]

    Lim D, Ro H, Seo Y, Seo Y, Lee J, Kim S, Lee J , Lee H 2017 J. Chem. Thermodyn. 106 16Google Scholar

    [8]

    Partoon B, Nashed O, Kassim Z, Sabil K, Sangwai J, Lal B 2016 Process Eng. 148 1220

    [9]

    Zhang L X, Yang L, Wang J Q, Zhao J F, Dong H S, Yang M J, Liu Y, Song Y C 2017 Chem. Eng. J. 308 40Google Scholar

    [10]

    Castellani B, Rossetti G, Tupsakhare S, Rossi F, Nicolini A, Castaldi M 2016 J. Petrol. Sci. Eng. 147 515Google Scholar

    [11]

    Michael K, Roland B, Edward B, Lewis R N 2004 J. Am. Chem. Soc. 126 9407Google Scholar

    [12]

    Xu C G, Cai J, Lin F H, Chen Z Y, Li X S 2015 Energy 79 111Google Scholar

    [13]

    Naeiji P, Varaminian F, Rahmati M 2017 J. Nat. Gas Sci. Eng. 44 122Google Scholar

    [14]

    朱金龙, 赵予生, 靳常青 2019 物理学报 68 018203

    Zhu J L, Zhao Y S, Jin C Q 2019 Acta Phys. Sin. 68 018203

    [15]

    Choudhary N, Kushwaha O S, Bhattacharjee G, Chakrabarty S, Kumara R 2017 Energt Procedia 105 5026Google Scholar

    [16]

    Alavi S, Ohmura R 2016 J. Chem. Phys. 145 154708Google Scholar

    [17]

    Huang Y, Yuan L, Su Y, Zhao J 2015 Mol. Simul. 41 1086Google Scholar

    [18]

    Nguyen A H, Molinero V 2014 J. Chem. Phys. 140 3440

    [19]

    Yagasaki T, Matsumoto M, Andoh Y, Okazaki S, Tanaka H 2014 J. Phys. Chem. B 118 11797

    [20]

    Druart M L, Michel L, Van D H 2014 Fluid Phase Equilib. 381 108Google Scholar

    [21]

    耿春宇, 丁丽颖, 韩清珍, 温浩 2008 物理化学学报 24 595Google Scholar

    Geng C Y, Ding L Y, Han Q Z, Wen H 2008 Acta Phys. Chim. Sin. 24 595Google Scholar

    [22]

    Lasich M, Ramjugernath D 2016 Philos. Mag. 96 15Google Scholar

    [23]

    Zong X, Cheng G, Qiu N, Huang Q, He J, Du S, Li Y 2017 Chem. Lett. 46 1141Google Scholar

    [24]

    Liu Y, Ojamäe L 2014 J. Phys. Chem. A 118 51

    [25]

    Zhang X, Qiu N, Huang Q, Zha X, He J, Li Y, Du S 2017 J. Mol. Struct. 1153 292

    [26]

    Waldron C J, English N J 2017 J. Chem. Phys. 147 024506Google Scholar

    [27]

    Zquierdo-Ruiz F, Otero-de-la-Roza A, Contreras-García J, Menéndez B, Recio J M 2015 High Pressure Res. 35 49Google Scholar

    [28]

    Liu C, Zhang Z, Guo G J 2016 RSC Adv. 6 106443Google Scholar

    [29]

    Tang L, Shi R, Su Y, Zhao J 2015 J. Phys. Chem. A 119 10971Google Scholar

    [30]

    郑朝阳, 赵纪军 2012 物理化学学报 28 1809Google Scholar

    Zheng C Y, Zhao J J 2012 Acta Phys. Chim. Sin. 28 1809Google Scholar

    [31]

    Huang Y, Zhu C, Wang L, Cao X, Su Y, Jiang X, Meng S, Zhao J, Zeng X C 2016 Sci. Adv. 2 e1501010Google Scholar

    [32]

    曹潇潇 2016 博士学位论文(大连: 大连理工大学)

    Cao X X 2016 Ph. D. Dissertation (Dalian: Dalian University of Technology University) (in Chinese)

    [33]

    曹潇潇, 苏艳, 赵纪军, 刘昌岭, 周潘旺 2014 物理化学学报 30 1437Google Scholar

    Cao X X, Su Y, Zhao J J, Liu C L, Zhou D W 2014 Acta Phys. Chim. Sin. 30 1437Google Scholar

    [34]

    Nityananda R, Hohenberg P, Kohn W 2017 Resonance 22 809Google Scholar

    [35]

    Kohn W, Sham L J 1965 Phys. Rev. 140 1133Google Scholar

    [36]

    Huang Y Y, Zhu C Q, Wang L 2017 Chem. Phys. Lett. 2017 186Google Scholar

    [37]

    Blöchl E 1994 Phys. Rev. B 50 17953Google Scholar

    [38]

    Kresse G, Joubert D 1999 Phys. Rev. B 59 1758

    [39]

    罗强, 唐斌, 张智, 冉曾令 2013 物理学报 62 077101Google Scholar

    Luo Q, Tang B, Zhang Z, Ran Z L 2013 Acta Phys. Sin. 62 077101Google Scholar

    [40]

    Wang Z H, Li Y, Meng W J, Guo P, Luo Q, Ran Z L 2017 Appl. Ecol. Environ. Res. 15 861Google Scholar

    [41]

    Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865Google Scholar

    [42]

    Monkhorst H J, Pack J D 1976 Phys. Rev. B 13 5188Google Scholar

    [43]

    Broyden C G 1970 J. Inst. Math. Applic. 6 76Google Scholar

    [44]

    Fletcher R 1970 Comput. J. 13 317Google Scholar

    [45]

    Goldfarb D 1970 Math. Comput. 24 23Google Scholar

    [46]

    Shanno D F 1970 Math. Comput. 24 647Google Scholar

    [47]

    Guo P, Qiu Y L, Li L L, Luo Q, Zhao J F, Pan Y K 2018 Chin. Phys. B 27 276

    [48]

    丁家祥, 史伶俐, 申小冬, 梁德青 2017 化工学报 68 4802

    Ding J X, Shi L L, Sheng X D, Liang D Q 2017 J. Chem. Ind. Eng. 68 4802

    期刊类型引用(4)

    1. 任超,黄益旺. 海洋环境噪声空间相关系数的统计特性. 哈尔滨工程大学学报. 2024(11): 2160-2167 . 百度学术
    2. 任超,黄益旺,夏峙. 宽频带海洋环境噪声矢量场空间相关特性建模. 物理学报. 2022(02): 131-141 . 百度学术
    3. 李明杨,赵航芳,孙超. 风成噪声背景下垂直阵阵列信噪比随声源深度的变化规律. 物理学报. 2022(04): 141-152 . 百度学术
    4. 蒋光禹,孙超,李沁然. 涡旋对深海风成噪声垂直空间特性的影响. 物理学报. 2020(14): 111-125 . 百度学术

    其他类型引用(4)

  • 图 1  模型示意图 (a) N型甲烷水合物骨架模型; (b) 优化后的结构

    Figure 1.  Model diagram: (a) N-methane hydrate skeleton model; (b) optimized structure of N-methane hydrate.

    图 2  N型甲烷水合物结构X射线衍射图

    Figure 2.  X-ray diffraction of N- methane hydrate structure.

    图 3  甲烷分子态密度和分波态密度

    Figure 3.  Density of states and partial density of states in methane molecules.

    图 4  水笼子的态密度及分波态密度

    Figure 4.  Density of state and partial density of state in water cage.

    图 5  能带结构

    Figure 5.  Energy band structures.

    表 1  平均晶格参数

    Table 1.  Average lattice parameters.

    晶格参数GGA原始值I型甲烷水合物[47]
    六边形边长/Å2.72312.68042.7152
    四边形边长/Å2.72332.7656/
    水分子键长/Å1.00560.9940.993
    水分子键角/(°)107.738109.406106.629
    甲烷分子键长/Å1.09731.11781.0917
    甲烷分子键角/(°)109.471109.454109.471
    甲烷分子到水分子距离/Å4.28314.2783.78115
    氢键键长/Å1.72931.70251.7128
    DownLoad: CSV

    表 2  GGA近似下的形成能

    Table 2.  Formation energy of GGA approximation

    能量GGAI型甲烷水合物[47]
    Eform/eV–0.247–0.581
    Etotal/eV–6074.782–21186.729
    Ecage/eV–5634.720–20976.902
    ECH4/eV–219.784–209.246
    DownLoad: CSV
  • [1]

    Holder G, Kamath V, Godbole S 1984 Annu. Rev. Energy 9 427

    [2]

    Max M, Lowrie A 1996 J. Petrol. Geol. 19 41Google Scholar

    [3]

    Collett T S 2002 AAPG Bull. 86 1971

    [4]

    Kvamme B, Kuznetsova T, Sapate A, Qorbani K 2016 J. Nat. Gas Sci. Eng. 35 1594Google Scholar

    [5]

    Zhao J, Chen X, Song Y, Zhu Z, Yang L, Tian Y, Wang J, Yang M, Zhang Y 2014 Energy Procedia 61 75Google Scholar

    [6]

    Khlebnikov V, Antonov S, Mishin A, Bakulin D, Khamidullina I, 梁萌, Vinokurov V, Gushchin P A 2016 天然气工业 36 40Google Scholar

    Khlebnikov V, Antonov S, Mishin A, Bakulin D, Khamidullina I, Liang M, Vinokurov V, Gushchin P A 2016 Nat. Gas Ind. 36 40Google Scholar

    [7]

    Lim D, Ro H, Seo Y, Seo Y, Lee J, Kim S, Lee J , Lee H 2017 J. Chem. Thermodyn. 106 16Google Scholar

    [8]

    Partoon B, Nashed O, Kassim Z, Sabil K, Sangwai J, Lal B 2016 Process Eng. 148 1220

    [9]

    Zhang L X, Yang L, Wang J Q, Zhao J F, Dong H S, Yang M J, Liu Y, Song Y C 2017 Chem. Eng. J. 308 40Google Scholar

    [10]

    Castellani B, Rossetti G, Tupsakhare S, Rossi F, Nicolini A, Castaldi M 2016 J. Petrol. Sci. Eng. 147 515Google Scholar

    [11]

    Michael K, Roland B, Edward B, Lewis R N 2004 J. Am. Chem. Soc. 126 9407Google Scholar

    [12]

    Xu C G, Cai J, Lin F H, Chen Z Y, Li X S 2015 Energy 79 111Google Scholar

    [13]

    Naeiji P, Varaminian F, Rahmati M 2017 J. Nat. Gas Sci. Eng. 44 122Google Scholar

    [14]

    朱金龙, 赵予生, 靳常青 2019 物理学报 68 018203

    Zhu J L, Zhao Y S, Jin C Q 2019 Acta Phys. Sin. 68 018203

    [15]

    Choudhary N, Kushwaha O S, Bhattacharjee G, Chakrabarty S, Kumara R 2017 Energt Procedia 105 5026Google Scholar

    [16]

    Alavi S, Ohmura R 2016 J. Chem. Phys. 145 154708Google Scholar

    [17]

    Huang Y, Yuan L, Su Y, Zhao J 2015 Mol. Simul. 41 1086Google Scholar

    [18]

    Nguyen A H, Molinero V 2014 J. Chem. Phys. 140 3440

    [19]

    Yagasaki T, Matsumoto M, Andoh Y, Okazaki S, Tanaka H 2014 J. Phys. Chem. B 118 11797

    [20]

    Druart M L, Michel L, Van D H 2014 Fluid Phase Equilib. 381 108Google Scholar

    [21]

    耿春宇, 丁丽颖, 韩清珍, 温浩 2008 物理化学学报 24 595Google Scholar

    Geng C Y, Ding L Y, Han Q Z, Wen H 2008 Acta Phys. Chim. Sin. 24 595Google Scholar

    [22]

    Lasich M, Ramjugernath D 2016 Philos. Mag. 96 15Google Scholar

    [23]

    Zong X, Cheng G, Qiu N, Huang Q, He J, Du S, Li Y 2017 Chem. Lett. 46 1141Google Scholar

    [24]

    Liu Y, Ojamäe L 2014 J. Phys. Chem. A 118 51

    [25]

    Zhang X, Qiu N, Huang Q, Zha X, He J, Li Y, Du S 2017 J. Mol. Struct. 1153 292

    [26]

    Waldron C J, English N J 2017 J. Chem. Phys. 147 024506Google Scholar

    [27]

    Zquierdo-Ruiz F, Otero-de-la-Roza A, Contreras-García J, Menéndez B, Recio J M 2015 High Pressure Res. 35 49Google Scholar

    [28]

    Liu C, Zhang Z, Guo G J 2016 RSC Adv. 6 106443Google Scholar

    [29]

    Tang L, Shi R, Su Y, Zhao J 2015 J. Phys. Chem. A 119 10971Google Scholar

    [30]

    郑朝阳, 赵纪军 2012 物理化学学报 28 1809Google Scholar

    Zheng C Y, Zhao J J 2012 Acta Phys. Chim. Sin. 28 1809Google Scholar

    [31]

    Huang Y, Zhu C, Wang L, Cao X, Su Y, Jiang X, Meng S, Zhao J, Zeng X C 2016 Sci. Adv. 2 e1501010Google Scholar

    [32]

    曹潇潇 2016 博士学位论文(大连: 大连理工大学)

    Cao X X 2016 Ph. D. Dissertation (Dalian: Dalian University of Technology University) (in Chinese)

    [33]

    曹潇潇, 苏艳, 赵纪军, 刘昌岭, 周潘旺 2014 物理化学学报 30 1437Google Scholar

    Cao X X, Su Y, Zhao J J, Liu C L, Zhou D W 2014 Acta Phys. Chim. Sin. 30 1437Google Scholar

    [34]

    Nityananda R, Hohenberg P, Kohn W 2017 Resonance 22 809Google Scholar

    [35]

    Kohn W, Sham L J 1965 Phys. Rev. 140 1133Google Scholar

    [36]

    Huang Y Y, Zhu C Q, Wang L 2017 Chem. Phys. Lett. 2017 186Google Scholar

    [37]

    Blöchl E 1994 Phys. Rev. B 50 17953Google Scholar

    [38]

    Kresse G, Joubert D 1999 Phys. Rev. B 59 1758

    [39]

    罗强, 唐斌, 张智, 冉曾令 2013 物理学报 62 077101Google Scholar

    Luo Q, Tang B, Zhang Z, Ran Z L 2013 Acta Phys. Sin. 62 077101Google Scholar

    [40]

    Wang Z H, Li Y, Meng W J, Guo P, Luo Q, Ran Z L 2017 Appl. Ecol. Environ. Res. 15 861Google Scholar

    [41]

    Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865Google Scholar

    [42]

    Monkhorst H J, Pack J D 1976 Phys. Rev. B 13 5188Google Scholar

    [43]

    Broyden C G 1970 J. Inst. Math. Applic. 6 76Google Scholar

    [44]

    Fletcher R 1970 Comput. J. 13 317Google Scholar

    [45]

    Goldfarb D 1970 Math. Comput. 24 23Google Scholar

    [46]

    Shanno D F 1970 Math. Comput. 24 647Google Scholar

    [47]

    Guo P, Qiu Y L, Li L L, Luo Q, Zhao J F, Pan Y K 2018 Chin. Phys. B 27 276

    [48]

    丁家祥, 史伶俐, 申小冬, 梁德青 2017 化工学报 68 4802

    Ding J X, Shi L L, Sheng X D, Liang D Q 2017 J. Chem. Ind. Eng. 68 4802

  • [1] Song Qing-Gong, Zhao Jun-Pu, Gu Wei-Feng, Zhen Dan-Dan, Guo Yan-Rui, Li Ze-Peng. Ductile and electronic properties of La-doped gamma-TiAl systems based on density functional theory. Acta Physica Sinica, 2017, 66(6): 066103. doi: 10.7498/aps.66.066103
    [2] Wen Jun-Qing, Zhang Jian-Min, Yao Pan, Zhou Hong, Wang Jun-Fei. A density functional theory study of small bimetallic PdnAl (n =18) clusters. Acta Physica Sinica, 2014, 63(11): 113101. doi: 10.7498/aps.63.113101
    [3] Wen Jun-Qing, Xia Tao, Wang Jun-Fei. A density functional theory study of small bimetallic PtnAl (n=18) clusters. Acta Physica Sinica, 2014, 63(2): 023103. doi: 10.7498/aps.63.023103
    [4] Ruan Wen, Xie An-Dong, Yu Xiao-Guang, Wu Dong-Lan. Geometric structure and electronic characteristics of NaBn (n=19) clusters. Acta Physica Sinica, 2012, 61(4): 043102. doi: 10.7498/aps.61.043102
    [5] Ge Gui-Xian, Yan Hong-Xia, Jing Qun, Zhang Jian-Jun. Density functional theory study on the structure and electronicproperties of Aun Sc3 (n =1—7) clusters. Acta Physica Sinica, 2011, 60(3): 033101. doi: 10.7498/aps.60.033101
    [6] Zhang Zhi-Long, Chen Yu-Hong, Ren Bao-Xing, Zhang Cai-Rong, Du Rui, Wang Wei-Chao. Density functional theory study on the structure and properties of (HMgN3)n(n=15) clusters. Acta Physica Sinica, 2011, 60(12): 123601. doi: 10.7498/aps.60.123601
    [7] Ge Gui-Xian, Jing Qun, Cao Hai-Bin, Yang Zeng-Qiang, Tang Guang-Hui, Yan Hong-Xia. The study on structures and properties of Run and Run Au(n=112) clusters by density functional theory. Acta Physica Sinica, 2011, 60(10): 103102. doi: 10.7498/aps.60.103102
    [8] Zhang Xiu-Rong, Li Yang, Yang Xing. Theoretical study on structural and electronic properties of WnNim(n+m=8) clusters. Acta Physica Sinica, 2011, 60(10): 103601. doi: 10.7498/aps.60.103601
    [9] Jin Rong, Chen Xiao-Hong. Structure and properties of ZrnPd clusters by density-functional theory. Acta Physica Sinica, 2010, 59(10): 6955-6962. doi: 10.7498/aps.59.6955
    [10] Sun Jian-Min, Zhao Gao-Feng, Wang Xian-Wei, Yang Wen, Liu Yan, Wang Yuan-Xu. Study of structural and electronic properties of Cu-adsorbed (SiO2)n(n=1—8) clusters with the DFT. Acta Physica Sinica, 2010, 59(11): 7830-7837. doi: 10.7498/aps.59.7830
    [11] Gao Hong, Zhu Wei-Hua, Tang Chun-Mei, Geng Fang-Fang, Yao Chang-Da, Xu Yun-Ling, Deng Kai-Ming. Density functional calculation on the geometric structure and electronic properties of the endohedral fullerene N2@C60. Acta Physica Sinica, 2010, 59(3): 1707-1711. doi: 10.7498/aps.59.1707
    [12] Liu Li-Ren, Lei Xue-Ling, Chen Hang, Zhu Heng-Jiang. Geometry and electronic properties of Bn(n=2—15) clusters. Acta Physica Sinica, 2009, 58(8): 5355-5361. doi: 10.7498/aps.58.5355
    [13] Li Xi-Bo, Wang Hong-Yan, Luo Jiang-Shan, Wu Wei-Dong, Tang Yong-Jian. Density functional theory study of the geometry, stability and electronic properties of ScnO(n=1—9) clusters. Acta Physica Sinica, 2009, 58(9): 6134-6140. doi: 10.7498/aps.58.6134
    [14] Chen Yu-Hong, Kang Long, Zhang Cai-Rong, Luo Yong-Chun, Pu Zhong-Sheng. Density functional theory study of the structures and properties of (Li3N)n(n=1—5) clusters. Acta Physica Sinica, 2008, 57(7): 4174-4181. doi: 10.7498/aps.57.4174
    [15] Chen Yu-Hong, Kang Long, Zhang Cai-Rong, Luo Yong-Chun, Yuan Li-Hua, Li Yan-Long. Density functional theory study on the structures and properties of (Ca3N2)n(n=1—4) clusters. Acta Physica Sinica, 2008, 57(10): 6265-6270. doi: 10.7498/aps.57.6265
    [16] Ge Gui-Xian, Luo You-Hua. Density functional theory study of the structure and electronic properties of MgnOn(n=2—8) clusters. Acta Physica Sinica, 2008, 57(8): 4851-4856. doi: 10.7498/aps.57.4851
    [17] Li Xi-Bo, Luo Jiang-Shan, Guo Yun-Dong, Wu Wei-Dong, Wang Hong-Yan, Tang Yong-Jian. Density functional theory study of the stability, electronic and magnetic properties of Scn, Yn and Lan (n=2—10) clusters. Acta Physica Sinica, 2008, 57(8): 4857-4865. doi: 10.7498/aps.57.4857
    [18] Wang Qing-Lin, Ge Gui-Xian, Zhao Wen-Jie, Lei Xue-Ling, Yan Yu-Li, Yang Zhi, Luo You-Hua. Density functional theory study on the structure and properties of CoBen(n=1—12) clusters. Acta Physica Sinica, 2007, 56(6): 3219-3226. doi: 10.7498/aps.56.3219
    [19] Chen Yu-Hong, Zhang Cai-Rong, Ma Jun. Density functional theory study on the structure and properties of MgmBn(m=1,2;n=1—4) clusters. Acta Physica Sinica, 2006, 55(1): 171-178. doi: 10.7498/aps.55.171
    [20] Tan Ming-Qiu, Tao Xiang-Ming, Xu Xiao-Jun, Cai Jian-Qiu. Density functional theory study on the electronic structure of UAl3 a nd USn3. Acta Physica Sinica, 2003, 52(12): 3142-3149. doi: 10.7498/aps.52.3142
  • 期刊类型引用(4)

    1. 任超,黄益旺. 海洋环境噪声空间相关系数的统计特性. 哈尔滨工程大学学报. 2024(11): 2160-2167 . 百度学术
    2. 任超,黄益旺,夏峙. 宽频带海洋环境噪声矢量场空间相关特性建模. 物理学报. 2022(02): 131-141 . 百度学术
    3. 李明杨,赵航芳,孙超. 风成噪声背景下垂直阵阵列信噪比随声源深度的变化规律. 物理学报. 2022(04): 141-152 . 百度学术
    4. 蒋光禹,孙超,李沁然. 涡旋对深海风成噪声垂直空间特性的影响. 物理学报. 2020(14): 111-125 . 百度学术

    其他类型引用(4)

Metrics
  • Abstract views:  11564
  • PDF Downloads:  109
  • Cited By: 8
Publishing process
  • Received Date:  19 December 2018
  • Accepted Date:  14 May 2019
  • Available Online:  01 August 2019
  • Published Online:  20 August 2019

/

返回文章
返回