Ameliorative effect of oyster enzymatic products on glucocorticoid-induced osteoporosis in rats
-
摘要:
牡蛎具有多种生物活性，随着对其活性开发的深入研究，发现其在骨质疏松的改善方面有巨大潜力。以牡蛎酶解产物 (Oyster enzymatic products, OP) 为研究对象，探讨其对糖皮质激素性骨质疏松症 (Glucocorticoids osteoporosis, GIOP) 的缓解作用。通过建立地塞米松诱导的GIOP模型，使用不同剂量的OP进行干预，从血清生化指标、微型CT (Micro-CT)、三维重建及组织病理切片等方面评估其对大鼠骨质疏松症的缓解作用。结果发现，与模型组对比，OP低、中、高剂量组的血清中成骨细胞标志物骨保护素、I型前胶原氨基端肽和睾酮平均提高了13.67%、11.51%和17.80% (p<0.05)；破骨细胞标志物基质金属蛋白酶、组织蛋白酶K和抗酒石酸酸性磷酸酶活性平均下降了15.44%、31.20%和14.43% (p<0.05)；骨密度测定、Micro-CT和股骨组织病理分析结果显示，OP高剂量组的骨密度、骨小梁面积占比和骨盐沉积面积分别比模型组增加了56.31%、18.59%和22.51% (p<0.05)。综上，OP可通过提高成骨细胞合成骨胶原能力来加速新骨的形成，降低骨基质降解速度和骨代谢的速度，减少骨量流失，并提高性激素水平从而增加大鼠的骨密度、骨体积分数，增强骨的钙盐沉积能力，改善地塞米松诱导的糖皮质激素性骨质疏松症状。
Abstract:Oyster has a variety of biological activities, and with the intensive research on the development of its activities, the great potential for the improvement of osteoporosis has been found. In this paper, the oyster enzymatic products (OP) were used to investigate their mitigating effects on glucocorticoid-type osteoporosis. A dexamethasone-induced glucocorticoids osteoporosis (GIOP) model was established, and different doses of OP were used to evaluate the mitigating effects of OP on the osteoporosis in rats in terms of serum biochemical indexes, Micro-CT, three-dimensional reconstruction, and histopathological sections. The results show that compared with the model group, the serum content of osteoprotegerin, Type I procollagen amino-terminal peptide and testosterone, which are the osteoclast markers, increased by an average of 13.67%, 11.51% and 17.80% in the low, medium and high dose groups, respectively (p<0.05); but the serum content of matrix metalloproteinases, histone K and anti-tartrate acid phosphatase, which are the osteoclast markers, decreased by 15.44%, 31.20% and 14.43%, respectively (p<0.05). The results of bone densitometry and histopathological analysis of Micro-CT and femur indicate that compared with those of the model group, the bone mineral density, the percentage of bone trabecular area and the area of bone salt deposition in the high-dose group increased by 56.31%, 18.59% and 22.51%, respectively (p<0.05). In conclusion, OP can accelerate the formation of new bone by increasing the ability of osteoblasts to synthesize collagen, decrease the rate of bone matrix degradation and bone metabolism, decrease the loss of bone mass, increase the level of sex hormones, increase the BMD, bone volume fraction, enhance the ability of bone to deposit bone calcium salts in rats, thus improve the symptoms of glucocorticoid-induced osteoporosis induced by dexamethasone.
-
Keywords:
- Oysters /
- Enzymatic hydrolysate /
- Osteoporosis /
- Glucocorticoids /
- Bone mineral density
-
台风是最严重的自然灾害之一，对渔业的影响尤为显著。中国南海台风频发[1]，过去几十年来，台风过境严重威胁网箱养殖设施的安全[2-3]，给水产养殖业造成了巨大损失。海南省马袅湾网箱养殖区共包含72口HDPE (高密度聚乙烯) 圆形深水网箱[4-5]，总养殖容积约12 943 m3[6]，是海南省重要的网箱养殖区。马袅湾为“U”字形半封闭海湾，位于海南省澄迈县西北部，毗邻琼州海峡西南侧。由于琼州海峡是台风过境海南的高频区域，加之马袅湾网箱养殖区的特殊地理位置，已成为海南省内受台风影响最为严重的网箱养殖区之一。因此，分析不同路径和强度的台风对马袅湾网箱养殖区的影响，并建立相应的预警体系具有重要的现实意义。
目前有关台风对网箱养殖区影响的研究多集中于调查分析和数据统计方面。姜雨青[7]基于1989—2018年西太平洋台风路径、设施渔业调研数据对我国沿海省市台风与设施渔业的关联性进行了评估；黄滨等[8]通过实地调查，分析了台风“米雷”在山东省网箱养殖业的灾害影响程度；张润川[9]指出台风浪导致的网箱破坏是福建省养殖渔业的重要致灾因子。但目前针对网箱养殖区的台风影响分析和灾害预警研究还极为匮乏。
由于过境台风的路径和强度不同，网箱养殖区的台风浪与流场特征存在明显差异[10-12]。其中，台风浪是评估网箱养殖区是否受到影响的关键因素。研究表明，台风浪可能导致网箱结构变形、系泊系统失效及网衣损坏等灾害[13]；Zhang等[14]通过数值模拟和现场观测发现，循环波浪载荷引起的疲劳断裂是网箱受损的重要原因；Li等[15]分析了有效波高对系泊绳张力的影响，发现最大系泊绳张力随有效波高的增加而迅速增大；Yin等[16]指出，有效波高超过4 m的海浪已属于灾害性海浪，对网箱结构造成严重威胁；孙熊雄[17]基于BP神经网络建立了网箱致灾因子与破坏程度之间的定量关系，认为波高是导致锚绳断裂、浮架应力增大的主要因素，而流速的影响相对较小。但现有研究结果也表明，较大的流速也可能导致网箱失效或破坏[18-19]。基于以上研究，本研究选择有效波高作为马袅湾网箱养殖区受影响的关键评估指标，同时考虑最大流速的影响。在波高分级方面，参考已有研究[20]，将网箱养殖区受影响程度分为3个等级：当养殖区内有效波高最大值HS_max≥5.5 m时，网箱养殖区受到严重影响，需立即预警；当4 m≤HS_max<5.5 m时，网箱养殖区受到中等程度影响，需快速预警；当2.5 m≤HS_max<4 m时，网箱养殖区受到轻微影响，需适当预警。在流速方面，结合相关研究[21]、抗风浪深水网箱养殖技术规程及马袅湾网箱设计标准，当养殖区内最大流速Vmax≥1 m·s−1时，视为网箱养殖区受到严重影响；当Vmax<1 m·s−1时，以波高的影响为主，网箱养殖区相对安全。本研究以马袅湾大型网箱养殖区为主要研究区域，考虑到台风过境过程中风暴潮与波浪的较强相互作用[22-25]，基于Holland台风模型和Delft3D软件构建了适用于该海域的波流耦合模型，并对模型进行验证。通过数值模拟方法分析了不同路径和强度台风过程对马袅湾网箱养殖区的有效波高及流速的影响规律，确定了不同强度台风过境时的预警路径范围，以期为该区域网箱养殖业的防灾减灾工作提供科学依据和预警参考。
1. 数值模型的建立与验证
1.1 台风模型
风场数据的准确性是决定台风浪模拟精度的关键因素。由于台风过境过程中天气条件恶劣，实测数据获取困难，而欧洲中期天气预报中心再分析数据 (ERA)、交叉校准多平台数据集 (CCMP) 和美国国家环境预报中心 (NCEP) 等再分析数据普遍存在时空分辨率较低的问题，尤其在台风中心附近区域，其最大风速往往显著小于实际观测值。为此，目前普遍采用参数化台风模型[26-28]。本研究选用国际通用的Holland[29]台风模型来计算风场和气压场。其台风气压场分布公式为：
$$ P\left( r \right) {\mathrm{=}} {P_c} {\mathrm{+}} \left( {{P_n} {\text{−}} {P_c}} \right) \cdot \exp {\left( { {\text{−}} \dfrac{{{R_{\max }}}}{r}} \right)^B} $$ (1) 台风风场分布公式为：
$$ \begin{array}{c} {V_g}\left( r \right) {\mathrm{=}} \\ \sqrt {\left( {{P_n} {\text{−}} {P_c}} \right)\dfrac{B}{{{\rho _A}}}{{\left( {\dfrac{{{R_{\max }}}}{r}} \right)}^B} \cdot \exp {{\left( { {\text{−}} \dfrac{{{R_{\max }}}}{r}} \right)}^B} {\mathrm{+}} {{\left( {\dfrac{{rf}}{2}} \right)}^2}} {\text{−}} \dfrac{{rf}}{2} \end{array} $$ (2) 式中：r表示计算点与台风中心的距离；Pn表示台风外围气压；Pc表示台风中心气压；B表示台风形状参数；Rmax表示台风最大风速半径；f表示柯氏系数；ρA表示空气密度。
参数B采用Vickery[30]的参数经验公式，计算公式为：
$$ B {\mathrm{=}} 1.881 {\text{−}} 0.005\;57{R_{\max }} {\text{−}} 0.012\;95\varphi $$ (3) 最大风速半径采用Willoughby[31]提出的风速半径公式，计算公式为：
$$ {R_{\max }} {\mathrm{=}} 51.6\exp \left( { {\text{−}} 0.022\;3{V_{\max }} {\mathrm{+}} 0.028\;1\varphi } \right) $$ (4) 式中：φ表示地理纬度；Vmax表示台风中心最大风速。
台风中心位置、中心气压、最大风速等信息采用中国气象局热带气旋资料中心所提供的CMA最佳路径集[32-33]。
1.2 水动力模型
在正交曲线坐标系下，水动力基本方程[34]为：
连续性方程：
$$\begin{array}{c} \dfrac{1}{{\sqrt {{G_{\xi \xi }}} \sqrt {{G_{\eta \eta }}} }}\left( {\dfrac{{\partial \left[ {\left( {d {\mathrm{+}} \zeta } \right)U\sqrt {{G_{\eta \eta }}} } \right]}}{{\partial \xi }}} { {\mathrm{+}} \dfrac{{\partial \left[ {\left( {d {\mathrm{+}} \zeta } \right)V\sqrt {{G_{\xi \xi }}} } \right]}}{{\partial \eta }}} \right) {\mathrm{=}}\\ ({\text{d {\mathrm{+}} }}\zeta )Q {\text{−}} \dfrac{{\partial \zeta }}{{\partial t}} \end{array} $$ (5) $$ Q {\mathrm{=}} \left( {d {\mathrm{+}} \zeta } \right)\int_{ {\text{−}} 1}^0 {\left( {{q_{{\text{in}}}} {\text{−}} {q_{{\text{out}}}}} \right)d\sigma {\mathrm{+}} P {\text{−}} E} $$ (6) 式中：$ \sqrt{{G}_{\xi \xi }} $、$ \sqrt{{G}_{\eta \eta }} $分别表示ξ、η上的坐标转换系数；U、V分别表示ξ、η方向上的平均速度；ξ表示垂向上的自由面坐标；d表示垂向上的水底坐标 (即空间参考平面以下的水深)；Q表示单位面积流量；P和E分别表示降水量和蒸发量，由于研究的是海湾水动力问题，不考虑降水量和蒸发量，故此处为0。
动量方程：
$$\begin{array}{c} \dfrac{{\partial u}}{{\partial t}} {\mathrm{+}} \dfrac{u}{{\sqrt {{G_{\xi \xi }}} }}\dfrac{{\partial u}}{{\partial \xi }} {\mathrm{+}} \dfrac{v}{{\sqrt {{G_{\eta \eta }}} }}\dfrac{{\partial u}}{{\partial \eta }}{\text{ {\mathrm{+}} }}\dfrac{w}{{d {\mathrm{+}} \zeta }}\dfrac{{\partial u}}{{\partial \sigma }} {\text{−}} \\ \dfrac{{vv}}{{\sqrt {{G_{\xi \xi }}} \sqrt {{G_{\eta \eta }}} }}\dfrac{{\partial \sqrt {{G_{\eta \eta }}} }}{{\partial \xi }} {\mathrm{+}} \dfrac{{uv}}{{\sqrt {{G_{\xi \xi }}} \sqrt {{G_{\eta \eta }}} }}\dfrac{{\partial \sqrt {{G_{\xi \xi }}} }}{{\partial \eta }} = \\ {\text{−}} \dfrac{1}{{{\rho _0}\sqrt {{G_{\xi \xi }}} }}{P_\xi } {\mathrm{+}} {F_\xi } {\mathrm{+}} \dfrac{1}{{{{\left( {d {\mathrm{+}} \zeta } \right)}^2}}}\dfrac{\partial }{{\partial \sigma }}\left( {{\nu _{\text{V}}}\dfrac{{\partial u}}{{\partial \sigma }}} \right) {\mathrm{+}} {M_\xi } {\mathrm{+}} fv \\ \end{array} $$ (7) $$\begin{array}{c} \dfrac{{\partial v}}{{\partial t}} {\mathrm{+}} \dfrac{u}{{\sqrt {{G_{\xi \xi }}} }}\dfrac{{\partial v}}{{\partial \xi }} {\mathrm{+}} \dfrac{v}{{\sqrt {{G_{\eta \eta }}} }}\dfrac{{\partial v}}{{\partial \eta }}{\text{ {\mathrm{+}} }}\dfrac{w}{{d {\mathrm{+}} \zeta }}\dfrac{{\partial v}}{{\partial \sigma }} {\mathrm{+}} \\ \dfrac{{uv}}{{\sqrt {{G_{\xi \xi }}} \sqrt {{G_{\eta \eta }}} }}\dfrac{{\partial \sqrt {{G_{\eta \eta }}} }}{{\partial \xi }} {\text{−}} \dfrac{{uu}}{{\sqrt {{G_{\xi \xi }}} \sqrt {{G_{\eta \eta }}} }}\dfrac{{\partial \sqrt {{G_{\xi \xi }}} }}{{\partial \eta }} = \\ {\text{−}} \dfrac{1}{{\rho _0}\sqrt {G_{\eta \eta }}}{P_{\eta} } {\mathrm{+}} {F_{\eta} } {\mathrm{+}} \dfrac{1}{{{{\left( {d {\mathrm{+}} \zeta } \right)}^2}}}\dfrac{\partial }{{\partial \sigma }}\left( {{\nu _{{{\mathrm{V}}}}}\dfrac{{\partial v}}{{\partial \sigma }}} \right) {\mathrm{+}} \\ {M_{\eta} } {\text{−}} fu \end{array}$$ (8) $$ \begin{array}{c} \dfrac{1}{{\sqrt {{G_{\xi \xi }}} \sqrt {{G_{\eta \eta }}} }}\left( {\dfrac{{\partial \left[ {\left( {d {\mathrm{+}} \zeta } \right)u\sqrt {{G_{\eta \eta }}} } \right]}}{{\partial \zeta }}} \right.\left. { {\mathrm{+}} \dfrac{{\partial \left[ {\left( {d {\mathrm{+}} \zeta } \right)v\sqrt {{G_{\xi \xi }}} } \right]}}{{\partial \eta }}} \right) {\mathrm{=}}\\ \left( {d {\mathrm{+}} \zeta } \right)\left( {{q_{{\text{in}}}} - {q_{{\text{out}}}}} \right) - \dfrac{{\partial w}}{{\partial \sigma }} - \dfrac{{\partial \zeta }}{{\partial t}} \end{array} $$ (9) 式中：v、u、w分别表示η、ξ、σ方向上的流速；Pξ
、Fξ分别表示ξ方向上的静压梯度、雷诺应力梯度；Fη 、Pη分别表示η方向上的雷诺应力梯度、静压梯度；Mξ、Mη表示外部动量的源或汇；f表示科氏力系数；vv表示垂向涡流黏度系数；qin和qout分别表示单位体积的源项和汇项。 水动力模型的潮位边界采用Tide Model Driver软件获取M2、S2、K1、O1、N2、P1、K2、Q1共8个主要潮汐分潮的调和常数进行驱动。
1.3 波浪模型
Delft3D-WAVE模块采用第三代波浪模型SWAN[35]进行波浪过程模拟，它考虑了近岸波浪传播的物理过程，包括波浪浅化、绕射、折射和波浪破碎等。其动谱平衡方程为：
$$ \dfrac{\partial }{{\partial t}}N {\mathrm{+}} \dfrac{\partial }{{\partial x}}{c_x}N {\mathrm{+}} \dfrac{\partial }{{\partial y}}{c_y}N {\mathrm{+}} \dfrac{\partial }{{\partial \sigma }}{c_\sigma }N {\mathrm{+}} \dfrac{\partial }{{\partial \theta }}{c_\theta }N = \dfrac{S}{\sigma } $$ (10) 式中：θ表示波浪传播方向；σ表示波浪相对频率；cx、cy、cθ和cσ分别表示在x、y、θ和σ空间中的传播速度；S表示能量作用量的源汇项，包括风能输入、波浪破碎、白浪耗散、底部摩擦引起的能量损耗等。
1.4 耦合模型及研究区域网格设置
马袅湾海域及部分琼州海峡区域的地形数据取自航保部电子海图C151586和C1515819，其余地形采用全球水深数据ETOPO1进行补充，海南省马袅湾地理位置、网箱分布如图1所示。
![]()
图 1 海南省马袅湾网箱位置示意图Figure 1. Geographic location, distribution of cages in Maniao Bay, Hainan Province台风过程的模拟需要考虑波流相互作用，尤其是在浅水中[36]。本研究采用Delft3D FLOW &WAVE模型耦合，波对流的作用主要通过在水动力模型中加入由波浪作用所产生的辐射应力[37]和波浪破碎[38]产生的动量源项，流对波的作用则体现为流对波浪的折射及频移等效应[39-40]。
水动力模型和波浪模型共用一套结构化网格，台风过程中影响海域的范围很大，使用二重嵌套网格，大模型网格计算范围为105.5°E—118°E、15°N—24°N，覆盖了大部分南海区域，网格分辨率设置为4 000 m；小模型网格计算范围为109.6°E—110.8°E、19.9°N—21.4°N，主要区域为琼州海峡附近，并在马袅湾海域进行局部加密，最小网格分辨率约为100 m，如图2所示。
1.5 模型验证
选取秀英站和湛江站天文潮的模拟结果与国家海洋局出版的潮汐表数据进行对比，结果如图3所示。模拟结果与实测数据之间的平均绝对误差 (MAE)、均方根误差 (RMSE) 和相关系数 (CC)，结果如表1所示。可以看出，天文潮的模拟结果具有较好的一致性，平均绝对误差小于0.14 m，均方根误差不超过0.17 m，相关系数均在0.97以上。
表 1 天文潮模拟精度Table 1. Simulation accuracy of astronomical tide观测站
Observation station平均绝对误差
MAE/m均方根误差
RMSE/m相关系数
CC秀英 Xiuying 0.134 0.166 0.978 湛江 Zhanjiang 0.109 0.145 0.990 为了验证风暴潮增水的模拟效果，将秀英站的计算结果与观测值[41]进行对比验证。结果见图4。其MAE为0.139 m，RMSE为0.180 m，相关系数为0.945。
选择“威马逊”、“海鸥”和“启德”台风期间的秀英、东方站和浮标测点Buoy2 [42]实测浪高对模型进行验证，结果如图5所示。模拟和实测误差的RMSE、MAE和平均绝对百分比误差 (MAPE)见表2。
![]()
图 5 “威马逊”“启德”和“海鸥”台风有效波高结果比较Figure 5. Verification of significant wave height results for typhoons Rammasun, Kaitai and Kalmaegi表 2 台风浪模拟精度Table 2. Significant wave height accuracy观测站
Observation station平均绝对
误差
MAE/m均方根
误差
RMSE/m平均绝对
百分比误差
MAPE/%相关系数
Correlation
coefficient秀英 Xiuying 0.179 0.255 7.4 0.992 浮标测点Buoy2 0.425 0.520 18.3 0.980 东方 Dongfang 0.238 0.335 14.6 0.934 模型在天文潮、风暴潮增水以及台风浪的模拟中，模拟结果均与实测值吻合较好，模型具有较高的准确性。
2. 不同台风对网箱养殖区的影响
2.1 台风路径分类
依据马袅湾的地理位置，参考先前的研究[43]，将影响该区域的台风路径划分为4种类型 (图6)：第I类台风路径从海南东北部或雷州半岛斜穿后向西北移动，进入北部湾；第II类台风路径自东向西穿过海南岛后进入北部湾；第III类台风路径沿海南岛东侧向北移动，最终在广东省沿岸登陆；第IV类台风路径从海南岛南部海域绕岛向北发展，进入北部湾。统计1995—2022年间对马袅湾产生显著影响的52场台风，各类台风路径的发生频次及强度如图7所示。结果表明，第I和第II类台风发生频次较高，均为18场 (34.6%)；第III类台风次之，11场 (21.2%)；第IV类台风最少，仅5场 (9.6%)。在52场台风中，“强台风”和“超强台风”共10场，其中第I、第II类台风各4场；第III类2场；第IV类路径台风等级较低。由此可见，第I、第II类台风对马袅湾海域构成的潜在风险更大。
![]()
图 7 不同台风路径下各等级台风发生次数Figure 7. Number of various levels of typhoons under different tracks2.2 第I类台风对网箱养殖区的影响
本研究选取2014年第
1409 号台风“威马逊”(Rammasun) 作为第I类路径的典型案例进行分析 (图8)。该台风于2014年7月10日在西北太平洋海域生成，7月15日首次在菲律宾登陆，随后于7月18日以中心气压910 hpa、最大风速60 m·s−1的强度在中国海南省文昌市二次登陆。台风穿越雷州半岛进入北部湾后，强度减弱至中心气压940 hpa、最大风速50 m·s−1。对“威马逊”台风过境期间的波流过程进行了模拟，获得了各网格点的有效波高和流场变化过程。进一步统计台风过程中各网格点有效波高的最大值，得到了研究区域有效波高最大值分布及养殖区内流速最大时刻流场 (图8)。养殖区内有效波高范围为3.84~5.94 m，流速范围为0.62 ~1.42 m·s−1。![]()
图 8 “威马逊”台风下养殖区有效波高最大值分布及最大流速时刻流场注：a. 养殖区有效波高最大值分布；b. 养殖区内最大流速时刻流场。Figure 8. Distribution of maximum significant wave height and maximum flow velocity moment flow field in aquaculture area under typhoon RammasunNote: a. Distribution of maximum significant wave height in the aquaculture area; b. Maximum flow velocity moment flow field in the aquaculture area.以“威马逊”台风为模板，将原台风路径沿经度方向南北平移，平移间距为0.3°，构造出R1—R8共8场台风 (图9)。以构造的路径为基础，将每条路径下台风的强度等级按GB/T 19201—2006《热带气旋等级》依次设置为“超强台风”、“强台风”和“台风”，不同强度台风对应的参数统一取值如表3所示，共构造出27场设计台风。图10为R1—R8路径不同强度作用下马袅湾养殖区内的有效波高及流速最大值。
表 3 台风等级相关参数信息Table 3. Information on parameters related to typhoon levels台风等级
Typhoon intensity
level最大风速
Maximum wind
speed/(m·s−1)中心气压
Central
pressure/hpa超强台风 Super Typhoon 51.0 944 强台风 Severe Typhoon 41.5 962 台风 Typhoon 32.7 977 ![]()
图 10 马袅湾养殖区内有效波高及流速最大值Figure 10. Maximum significant wave height and current velocity in aquaculture area of Maniao Bay总体而言，随着台风路径与养殖区距离的减小，各强度台风在养殖区内引发的HS_max与Vmax均呈递增趋势。R4路径台风导致的HS_max与Vmax均为最大，对养殖区的威胁最显著。这是由于R4台风路径位于马袅湾东北方向，台风中心距离较近，其逆时针旋转风场驱动波浪自湾口向湾内传播所致。在相同路径条件下，台风强度是影响网箱养殖区的关键因素，当台风强度减弱时，养殖区HS_max与Vmax均相应降低，威胁程度减小。在“台风”、“强台风”和“超强台风”等级下，路径每平移0.3°，HS_max平均变化幅度为12.77%~15.57%，Vmax平均变化幅度为12.33%~14.45%。
当台风等级为“超强台风”时，R3—R8路径下养殖区内Vmax≥1 m·s−1，为严重影响区域。其中，R4路径台风引发的HS_max>5.5 m且Vmax>1.5 m·s−1，是影响程度最严重的路径。对于“强台风” 等级，R3—R5路径下Vmax≥1 m·s−1，养殖区受严重影响；R6路径下网箱养殖区内Vmax<1 m·s−1，但由于4 m≤HS_max<5.5 m，养殖区受中等程度影响；而在R3以北和R6以南路径下，养殖区内Vmax<1 m·s−1且2.5 m≤HS_max<4 m，网箱养殖区受轻微影响。在“台风”等级情况下，“威马逊”至R6路径之间的范围均会导致养殖区内Vmax≥1 m·s−1，网箱养殖区受严重影响；而“威马逊”以北和R6以南路径下，Vmax<1 m·s−1且2.5 m≤HS_max<4 m，养殖区仅受轻微影响。
2.3 第II类台风对网箱养殖区的影响
选取2021年第
2118 号台风“圆规”(Kompasu) 作为第II类台风路径的典型代表。图11为该台风下养殖区有效波高最大值分布及最大流速时刻流场，马袅湾内部有效波高较小，均在4 m以下，特别是在网箱养殖区，波高介于2.5~3.3 m，流速介于0.33 ~0.65 m·s−1。![]()
图 11 “圆规”台风下养殖区有效波高最大值分布及最大流速时刻流场注：a. 养殖区有效波高最大值分布；b. 养殖区内最大流速时刻流场。Figure 11. Distribution of maximum significant wave height and maximum flow velocity moment flow field in aquaculture area under typhoon KompasuNote: a. Distribution of maximum significant wave height in the aquaculture area; b. Maximum flow velocity moment flow field in the aquaculture area.将次平移的间距设定为0.3°，构建了K1、K2、K3、K4和K5共5条设计台风路径，如图12所示。将每条路径的台风强度分别设定为“超强台风”、“强台风”和“台风”，共构建出18场设计台风。图13为不同强度和路径的台风作用下网箱养殖区内有效波高及流速最大值。整体上看，随着台风路径从南向北移动，HS_max与Vmax均随路径的北移呈先增后减的趋势。“台风”、“强台风”、“超强台风”等级下，每平移0.3°，HS_max平均变化幅度为13.91%~16.90%，Vmax平均变化幅度为5.57%~14.55%。
![]()
图 13 马袅湾养殖区内有效波高及流速最大值Figure 13. Maximum significant wave height and current velocity in aquaculture area of Maniao Ba当台风等级为“超强台风”时， K2和K3路径下的网箱养殖区Vmax≥1 m·s−1，属于严重影响区域；K1和K4路径下的养殖区Vmax<1 m·s−1且4 m≤HS_max<5.5 m，属于中等影响区域；而“圆规”和K5路径下的Vmax<1 m·s−1且2.5 m≤HS_max<4 m，属于轻微影响区域。当台风等级为“强台风”时，各路径下的养殖区Vmax均小于1 m·s−1，其中K2—K4路径下4 m≤HS_max<5.5 m，网箱养殖区会受中等程度影响，“圆规”、K1和K5路径下2.5 m≤HS_max<4 m，网箱养殖区受轻微影响。台风等级为“台风”时，各路径下的养殖区Vmax均小于1 m·s−1，K1—K5路径下2.5 m≤HS_max<4 m，网箱养殖区受轻微影响， “圆规”路径下HS_max<2.5 m，网箱养殖区不受影响。
2.4 第III类台风对网箱养殖区的影响
选取2002年第0214号强台风“黄蜂”(VONGFONG) 作为第III类路径的典型代表。“黄蜂”过境过程中在马袅湾引起的有效波高最大值分布及在养殖区内引起最大流速时刻的马袅湾流场如图14所示。可以看出，在网箱养殖区内，波高主要集中在3.0~3.3 m之间，流速介于0.33~0.71 m·s−1。
![]()
图 14 “黄蜂”台风下养殖区有效波高最大值分布及最大流速时刻流场注：a. 养殖区有效波高最大值分布；b. 养殖区内最大流速时刻流场。Figure 14. Distribution of maximum significant wave height and maximum flow velocity moment flow field in aquaculture area under typhoon VongfongNote: a. Distribution of maximum significant wave height in the aquaculture area; b. Maximum flow velocity moment flow field in the aquaculture area.依据路径类型定义范围，设置平移间距为0.6°，将原台风路径南向平移构建出V1、V2和V3共3条新的设计台风路径 (图15)。在此基础上，将等级分别设置为“超强台风”、“强台风”和“台风”，共构建出12场不同的设计台风。图16为不同强度和路径的台风作用下网箱养殖区内HS_max和Vmax。总体上看，随着台风路径向南移动，HS_max和Vmax逐渐增大。“台风”、“强台风”和“超强台风”等级下，每平移0.6°，HS_max平均变化幅度为9.17%~11.76%，Vmax平均变化幅度为4.01%~7.57%。
![]()
图 16 马袅湾养殖区内有效波高及流速最大值Figure 16. Maximum significant wave height and currentvelocity in aquaculture area of Maniao Bay当台风等级为“超强台风”时，V2和V3路径下的养殖区Vmax≥1 m·s−1，属于严重影响区域；V1路径下的养殖区Vmax<1 m·s−1且4 m≤HS_max<5.5 m，属于中等影响区域；而“黄蜂”路径下的养殖区2.5 m≤ HS_max<4 m且Vmax<1 m·s−1，属于轻微影响区域。当台风等级为“强台风”时，V3路径下的养殖区Vmax≥1 m·s−1，受严重影响；在“黄蜂”至V2路径范围内，养殖区Vmax均小于1 m·s−1，其中V2路径下4 m≤HS_max<5.5 m，网箱养殖区受中等程度影响；而“黄蜂”和V1路径下2.5 m≤HS_max<4 m，网箱养殖区受轻微影响。在台风等级为“台风”时，各路径下的养殖区Vmax均小于1 m·s−1。其中，V2和V3路径下网箱养殖区4 m≤HS_max<5.5 m，受中等程度影响；而“黄蜂”和V1路径下2.5 m≤HS_max<4 m，均受轻微影响。
2.5 第IV类台风对网箱养殖区的影响
选取1995年第
9516 号台风“泰德”(Ted) 作为第IV类台风路径的典型代表案例 。图17为“泰德”台风下养殖区的有效波高最大值分布及最大流速时刻的流场。可以看出，该台风在养殖区引发的最大波高为0.8~0.9 m，流速范围为0.17~0.55 m·s−1。![]()
图 17 “泰德”台风下养殖区有效波高最大值分布及最大流速时刻流场注：a. 养殖区有效波高最大值分布；b. 养殖区内最大流速时刻流场。Figure 17. Distribution of maximum significant wave height and maximum flow velocity moment flow field in aquaculture area under typhoon TedNote: a. Distribution of maximum significant wave height in the aquaculture area; b. Maximum flow velocity moment flow field in the aquaculture area.结合路径类型的定义范围，将原路径向东以0.3° 为间距平移，构造出T1、T2、T3和T4共4条设计台风路径 (图18)。将每条路径的台风强度分别设定为“超强台风”、“强台风”和“台风”，共构建出15场设计台风。图19为不同强度与路径台风作用下网箱养殖区的HS_max和Vmax。总体上看，随着台风路径向东移动，HS_max和Vmax逐渐增大。在“台风”、“强台风”和“超强台风”等级下，每东移0.3°，HS_max平均变化幅度为22.32%~23.88%，Vmax平均变化幅度为13.11%~19.02%。“超强台风”等级下，距离养殖区最近的T4路径引发的HS_max最大 (2.49 m)。仅在“超强台风”和“强台风”等级的T4路径下，养殖区出现Vmax>1 m·s−1。本研究表明，第IV类台风对网箱养殖区安全构成的威胁相对较小。
![]()
图 19 马袅湾养殖区内有效波高及流速最大值Figure 19. Maximum significant wave height and currentvelocity in aquaculture area of Maniao Bay3. 结论
本研究基于Holland风场模型和Delft3D软件建立了适用于海南省马袅湾的台风波流耦合模型，深入探究了台风路径及其强度对网箱养殖区的影响，得到以下结论：
1) 第I类台风路径发生频率高，覆盖范围广，影响程度大。“超强台风”等级时，路径每平移0.3°，HS_max平均变化15.57%，Vmax平均变化13.15%。R3—R8路径均会使马袅湾养殖区受到严重影响，其中R4路径下，养殖区HS_max>5.5 m，Vmax>1.5 m·s−1，遇到类似台风时应立即预警并做好周全保护措施，而R3以北路径台风对网箱养殖区影响轻微；当等级为“强台风”时，HS_max平均变化14.45%，Vmax平均变化14.84%，在R3—R5路径下，网箱养殖区受严重影响，应立即预警，R6路径下，养殖区受中等程度影响，R3以北和R6以南路径的台风，对网箱养殖区的影响轻微；当等级降至“台风”时，HS_max平均变化12.77%，Vmax平均变化12.33%，“威马逊”至R6路径的范围会使网箱养殖区受严重影响，均应立即预警，而“威马逊”以北和R6以南路径台风的影响程度轻微，可适度预警。
2) 第II类台风作用下，“台风”、“强台风”、“超强台风”等级下，每平移0.3°，HS_max平均变化幅度为13.91%~16.90%，Vmax平均变化幅度为5.57%~14.55%。“超强台风”等级时K2—K3路径下，网箱养殖区受严重影响；“超强台风”等级时K1与K4路径、“强台风”等级时K2—K4路径，网箱养殖区受中等程度影响，应快速预警；其余情况下，网箱养殖区受轻微及以下影响，可适度预警并采取设施检查等措施。
3) 第III类台风作用下，“台风”、“强台风”和“超强台风”等级下，每平移0.6°，HS_max平均变化幅度为9.17%~11.76%，Vmax平均变化幅度为4.01%~7.57%。等级为“超强台风”时的V2与V3路径下、“强台风”等级的V3路径下，网箱养殖区受严重影响；等级为“超强台风”时的V1路径、“强台风”和“台风”等级的V2路径、“台风”等级的V3路径下，网箱养殖区受中等程度影响，应快速预警，可采取设施加固等措施；其余情况下，网箱养殖区受轻微影响，应适度预警。
4) 第IV类台风作用下，“台风”、“强台风”和“超强台风”等级下，每平移0.3°，HS_max平均变化幅度为22.32%~23.88%，Vmax平均变化幅度为13.11%~19.02%。第IV类台风中心距马袅湾较远，仅在“超强台风”和“强台风”等级的T4路径下，养殖区Vmax>1 m·s−1，对养殖区影响严重。其余情况下，养殖区处于安全状态，无需预警。
本研究采用有效波高最大值和流速最大值作为预警依据，所得结果可为后续建立马袅湾网箱养殖区的台风防灾减灾快速预警系统提供依据。然而，网箱破坏是一个复杂的过程，后续将依据马袅湾不同的波流工况对网箱承载能力进行研究，以获得更科学的预警判据，提升预警的精度。
-
![]()
图 1 各组血清生化指标检测
注：不同字母表示差异显著 (p<0.05)。
Figure 1. Serum biochemical indexes detection in each group
Note: Different letters represent significant differences (p<0.05).
![]()
图 2 各组骨小梁指标检测
注：不同字母表示差异显著 (p<0.05)。
Figure 2. Bone trabecular index test for each group
Note: Different letters represent significant differences (p<0.05).
![]()
图 3 各组股骨干骺端Micro-CT扫描图像
注：箭头所指为骨髓腔中骨小梁。a. 空白组；b. 模型组；c. 阳性组；d. 牡蛎酶解产物低剂量组；e. 牡蛎酶解产物中剂量组；f. 牡蛎酶解产物高剂量组。
Figure 3. Micro-CT scans of the femoral epiphysis of each group
Note: Arrows point to trabeculae in the marrow cavity. a. CK group; b. DEX group; c. AL group; d. OP-L group; e. OP-M group; f. OP-H group.
![]()
图 4 各组Micro-CT 3D重建图像
注：上为3D重建正视图，下为3D重建俯视图。a. 空白组；b. 模型组；c. 阳性组；d. 牡蛎酶解产物低剂量组；e. 牡蛎酶解产物中剂量组；f. 牡蛎酶解产物高剂量组。
Figure 4. Micro-CT 3D reconstructed images of each group
Note: The top images are the 3D reconstruction front view, and the bottom images are the 3D reconstruction top view. a. CK group; b. DEX group; c. AL group; d. OP-L group; e. OP-M group; f. OP-H group.
![]()
图 5 各组股骨HE染色结果 (200×)
注：箭头所指为骨小梁。a. 空白组；b. 模型组；c. 阳性组；d. 牡蛎酶解产物低剂量组；e. 牡蛎酶解产物中剂量组；f. 牡蛎酶解产物高剂量组。
Figure 5. HE staining results of femur in each group (200×)
Note: Arrows point to trabeculae. a. CK group; b. DEX group; c. AL group; d. OP-L group; e. OP-M group; f. OP-H group.
![]()
图 6 各组股骨钙结节指染色结果
注：箭头所指为钙盐沉积部分。a. 空白组；b. 模型组；c. 阳性组；d. 牡蛎酶解产物低剂量组；e. 牡蛎酶解产物中剂量组；f. 牡蛎酶解产物高剂量组。
Figure 6. Finger staining results of femoral calcaneal nodes in each group
Note: Arrows point to the portion of calcium salt deposits. a. CK group; b. DEX group; c. AL group; d. OP-L group; e. OP-M group; f. OP-H group.
表 1 动物模型造模及给药
Table 1 Animal model modeling and drug administration
组别
Group肌肉注射
Intramuscular injection灌胃
Gastric lavage对照组 CK 0.9% (w) 氯化钠注射，2.5 mg·kg−1，每周2次 蒸馏水，每天1次 模型组 DEX 地塞米松磷酸钠注射液，2.5 mg·kg−1，每周2次 蒸馏水，每天1次 阳性药物组 AL 地塞米松磷酸钠注射液，2.5 mg·kg−1，每周2次 阿仑膦酸钠片，6.3 mg·kg−1，每周1次 牡蛎酶解产物给药低剂量组 OP-L 地塞米松磷酸钠注射液，2.5 mg·kg−1，每周2次 牡蛎酶解产物，300 mg·kg−1，每天1次 牡蛎酶解产物给药中剂量组 OP-M 地塞米松磷酸钠注射液，2.5 mg·kg−1，每周2次 牡蛎酶解产物，900 mg·kg−1，每天1次 牡蛎酶解产物给药高剂量组 OP-H 地塞米松磷酸钠注射液，2.5 mg·kg−1，每周2次 牡蛎酶解产物，2 700 mg·kg−1，每天1次 
下载: 导出CSV
-
[1] JULIET C. Glucocorticoid-induced osteoporosis: update[J]. Endocrine, 2018, 61(1): 7-16. doi: 10.1007/s12020-018-1588-2
[2] JHA S S. Correction: glucocorticoid-induced osteoporosis (GIOP)[J]. Ind J Orthop, 2024, 58(2): 222. doi: 10.1007/s43465-023-01078-z
[3] 中华医学会骨质疏松和骨矿盐疾病分会. 中国骨质疏松症流行病学调查及“健康骨骼”专项行动结果发布[J]. 中华骨质疏松和骨矿盐疾病杂志, 2019, 12(4): 317-318. doi: 10.3969/j.issn.1674-2591.2019.04.001 [4] 中华医学会物理医学与康复学分会. 骨质疏松症康复治疗指南 (2024版)[J]. 中国循证医学杂志, 2024, 24(6): 626-636. [5] 王佳佳, 张静琳. 糖皮质激素性高眼压小梁网病理改变及相关分子机制进展[J]. 国际眼科杂志, 2024, 24(7): 1088-1092. doi: 10.3980/j.issn.1672-5123.2024.7.15 [6] ZHANG L, LI J H, XIE R S, et al. Osteoporosis guidelines on TCM drug therapies: a systematic quality evaluation and content analysis[J]. Front Endocrinol, 2024, 14: 1276631. doi: 10.3389/fendo.2023.1276631
[7] ZHUO Y, LI M, JIANG Q Y, et al. Evolving roles of natural terpenoids from traditional chinese medicine in the treatment of osteoporosis[J]. Front Endocrinol, 2022, 13: 901545. doi: 10.3389/fendo.2022.901545
[8] DU J Z, WANG Y C, WU C L, et al. Targeting bone homeostasis regulation: potential of traditional Chinese medicine flavonoids in the treatment of osteoporosis[J]. Front Pharmacol, 2024, 15: 1361864. doi: 10.3389/fphar.2024.1361864
[9] GU Z X, ZHOU G H, ZHANG X Q, et al. Research progress of plant medicine and Chinese herbal compounds in the treatment of rheumatoid arthritis combined with osteoporosis[J]. Front Med, 2024, 10: 1288591. doi: 10.3389/fmed.2023.1288591
[10] 胡子聪, 杨可君, 房翠兰. 南极磷虾油的营养特性及其改善骨质疏松症的研究进展[J]. 中国油脂, 2024, 49(4): 32-39. [11] 梅刚. 鲑鱼降钙素注射液对合并骨关节炎绝经后骨质疏松症患者疼痛、骨密度及膝关节功能的影响分析[J]. 中国处方药, 2021, 19(2): 93-94. doi: 10.3969/j.issn.1671-945X.2021.02.048 [12] 马婷, 吴谦, 申铉日. 珍珠贝外套膜胶原蛋白肽及其锌螯合物的体外抑制骨质疏松作用[J]. 肉类研究, 2018, 32(2): 36-42, 9. [13] 逯冠宏. 牡蛎多肽提取工艺的建立及其对RAW264.7巨噬细胞免疫调节作用的研究[D]. 长春: 吉林大学, 2021: 47-60. [14] 钟凯翠, 邹雯清, 陈忠琴, 等. 海红米花色苷-牡蛎肽新型复合饮料制备工艺及降血糖作用研究[J/OL]. 食品与发酵工业, 1-13. https://doi.org/10.13995/j.cnki.11-1802/ts.039921. [15] 蔡雯雯, 葛小东, 李娜, 等. 牡蛎蛋白酶解物对小鼠慢性酒精性肝损伤的保护作用[J]. 食品科学, 2022, 43(17): 147-155. doi: 10.7506/spkx1002-6630-20210720-228 [16] 谭强来, 陆馨敏, 吴兰兰, 等. 牡蛎多糖缓解体力疲劳及其对力竭小鼠肠道菌群的影响[J]. 中国微生态学杂志, 2024, 36(4): 390-400. [17] 黄倩倩, 吴海英, 秦小明. 牡蛎肽对帕罗西汀致雄性小鼠性功能障碍的作用效果及潜在机制[J]. 南方水产科学, 2024, 20(1): 161-172. doi: 10.12131/20230153 [18] HU Y Y, HOU Z X, LIU Z Q, et al. Oyster mantle-derived exosomes alleviate osteoporosis by regulating bone homeostasis[J]. Biomaterials, 2024, 311: 122648.
[19] IHN H J, KIM J A, LIM S, et al. Fermented oyster extract prevents ovariectomy-induced bone loss and suppresses Osteoclastogenesis[J]. Nutrients, 2019, 11(6): 1392. doi: 10.3390/nu11061392
[20] MOLAGODA I M N, ATHAPATHTHU A M G K, PARK E K, et al. Fermented oyster (Crassostrea gigas) extract cures and prevents prednisolone-induced bone resorption by activating osteoblast differentiation[J]. Foods, 2022, 11(5): 678. doi: 10.3390/foods11050678
[21] CHEN H, SHI P J, XU Z, et al. Oral administration of oyster peptide prevents bone loss in ovariectomized mice[J]. eFood, 2020, 1(4): 298-309. doi: 10.2991/efood.k.200812.001
[22] JEONG J W, CHOI S H, HAN M H, et al. Protective effects of fermented oyster extract against RANKL-induced Osteoclastogenesis through scavenging ROS generation in RAW 264.7 cells[J]. Int J Mol Sci, 2019, 20(6): 1439. doi: 10.3390/ijms20061439
[23] ASHA K K, KUMARI K R R, KUMAR K A, et al. Sequence determination of an antioxidant peptide obtained by enzymatic hydrolysis of oyster Crassostrea madrasensis (Preston)[J]. Int J Pept Res Ther, 2016, 22: 421-433. doi: 10.1007/s10989-016-9521-0
[24] 程亚军, 陈绍丰, 李明, 等. 骨转换标志物在绝经后骨质疏松症中的临床应用[J]. 海军军医大学学报, 2022, 43(2): 182-187. [25] 《中国骨质疏松杂志》骨代谢专家组. 骨代谢生化指标实验规范 (2024)[J]. 中国骨质疏松杂志, 2024, 30(8): 1093-1101. doi: 10.3969/j.issn.1006-7108.2024.08.001 [26] 冯佳, 向阳, 夏燕, 等. 三种大鼠骨质疏松模型的比较研究[J]. 实验动物与比较医学, 2015, 35(1): 52-55. doi: 10.3969/j.issn.1674-5817.2015.01.011 [27] da COSTA N D S D, LIMA S L, GALICIOLLI M E A, et al. Drug-induced osteoporosis and mechanisms of bone tissue regeneration through trace elements[J]. J Trace Elem Med Biol, 2024, 84: 127446. doi: 10.1016/j.jtemb.2024.127446
[28] BACHMANN S, DUNMORE C J, SKINNER M M, et al. A computational framework for canonical holistic morphometric analysis of trabecular bone[J]. Sci Rep, 2022, 12(1): 5187. doi: 10.1038/s41598-022-09063-6
[29] MÜLLER R. Hierarchicl microimaging of bone structure and function[J]. Nat Rev Rheumatol, 2009, 5(7): 373-381. doi: 10.1038/nrrheum.2009.107
[30] 苏开鑫, 谢华, 王宏芬, 等. 牡蛎肉提取物对类固醇性骨质疏松大鼠骨代谢的影响[J]. 中国自然医学杂志, 2009, 11(2): 97-99. [31] WAWRZYNIAK A, BALAWENDER K. Structural and metabolic changes in bone[J]. Animals, 2022, 12(15): 1946. doi: 10.3390/ani12151946
[32] ZIMMERMANN E A, RIEDEL C, SCHMIDT F N, et al. Mechanical competence and bone quality develop during skeletal growth[J]. J Bone Miner Res, 2019, 34(8): 1461-1472.
[33] CIOSEK Ż, KOT K, KOSIK-BOGACKA D, et al. The effects of calcium, magnesium, phosphorus, fluoride, and lead on bone tissue[J]. Biomolecules, 2021, 11(4): 506.
计量
- 文章访问数: 0
- HTML全文浏览量: 0
- PDF下载量: 0
粤公网安备 44010502001741号
