Page 162 - 《水产学报》2026年第01期

1 期 水 产 学 报 50 卷

表 1 DEB 模型中所用的主要关系式
Tab. 1 Main function in the DEB model
序号 方程 定义
number equation definition
( ) [ ( ) ( )] −1
T A T A T AL T AL T AH T AH
1 k(T) = k 1 exp − · 1+exp − +exp − 温度依赖关系 temperature dependence
T 1 T T T L T H T
 { } 2/3 
[E]     [E G ]· ˙ P Am V [ ]    
2 ˙ P = k(T) 代谢率 catabolic rate
[E G ]+κ ·[E] ·    [E M ] + ˙ P M ·V  
E
3 [E] = 能量密度 energy density
V
{
}
4 ˙ P A = k(T)· f · ˙ P Am ·V 2/3 吸收率 assimilation rate
F
5 f = 功能性反应 functional response
F + F H
[ ]
6 ˙ P M = k(T)· ˙ P M ·V 维持率 maintenance rate
( )
[ ] 1−κ
7 ˙ P J = k(T)·min(V,V P )· ˙ P M · 繁育维持率 maturity maintenance rate
κ
dE
8 = ˙ P A − ˙ P C 储能动力学 reserve dynamic
dt
dE R
9 = (1−κ)· ˙ P C − ˙ P J 繁育储能变化率 reproductive reserve dynamic
dt
dV κ · ˙ P C − ˙ P M
10 = 生物体积生长 bio-volume growth
dt [E G ]
E k R + E R
11 DW = + +V ·ρ 组织干重 dry tissue weight
µ E µ E
E 为分配到繁育中的能量 (J)。 2 结果
R
模型运行与验证 使用 R 软件进行模型模
拟。状态变量的初始值由第 1 次采样的测量值估 2.1 模型参数获取
算，繁育储能 E 的初始值设为 0。水温 (℃) 和
R
δ 对咬齿牡蛎软体部湿重的立方根和壳
m
Chl.a 浓度 (μg/L) 是模型运行的强制函数。大亚湾
长进行线性回归，得到的斜率即为 δ ，从图 1 中
m
咬齿牡蛎的生长模拟期为 2022 年 6 月—2024 年
可知咬齿牡蛎的 δ 值为 0.270。
m
4 月，共运行 700 d。以每 6 个月获得的咬齿牡蛎
T 在 15~30 ℃，咬齿牡蛎的耗氧率随水
生长数据 (壳长和软体部干重) 作为模型的验证。 A
温升高而增加，在 30 ℃ 达到峰值后，随水温升
咬齿牡蛎生长限制性分析 根据表 1 中
高而降低 (图 2)，其耗氧率与水温的热力学温度倒
公式 5，在食物半饱和常数一定的情况下，食物
数的线性回归关系 (图 3)：ln R=−4 900.10T + 15.73
−1
的功能性反应 f 值依据环境中食物浓度而变化，
2
(R =0.79)，根据公式获得咬齿牡蛎 T 为 4 900.10 K。
A
可反映环境中食物对牡蛎能量获得及生长的影响
[E ]、[E ]、[Ṗ ] 持续饥饿后，咬齿牡蛎
M
M
G
情况。f 值介于 0~1，值越小，说明对牡蛎获取能
软体部干重最终维持在 (0.43±0.06) g。初始软体部
量的限制性越大。表 1 中公式 1 为温度依赖关系
有机物含量为 89.80%±3.67%，实验结束时，软体
k(T)，可反映随环境水温对贝类生理反应速率的
部有机物含量下降至 65.46%±2.80% (图 4)。饥饿
影响情况。k(T) 的值>0，值越小，说明水温对牡
实验进行20 d 时，耗氧率趋于稳定，咬齿牡蛎实
蛎生长的限制性越大，当外界水温为最适温度时，
验结束时稳定在 0.18 mg/(g DW·h)，相较于初始耗
k(T) 为最大值。
氧率下降了 50.00% (图 5)。根据能量转换 (1 mL
1.3 数据分析
O =20.3 J) 及 δ ，获得咬齿牡蛎维持生命所需能量
m
2
3
实验数据采用 Excel 2010 软件进行整理、线 [Ṗ ] 为 26 J/(cm ·d)。根据有机物能值、δ 和反应
m
M
性回归拟合并作图，结果以平均值±标准差表示。 生长效率系数 40%，获得单位体积结构物质所需
https://www.china-fishery.cn 中国水产学会主办 sponsored by China Society of Fisheries
4