1990—2015年中国沿海低地地震高危险性地区的人口暴露研究

冯思琦; 吴清; 沈鎏澄; 王惠铎; 温家洪

doi:10.11899/zzfy20220412

1990—2015年中国沿海低地地震高危险性地区的人口暴露研究

doi: 10.11899/zzfy20220412

1.
上海师范大学, 环境与地理科学学院, 上海 200234
2.
中国地震局地球物理研究所, 北京 100081

基金项目: 国家重点研发计划（2017YFC1503001）；国家自然科学基金（42171080）；中国地震局重大政策理论与实践问题研究课题（CEAZY2019JZ10）

详细信息

作者简介:
冯思琦，女，生于1994年。硕士研究生。主要从事沿海灾害风险管理研究。E-mail：f19821824845@163.com

通讯作者:
温家洪，男，生于1966年。教授，博士生导师。主要从事气候变化与极端事件、自然灾害风险分析与管理研究。E-mail： jhwen@shnu.edu.cn

2 https://go.climatecentral.org/coastaldem/
3 http://www.worldpop.org
计量
- 文章访问数: 318
- HTML全文浏览量: 19
- PDF下载量: 16
- 被引次数: 0
出版历程
- 收稿日期: 2022-01-13
- 刊出日期: 2022-12-31

Population Expousure in the High Seismic Hazard AREA of Low-elevation Coastal Zone in China from 1990 to 2015

1.
School of Environmental and Geographical Sciences, Shanghai Normal University, Shanghai 200234, China
2.
Institute of Geophysics, China Earthquake Administration, Beijing 100081, China

摘要

摘要: 中国沿海低地人口稠密、经济发达，是包括地震在内的自然灾害频发地区。由于已有研究成果缺乏对地震灾害暴露和风险的研究，本文利用CoastalDEM、Worldpop人口和中国地震动参数区划等数据，分析1990—2015年中国沿海低地地震高危险性地区人口暴露时空变化特征。研究结果表明，中国沿海低地地震高危险性地区面积为15.1×10⁴ km²，约占全国沿海低地总面积的70.6%；1990—2015年，沿海低地地震高危险性地区暴露人口不断增加，2015年暴露人口达16869万人，与此同时，城市人口快速增长，而农村人口大幅下降。
- 沿海低地 /
- 地震高危险性地区 /
- 暴露人口 /
- 人口时空变化
Abstract: The low-elevation coastal zone in china is densely populated and economically developed. It is also prone to natural disasters, including earthquakes. However, previous studies had a research gap on the exposure and risks of earthquakes. In this work, we analyzed the spatiotemporal patterns of population exposure in the high seismic hazard area of low-elevation coastal zone in China from 1990 to 2015, based on datasets of CoastalDEM, Worldpop population data and seismic ground motion parameters zonation map of China, etc. Our results show that, the high seismic hazard area of low-elevation coastal zone is 0.151 million km², accounting for 70.6% of the low-elevation coastal zones in China. The population in the high seismic hazard area of low-elevation coastal zones in 2015 was approximated 168.69 million. The exposured population of the high seismic hazard area in low-elevation coastal zones increased continuously, while the urban population grew rapidly and the rural population decreased significantly over the period.
- Low-elevation coastal zone /
- High seismic hazard area /
- Exposure /
- Spatiotemporal change of population
注释:

1) 2 https://go.climatecentral.org/coastaldem/

2) 3 http://www.worldpop.org

HTML全文

引言

近年来，面对日渐严峻的环境与能源问题，新能源的发展和应用受到关注。风能是目前发展迅速、产业前景可观的可再生能源之一，开发利用海上风能是重要的能源选择，研究探索海上风电技术逐渐成为能源发展重心。海上风机（Offshore Wind Turbines，OWTs）在设计时，系统基频应避免与荷载激励频率重叠而发生共振。目前国际上多采用“软-刚”模式进行风机支撑结构设计（Kühn等，1998），即风机基频介于涡轮机的转动频率（1 P频率）和叶轮扫掠频率（2 P/3 P频率）之间。针对单桩式海上风机结构系统自振频率的求解，关键在于塔筒变截面特性以及桩土模式的选择。为克服上述关键难点，各国学者采用了多种方式与方法对其开展了相关研究。对于变截面塔筒，Adhikari等（2012）将塔筒按等截面方式考虑，忽略了风机塔筒弯曲刚度和质量密度的变换特征，将其简化为常数，截面属性不随高度变换，在考虑桩-土相互作用的前提下对风机进行了自振频率求解。杨春宝等（2018）考虑风机塔筒变截面特性和桩-土相互作用，利用有限元法建立了单桩式海上风机结构系统自振频率求解方法。风机下部结构为大直径单桩，但在前期的研究中，为了简化计算，Uściłowska等（1998）忽略了桩-土相互作用，利用底部固定的方式，并采用顶部带有集中质量的Euler-Bernoulli梁模型对风塔等部分浸没的海上结构物进行模拟，推导了结构各阶自振频率的解析解。值得注意的是，风机运行期间，基础作为嵌入土体的一部分，难以避免地受到桩侧土体的作用，因此风机动力特性也会受到影响。Barltrop等（1991）在考虑土体条件的情况下，认为泥面以下管桩等效桩长可由3.5～8倍桩径择定，桩端采用固接方式推导相应的刚度矩阵，进而获得结构风机系统基频；Zuo等（2018）利用非线性Winkler地基模型模拟桩-土作用，发现桩-土相互作用对风机自振特性的影响不可忽略。

塔筒变截面控制方程为四阶变系数偏微分方程，一般其解析解不易得到。李大潜等（1974）利用有限元素法的思想将变截面直梁近似地看成若干等截面直梁的组合体，利用余量法求解等截面直梁的固有频率。Van等（2002）基于Euler-Bernoulli梁理论，将风机塔筒简化为等截面梁，顶部简化为集中质量，得到了风机自振频率。上述研究多通过等效截面的方式对变截面梁进行简单模拟。因此，在海上风机前期设计中，风机系统基频的准确预估仍需可全面考虑塔筒变截面特性的计算方法。微分变换法（Different Transform Method，DTM）是基于Taylor级数求解微分方程的半解析法（赵家奎等，1988），利用充分可微的多项式作为精确解的近似形式，将边界条件和微分方程转化为易于进行计算机编程且多次迭代的代数方程。该方法方便求解，原理简单，目前逐渐应用于结构动力学分析中（Çatal，2008；Balkaya等，2009；Torabi等，2012）。

本文基于Euler-Bernoulli梁（E-B）和Timoshenko梁（T-B）理论，采用微分变换法考虑近海风机塔筒变截面特性，通过三弹簧模型模拟桩-土相互作用，以提供简便可行、结果精准的单桩式海上风机系统频率预估方法，进而为风机基频研究提供参考。

1. 微分变换法（DTM）

设解析函数f（x）的定义域为D，x = x₀表示自变量x在定义域D中的任意一点，对函数f（x）关于自变量x求k阶导数，可得到原函数f（x）在x = x₀处变换后的变换式F（k）为：

$$ F\left( k \right) = \frac{1}{{k!}}{\left( {\frac{{{{\mathrm{d}}^k}f\left( x \right)}}{{{\mathrm{d}}{x^k}}}} \right)_{x = {x_0}}} $$

(1)

根据Taylor级数收敛性质，对原函数f（x）进行Taylor级数展开，可得到F（k）逆变换式f（x）为：

$$ f\left( x \right) = \sum\limits_{k = 0}^\infty {\frac{{{{(x - {x_0})}^k}}}{{k!}}} {\left( {\frac{{{{\mathrm{d}}^k}f\left( x \right)}}{{{\mathrm{d}}{x^k}}}} \right)_{x = {x_0}}} $$

(2)

将式（1）代入式（2）可得:

$$ f\left( x \right) = \sum\limits_{k = 0}^\infty {F\left( k \right){{(x - {x_0})}^k}} $$

(3)

在实际问题中，仅需保证满足精度要求，函数便可采用有限多项式进行求和，即：

$$ f\left( x \right) = \sum\limits_{k = 0}^N {F\left( k \right){{(x - {x_0})}^k}} $$

(4)

式中，N为计算项数，主要取决于所求特征值的精度要求。

微分变换法的基本变换法则如表1所示，依据基本变换法则对动力方程进行DTM转换，原函数中各物理量用小写字母表示，转换后各物理量用与之相对应的大写字母表示。

表 1 DTM基本变换法则

Table 1. Fundamental transformation theorems of DTM

原函数	转换函数
$ f(x) = g(x) \pm h(x) $	$ F(k) = G(k) \pm H(k) $
$ f(x) = c \cdot g(x) $	$ F(k) = c \cdot G(k) $
$ f(x) = g(x) \cdot h(x) $	$ F(k) = \displaystyle\sum\limits_{l = 0}^k {G(l) \cdot H(k - l)} $
$ f(x) = \dfrac{{{{\mathrm{d}}^n}g(x)}}{{{\mathrm{d}}{x^n}}} $	$ F(k) = (k + 1)(k + 2)\cdots(k + n)G(k + n) $
$ f(x) = {x^n} $	$ F(k) = \delta (k - n) = x ,\; \left\{ {x = 1,k = n;x = 0,k \ne n} \right\} $

下载: 导出CSV

| 显示表格

2. 方程建立与求解

单桩式海上风机支撑结构主要包括叶轮-机舱组合件（Rotor Nacelle Assembly，RNA）、塔筒、连接段及单桩。叶轮-机舱组合件由风机叶片、轮毂和机舱组合件组成，其与风机平台之间为变截面塔筒，塔筒顶部直径最小，随着高度降低在平台处直径最大。单桩位于风机底部，与塔筒之间采用连接段和刚性法兰连接，如图1（a）所示。

图 1 海上风机计算模型

Figure 1. Calculation models of OWTs

下载: 全尺寸图片幻灯片

2.1 计算模型建立与简化

近海风机体量相对较小，其连接段长度较塔筒高度而言差距较大，且连接段直径与塔筒底部直径相近。为简化计算，可近似考虑塔筒底部与泥面相接，机舱至泥面高度为变截面梁高，即塔筒总体高度。因此，将风机泥面以上的部分简化为与泥面相接的变截面梁，其中塔顶机舱、轮毂和叶片等效为集中质量m_RNA。泥面以下根据是否考虑桩-土相互作用分为风机底部弹簧模型和风机底部固接模型。风机底部弹簧模型采用三弹簧模型考虑桩-土相互作用，如图1（b）所示。K_L为侧向弹簧刚度，即桩头单位侧移所需的力；K_R为旋转弹簧刚度，即桩头单位旋转所需的力矩；K_LR为耦合弹簧刚度，即桩头力和力矩之间的耦合关系，单位为N。风机底部固接模型忽略桩-土相互作用，按照底部固接的方式考虑，其计算模型如图1（c）所示。基本假设如下：

（1）采用集中质量m_RNA的方式对风机顶部叶轮-机舱组合件进行简化；

（2）将塔筒与连接段并入一起考虑，组合为变截面梁，其中塔筒高度为L_T，连接段高度为L_S。合并之后的变截面梁整体高度为L_T+L_S，并用L表示；

（3）假定连接段壁厚与风机塔筒壁厚一致，用t_T表示整体变截面梁壁厚；

（4）假定连接段截面直径与风机塔筒底部直径相同，用D（x）表示变截面梁任一截面处的直径。

2.2 变截面梁参数确定

变截面梁如图2所示，其中x为梁沿着轴线的坐标，L为梁高，D_t为梁顶部直径，D_b为梁底部直径，t_T为梁厚，A（x）为梁x位置的圆环截面积，I（x）为圆环截面惯性矩。

图 2 变截面梁示意

Figure 2. Schematic diagram of variable section beam

下载: 全尺寸图片幻灯片

根据图2引入变截面梁锥度比$ \beta $：

$$ \beta = \frac{{{D_{{\mathrm{t}}}}}}{{{D_{{\mathrm{b}}}}}} $$

(5)

在x处截面直径为D（x），由相似关系可得：

$$ D\left( x \right) = {D_{{\mathrm{b}}}}\left[ {1 + (\beta - 1)\frac{x}{L}} \right] $$

(6)

变截面锥形梁任一截面处的A（x）与I（x）为：

$$ A(x) = {\text{π} }[D(x){t_{{\mathrm{T}}}} - t_{{\mathrm{T}}}^{2}] $$

(7)

$$ I(x) = \frac{{\text{π} }}{2}D(x)t_{{\mathrm{T}}}^{3} + \frac{{\text{π} }}{8}{D^3(x)}{t_{{\mathrm{T}}}} - \frac{3}{8}{\text{π} }{D^2(x)}t_{{\mathrm{T}}}^{2} - \frac{{\text{π} }}{4}t_{{\mathrm{T}}}^{4} $$

(8)

根据表1分别对式（7）和式（8）进行DTM转换，简化可得：

$$ A(k) = {\text{π} }{t_{{\mathrm{T}}}}\left\{ {{D_{{\mathrm{b}}}}\left[ {(\beta - 1)\delta (k - 1) + \delta (k)} \right] - {t_{{\mathrm{T}}}}\delta (k)} \right\} $$

(9)

$$\begin{aligned} & I(k) = \frac{{\text{π} }}{2}{D_{{\mathrm{b}}}}t_{{\mathrm{T}}}^{3}[(\beta - 1)\delta (k - 1) + \delta (k)] + \frac{{\text{π} }}{8}{D_{{\mathrm{b}}}^{3}}{t_{{\mathrm{T}}}}[{(\beta - 1)^3}\delta (k - 3) + \\ &\qquad 3{(\beta - 1)^2}\delta (k - 2) + 3(\beta - 1)\delta (k - 1) + \delta (k)] - \frac{3}{8}{\text{π} }t_{{\mathrm{T}}}^{2}{D_{{\mathrm{b}}}^{2}}[{(\beta - 1)^2}\delta (k - 2) + \\ &\qquad 2(\beta - 1)\delta (k - 1) + \delta (k)] - \frac{{\text{π} }}{4}t_{{\mathrm{T}}}^{4}\delta (k) \end{aligned} $$

(10)

2.3 基于Euler-Bernoulli梁理论的基本计算

基于E-B理论，建立变截面梁横向振动方程为：

$$ \frac{{{\partial ^2}}}{{\partial {x^2}}}\left[ {EI(x)\frac{{{\partial ^2}y(x,t)}}{{\partial {x^2}}}} \right] + \rho A(x)\frac{{{\partial ^2}y(x,t)}}{{\partial {t^2}}} = 0 $$

(11)

式中，x为空间坐标中的任意一点；t为时间；$ y\left( {x,t} \right) $为梁横向位移；E为梁杨氏模量；$ \rho $为材料密度，数值由材料本身决定。

根据周叮（1996）的研究，当梁发生横向振动时，式（11）的简谐振动反应可表示为：

$$ y\left( {x,t} \right){\text{ = }}Y\left( x \right){{\mathrm{e}}^{{{\mathrm{i}}}\omega t}} $$

(12)

式中，$ \omega $为频率；$ Y\left( x \right) $为梁横向位移形函数，描述了梁空间位移变化。

引入无量纲量$ \xi = x/L $，将式（12）代入式（11）并进行无量纲处理，简化可得：

$$ I(\xi ){Y{''''}}(\xi ) + 2{I{'}}(\xi ){Y{'''}}(\xi ) + {I{''}}(\xi ){Y{''}}(\xi ) - {\varOmega ^2}A(\xi )Y(\xi ) = 0 $$

(13)

式中，$ {\varOmega ^2} = \dfrac{{\rho {\omega ^2}{L^4}}}{{E}} $，下同。

对式（13）进行DTM转换可得：

$$ \begin{aligned} & \sum\limits_{l = 0}^k {I(k - l)} Z(l + 4)(l + 1)(l + 2)(l + 3)(l + 4) + 2\sum\limits_{l = 0}^k {I(k - l + 1)Z(l + 3)(k - l + 1)(l + 1)(l + 2)(l + 3)} + \\ &\qquad \sum\limits_{l = 0}^k {I(k - l + 2)} Z(l + 2)(k - l + 1)(k - l + 2)(l + 1)(l + 2) + \sum\limits_{l = 0}^k {A(k - l)} Z(l){\varOmega ^2}{\text{ = }}0 \end{aligned}$$

(14)

式中，$ Z(l) $为$ Y(\xi ) $的微分变换形式。

由于递推式的周期性，梁横向位移的DTM转换通项表达式$ Z(k + 4) $由式（14）可得：

$$ \begin{aligned} & Z(k{\text{ + 4}}) = \{ {\varOmega ^2}A(1)Z{\text{(}}k{\text{ - 1)}} + {\varOmega ^2}A(0)Z{\text{(}}k{\text{)}} -I(3)[(k - 2)(k - 1)k(k + 1) + \\ &\qquad 6(k - 1)k(k + 1) + 6k(k + 1)]Z(k + 1) - I(2)[(k - 1)k(k + 1)(k + 2) + \\ &\qquad 4k(k + 1)(k + 2) + 2(k + 1)(k + 2)]Z(k + 2) -I(1)[k(k + 1)(k + 2)(k + 3) + \\ &\qquad 2(k + 1)(k + 2)(k + 3)]Z(k + 3)\} /(k + 1)(k + 2)(k + 3)(k + 4)I(0) \end{aligned}$$

(15)

2.4 基于Timoshenko梁理论的基本计算

基于T-B理论，建立变截面梁横向振动方程为：

$$ \frac{\partial }{{\partial x}}\left[ {\kappa A(x)G\left( {\frac{{\partial u\left( {x,t} \right)}}{{\partial x}} - \theta \left( {x,t} \right)} \right)} \right] - \rho A\left( x \right)\frac{{{\partial ^2}u\left( {x,t} \right)}}{{\partial {t^2}}} = 0 $$

(16)

$$ \frac{\partial }{{\partial x}}\left[ {EI(x)\frac{{\partial \theta \left( {x,t} \right)}}{{\partial x}}} \right] + \kappa A\left( x \right)G\left[ {\frac{{\partial u\left( {x,t} \right)}}{{\partial x}} - \theta \left( {x,t} \right)} \right] - \rho I(x)\frac{{{\partial ^2}\theta \left( {x,t} \right)}}{{\partial {t^2}}} = 0 $$

(17)

式中，$ u\left( {x,t} \right) $和$ \theta \left( {x,t} \right) $分别为梁横向位移和转角；G为梁剪切模量；$ \kappa $为截面有效剪切系数。其中，截面有效剪切系数$ \kappa $可由下式确定（Cowper，1966）：

$$ \kappa = \frac{{2(1 + \upsilon )}}{{4 + 3\upsilon }} $$

(18)

式中，$ \upsilon $为泊松比，钢材取0.3。

类似的，式（16）和式（17）变量分离解为：

$$ u(x,t) = U(x){{\mathrm{e}}^{{{\mathrm{i}}}\omega t}},\theta \left( {x,t} \right){\text{ = }}\varphi \left( x \right){{\mathrm{e}}^{{{\mathrm{i}}}\omega t}} $$

(19)

式中，$ U\left( x \right) $和$\varphi \left( x \right)$分别为梁横向位移形函数和转角形函数。

将式（19）代入式（16）、式（17）并进行无量纲处理，简化可得：

$$ A(\xi )\left[ {{U{''}}(\xi ) - {\varphi {'}}\left( \xi \right)} \right] + {A{'}}(\xi )\left[ {{U{'}}(\xi ) - \varphi \left( \xi \right)} \right] + {\alpha _1}A(\xi )U(\xi ) = 0 $$

(20)

$$ I(\xi ){\varphi {''}}\left( \xi \right) + {I{'}}(\xi ){\varphi {'}}\left( \xi \right) + {\alpha _2}A(\xi )\left[ {{U{'}}(\xi ) - \varphi \left( \xi \right)} \right] + {\alpha _3}I(\xi )\varphi \left( \xi \right) = 0 $$

(21)

式中，$ {\alpha _1} = \dfrac{{\rho {\omega ^2}{L^2}}}{{\kappa G}} $，$ {\alpha _2} = \dfrac{{\kappa G{L^2}}}{E} $，$ {\alpha _3} = \dfrac{{\rho {\omega ^2}{L^2}}}{E} $，下同。

对式（20）和式（21）进行DTM转换，简化可得：

$$ \begin{aligned} & \sum\limits_{l = 0}^k {A(l)} (k - l + 1)(k - l + 2)\nu (k - l + 2) + \sum\limits_{l = 0}^k {(l + 1)A(l + 1)(k - l + 1)\nu (k - l + 1)} + \\ &\qquad {\alpha _1}\sum\limits_{l = 0}^k {A(l)} \nu (k - l) - \sum\limits_{l = 0}^k {A(l)} (k - l + 1)\varPhi \left( {k - l + 1} \right) - \sum\limits_{l = 0}^k {(l + 1)A(l + 1)} \varPhi \left( {k - l} \right){\text{ = }}0 \end{aligned}$$

(22)

$$ \begin{aligned} & \sum\limits_{l = 0}^k {I(l)} (k - l + 1)(k - l + 2)\varPhi (k - l + 2) + \sum\limits_{l = 0}^k {(l + 1)I(l + 1)(k - l + 1)\varPhi (k - l + 1)} + \\ &\qquad{\alpha _2}\sum\limits_{l = 0}^k {A(l)} (k - l + 1)\nu (k - l + 1) - {\alpha _2}\sum\limits_{l = 0}^k {A(l)} \varPhi \left( {k - l} \right) + {\alpha _3}\sum\limits_{l = 0}^k {I(l)} \varPhi \left( {k - l} \right){\text{ = }}0 \end{aligned} $$

(23)

式中，$ \nu (k - l) $和$ \varPhi (k - l) $分别为$ U(\xi ) $和$ \varphi (\xi ) $的微分变换形式。

梁横向位移和转角DTM转换的通项表达式$ \nu (k + 2) $和$ \varPhi (k + 2) $由式（22）和式（23）可得：

$$ \begin{aligned} & \nu (k + 2) = [(k + 1)A(0)\varPhi (k + 1) + (k + 1)A(1)\varPhi (k) - {(k + 1)^2}A(1)\nu (k + 1) - {\alpha _1}A(0)\nu (k) - \\ &\qquad {\alpha _1}A(1)\nu (k - 1)]/(k + 1)(k + 2)A(0) \end{aligned} $$

(24)

$$ \begin{aligned} & \varPhi (k + 2) = \{ - {\alpha _2}A(1)k\nu (k) - {\alpha _2}A(0)(k + 1)\nu (k + 1) - {\alpha _3}I(3)\varPhi {{({{k}} - 3)}} - {\alpha _3}{\text{I(2)}}\varPhi {{({{k}} - 2)}} - \\ &\qquad \left[ {(k - 1)(k + 1)I(3) + {\alpha _3}I(1) - {\alpha _3}A(1)} \right]\varPhi (k - 1) - \\ &\qquad \left[ {k(k + 1)I(2) + {\alpha _3}I(0) - {\alpha _2}A(0)} \right]\varPhi (k) - \\ &\qquad {(k + 1)^2}I(1)\varPhi (k + 1)\} /(k + 1)(k + 2)I(0) \end{aligned}$$

(25)

2.5 边界条件与基频求解

底部弹簧模型与固接模型在风机顶部的边界条件相同，而在风机底部存在明显差异。底部弹簧模型采用三弹簧模型模拟桩-土相互作用，Poulos等（1980）利用桩长、桩径和地基水平抗力系数进行了弹簧刚度求解。弹簧刚度与单桩桩头位移之间的关系为（Arany等，2015）：

$$ \left[ \begin{gathered} M \\ Q \\ \end{gathered} \right] = \left[ {\begin{array}{*{20}{c}} {{K_{{\mathrm{R}}}}}&{{K_{{{\mathrm{LR}}}}}} \\ {{K_{{{\mathrm{LR}}}}}}&{{K_{{\mathrm{L}}}}} \end{array}} \right]\left[ \begin{gathered} \theta \\ u \\ \end{gathered} \right] $$

(26)

式中，$ M、Q $分别为单桩桩头弯矩和剪力；$ u、\theta $分别为桩头横向位移和转角；K_L为侧向弹簧刚度；K_R为旋转弹簧刚度；K_LR为耦合弹簧刚度。

当忽略了桩-土相互作用的影响，即风机底部按固接考虑时，其在端点位置（ξ=0）线位移和角位移均为0。当ξ=0，1时，T-B与E-B理论底部弹簧模型和固接模型顶部边界条件、底部边界条件和DTM转换结果如表2所示。以E-B理论底部弹簧模型为例，当ξ=0时，根据式（26）可得到桩底边界条件为：

表 2 边界条件

Table 2. Boundary conditions table

边界条件	物理量	T-B理论	E-B理论
底部弹簧模型	弯矩	$ \dfrac{{E{I_{{\mathrm{b}}}}}}{L}{\varphi {'}}(0) - {K_{\text{R}}}\varphi (0) + {K_{{\text{LR}}}}LU(0) = 0 $	$ \dfrac{{E{I_{{\mathrm{b}}}}}}{{{L^{2}}}}{Y{'}}{'}(0) - \dfrac{{{K_{\text{R}}}}}{L}{Y{'}}(0) + {K_{{\text{LR}}}}Y(0) = 0 $
底部弹簧模型	剪力	$ \kappa {A_{{\mathrm{b}}}}G\left[ {{U{'}}(0) - \varphi (0)} \right] + {K_{{\text{LR}}}}\varphi (0) - {K_{\text{L}}}LU(0) = 0 $	$ \dfrac{{E{I_{{\mathrm{b}}}}}}{{{L^{3}}}}{Y{'''}}(0) + \dfrac{{{K_{{\text{LR}}}}}}{L}{Y{'}}(0) - {K_{\text{L}}}Y(0) = 0 $
底部固接模型	位移	$ U\left( 0 \right) = 0 $	$ Y(0) = 0 $
底部固接模型	转角	$ \varphi \left( 0 \right) = 0 $	$ {Y{'}}(0) = 0 $
风机顶部	弯矩	$ E{I_{{\mathrm{t}}}}{\varphi {'}}(1) = 0 $	$ E{I_{{\mathrm{t}}}}{Y{''}}(1) = 0 $
风机顶部	剪力	$ \kappa {A_{t}}G\left[ {{U{'}}(1) - \varphi (1)} \right] + m{\omega ^2}U(1) = 0 $	$ E{I_{{\mathrm{t}}}}{Y{'''}}(1) + m{\omega ^2}Y(1) = 0 $
DTM底部弹簧模型	弯矩	$ \varPhi (1) - {\eta _{\text{R}}}\varPhi (0) + {\eta _{{\text{LR}}}}\nu (0) = 0 ,\;{\eta _{\text{R}}} = \dfrac{{{K_{\text{R}}}L}}{{E{I_{{\mathrm{b}}}}}} ,\; {\eta _{{\text{LR}}}} = \dfrac{{{K_{{\text{LR}}}}{L^2}}}{{E{I_{{\mathrm{b}}}}}}$	$ 2 Z(2) - {\eta _{\text{R}}}Z(1) + {\eta _{{\text{LR}}}}Z(0) = 0 ,\;{\eta _{\text{R}}} = \dfrac{{{K_{\text{R}}}L}}{{E{I_{{\mathrm{b}}}}}},\;{\eta _{{\text{LR}}}} = \dfrac{{{K_{{\text{LR}}}}{L^2}}}{{E{I_{{\mathrm{b}}}}}}$
DTM底部弹簧模型	剪力	$ \nu (1) - \varPhi (0) + {\psi _{{\text{LR}}}}\varPhi (0) - {\psi _{\text{L}}}\nu (0) = 0,\;{\psi _{{\text{LR}}}} = \dfrac{{{K_{{\text{LR}}}}}}{{\kappa {A_{b}}G}} ,\; {\psi _{\text{L}}} = \dfrac{{{K_{\text{L}}}L}}{{\kappa {A_{b}}G}}$	$ 6 Z(3) + {\eta _{{\text{LR}}}}Z(1) - {\eta _{\text{L}}}Z(0) = 0 ,\;{\eta _{\text{L}}} = \dfrac{{{K_{\text{L}}}{L^3}}}{{E{I_{{\mathrm{b}}}}}} ,\;{\eta _{{\text{LR}}}} = \dfrac{{{K_{{\text{LR}}}}{L^2}}}{{E{I_{{\mathrm{b}}}}}} $
DTM底部固接模型	位移	$ \nu \left( 0 \right) = 0 $	$ Z(0) = 0 $
DTM底部固接模型	转角	$ \varPhi \left( 0 \right) = 0 $	$ Z(1) = 0 $
DTM风机顶部	弯矩	$ \displaystyle\sum\limits_{k = 0}^N {(k + 1)\varphi (k + 1)} = 0 $	$ \displaystyle\sum\limits_{k = 0}^N {(k + 2)(k + 1)Z(k + 2)} = 0 $
DTM风机顶部	剪力	$ \displaystyle\sum\limits_{k = 0}^N {(k + 1)U(k + 1)} - \sum\limits_{k = 0}^N {\varphi (k)} + {\alpha _n}\sum\limits_{k = 0}^N {U(k)} = 0,\;{\alpha _n} = \dfrac{{m{\omega ^2}L}}{{\kappa {A_{{\mathrm{t}}}}G}} $	$ \displaystyle\sum\limits_{k = 0}^N {(k + 3)(k + 2)(k + 1)Z(k + 3)} + {\alpha _m}\sum\limits_{k = 0}^N {Z(k)} = 0 ,\;{\alpha _m} = \dfrac{{m{\omega ^2}{L^{3}}}}{{E{I_{{\mathrm{b}}}}}} $

下载: 导出CSV

| 显示表格

$$ \frac{{E{I_{{\mathrm{b}}}}}}{{{L^{2}}}}{Y{'}}{'}(0) - \frac{{{K_{\text{R}}}}}{L}{Y{'}}(0) + {K_{{\text{LR}}}}Y(0) = 0 $$

(27)

$$ \frac{{E{I_{{\mathrm{b}}}}}}{{{L^{3}}}}{Y{'''}}(0) + \frac{{{K_{{\text{LR}}}}}}{L}{Y{'}}(0) - {K_{\text{L}}}Y(0) = 0 $$

(28)

对边界条件进行DTM转换可得：

$$ 2Z(2) - {\eta _{\text{R}}}Z(1) + {\eta _{{\text{LR}}}}Z(0) = 0 $$

(29)

$$ 6Z(3) + {\eta _{{\text{LR}}}}Z(1) - {\eta _{\text{L}}}Z(0) = 0 $$

(30)

式中，$ {\eta _{\text{R}}} = \dfrac{{{K_{\text{R}}}L}}{{E{I_{{\mathrm{b}}}}}} $，$ {\eta _{{\text{LR}}}} = \dfrac{{{K_{{\text{LR}}}}{L^2}}}{{E{I_{{\mathrm{b}}}}}} $，$ {\eta _{\text{L}}} = \dfrac{{{K_{\text{L}}}{L^3}}}{{E{I_{{\mathrm{b}}}}}} $。

令$ Z\left( 0 \right) = j $，$ Z\left( 1 \right) = f $，其中j和f为不为0的常数。故$ Z\left( 2 \right) = \left[ {{\eta _{\text{R}}}f - {\eta _{{\text{LR}}}}j} \right]/2,Z(3) = \left[ {{\eta _{\text{L}}}j - {\eta _{{\text{LR}}}}f} \right]/6 $。因此关于$ Z(k) $的各项系数均可由含j和f的公式表示。

当ξ=1时，根据风机顶部简化方式可得边界条件为：

$$ E{I_{{\mathrm{t}}}}{Y{''}}(1) = 0 $$

(31)

$$ E{I_{{\mathrm{t}}}}{Y{'''}}(1) + m{\omega ^2}Y(1) = 0 $$

(32)

式中，m为风机顶部集中质量。

对风机顶部边界条件进行DTM转换，再将梁横向位移函数DTM转换式汇总并代入顶部边界条件DTM转换式中可得：

$$ \sum\limits_{l = 0}^N {(k + 2)(k + 1)Z(k + 2)} = X_{11}^Nj + X_{12}^Nf{\text{ = }}0 $$

(33)

$$ \sum\limits_{k = 0}^N {(k + 3)(k + 2)(k + 1)Z(k + 3)} + {\alpha _m}\sum\limits_{k = 0}^N {Z(k)} = X_{21}^Nj + X_{22}^Nf = 0 $$

(34)

式中，$ {\alpha _m} = \dfrac{{m{\omega ^2}{L^{3}}}}{{E{I_{{\mathrm{b}}}}}} $。

变截面塔筒频率特征方程为：

$$ \left[ {\begin{array}{*{20}{c}} {X_{11}^N}&{X_{12}^N} \\ {X_{21}^N}&{X_{22}^N} \end{array}} \right]\left\{ {\begin{array}{*{20}{c}} j \\ f \end{array}} \right\} = \left\{ {\begin{array}{*{20}{c}} 0 \\ 0 \end{array}} \right\} $$

(35)

式中，$ {X}_{11}^{N}、{X}_{12}^{N}、{X}_{21}^{N}、{X}_{22}^{N} $为展开项数为N时求得的关于结构自振频率$ \omega $的多项式函数，形如通项$Z(k + 4) $。

为使式（35）有非零解，其系数矩阵行列式必定为零，即：

$$ \left| {\begin{array}{*{20}{c}} {X_{11}^N}&{X_{12}^N} \\ {X_{21}^N}&{X_{22}^N} \end{array}} \right|{\text{ = }}0 $$

(36)

3. 方法可行性验证

为验证微分变换法对风机系统基频求解的精确性与有效性，现选取Arany等（2016）中风机Irene Vorrink23、Irene Vorrink28和LelyA2作为算例，其规格参数如表3所示。

表 3 风机规格参数

Table 3. Specifications of Irene Vorrink23, 28 and LelyA2

规格参数	Irene Vorrink23	Irene Vorrink28	LelyA2
材料密度/（kg·m⁻³）	7 850	7 850	7 850
RNA质量/t	35.7	35.7	32.0
塔筒高/m	44.5	44.5	41.5
塔筒底部直径/m	3.5	3.5	3.2
塔筒顶部直径/m	1.7	1.7	1.9
塔筒壁厚/mm	13	13	13
材料杨氏模量/GPa	210	210	210
平台高/m	6.0	5.8	4.6

下载: 导出CSV

| 显示表格

将相应风机规格参数代入本方法中，利用MATLAB软件求解风机系统基频，并将计算结果与相应的实测频率和固接频率进行对比，分析不同理论与计算模型对风机系统基频的影响，如表4所示。

表 4 风机系统基频结果验证

Table 4. Verification of OWTs fundamental frequency results

风机	基频/Hz						相对误差/%
	文献实测频率	文献固接频率	T-B理论		E-B理论		T-B理论		E-B理论
	文献实测频率	文献固接频率	弹簧	固接	弹簧	固接	弹簧	固接	弹簧	固接
LelyA2	0.634	0.713	0.642	0.691	0.644	0.688	1.26	3.09	1.57	3.51
Irene Vorrink23	0.563	0.590	0.582	0.616	0.584	0.618	3.37	4.41	3.73	4.75
Irene Vorrink28	0.560	0.586	0.580	0.614	0.581	0.615	3.57	4.79	3.75	4.95

下载: 导出CSV

| 显示表格

首先，基于T-B理论，当考虑桩-土相互作用时，LelyA2风机所得基频的相对误差最小，结果为0.642 Hz，与实测频率0.634 Hz相比较相对误差为1.26%；风机底部按固接计算时所得基频为0.691 Hz，与文献固接频率0.713 Hz之间的相对误差为3.09%。对于Irene Vorrink28风机，当考虑弹簧地基系统时所得相对误差为3.57%，与按底部固接所得相对误差4.79%相比更精准。类似地，LelyA2风机利用E-B理论底部弹簧模型算得基频相对误差为1.57%，底部固接模型相对误差为3.51%，前者所得结果更接近风机实测频率，说明风机底部按考虑桩-土相互作用计算系统基频更精确。

其次，以E-B为例，对Irene Vorrink23与Irene Vorrink28风机基频计算结果进行分析。Irene Vorrink23风机连接段尺寸为5.8 m，相比Irene Vorrink28风机6.0 m的连接段略短。风机底部固接时，Irene Vorrink23与Irene Vorrink28风机算得系统基频的相对误差分别为4.75%和4.95%，当考虑桩-土相互作用时，相对误差分别为3.73%和3.75%，可知利用本方法算得Irene Vorrink23风机系统基频始终较Irene Vorrink28风机更精准。同时，基于T-B理论进行结果验证分析也得到相同的结论。说明利用本方法计算连接段相较于塔筒更短小的风机系统基频更有优势。

最后，对比T-B与E-B理论，当均考虑桩-土相互作用时，T-B与E-B理论计算得到的LelyA2风机基频分别为0.642、0.644 Hz，与实测频率0.634 Hz相比，相对误差分别为1.26%和1.57%；如果风机底部固接，T-B与E-B理论计算得到的LelyA2风机基频分别为0.691、0.688 Hz，相对误差分别为3.09%和3.51%。说明基于T-B与E-B理论计算得到的风机系统基频结果较接近，且基于T-B理论求解结果更精准。Irene Vorrink23、Irene Vorrink28风机较好地验证了该结果。

以上分析表明，利用本方法对近海风机系统基频进行求解，其结果与现场实测频率基本一致，误差控制在合理范围之内，证明微分变换法对于风机系统基频的求解是有效且准确的。基于T-B理论算得基频较E-B理论更接近风机实际频率。风机底部在考虑桩-土相互作用后，算得风机系统基频结果更精准。因此，风机进行前期设计时，建议采用T-B理论并考虑桩-土相互作用预估基频，这对于提高风机系统基频计算精度具有重要意义。

4. 参数影响分析

风机系统基频限制要求高，风机规格参数和弹簧模量等因素会对其产生一定影响。以LelyA2海上风机为例，量化各因素对风机基频的影响。

4.1 计算项数及塔筒高度的影响

采用底部弹簧模型对比分析T-B与E-B理论计算得到微分变换法中计算项数N对风机基频的影响规律，如图3所示。由图3可知，利用微分变换法得到的DTM解随计算项数N的增加均有不同程度的变换。图3（a）中，当塔筒高度为36.5 m，计算项数N为4时，基于T-B理论得到的DTM解为0.780 Hz，与N为10时得到的基频0.769 Hz相比降低了1.41%。而E-B理论计算结果不同，如图3（b）所示，当计算项数N为4时，所得风机基频最小，仅为0.741 Hz；随后基频随计算项数N的增加而增大，当计算项数N为10时达到0.765 Hz，提升了2.68%。T-B和E-B理论计算得到的计算项数N＜10时风机基频有不同程度的上升或下降趋势。另外，当计算项数N为10时，基于T-B与E-B理论算得的DTM解相差0.52%，较接近。当计算项数N增至10后，T-B和E-B理论计算结果变化不明显，曲线逐渐趋于水平，达到收敛点。

图 3 计算项数与塔筒高度对风机基频的影响

Figure 3. Influence of value of N and tower height on fundamental frequency

下载: 全尺寸图片幻灯片

风机基频对塔筒高度变化敏感性较强，由图3（a）可知，当计算项数N为14时，以LelyA2风机塔筒为基础，塔筒高度变化分别为−10、−5、5、10 m，即塔筒高度L_T分别为31.5、36.5、46.5、51.5 m时，由T-B理论算得风机基频分别为0.936、0.768、0.550、0.476 Hz，与LelyA2风机塔筒高度为41.5 m的实测频率0.634 Hz相较分别变化了47.63%、21.14%、−13.25%、−24.92%。同样地，对应于E-B理论，其DTM解在计算项数N为14且塔筒高度分别为31.5、36.5、46.5、51.5 m时，风机基频分别变化了47.00%、20.66%、−13.56%、−25.24%。这表明随着塔筒高度的增加，风机基频逐渐降低。对应于塔筒高度变化相同幅度的情况，塔筒高度减小对基频的影响更大。

分析表明，利用微分变换法进行基频求解时，DTM解在计算项数N为10及以上时趋于水平收敛，所得DTM解已基本满足精度要求且T-B与E-B理论计算的基频结果较接近。由于风机基频会随着塔筒高度的增加而减小，在满足安全性和经济性的前提下可适当通过调整塔筒高度实现对风机系统基频的调控。

4.2 锥度比的影响

T-B与E-B理论计算得到的变截面梁锥度比与塔筒直径变换的基频变换曲线如图4所示。由图4可知，T-B与E-B理论计算结果均表明风机基频随着锥度比的增加逐渐降低。T-B理论计算结果表明，当塔筒顶部直径D_T为2 m时，随着锥度比的增加，塔筒底部直径逐渐减小，风机系统基频明显降低。本风机实际锥度比为0.6，当锥度比分别为0.4和0.8时，算得基频分别为1.072、0.519 Hz，与原基频相比分别变换了69.09%和−18.14%。说明塔筒顶部直径不变时，锥度比的降低对风机基频有更大的提升作用。另一方面，在塔筒顶部直径一定的前提下，锥度比的增加意味着塔筒底部直径减小，说明风机系统基频随着塔筒底部直径的减小而降低。E-B理论计算结果表明，当锥度比为0.4时，基于E-B理论算得风机基频为1.037 Hz，小于T-B理论计算结果。当锥度比＞0.6时，随着锥度比的增加，T-B与E-B理论计算得到的变换曲线趋于重合，风机基频基本相同。同时，基于T-B理论得到基频随锥度比的增加而降低的幅度相对较大。

图 4 锥度比对风机基频的影响

Figure 4. Influence of taper ratio on fundamental frequency

下载: 全尺寸图片幻灯片

分析表明，锥度比与塔筒底部直径会影响风机基频。随着塔筒锥度比的增大，T-B与E-B理论算得风机基频均有不同程度的降低，且前者降低幅度较大。基于锥度比变换关系，当塔筒顶部直径一定时，塔筒底部直径增加使风机基频明显提升。因此，对风机进行前期设计时，可通过适当调整变截面梁锥度比和塔筒底部直径调控风机系统基频。

4.3 弹簧刚度的影响

三弹簧刚度与风机基频关系曲线如图5所示。为充分考虑弹簧刚度对风机基频的影响，首先分别取3个弹簧刚度的较小值，选定其中2个弹簧刚度不变，研究另一弹簧刚度变化对风机基频的影响程度；然后通过比较选定其中影响程度较大的弹簧，对除该弹簧以外的2个弹簧取刚度较大值进行研究，得到关于影响程度最大的弹簧在整体范围内对风机基频的影响规律。

图 5 三弹簧刚度对基频的影响

Figure 5. Effect of three spring stiffness on fundamental frequency

下载: 全尺寸图片幻灯片

由Arany等（2016）可知，旋转弹簧刚度K_R为20～210 ×10⁹ N∙m/rad，侧向弹簧刚度K_L为0.5～2.0 ×10⁹ N/m，耦合弹簧刚度K_LR为−14～−4 ×10⁹ N。因此，旋转弹簧刚度K_R、侧向弹簧刚度K_L、耦合弹簧刚度K_LR较小值分别为20.0 ×10⁹ N∙m/rad、0.5 ×10⁹ N/m、−14.0 ×10⁹ N。

由图5（a）可知，T-B与E-B理论计算得到的旋转弹簧刚度与风机基频关系曲线变化趋势基本一致。以E-B理论计算结果为例，当K_LR为−14.0 ×10⁹ N，K_L为0.5 ×10⁹ N/m时，风机基频随旋转弹簧刚度K_R的增加逐渐增大，但增幅较小，增大了约0.001 Hz，变化了0.01%。说明旋转弹簧刚度K_R对风机基频的影响较小。

由图5（b）可知，随着侧向弹簧刚度的增加，风机基频不断变大，变化幅度相较于旋转弹簧刚度稍大。基于E-B理论算得基频结果变换了0.006 Hz，E-B与T-B理论计算得到的基频变换曲线变化趋势较相近。同样地，耦合弹簧刚度增加也会使风机基频变大，其对风机基频的影响程度较前二者更明显。

由图5（c）可知，T-B与E-B理论计算得到的风机基频在耦合弹簧刚度增加的情况下其增长幅度较大，分别变化了5.96%和4.88%。分析表明，耦合弹簧刚度变化对风机基频的影响最突出。

侧向弹簧刚度和旋转弹簧刚度最大值分别为2.0 ×10⁹ N/m和210.0 ×10⁹ N∙m/rad，由图5（d）可知，随着耦合弹簧刚度的增加风机基频逐渐增加，但整体变化幅度较小，E-B理论计算得到风机基频提高了约0.002 Hz，T-B理论计算得到风机基频增长了0.001 Hz。

分析表明，任一弹簧刚度增加均会使风机基频增加，其中耦合弹簧刚度的影响最突出，且在侧向弹簧刚度和旋转弹簧刚度取较小值时，影响更显著。

4.4 风机基频参数量化分析

综上所述，塔筒高度、塔筒直径、弹簧刚度均会对风机系统基频产生不同程度的影响。进一步分析参数类型对风机基频影响程度，结果如图6所示。由图6可知，以LelyA2风机塔筒高度、塔筒底部直径、耦合弹簧刚度为研究对象，分别对该风机相应参数进行20%的增减，将所得结果与LelyA2风机实测频率0.634 Hz进行比较，得出相对误差，进而判定对风机基频影响程度最大的参数类型。由图6可知，风机基频受塔筒高度变化的敏感性最强，受耦合弹簧刚度变化的敏感性最弱，塔筒底部直径介于二者之间。当塔筒高度降低20%后，相对误差为42.29%，对基频影响程度最大；塔筒底部直径提升20%对基频的影响略大于降低的情况；将耦合弹簧刚度提升20%后算得相对误差为3.74%，对基频影响最小。

图 6 风机基频参数量化分析

Figure 6. Analysis of fundamental frequency parameter results

下载: 全尺寸图片幻灯片

5. 结语

本文基于Euler–Bernoulli梁和Timoshenko梁理论，将风机塔筒和连接段组合为变截面梁，推导了变截面梁横向振动方程。风机底部根据是否考虑桩-土相互作用建立了底部弹簧模型和底部固接模型，并利用微分变换法对变截面梁横向振动方程进行求解，得到了单桩式海上风机与弹簧地基系统相连及底部固接情况下系统基频简化求解方法。通过实际风机实测数据结果对本文预估频率方法进行有效性验证。

（1）与底部固接模型相比，考虑风机塔筒的变截面特性和桩-土相互作用影响算得风机系统基频更接近于实测数据。同时，基于Timoshenko梁对风机系统基频预估更精准。

（2）微分变换法可用于变系数微分方程计算，与其他方法相比，本方法耗费时间和精力较少，更适用于连接段相较塔筒而言较小的近海风机结构。当计算项数N为10及以上时，所得DTM解趋于水平收敛，DTM解已基本满足精度要求。

（3）风机系统基频对塔筒高度减小的敏感性强。随着塔筒高度的减小，基频逐渐增加。

（4）基于变截面梁锥度比的变化关系，塔筒底部直径减小时，风机系统基频明显减小。

（5）在三弹簧地基系统中，弹簧刚度的增加对风机系统基频有提升作用，其中耦合弹簧刚度对风机基频的影响最突出。

图 1 中国沿海低地地震高危险性地区空间分布

Figure 1. Spatial distribution in the high seismic hazard area of low-elevation coastal zones in China

下载: 全尺寸图片幻灯片

图 2 1990—2015年中国沿海低地地震高危险性地区人口时间变化

Figure 2. Population changes in the high seismic hazard area of low-elevation coastal zones in China from 1990 to 2015

下载: 全尺寸图片幻灯片

图 3 中国沿海低地地震高危险性地区人口密度时空变化

Figure 3. Spatiotemporal changes of population density in the high seismic hazard area of low-elevation coastal zones in China

下载: 全尺寸图片幻灯片

表 1 地震动峰值加速度与地震烈度的对比

Table 1. Comparison table between peak acceleration of ground motion and seismic intensity

地震烈度/度	地震动峰值加速度/g
Ⅶ	0.1、0.15
Ⅷ	0.2、0.3
Ⅸ	0.4

下载: 导出CSV

表 2 1990—2015年沿海低地地震高危险性地区不同烈度区的暴露人口

Table 2. Population exposed to different intensity zones in the high seismic hazard area of low-elevation coastal zones in China from 1990 to 2015

时间/年	烈度区	暴露人口/万人	占比/%	人口密度/（人·km⁻²）
1990	Ⅶ度区	8 838	84.2	698
	Ⅷ度区	1 641	15.6	868
	Ⅸ度区	20.8	0.2	42
2000	Ⅶ度区	10 990	83.4	867
	Ⅷ度区	2 138	16.2	1 131
	Ⅸ度区	47.3	0.4	95
2010	Ⅶ度区	13 543	84.7	1 069
	Ⅷ度区	2 406	15.0	1 273
	Ⅸ度区	47.3	0.3	95
2015	Ⅶ度区	14 265	84.6	1 126
	Ⅷ度区	2 557	15.2	1 353
	Ⅸ度区	46.5	0.3	93

下载: 导出CSV

表 3 1990—2015年中国沿海低地地震高危险性地区不同烈度区的城乡人口变化

Table 3. Urban and rural population changes in different intensity zones of the high seismic hazard area in low-elevation coastal zones, China from 1990 to 2015

烈度区	总人口变化量/万人	城市人口变化量/万人	农村人口变化量/万人
Ⅶ度区	5 427（增幅61.4%）	8 143（增幅568%）	−2 716（减幅36.7%）
Ⅷ度区	916（增幅55.8%）	1 320（增幅259.8%）	−404（减幅35.7%）
Ⅸ度区	25.7（增幅124.2%）	29.2（增幅695.2%）	−3.5（减幅20.7%）

下载: 导出CSV

参考文献(38)

段光贤, 曲则生, 余安东, 1980. 1979年溧阳6级地震震害调查. 上海地质, （1）: 70—82.

高战武, 缑亚森, 钟慧等, 2021. 中国东部海域断裂构造格架与地震活动研究. 震灾防御技术, 16（1）: 11—18 doi: 10.11899/zzfy20210102

Gao Z. W. , Gou Y. S. , Zhong H. , et al. , 2021. Fault structure frame and seismicity in the sea on the EastSide of Chinese mainland. Technology for Earthquake Disaster Prevention, 16(1): 11—18. (in Chinese) doi: 10.11899/zzfy20210102

李巧萍, 赵敏, 2020. 《中国地震动参数区划图》施行4周年. 防灾博览, （3）: 19—21. doi: 10.3969/j.issn.1671-6310.2020.03.005

李小军, 李娜, 陈苏, 2021. 中国海域地震区划及关键问题研究. 震灾防御技术, 16（1）: 1—10 doi: 10.11899/zzfy20210101

Li X. J. , Li N. , Chen S. , 2021. Study on seismic zoning in China sea area and its key issues. Technology for Earthquake Disaster Prevention, 16(1): 1—10. (in Chinese) doi: 10.11899/zzfy20210101

梁亚婷, 温家洪, 杜士强等, 2015. 人口的时空分布模拟及其在灾害与风险管理中的应用. 灾害学, 30（4）: 220—228

Liang Y. T. , Wen J. H. , Du S. Q. , et al. , 2015. Spatial-temporal distribution modeling of population and its applications in disaster and risk management. Journal of Catastrophology, 30(4): 220—228. (in Chinese)

刘光鼎, 1992. 中国海地球物理场和地球动力学特征. 地质学报, 66（4）: 300—314

Liu G. D. , 1992. Features of geophysical fields and geodynamics of the sea areas of China. Acta Geologica Sinica, 66(4): 300—314. (in Chinese)

潘建雄, 黄日恒, 1994.1918年南澳地震的破坏强度及震中位置. 华南地震, 14（2）: 17—23

Pan J. X. , Huang R. H. , 1994. To discuss the damage intensity and the epicentral location of nanao earthquake in 1918. South China Journal of Seismology, 14(2): 17—23. (in Chinese)

潘顺, 杜士强, 徐慧等, 2016. 长三角地区沿海低地及其人口的时空变化分析. 地域研究与开发, 35（4）: 161—165 doi: 10.3969/j.issn.1003-2363.2016.04.030

Pan S. , Du S. Q. , Xu H. , et al. , 2016. Spatio-Temporal analysis on low-elevation low-elevation coastal zone and its population in the Yangtze River Delta Region. Areal Research and Development, 35(4): 161—165. (in Chinese) doi: 10.3969/j.issn.1003-2363.2016.04.030

彭艳菊, 孟小红, 吕悦军等, 2008. 我国近海地震活动特征及其与地球物理场的关系. 地球物理学进展, 23（5）: 1377—1388

Peng Y. J. , Meng X. H. , Lü Y. J. , et al. , 2008. The seismicity of China offshore seas and its relationship with geophysical fields. Progress in Geophysics, 23(5): 1377—1388. (in Chinese)

温家洪, 袁穗萍, 李大力等, 2018. 海平面上升及其风险管理. 地球科学进展, 33（4）: 350—360

Wen J. H. , Yuan S. P. , Li D. L. , et al. , 2018. Sea level rise and its risk management. Advances in Earth Science, 33(4): 350—360. (in Chinese)

谢卓娟, 李山有, 吕悦军等, 2020. 中国海域及邻区统一地震目录及其完整性分析. 地震地质, 42（4）: 993—1019 doi: 10.3969/j.issn.0253-4967.2020.04.015

Xie Z. J. , Li S. Y. , Lü Y. J. , et al. , 2020. Unified earthquake catalog for China's seas and adjacent regions and its completeness analysis. Seismology and Geology, 42(4): 993—1019. (in Chinese) doi: 10.3969/j.issn.0253-4967.2020.04.015

赵琳, 韩群, 孙俊仁, 1990. 江苏常熟-太仓5.1级地震经济损失分析. 国际地震动态, （7）: 6—10

Zhao L. , Han Q. , Sun J. R. , 1990. Analysis of economic losses caused by the Changshu-Taicang earthquake of M5.1 in Jiangsu Province. Progress in Earthquake Sciences, (7): 6—10. (in Chinese)

中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会, 2016. GB 18306—2015 中国地震动参数区划图. 北京: 中国标准出版社.

General Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China, Standardization Administration of the People's Republic of China, 2016. GB 18306—2015 Seismic ground motion parameters zonation map of China. Beijing: Standards Press of China. (in Chinese)

Brown S., Nicholls R. J., Woodroffe C. D., et al., 2013. Sea-level rise impacts and responses: a global perspective. In: Finkl C. W., ed., Coastal Hazards. Dordrecht: Springer, 117—149.

Dou Y. Y. , Huang Q. X. , He C. Y. , et al. , 2018. Rapid population growth throughout Asia’s earthquake-prone areas: a multiscale analysis. International Journal of Environmental Research and Public Health, 15(9): 1893. doi: 10.3390/ijerph15091893

Edmonds D. A. , Caldwell R. L. , Brondizio E. S. , et al. , 2020. Coastal flooding will disproportionately impact people on river deltas. Nature Communications, 11(1): 4741. doi: 10.1038/s41467-020-18531-4

Fang Y. Q. , Du S. Q. , Scussolini P. , et al. , 2018. Rapid population growth in Chinese floodplains from 1990 to 2015. International Journal of Environmental Research and Public Health, 15(8): 1602. doi: 10.3390/ijerph15081602

Freire S. , Aubrecht C. , 2012. Integrating population dynamics into mapping human exposure to seismic hazard. Natural Hazards and Earth System Sciences, 12(11): 3533—3543. doi: 10.5194/nhess-12-3533-2012

He C. Y. , Huang Q. X. , Dou Y. Y. , et al. , 2016. The population in China's earthquake-prone areas has increased by over 32 million along with rapid urbanization. Environmental Research Letters, 11(7): 074028. doi: 10.1088/1748-9326/11/7/074028

He C. Y. , Huang Q. X. , Bai X. M. , et al. , 2021. A global analysis of the relationship between urbanization and fatalities in earthquake-prone areas. International Journal of Disaster Risk Science, 12(6): 805—820. doi: 10.1007/s13753-021-00385-z

Huang Q. X. , Meng S. T. , He C. Y. , et al. , 2019. Rapid urban land expansion in earthquake-prone areas of China. International Journal of Disaster Risk Science, 10(1): 43—56. doi: 10.1007/s13753-018-0207-4

Kulp S. A. , Strauss B. H. , 2018. CoastalDEM: a global coastal digital elevation model improved from SRTM using a neural network. Remote Sensing of Environment, 206: 231—239. doi: 10.1016/j.rse.2017.12.026

Lang C. , Gao M. T. , Wu G. C. , et al. , 2019. The concentration of population and GDP in high earthquake risk regions in China: temporal–spatial distributions and regional comparisons from 2000 to 2010. Pure and Applied Geophysics, 176(10): 4161—4175. doi: 10.1007/s00024-019-02126-2

Li Y. X. , Gong P. , Zhou Y. Y. , et al. , 2020. Mapping global urban boundaries from the global artificial impervious area (GAIA) data. Environmental Research Letters, 15(9): 094044. doi: 10.1088/1748-9326/ab9be3

Lichter M. , Vafeidis A. T. , Nicholls R. J. , et al. , 2011. Exploring data-related uncertainties in analyses of land area and population in the “low-elevation coastal zone” (LECZ). Journal of Coastal Research, 27(4): 757—768.

Liu J. L. , Wen J. H. , Huang Y. Q. , et al. , 2015. Human settlement and regional development in the context of climate change: a spatial analysis of low-elevation low-elevation coastal zones in China. Mitigation and Adaptation Strategies for Global Change, 20(4): 527—546. doi: 10.1007/s11027-013-9506-7

Lv Y. M. , Li W. J. , Wen J. H. , et al. , 2021. Population pattern and exposure under sea level rise: low-elevation low-elevation coastal zone in the Yangtze River Delta, 1990-2100. Climate Risk Management, 33: 100348. doi: 10.1016/j.crm.2021.100348

McGranahan G. , Balk D. , Anderson B. , 2007. The rising tide: assessing the risks of climate change and human settlements in low-elevation low-elevation coastal zones. Environment and Urbanization, 19(1): 17—37. doi: 10.1177/0956247807076960

Micheal O., Glavovic B., Hinkel J., et al., 2019. Sea level rise and implications for low lying islands, coasts and communities. In: IPCC, ed., IPCC Special Report on the Ocean and Cryosphere in a Changing Climate. Cambridge: Cambridge University Press.

Neumann B. , Vafeidis A. T. , Zimmermann J. , et al. , 2015. Future coastal population growth and exposure to sea-level rise and coastal flooding - a global assessment. PLoS One, 10(3): e0118571. doi: 10.1371/journal.pone.0118571

Nicholls R. J. , Cazenave A. , 2010. Sea-level rise and its impact on coastal zones. Science, 328(5985): 1517-1520. doi: 10.1126/science.1185782

Pesaresi M., Ehrlich D., Kemper T., et al., 2017. Atlas of the human planet 2017: global exposure to natural hazards. EUR 28556 EN. Luxembourg: Publications Office of the European Union.

Poljanšek K., Ferrer M. M., De Groeve T., et al., 2017. Science for disaster risk management 2017: knowing better and losing less. Luxembourg: Publications Office of the European Union.

Shi P. J. , 2016. Natural disasters in China. Berlin, Heidelberg: Springer.

Stevens F. R. , Gaughan A. E. , Linard C. , et al. , 2015. Disaggregating census data for population mapping using random forests with remotely-sensed and ancillary data. PLoS One, 10(2): e0107042. doi: 10.1371/journal.pone.0107042

UNISDR, 2015. Making development sustainable: the future of disaster risk management. global assessment report on disaster risk reduction. Geneva, Switzerland: United Nations Office for Disaster Risk Reduction (UNISDR).

Vafeidis A., Neumann B., Zimmerman J., et al., 2011. MR9: analysis of land area and population in the low-elevation coastal zone (LECZ). London: Government Office for Science.

Yang X. C. , Yao C. M. , Chen Q. , et al. , 2019. Improved estimates of population exposure in Low-Elevation Coastal Zones of China. International Journal of Environmental Research and Public Health, 16(20): 4012. doi: 10.3390/ijerph16204012

施引文献

资源附件(0)

访问统计