DreamBlooms/HansPub-15000-to-45000 huggingface dataset:
system: [SYSTEM]

user: 2. 临床资料 2.1. 一般资料 临床资料 对照组共98例，男52例，女46例；年龄58~89岁，平均67.8岁；急性发病2~6天；治疗组共100例，男53例，女47例；年龄55~92岁，平均68.8岁；急性发病2~7天；二组资料经检验，P > 0.05，差异无统计学意义。 2.2. 纳入标准 ① 术前明确诊断：胆囊结石伴慢性胆囊炎急性发作或急性非结石性胆囊炎；② 术中见胆囊结石嵌顿伴急性化脓性胆囊炎或伴坏疽性胆囊炎；急性非结石性胆囊炎胆囊坏疽或胆囊穿孔；③ 排除胆管结石。 2.3. 治疗方法 ① 二组所有病人均同一方案常规做好术前准备，采用全麻、造CO 2 气腹12~14 mmHg，头高足低，右侧抬高位，吸净胆囊周围积液，分离胆囊周围粘连，仔细分辨胆囊三角解剖。② 对照组：常规解剖胆囊三角，辨别胆总管，肝总管及胆囊管，结扎胆囊管，逐一顺行或逆行切除胆囊，烧灼胆囊床止血，常规置引流管一根于温氏孔，观察出血及胆瘘情况。胆囊张力过高不易抓持胆囊则于胆囊底部减压，吸净胆汁，取出主要结石，再切除胆囊。若胆囊三角难以解剖，胆囊三角难以辨别则逆行胆囊切除，实在无法完成则行胆囊大部切除术，残余胆囊结石保证取净，有胆汁流出，残端上生物夹，不能上夹子则缝合封闭残端。③ 治疗组：暴露胆囊后首先胆囊体部切开减压，取出主要结石，使胆囊壶腹部及胆囊颈部空虚；选合适位置切开胆囊壁浆膜及浆膜下疏松结缔组织层到达肌层，保留粘膜层及肌层；用吸引器沿解剖层次逐一分离胆囊床部分的胆囊壁；沿胆囊床边缘切除胆囊；沿解剖层次分离出胆囊管并上生物夹，离断胆囊管；烧灼处理出血部位，常规置引流管一根于温氏孔，观察出血及胆瘘情况，结束手术；对于非结石性胆囊炎直接从切开胆囊壁浆膜及浆膜下层至肌层，保留粘膜层及粘膜下肌层钝性分离开始，其余步骤都一样。 2.4. 统计指标 ① 手术时间：从入腹开始到置引流管结束的时间。② 出血量：术中引流总量减胆汁量及腹腔积液量。③ 引流管放置时间：引流量小于5 ml，彩超检查无明显腹腔残余液体即拔管，记录。④ 术后住院时间。⑤ 中转开腹：用以上方法不能完成手术，或有胆管损伤可能，决定中转开腹。 2.5. 统计处理 本资料计量资料数据用均数 ± 标准差(χ ± s)表示，采用t检验；计数资料用χ 2 检验。所有数据采用SPSS11.5软件进行统计学分析。P < 0.05为差异有统计学意义。
system: [SYSTEM]

user: 3. 结果 见表1。
system: [SYSTEM]

user: 4. 讨论 严重急性胆囊炎至今仍是胆囊切除手术难点之一，也是很多大医院LC术中转开腹的指佂。正常胆囊壁解剖：胆囊壁由粘膜层、肌层和外膜组成，粘膜层和肌层均较薄，外膜层较厚，由大量疏松结缔组 Table 1 The treatment results of two group 观察指标 对照组(n = 98) 治疗组(n = 100) t値或χ 2 値 P值 手术时间(min) 术中出血量(ml) 引流管保留时间(d) 术后住院时间(d) 中转开腹(n) 28 ± 46.5 65 ± 36.5 3.0 ± 0.8 5.0 ± 1.8 54.8% (9/98) 85 ± 32.3 42 ± 23.6 2.0 ± 0.5 4.0 ± 0.7 30.0% (2/100) 7.57 5.28 10.57 5.17 4.87 P < 0.001 P < 0.001 P < 0.001 P < 0.001 P < 0.05 表1. 两组治疗结果 织组成，内含血管、淋巴管和神经，浆膜层由肝包膜移行而来。我院对严重急性胆囊炎病理解剖研究发现：① 胆囊壁增厚，较慢性胆囊炎增厚5~10倍；② 胆囊壁组织严重水肿，粘膜层及肌层变化不大，外膜水肿明显，分层明显，切开浆膜层，到达肌层，在肌层与外膜层间水肿带沿解剖层次钝性分离也很容易；③ 结石嵌顿于胆囊颈部后首先出现胆囊颈肌肉痉挛，有限时间内未解除嵌顿，胆囊静脉血管闭塞，胆囊液回流障碍，组织水肿，胆囊压力增高，压迫胆囊动脉致胆囊动脉完全闭塞，出现胆囊壁坏疽，甚至穿孔；非结石性胆囊炎则主要是糖尿病、高血压等病人出现胆囊动脉僵硬、壁增厚、管腔狭窄，特定条件下血管完全闭塞导致胆囊坏疽，甚至穿孔 [ 3 ] ；由于胆囊血管闭塞，此时在胆囊壁粘膜下组织间钝性分离反而较轻型急性胆囊炎容易，钝性分离，挑断血管也不易出血；④ 胆囊炎症水肿严重，波及胆囊床周围肝脏组织，导致肝脏炎症水肿严重，触之易出血，出血后难以烧灼止血，缝合也困难，由于肝包膜已破坏，术后渗血、渗液也较多。因此，保证肝包膜的完整性是决定手术是否顺利的关键。对于以上病理解剖特征，我们总结出了一套经验，取名“浆膜下剥离术”。 与传统腹腔镜胆囊切除术对比发现，“浆膜下剥离术”术中出血少，手术时间短，术后渗液、渗血较少，术后引流管安置时间短，与对照组比，差异有统计学意义，P < 0.001。可能原因是：① “浆膜下剥离术”是钝性分离，血管闭塞，在胆囊壁浆膜下疏松结缔组织与肌层间钝性分离较易，分离过程中几乎不出血，钝性挑断血管也不出血；② 肝脏浆膜层完整，肝脏无出血，免去了肝脏止血时间，减少了肝脏裸露面的渗液；③ 在浆肌层间分离直达胆囊管，胆囊管暴露较容易，无需看清“三管”即可做出判断，结扎胆囊管，胆囊管较细上生物夹容易。 “浆膜下剥离术”胆瘘发生率较常规LC低，可能原因：① 前者胆囊床完整，残余胆囊壁烧灼不出血，迷走胆管瘘不易发生；② 前者胆囊管的暴露在粘膜下层逐层分离即可找到，无需在整个炎症水肿的一片组织中寻找，故胆囊管不易损伤，胆总管也没波及到，难以损伤，中转开腹大多因胆管损伤所致，因此中转开腹也低于对照组，差异有统计学意义，P < 0.05。 胆囊三角的解剖，胆囊周围粘连程度和全身情况是LC是否成功的关键影响因素 [ 2 ] 。“浆膜下剥离术”在浆膜下逆行剥离胆囊管，无需暴露胆总管和肝总管，在胆囊三角致密粘连呈“冰冻状”，难以显露的情况下 [ 4 ] ，胆囊管也能很好结扎，降低了中转开腹率。治疗组2例中转开腹均是由于术前未能及时明确诊断，手术时距离发病时间达7天，腹腔肠管严重水肿，与胆囊、胆管、网膜粘连成片状难以分离，无法暴露胆囊，刻意去分离容易分破肠管，而选择了中转开腹。差异有统计学意义，P < 0.05。 总之，“浆膜下剥离术”充分利用严重急性胆囊炎病理解剖特征，在浆膜下疏松结缔组织与肌层间钝性分离解剖，由于血管闭塞后离断不出血，又充分保留了肝脏包膜，减少出血、炎性渗出机会，出血少，手术时间短，术后渗出少，引流管留置时间短，严重急性胆囊炎是较好的适应症，临床疗效显著，且较易推广。
system: [SYSTEM]

user: 摘 要 方向估计的主要目的是计算出图像等多维信号各点的方向信息，在图像处理和机器视觉的底层处理中具有广泛的应用。在总结现有的方向估计方法的基础上，对现有方法进行了误差比较，从而有利于对这类方法进行深入研究或设计更准确的方向估计方法。 关键词 :方向估计，图像处理，误差比较
system: [SYSTEM]

user: 1. 引言 方向估计可以计算出图像中各点的方向信息，在边缘提取、图像分割、纹理分析、自适应滤波、立体视差估计、图像增强和运动分析中都有直接应用，对位移和光流场的计算即对信号在三维空间中进行方向估计。如果把图像处理分为三个层次：底层视觉、中层视觉和高层视觉，那么底层视觉的主要任务就是提取不同尺度下的边缘及其方向信息。 本文的主要工作是对现有的方向估计方法进行分类总结，并进行误差比较，从而有利于对这类方法进行进一步研究。
system: [SYSTEM]

user: 3. 结束语 本文对现有的方向估计方法进行了分类总结和实验分析，从而可以看出各种方法对噪声的鲁棒性程 度。在此基础上，本文进一步对其中性能较好的两种方法：基于自相关函数的方法和基于正交滤波的方法进行了计算复杂度分析，从而可以看出基于正交滤波的方向估计方法的优越性。我们下一步的工作将对基于正交滤波的方法进行进一步的改进。
system: [SYSTEM]

user: 摘 要 大力发展高等职业技术教育，应有一支师德高尚、业务精湛、结构合理、综合素质高、理论知识丰富、实践动手能力强的“双师型”师资队伍。文章提出了目前高职院校数学专业“双师型”师资队伍建设的现状与存在的主要问题；指出了加强数学专业“双师型”师资队伍建设，首先要稳定队伍建设机制，加大人才培养力度，给予“双师型”人才良好的学习和生活空间；其次要积极推进实践教学改革，提升“双师型”师资队伍建设整体水平。 关键词 :高职院校，数学，“双师型”，师资队伍建设
system: [SYSTEM]

user: 1. 引言 随着高等职业技术教育办学实践的特点和规律，“双师型”师资队伍建设越来越显得重要。任何一所高职院校，没有一支高水平、高技能的师资队伍，就培养不出高素质的人才，也不会提升教育教学质量，相反，一所学校有了高素质的师资队伍和充足的“双师型”人才，那么，这所院校的办学质量和教学水平就一定能够提高 [ 1 ] 。衡量一所院校办学质量和办学水平的高低，不只是评判仪器设备台件数的多少，总额多少，还要看师资队伍建设情况。高职院校近年来办学实践表明，要办好高职教育，必须突出高职特色，办出自己的亮点，依据人才培养总体目标要求，为地方经济建设培养生产、服务和管理一线高尖端、专业技能型人才，要实现这个目标首先要加强“双师型”师资队伍建设，特别是数学专业的教师，光凭着“黑板上画图，只讲X + Y = Z”已经不适应时代发展的要求，必须理论与实践相结合，做到能文能武，上课能讲理论知识，下课能动手实践。假如：一所院校有再多高尖端的仪器设备，无人会用，坏了又无人会修，把仪器当摆设放在那里，试问再多仪器设备有何用？因此，高职院校只有加强“双师型”师资队伍建设，落实人才强校的发展战略，努力打造一支数量适宜、结构合理、综合素质优良的“双师型”师资队伍，才能不断提高高职院校教育教学质量，才能确保高职教育恒久的生命力和可持续发展的再生能力，才能办出自己的特色和水平，同时也是高职院校安身立命的根本所在 [ 2 ] 。
system: [SYSTEM]

user: 2. 高职院校数学专业“双师型”师资队伍建设现状及存在问题 2.1. “双师型”师资队伍总量匮乏 从目前情况看，持有“双证”的专任教师数量不足，多数教师只具有教师资格证书，缺少具有本专业实际工作的职业资格证书。所以，有些教师没有企业、行业实践工作的阅历，专业实践能力不强，社会实践技能不高。 2.2. “双师型”师资队伍硕博比偏低 高职院校“双师型”专任教师的硕博比偏低，多数学校达不到教育部规定要求，一些高学历人才毕业后，不愿跨入职业院校的门槛，而是愿意走向重点院校或科研机构；教育部在职业院校教学水平评估等级中规定标准，专任教师硕博比的A级标准是达到或超过50%，C级标准是达到30%~40%，虽然有些学校达到了C级标准，但距离A级标准还有一定的差距。 2.3. “双师型”师资队伍不稳定 由于不同地域经济发展不平衡，校际间人才竞争日趋激烈，待遇相差较大，在加上高职院校评职较困难、研究场所不足、学校福利偏低等原因，造成高职院校“双师型”师资队伍不稳定。
system: [SYSTEM]

user: 摘 要 本文建立了分光光度法测量海洋沉积物碘的方法。对采集自东太平洋、东海、胶州湾，和福建沿海沉积物岩芯的碘进行了研究。处理样品时，将沉积物样品与KOH混合灼烧，之后用水沥取样品中的碘，在酸性条件下，用饱和溴水将I − 氧化成IO 3− ，加入淀粉和KI混合溶液显色剂，IO 3− 与显色剂中过量的I − 反应生成I 2 ，最后用紫外可见分光光度计进行测量。研究发现用KOH作为溶剂处理样品，沉积物与KOH的质量比以1:3至1:5为宜，灼烧温度和灼烧时间分别为500℃~800℃和1~4 h为佳。测量时最大光吸收波长为580 nm。对碘含量为2~100 mg/kg的沉积物样品，以2 g样品为例，沥取后定容100 cm 3 ，加入6~12滴20 g/dm 3 淀粉和10 g/dm 3 KI混合溶液显色剂较为合适；上机测量的样品溶液pH值在1.8~2.4之间较为适宜。最后我们对测定得海洋沉积物岩芯中的碘含量水平与文献结果进行了比较。 关键词 :碘，分光光度法，海洋沉积物岩芯
system: [SYSTEM]

user: 1. 引言 碘是具有重要生物效应的微量元素之一，关于碘与人类健康的关系已经引起高度重视，所以，人们对环境各种介质和食品等中的碘进行了广泛研究 [ 1 ] 。 全球碘分布主要受海洋系统控制。海洋沉积物是自然界中碘含量最高的地质体 [ 2 ] ，其碘储量占地壳总碘储量的70% [ 3 ] 。碘在沉积物中的分布受到物源及各种物理化学过程及生物过程的制约，同时通过各种生物与地球化学过程，海洋沉积物中的碘又影响海水，从而影响全球碘的输运与分布。也有报道表明海洋边界层的碘影响对流层臭氧的浓度 [ 4 ] [ 5 ] 。所以海洋沉积物碘研究除具有重要的生物地球化学意义外，对全球变暖等目前极为关注的气候变化研究也具有重要的意义。但是人们对海洋沉积物中碘的地球化学研究并不多，国内这方面的研究就更少 [ 6 ] 。 有各种碘的测量方法。包括容量法、分光光度法、色谱法、原子光谱法、电化学方法和ICP-MS法等 [ 7 ] - [ 9 ] 。这些方法尽管都有着自身的优点，但也存在着不足之处，其中有的仪器设备昂贵，不易普及推广；有的则灵敏度较低，耗时长，在海洋沉积物中碘的分析测量中难以广泛开展，这也是至今海洋沉积物碘研究报道较少的原因之一。 用分光光度法测量环境样品中的碘，仪器费用低，数千元价格的分光光度计就可用于多种样品中碘的测量，而且样品前处理方法并不比其他方法复杂。本文建立了分光光度法测量沉积物中碘的方法。研究了测量过程和实验条件，并用所建立的方法测量了多个海域沉积物岩心中的碘含量，研究了碘在海洋沉积物中的分布 [ 10 ] - [ 13 ] 。
system: [SYSTEM]

user: 2. 方法 2.1. 方法步骤 处理样品时将试样与固体溶剂混合，在高温下加热，利用试样与溶剂的复分解反应，使样品中的碘易于被沥取。之后用水沥取样品中的碘。在酸性条件下，用饱和溴水将沥取液中的I − 氧化成 ，加入淀粉-KI显色剂， 与显色剂中过量的I − 反应生成I 2 。最后用紫外可见分光光度计进行测量。实验中所用化学试剂均为分析纯。样品的处理和测量步骤如下： 1) 将样品凉干、磨细。称量约2 g样品，烘干称重，加入铺有约2 g KOH的瓷坩埚中，上面再覆盖约4 g的KOH。 2) 将坩埚放入马弗炉灼烧，慢慢升温至300℃保温0.5 h。然后升温至600℃，保温2 h。 3) 冷却后将样品取出加水，在电热板上加热分散，如果溶液呈蓝绿色，加入数滴乙醇还原Mn 6+ 。 4) 冷却后，将分散物用去离子水转入100 cm 3 离心管中，离心分离。取上清液转入100 cm 3 容量瓶定容。 5) 取10 cm 3 上清液于25 cm 3 比色管中，加入0.2 cm 3 4.7 mol/L H 3 PO 4 ，摇匀后，加入过量饱和溴水，静置2 min；然后加入2滴碱性酚钠溶液使溴水的黄色刚好褪去；最后加入6滴20 g/dm 3 淀粉和10 g/dm 3 KI混合溶液显色剂，摇匀。 6) 15 min后在580 nm波长处测量试样溶液的吸光度，根据标准曲线计算沉积物中碘含量。 7) 同时以同样操作的试剂空白作参比。 2.2. 条件实验 对所建立方法的实验条件进行了研究，测量参数包括测量时所选最大吸收波长、显色剂用量、显色时间、酸度的影响、灼烧样品的KOH用量、灼烧温度、灼烧时间等几方面。实验结果随条件的变化如图1所示。其中，最大吸收波长、显色剂用量、显色时间和酸度的影响用KI标准溶液进行检验，灼烧样品的KOH用量、灼烧温度和灼烧时间用沉积物样品进行。另外还进行了回收率测量。KI标准溶液由分析纯KI试剂制备，制备时将KI 105℃烘干，准确称量0.13080 g，溶于水中，加入2 ml 0.1 mol/l NaOH溶液，纯净水100 ml定容，用时稀释至浓度为5 mg/dm 3 。 2.2.1. 最大吸收波长 用KI标准溶液作为测量样品，按2.1节方法(5、6)步制样，在550~610 nm范围内的选择10个波长，测量吸光值，结果如图1(a)所示，最大吸收波长为580 nm。 2.2.2. 显色剂用量 用KI标准溶液按2.1节方法(5、6)步制样，依次加入3、6、9、12、15、18、21滴淀粉和KI混合溶液显色剂，用分光光度法于波长580 nm测量吸光值，吸光值随显色剂用量变化测量结果如图1(b)所示。当显色剂为6~12滴时，吸光值最大，且变化平稳。 2.2.3. 显色时间 用标准溶液按按2.1节方法(5、6)步制样。加显色剂后，放置0~85 min测量吸光度，见图1(c)。结果表明，吸光度在0~15 min内随时间延长而增加，在15 min后基本保持不变，到60 min后有下降趋势。因此显色15 min后是测量的最佳时间段，且测量最好在60 min内完成。 2.2.4. 酸度实验 分别移取 4 cm 3 KI标准溶液置于8只25 cm 3 比色管中依次加入4.7 mol/L H 3 PO 4 1、2、4、6、8、10、12、14滴，其余按2.1节方法(5、6)步进行操作。结果见图1(d)。4.7 mol/L H 3 PO 4 用量在4~14滴(对应pH值范围约为1.80~2.40)之间吸光度变化不大。故测量溶液pH值在1.8~2.4之间(即加入4~14滴4.7 mol/L的H 3 PO 4 )较为适宜。 图1. 分光光度法测量沉积物中的碘的7种影响因素的条件实验曲线。(a) 最大吸收波长，(b) 显色剂用量，(c) 吸光度随显色时间变化，(d) 试液酸度影响，(e) KOH用量，(f) 灼烧温度影响，(g) 灼烧时间影响 2.2.5. KOH用量 样品量和溶剂量的比例及灼烧温度是样品能否达到良好分解的关键。同一个样品，取6份，每份2 g，分别与2、3、4、6、8、10 g KOH混合。KOH用量为2、3 g时，KOH直接覆盖在样品上面；其他量的KOH，先将约2 g KOH平铺在坩埚中，然后加入沉积物样品，再将其余的KOH覆盖在样品上。按2.1节方法(2~6)步进行样品处理和测量。结果表明(图1(e))，KOH用量在2~6 g范围内，随KOH用量的增加样品碘含量增加；在6~10 g内碘含量基本保持不变。这表明沉积物与溶剂量的重量比在1:3至1:5之间可达到全分解效果。在能达到样品全分解的前提下，应该使用较小的溶剂量，这样可减少浸取液中盐分量及试剂空白 [ 9 ] 。 2.2.6. 灼烧温度 同一个沉积物样品分别取9份按2.1节步骤(1)分别与6 g KOH混合，在200℃、300℃、400℃、450℃、500℃、550℃、600℃、650℃、700℃温度下灼烧样品。之后按测量步骤(3~6)测量样品碘含量，结果如图1(f)所示。在200℃~500℃内碘含量随灼烧温度升高而增加；在500℃~800℃内基本保持不变。这表明在500℃~800℃温度范围内，沉积物中的碘已完全释放被KOH吸收。Muramatsu等人采用高温分解方法处理样品时，发现当灼烧温度在700˚时，大部分碘已释放出来 [ 14 ] 。 2.2.7. 灼烧时间 分别称量2 g同一沉积物样品6份，按2.1节步骤(1)与6 g KOH混合后置于坩埚中，在600℃下分别加热0.5、1.0、2.0、3.0、4.0、5.0 h后，测量碘含量。结果表明(图1(g))，灼烧0.5~1 h测量得到的碘含量随灼烧时间增加；当灼烧时间为1~4 h时，得到的碘含量基本保持不变；当灼烧时间大于4 h后，得到的碘含量随时间逐渐降低。 2.2.8. 回收率 向沉积物样品中加入一定量KI试剂后按分析步骤测量，得到样品回收率为70.8%~83.2%，平均回收率为75.2%。Aumann等人(1981)报道在600℃下用碱熔融法处理沉积物时，其回收率80%左右 [ 15 ] 。本实验方法回收率稍低于Aumann等的回收率。 2.3. 样品测量与数据处理 每一批样品在测量前和测量后做标准工作曲线，测量样品时按以上方法步骤和条件实验选定的参数进行。样品碘含量误差由各测量量的误差按误差传递关系计算得到 [ 16 ] 。 2.4. 方法精密度和检出限 同一沉积物样品取7份，按测量方法测量其中的碘含量。得到的结果平均值相对标准偏差为3.6%。根据国际理论和应用化学联合会(IUPAC)规定，检测限的计算公式为： (1) 式(1)中，L是检测限，S是空白样品碘含量的标准偏差，空白样品的处理方法同样品分析一样。K为标准曲线的斜率。通过对8个空白样品的分析结果计算标准偏差，得到本文方法的最低检测限为0.006 mg/dm 3 。
system: [SYSTEM]

user: 3. 海洋沉积物岩芯中碘分布研究 本实验室从东太平洋、东海、胶州湾和福建近海采集沉积物岩心，测量了其中的碘含量，研究了岩芯中的碘分布。 3.1. 研究海域与岩芯情况介绍 东太平洋岩芯(8˚45'0N， 145˚22'30W)，站位水深5148 m，岩心长32 cm，表层和底层样品厚度为1 cm，其余样品厚度为2 cm，共17个样品。 东海岩心(26˚22'34N， 121˚18'17E)，站位水深79 m，岩心长度为157.4 cm，按1 cm厚度分割样品，共157个样品。 在胶州湾采集了两个岩心，一个长岩心，长271 cm，上层36 cm按1 cm间隔分割，其余按2 cm间隔分割，共155个子样；一个短岩心，长48 cm，按4 cm间隔分割，共12个样品；另外还测定了6个采集自胶州湾海域的表层沉积物和10个采集自胶州湾周边陆地的沉积物样品。 在福建近海三个海域采集了沉积物岩芯。在兴化湾一个岩心，采自湾顶部潮间带，岩心长93 cm，按4 cm间隔分割成23个样品。在东山湾一个岩心，采自西海岸潮间带，岩心长76 cm，按4 cm间隔分割为19个样品。在厦门海域采集了5个岩芯，一个采集自厦门西港海域，岩芯长50 cm，按2 cm间隔分样；一个采集于九龙江口鸡屿附近海域，岩芯长44 cm，按2 cm间隔分样；另外三个在潮间带采集，一个位于九龙江口的北港口，岩心长86 cm，按4 cm间隔分割为22个样品，另外两个采自九龙江口南港海区浮宫附近红树林保护区外的潮间带，岩芯长分别为54 cm和44 cm，按2 cm间隔分样。 3.2. 沉积物岩芯中的碘含量水平及其与其他环境介质碘含量水平的比较 表1给出了海水、岩石、海洋沉积物、海藻、福建土壤和中国海沉积物的碘含量水平。比较表1给出的数据可以看出，海藻的碘含量最高，在100~1000 mg/kg量级；海洋沉积物在1~100 mg/kg量级；海湾和河流沉积物的碘含量低于海洋沉积物；而土壤的碘含量水平更低一些。火成岩碘含量最低。沉积岩碘含量水平也低于海洋沉积物。本研究测定得全部样品碘含量为2.5~105.6 mg/kg，按岩芯平均为15.2~82.4 mg/kg。以东太平洋岩芯碘含量最高，胶州湾次之，福建沿海沉积物中的碘含量总体低于东太平洋、东海和胶州湾，稍高于福建土壤的碘含量。我们的结果总体在文献报道的海洋沉积物碘含量水平范围内。 Table 1 The iodine contents in some environment material 介质 碘含量(mg/kg) 文献 含量范围 平均值 大气 10~20 ng/m 3 Whitehead， 1984 [ 17 ] 海水 0.06 Wong， 1991 [ 2 ] 火成岩 0.7~1.0 Whitehead， 1984 [ 17 ] 沉积岩 0.2~10.0 海洋沉积物 3~400 近代海洋沉积物 5~200 Fuge and Johnson， 1986 [ 18 ] 浅海富氧沉积物 20~220 Harvey， 1980 [ 19 ] 波罗的海沉积物 42~76 Aldahan et al.， 2007 [ 20 ] 密西西比河沉积物岩心 3~34 Oktay et al.， 2000 [ 21 ] 白令海沉积物 50.0~105.0 73.8 高爱国等，2003 [ 22 ] 楚可奇海沉积物 42.5~139.5 46.4 海藻 90~2500 Whitehead， 1984 [ 17 ] 渤海黄海东海南海沉积物 14~22 赵一阳等，1994 [ 6 ] 冲绳海槽沉积物 64 东太平洋 68.2~105.6 82.4 纪丽红，2011 [ 10 ] 东海 35.2~77.8 50.6 胶州湾 51.6~94.3 70.8 环胶州湾陆地沉积物 36.5~50.3 42.0 厦门西港 14.5~26.1 17.2 王荣元，2013 [ 11 ] 九龙江口鸡屿附近海域 13.8~25.7 16.9 九龙江口潮间带沉积物 17.3~39.6 23.8 纪丽红，2013 [ 12 ] ；王荣元，2013 [ 11 ] 兴化湾潮间带沉积物 2.5~43.7 15.2 纪丽红，20113 [ 12 ] 东山湾潮间带沉积物 11.3~29.1 19.5 福建土壤 0.9~14.0 魏复盛等，1990 [ 23 ] 表1. 一些环境介质中的碘含量 3.3. 沉积物岩心中碘分布 我们测定的9个沉积物岩芯，其中5个岩心的碘含量分布如图2所示。东太平洋岩芯用 230 Th ex 和 231 Pa ex 方法测年，其余岩芯用 210 Pb ex 方法测年 [ 24 ] 。岩芯中的碘含量分布分两种情况。在开阔海域采集的岩芯，东太平洋、东海和胶州湾，其中的碘含量随深度变化起伏较大，但总体变化趋势不明显，细分析可能发现存在分段随深度呈增加趋势。福建近海的5个岩芯，均呈现出碘含量随深度降低趋势，而且在某些深度出现极大值。
system: [SYSTEM]

user: 4. 结语 本文建立了分光光度法测量海湾潮间带沉积物碘的方法，并对实验条件进行了研究。用碱熔法处理样品，将试样与固体溶剂KOH混合，在高温下灼烧，使样品中的碘易于被沥取。之后用水沥取样品中的 图2. 五个不同海域沉积物岩心中的碘分布。(a) 东太平洋，(b) 东海，(c) 胶州湾，(d) 厦门鸡屿附近海域，(e) 九龙江口南港附近潮间带 碘。在酸性条件下，用饱和溴水将沥取液中的I − 氧化成 ，加入淀粉和KI混合溶液显色剂， 与显色剂中过量的I − 反应生成I 2 。最后用紫外可见分光光度计进行测量。该方法所用仪器价廉，流程简单，过程时间短，易普及推广。对9个海洋沉积物岩芯中碘的测量，证明了本方法的可用性。 本研究建立的方法平均回收率为75.2%，稍低于文献报道的结果。多次重复试验，并未找出关键的影响因素，推测可能是样品预处理过程中，当碘在无机溶液中处于单质态时，部分碘升华造成的。另外，整个过程还是无法避免使用塑料器皿，这可能是影响全程回收率的第二个主要因素。 我们用所建立的方法研究了多个海域包括东太平洋、东海、胶州湾、福建近海等的沉积物岩芯中的碘含量水平与分布。研究发现，开阔海域，沉积物岩芯中的碘含量随深度变化有起伏，但变化趋势不明显。近岸海域，本文中的福建近岸海域，沉积物岩芯中的碘含量随深度呈逐渐降低趋势。
system: [SYSTEM]

user: 摘 要 在中国经济进入新常态的大背景下，会计在企业中的地位愈发凸显，而会计在核算上存在的问题、风险也越发突出。本文采用逻辑推理、对比分析等方法，总结了税务经济环境与会计核算风险的关系与影响，重点剖析了税务经济环境影响下的会计核算风险表现形式及危害，并较有针对性地提出了有效防范和规避会计涉税风险的相关策略。 关键词 :税务经济环境，会计核算风险，风险防范
system: [SYSTEM]

user: 1. 引言 改革开放以来，我国经济体制逐步由计划经济向市场经济转变。自此，税收也就成为连接国家财政与企业之间最重要的宏观调控手段。在很大程度上讲，国家与企业的关系就是税收征纳关系。会计核算的重要性决定了其在会计职能中占有绝对的主导地位，会计核算在遵从会计准则的同时，必须严格遵守税收法规的规定。一旦会计核算与客观事实出现矛盾不符，便会遭到税务机关查处的可能性，即会计核算的税收风险。可见，在税务经济环境下针对会计核算风险进行研究是十分必要和迫切的。
system: [SYSTEM]

user: 2. 税务经济环境与会计核算风险的关系和特点 2.1. 税务经济环境的涵义 税收作为国家意志的财政聚财工具，有着绝对的强制性和单项性，在调节国家宏观经济的同时，也影响着市场经济的发展。税务经济可以概括成是在税务调控影响下发展国民经济，或者税务行为干预下的经济活动。就经济决定税收关系理论讲，税务经济环境应有两层含义：一是影响税收活动的各种经济因素的总称，二是干预和影响经济运行发展的各种税收体制、征管活动等因素的总称。本文特指后者，即在各税收要素影响下的经济发展环境。由于市场经济条件下，国家对企业由行政控制改为以税收为主导的宏观经济调控，促成了税务经济环境的形成。该环境内，经济领域内所有相关活动都必然受到税收牵连和冲击，而任何税收法律制度与征管体制的变化，都将对社会经济发展速度、结构和质量产生巨大影响和变化。会计作为经济核算中心，成为企业经济力攸关的风险纠结关键点。 2.2. 税务经济环境对会计核算的影响与要求 会计和税收是经济领域中具有关联性的两个重要分支，在市场经济的发展下，会税之间的分离关系愈发显现。会计核算必须遵循公认的会计原则，为管理者和决策者提供真实、完整的财务报表和经营成果及现金流量情况。但税收作为国家宏观调控的重要手段，具有单向性和强制性。虽然会计的产生早于税法，但自从国家产生，纳税就成为社会组织和个人必不可少的一项重要活动。会计作为纳税义务实现的财务工具，也是企业税收信息产生的唯一途径。根据法律至上原则，当会计准则同税收法律制度出现规定不一致时，税法的强制性和权威性要求会计核算应无条件地遵从税收法律制度。所以，税收与会计存在着制约与服务、强制与遵从的天然联系。
system: [SYSTEM]

user: 3. 税务经济环境下会计核算风险表现及危害 3.1. 税务经济环境下会计核算风险的表现 3.1.1. 会计操作风险 1) 技术性风险 实践中，多数会计都按照书本学习的专业套路和执业习惯“见票做账”，严重忽略了对税收的影响。例如，企业支援弱小企业发展过程中实施的无偿且垫付资金的代购商品物资业务，企业在到货后根据仓库转来《验收单》等票证就进行“原材料”和“库存商品”等业务核算，待受托企业交款提货时做出以货币资金直接冲减“原材料”和“库存商品”的账务处理。会计账结平了，风险也就积存下来。一经国税机关检查，企业将为此付出补缴增值税甚至支付税收滞纳金和偷税罚款的惨重代价。 2) 非技术性风险 会计人员遵照会计准则做账、报税，在不知缘由的情况下却遭到税务机关的处罚，这在很大程度上是由于会计人员缺乏风险规避意识。认为会计核算只是数字的叠加，却不重视各要素之间的勾稽关系和核算方式方法等问题。这些不良行为都是由于工作人员的会计核算风险意识薄弱、对税收法规没有清楚的认识和判断，虽然在主观上没有偷税的愿望，但在工作细节和纳税行为上没有按照会计操作规范去操作，遭到税务机关的查处。税法中规定，由于纳税人原因造成少纳税或未纳税的，补缴少纳或未纳税款并按少纳或未纳税额的50%以上、5倍以下处以罚款。可见，企业在无意识地情况下偷税、漏税行为，不仅给企业带来经济损失，更在日后的经营中有损声誉。一些企业通过“关系税”、“人情税”的方式解决纳税问题，但从长远角度看，这种行为不但不能为企业减轻税负，反而这样的违法行为又为企业构成了一项新的纳税风险，成为企业未来发展中的绊脚石。由此可见，在非会计核算技术方面疏漏而可能造成的损失并不亚于技术失误造成的会计风险 [ 1 ] 。 3.1.2. 税收政策风险 1) 政策变动风险 我国现有的税收法律、法规层次较多，部门规章和地方性法规众多，企业要应对国税和地税两套征纳系统，对涉税会计来讲已经是应接不暇，又因为国税和地税在同一方面的规定不尽相同或者模糊不清，导致纳税人无法正确划分业务种类或采用正确的计税方法。加之近期“营改增”政策的全面推行，更使得很多企业应对混乱，无从下手。很多政策通知时间晚、时效短，没有充分的学习理解时间，造成很多涉税会计在没有充分掌握政策的情况下就投入到实际操作当中，再加之会计人员专业素养各异、政策驾驭能力有别，就难免会有疏忽和失误的情况出现，企业很可能因此在计算税额上出现问题。 2) 政策执行风险 客观地讲，由于税收具有较强的国家意志和单向性，企业履行纳税义务就意味着经营成果的减少，因此没有任何一个纳税人心甘情愿地付税。而企业往往处于为了纳税而纳税的现实情况，对税收政策理解一知半解就盲目照搬之前经验等来执行税收政策，不可避免地出现大量失误和疏忽甚至出现多纳或少纳税的错误，如记反了借贷方向、纳税金额核算错误、拖延报税时间等等。对于多缴的税款不仅是对企业经营成果的侵蚀，在三年内若没有得以更正将造成永久性的经济损失；而对于少缴的税款，一旦被税务机关查实，纳税人必将补缴税款、加收滞纳金、罚款等，构成偷税的，还将有面临被追究刑事责任之忧。 3.1.3. 领导意志风险 会计信息作假是一种企业行为，并非只是会计人员作假。在利益驱使下，管理层只考虑自身的局部利益，并不顾全大局甚至违背国家利益，领导者在主观上就缺乏对税收的正确认识，致使出现管理层授意指使会计人员对会计信息进行美化和粉饰。如通过隐瞒、篡改、伪造纳税凭证等，隐瞒收入、虚列支出等作假帐的行为，已达到不交税，少交税的目的。虽然在财经法规中赋予会计人员一定的执行权利，对于领导授意致使的违法行为，会计应对以拒绝和举报，但若想严格遵守并不具有现实性。在纳税遵从和领导遵从之间，会计常两难选择。 3.2. 税务经济环境下会计核算风险的危害 税务经济环境下，会计核算的最普遍风险就是错计、错纳税的可能性，其中最重要的当属缺乏税收牵连的周密考虑而盲目核算造成偷税处罚等可能性。这种核算风险一经税务审计查实，企业遭受经济损失在所难免，甚至招致灭顶之灾。 例如，某公司在无需垫付资金条件下，代友好企业采购材料和产品化验设备共351万元，业务内容与会计核算流程如下： 1) 收到代购材料款 借：银行存款 3，510，000.00 贷：预收账款 3，510，000.00 2) 将款项汇给供应方账户 借：预付账款 3，510，000.00 贷：银行存款 3，510，000.00 3) 代购货物到达、库管验收后 借：原材料 2，340，000.00 在建工程 1，170，000.00 贷：预付账款 3，510，000.00 4) 委托方提回材料和设备仪器 借：预收账款 3，510，000.00 贷：原材料 2，340，000.00 在建工程 1，170，000.00 至此，会计对代购业务处理完毕，账面余额清零。后来，国税分局检查认定该业务系采用预收款结算方式的货物与固定资产销售而不反映收入，造成偷增值税51万元，要求企业补缴增值税，加收滞纳金22，950元(510，000 × 0.5‰ × 90天)，并处以应纳税额50%罚款25.1万元，合计支付783，950元，令企业叫苦不迭但也只能遵从处理。该业务中最突出的问题有两点：一是会计将收取代存款错记成反映交易内容的“预收账款”，给检察人员留下销售业务悬念；二是会计误将代购材料设备记入本单位资产账，继而在友好企业提走材料设备环节一错再错地用“资产”冲销“负债”，验证了检查人员对销售业务的判断。企业支援友好企业的一种“义举”，在会计的错误核算下，竟然承受了78万多元的“冤枉债”，最让人委屈的是还背负“偷税黑锅”。 企业在日常生产经营和经济核算中都离不开纳税款项的核算，而能否及时、准确的核算涉税金额更成为企业整体会计职能中重要的体现，尤其在税务经济环境下的会计核算更显得责任重大、使命突出。与此同时，由于会计核算差错而产生的错误的会计信息不仅会给市场经济环境带来不安定的因素，产生不良影响、招致税务机关查处，而且会给会计信息的使用者如债权人、投资者等产生误导，造成不同程度的损失。另外，在社会风气日下的环境下，也给贪污腐败、行贿受贿等违法行为留有可乘之机。可见，涉税会计核算的风险如不能得到有效控制，不仅企业和会计自身会遭受损失，国家也很可能付出昂贵的代价。
system: [SYSTEM]

user: 4. 税务经济环境下会计核算风险的防范策略 4.1. 端正领导层纳税态度 企业的领导层应正确看待纳税，税款毕竟是收入中的一小部分，是完全可以通过科学合理的管理筹划实现“增收减支”的，一味只追求纳税最少并非是对企业日后的经营管理最好的选择 [ 2 ] 。企业管理层应该在坚持依法纳税、合法经营的基础上，重视会计环境的建设，杜绝打税收“小算盘”的不良现象出现。领导者应将眼光多放在生产经营的各个环节而不是只盯在税负上，给予会计参与合同签订、采购、资产处置等管理方面的决策的权利，时刻以会计核算风险为导向，将会计的管理监督职能渗透到企业经营的各个环节方面。企业应不断加大会计队伍的人力资本投入，支持会计积极参加各种职业技能和财税业务培训活动，督促会计努力晋升职称，调动会计虚心学习研究业务和认真处理涉税经济业务的积极性，真正使会计成为企业经济利益的“保护神”。 4.2. 提高会计核算技术水平 企业的会计核算人员，应清楚地认识到自己所处的税务经济环境的严肃性，对所核算的经济业务要时刻保持小心谨慎的态度，时刻谨记自己所做的账目是给税务部门和审计机构查看的，必须从建账、记录到核算、编制报表都严格依照会计准则规范和税收法规规定执行，更不能有漏记少记的现象发生。实践中，会计未能按照权责发生制全面核算导致当期实际发生的支出少记，反之却设法虚列成本费用调节盈亏，比如期末结算前不清理外欠款业务发票来核算相关成本费用即违背了权责发生制，会计反而通过虚列成本费用调节盈亏造成偷逃所得税，这是很多企业都普遍存在的愚蠢会计行径。所以，财务会计对会计原理、准则的深刻理解和熟练运用至关重要，对纳税检查对象的账证表单和会计记录要十分认真百倍谨慎，防范掩埋会计核算风险。 4.3. 培养涉税会计人员的综合素质 从实际角度讲，一位综合素质过硬的涉税会计人员将是企业在瞬息万变的经济环境中平稳、安全的得以发展的必要保障。企业在注重专业素养的同时，必须关注会计的整体素养，一位只重视专业技术而缺乏基本职业道德的涉税会计可能造成的后果是不可想象的。涉税会计的综合素质是其做好纳税核算工作的基础，不仅有过硬的专业技术、政策驾驭能力、较高的财税法规认知实力，在税收法规与会计准则出现出入时，能够自觉、准确地进行调整。同时更要有一定的职业道德和良好的工作作风，不断提高自身的文化修养，定期对会计人员进行考核测试，从实际业务操作能力到核算税款的精准程度及与税务机关的沟通能力等等，在加强会计人员的紧迫感同时，也培养了高素质的综合型涉税会计人才。 4.4. 设立专门涉税会计岗位 与财务会计不同，涉税会计必须时刻关注国家税收有关政策的实施和变动，一旦没有及时掌握，都有可能对企业造成经济损失和商业失信。因此，在条件允许的企业应设立专职纳税会计岗或单独的涉税会计部门，专门负责税收相关工作，并赋予涉税会计高于财务会计的职业地位，在企业中享有直接向企业管理层汇报的权利。没有条件的企业必须做到涉税会计人员及时关注税收政策，及时与税务机关沟通，定期聘请税收专家顾问到企业来进行实操指导，在此基础上也可向中介机构进行咨询，做到问题及时解决、不留隐患 [ 2 ] 。
system: [SYSTEM]

user: 5. 结论 总之，对处于税务经济环境下企业而言，会计核算风险是企业面临的最主要的风险之一，企业应设立必要的纳税会计岗位的同时，赋予会计以更多的经营参与权和内部控制话语权，努力将会计理财思想和纳税筹划意识渗透于企业生产经营的方方面面以及微宏观决策的全部过程。会计也要珍惜自己的专业声望和职业荣誉，在不断潜心钻研本岗业务技能的同时，始终用风险意识武装自己的工作活动和行为，谨慎地对待每一笔经济业务，坚持冷静分析和准确无误地确认业务性质，确保企业经济和声誉上的共同发展。
system: [SYSTEM]

user: 摘 要 钙离子是细胞内广泛存在的第二信使，控制着几乎所有的生命过程。本文以钙信号的生物背景为基础，首先介绍了构建钙离子信号网络模型的数学方法，包括生化反应建模、参数拟合等；接着我们介绍了关于1，4，5-三磷酸肌醇受体(IP 3 R)通道和线粒体钙单向转运体(MCU)通道的各种模型；然后我们讨论了关于钙离子信号震荡的各种点源模型，包括内质网释放的钙信号震荡模型和内质网-线粒体相互作用调控的钙信号震荡模型；最后我们建议未来对钙信号相关的细胞信号网络建模展开进一步的研究。 关键词 :钙离子通道，钙信号，建模分析
system: [SYSTEM]

user: 2. 相关数学基础 数学建模就是将生物系统之间的相互作用用数学语言表现出来。下面就介绍几个建立生物模型的基本数学工具。 2.1. 微分方程建模 生物化学反应主要是通过常微分方程(Ordinary Differential Equations， ODEs)来表示，因为ODEs可以比较准确地描述系统动力学的变化，即各个参与反应物质的浓度随时间的变化 [ 38 ] 。ODEs构建的基本方法主要是质量作用定律。质量作用定律是19世纪G.M.古德贝格和P.瓦格提出的，主要内容是：化学反应速率与反应物的有效质量成正比 [ 39 ] 。其中有效质量是指反应物的有效浓度。 假设有两种反应物A，B相互作用，生成产物AB。其反应方程为： ， (2.1) 其中[A]，[B]分别代表反应物A，B的浓度，[AB]代表生成物AB的浓度；k 1 和k −1 分别为结合速率常数和解离速率常数。结合速率r 1 和解离速率r 2 分别表示为： ， (2.2) 。 (2.3) 对于反应物A，B的浓度变化主要是在正反应过程中减少，在逆反应过程中增多，而生成物AB的浓度变化则与A，B相反；用公式可以表示为： ， (2.4) 。 (2.5) 当反应达到平衡状态，即结合速率与解离速率相等，反应物与生成物地浓度都保持不变，可以得到： ， (2.6) 其中，K d 称为平衡常数(equilibrium constant)，又称为解离常数(dissociation constant)。K d 越大，说明反应物A，B的结合能力越弱，生成物AB越少，AB的解离速度更快。 2.2. 马尔科夫过程 马尔科夫过程是以俄国数学家安德烈·马尔可夫命名的随机过程，马尔科夫过程被认为是“无记忆”的，因为系统可以仅仅根据其目前的状态进程对未来作出预测，并不考虑之前的状态。 在生物系统建模中，微分方程可以用来表示生化反应过程。但是，对于随机性比较强的反应过程，微分方程就不能具体描述。因此，马尔科夫随机过程就成为解释生物过程中随机过程的一种重要的手段。对于钙信号模型，由于钙门控通道的随机性特点，其通道亚基的激活、失活和恢复过程就可以用马尔科夫过程来描述。 2.3. 参数优化 由于生物建模过程中有许多未定的参数，为了尽量增大模型的输出结果与实验结果的达到更好的拟合效果，尽量还原各种生物学反应。对于模型来说，参数优化是一个重要的课题。本节主要对模拟退火算法、遗传算法和免疫算法作简单的介绍。 2.3.1. 模拟退火算法 模拟退火算法成为了一种求解大规模组合优化的随机性方法。它是由美国物理学家N. Metropolis等人在1953年提出的 [ 40 ] ，1983年，S. Kirkpatrick等人将模拟退火思想引入到组合优化的领域中 [ 41 ] 。 模拟退火算法是基于蒙特卡洛迭代法，其思想来源于固体物质的退火过程。因为在自然规律中，物质总是趋向于最低能态，因为最低能态是最稳定的状态。 模拟退火算法的基本思想是先把固体加热至足够的温度，使固体中的所有粒子处于无序的状态，然后慢慢地将温度降低，使粒子活动逐渐有序，在常温时粒子达到最低能态，即寻找到最优解。为了防止陷入局部最优解，常常使用Metropolis准则。 Metropolis准则是根据概率的大小判断是否接受新的状态。在高温下，可以接受与当前状态能量差较大的新状态；而在低温下，只接受与当前状态能量差较小的新状态。 模拟退火算法具有可并行性、扩展性和通用性，它可以高效的解决几乎所有的组合优化过程。 2.3.2. 遗传算法 随着达尔文的进化论和孟德尔的遗传学说被广泛应用之后，在计算机领域，人们将这两种学说应用于算法中，得到了遗传算法。遗传算法是由美国Holland教授首先提出的 [ 42 ] ，是一种通过模拟自然进化过程随机搜索全局最优解的方法 [ 43 ] 。 遗传算法将问题的求解过程模拟为群体的生存过程，通过一代代的优胜劣汰的筛选进化，找出问题的最优解。它的基本思想是，首先随机产生初代种群，然后根据进化理论，逐代进行优胜劣汰的选择，产生越来越接近的最优解。在每一代中，根据适应度的大小来选择个体，并运用遗传算子对个体进行交叉和变异，改变个体的基因，从而产生新的解的解集。 由于遗传算法具有使用方便、可以并行计算、自适应强等优点，因而在函数优化、组合优化、人工智能等领域被广泛应用。但是遗传算法也存在一定的不足之处，比如遗传算子可能引起系统退化，遗传算子相对固定，缺少灵活性等。 2.3.3. 免疫算法 因为遗传算法存在着一定的不足，因此人们在保留遗传算法优点的基础上，为了避免其缺陷，引入人体的免疫系统机制，形成了免疫算法。免疫算法是人们模仿生物免疫系统基本机制，模拟免疫系统对外来病菌的识别能力而设计出来的多峰值搜索方法 [ 43 ] ，最早的免疫系统的模型是在1973年由Jerne提出的 [ 44 ] 。 免疫算法的基本思想是：对于具体要解决的问题(抗原)，首先是必须对问题进行具体分析，提取问题最基本的特征信息(疫苗)，对特性信息进行分析处理，得到待解决问题的一种方案(抗体)，然后对抗原和抗体进行亲和力分析，如果没有得到满意的结果，免疫算子进行具体的操作，形成新的抗体。免疫算子主要通过接种疫苗和免疫选择两个步骤完成，接种疫苗是为了确保提高亲和力，而免疫选择是为了防止群体退化 [ 45 ] 。
system: [SYSTEM]

user: 4. 内质网调控Ca 2+ 模型的研究 在本综述中，我们将忽略钙信号的空间扩散行为，而关注于钙信号的震荡动力学，所以我们将只讨论钙信号震荡的点源模型，不考虑空间扩散动力学，所以其数学描述仅仅是常微分方程。本章主要讨论内质网调控Ca 2+ 模型。 4.1. 全局钙震荡模型 在DeYong-Keizer模型和Li-Rinze模型中，细胞质中的钙浓度是随着时间震荡的，而IP 3 浓度是不变的。细胞质内的钙浓度都是由IP 3 R通道、漏通道和钙泵三项构成，其微分方程为： 图8. MCU的4态动力学模型 [ 48 ] ， (4.1) 其中，[Ca 2+ ]表示细胞质内自由的Ca 2+ 浓度，J 1 表示由于IP 3 R通道打开引起的Ca 2+ 的外流，J 2 表示由漏通道流出的Ca 2+ 。J 1 和J 2 项是表示Ca 2+ 的外流，引起细胞溶质中的[Ca 2+ ]升高。J 3 表示通过钙泵将Ca 2+ 逆浓度梯度从细胞质中输送至内质网中，是引起细胞溶质内[Ca 2+ ]减小的过程。 因此，J 1 项由方程可以表示为： ， (4.2) 其中，v 1 表示IP 3 R通道的最大流量。c 1 表示内质网与细胞质的体积比。 表示内质网中的[Ca 2+ ]。P open 则是前章所计算的通道开通几率。 J 2 项可以用方程表示为： ， (4.3) 其中，v 2 表示漏通道的最大流量。 内流的Ca 2+ 主要是通过钙泵(SERCA)将Ca 2+ 逆浓度梯度输送到内质网中，其表达式主要是通过希尔方程表示： 。 (4.4) 根据上述的数学模型，DeYong-Keizer模型和Li-Rinze模型都可以得到了钙的震荡曲线。图9是DeYong-Keizer模型得到的分岔图。 4.2. Rudiger-Shuai局域钙模型 2010年，Rudiger和Shuai在PRL上发表了一篇关于运用决定性和随机性模拟方法模拟单个集团内几个IP 3 Rs通道释放局域钙信号的网络模型 [ 56 ] 。在模型中，9个IP 3 R通道均匀地分布在立方体的一个表面上(如图10(a)所示)。由于IP 3 R通道集团相对较大，并且Ca 2+ 扩散速度相对较快，同一集团内的IP 3 R通道往往处于不同浓度的Ca 2+ 环境中，因此他们对于关闭和开放的IP 3 R通道引入了不同的[Ca 2+ ]。对于一个开放的IP 3 R通道，其通道口的[Ca 2+ ]可以达到100 μM甚至更高；而对于关闭的IP 3 R通道取决于其他打开的IP 3 R通道。 图10(b)可以看出，关闭通道口的平均[Ca 2+ ]与打开通道数成线性关系。因此，可以用如下公式表示： ， (4.5) 其中， 表示某一关闭通道得到其他n个打开通道扩散的[Ca 2+ ]的平均值。 根据上述公式，Rudiger-Shuai模型很好的模拟出了实验中所观测到的局域钙Puff信号幅度分布的结果，为之后的理论建模提供了新的思路。 4.3. Qi局域钙模型 2014年，Qi建立了一个双Ca 2+ 浓度的模型研究IP 3 R通道集团内的产生的局域钙blip (单个通道打开)和puff (多通道打开)信号动力学行为 [ 57 ] 。该模型由一个简单的钙动力学方程和一个马尔科夫过程决定。对于内置网上的Ca 2+ 通道，在该模型中只考虑IP 3 R通道对Ca 2+ 动力学的影响。简化模型的同时，还考虑了Ca 2+ 结合蛋白对Ca 2+ 动力学的影响。对于IP 3 R通道的开关动力学，主要是由修饰过的DeYong-Keizer模型来讨论，其中亚基的状态变换是由马尔科夫动力学决定。 图9. 分岔图。图中显示了稳定和不稳定状态，给出了在稳定轨道周期内细胞溶质钙浓度的最大值和最小值 [ 46 ] 图10. Ca 2+ 的释放过程。(a) 8*8*5 m 3 的立方体表示细胞质溶质，IP 3 R通道分布于立方体的某一表面上。(b) 所有关闭通道口的平均[Ca 2+ ]与打开通道数的关系 [ 56 ] 对于打开和关闭的IP 3 R通道周围的Ca 2+ 浓度分别采用不同的[Ca 2+ ]。对于开放的IP 3 R通道，它将产生一个很高的自抑制浓度，在模型中，设置[Ca 2+ ] = 120 μM。而对于关闭的IP 3 R通道，根据Rudiger-Shuai模型 [ 56 ] 的推导，在平衡状态下关闭通道平衡时周围的[Ca 2+ ]可以表示为： ， (4.6) 其中，[Ca 2+ ] Rest 为静息状态下细胞溶质中的[Ca 2+ ]，∆C为打开通道的Ca 2+ 扩散到邻近通道，引起[Ca 2+ ]上升的平均值，N Open 为打开的通道数。 当一个通道由打开状态转向关闭状态，达到平衡状态的[Ca 2+ ] Equil 之前，[Ca 2+ ]会受到Ca 2+ 的扩散和Ca 2+ 缓冲蛋白的影响。因此，对于关闭状态的[Ca 2+ ]，可以用微分方程表示为： ，(4.7) 其中，γ为Ca 2+ 的衰减率。 根据上述模型，可以得到集团数目和[IP 3 ]的增多会减少blip的比例和持续时间。相邻局域钙信号之间的时间间隔及第一次释放事件的响应延迟时间与集团大N的倒数成线性关系等结论。图11是上述模型得到的[IP 3 ]对blip比例和持续时间的影响。 4.4. Cao局域钙模型 2013年，Cao发表了一篇在单通道数据的基础上得到的puff信号的随机模型 [ 58 ] 。该模型的IP 3 R通道模型是引用了Siekmann的IP 3 R模型 [ 53 ] 。Cao的模型除了考虑IP 3 R通道对Ca 2+ 动力学影响，还考虑了漏通道、快速缓冲蛋白和慢速缓冲蛋白dye(fluo-4)对钙信号的影响。在模型中，主要的动力学微分方程可以表示为： ， (4.8) ， (4.9) 其中， 表示可以是细胞质内钙浓度升高的钙流。 表示使Ca 2+ 浓度减少的钙流，主要是通过Ca 2+ 的扩散和钙泵的影响。 表示内质网上漏通道的影响。 和 分别表示dye的浓度和与Ca 2+ 结合的dye的浓度。 表示打开的通道数。 Cao根据以上的模型，发现在小IP 3 R集团内，puff的平均幅值与集团大小成线性关系，而对于较大IP 3 R集团来说，呈非线性关系。Puff的衰减时间对集团大小敏感性比上升时间对其敏感性更强。 4.5. 全局钙和IP 3 震荡模型 前面介绍的模型中，IP 3 的浓度是不随时间改变的一个常数值。而对于有些类型的细胞，其IP 3 的浓度是由钙离子浓度所影响的，这样钙离子浓度的震荡将导致IP 3 震荡，这里我们也讨论其相关的建模。 Sneyd等人 [ 59 ] 在PANS上发表一篇关于IP 3 浓度影响钙震荡的模型。该模型主要包含了一个钙震荡方程和一个IP 3 震荡方程，表示如下： ， (4.10) ， (4.11) ，(4.12) 。 (4.13) 图11. (a) blip事件的比例和(b) blip持续时间在不同集团大小情况下随IP 3 分子浓度的变化 [ 57 ] 其中，[Ca 2+ ]表示全局钙浓度，[Ca 2+ ] ER 表示内质网中的钙浓度，[IP 3 ]表示IP 3 的浓度，n表示IP 3 受体的比例。通过该模型，他们很好的模拟出了在胰腺腺泡细胞中得到的IP 3 浓度影响钙震荡的实验结果(如图12所示)。 Shuai等人 [ 60 ] 在文章中用了另一种方程表示IP 3 分子的生成。可以用方程表示为： ， (4.14) 其中，IP 3 的降解J D 表示为： ， (4.15) IP 3 的生成可以表示为： 。 (4.16)