1、1基于核磁共振技术的牙膏含氟量测量方法 俞 捷 1，陈珊珊 1，李 然 1，汪红志 1,2,3，张学龙 1,2,3 （1. 上海理工大学医学影像工程研究所， 2. 上海医疗器械高等专科学校，上海 2000933. 上海纽迈-医专核磁共振分析测试中心，上海，200062）摘 要 核磁共振是一种物理现象，作为一种分析测试手段在许多领域都有着极为广泛的应用，如地质勘探、石油化工、食品农业等。本文应用核磁共振技术对牙膏中的含氟量进行测量方法研究，包括应用相对测量、实验原理方法以及相关细节。目的是运用一种新的方法来快速便捷的测量工业生产中牙膏的含氟量，并分析此方法的优缺点、精度和误差。关键词 核磁共振；
2、FID 牙膏；含氟量 中图分类号 O482.53；R318.6 文献标识码:A 文章编号: Measurement Method of Fluoride Content in the Toothpaste Based on NMRYu jie1，Chen Shanshan 1，Li Ran 1，Wang Hongzhi 1， 2， 3，Zhang Xuelong 1， 2， 3+ (1,Medical Instrumentation College, 2,University of Shanghai for Science and Technology, Shanghai 200093, P.
3、R.China； 3,NMR analyst center of niumag and Medical Instrumentation College, shanghai ,200062, P.R.China)Abstract：NMR is a physical phenomenon, being an analytical tool in many areas, it has a very wide application, such as food, geologic, industrial ,chemical and agricultural analysis. In this arti
4、cle, you can see an NMR-Based method to measure the fluoride content in the toothpaste. It consists of the application of the relative measurement, methods and details of the experimental principle. And the aim of the experiment is to use a new, fast and convenient approach measure the fluoride cont
5、ent in industrialized producing toothpaste, and to analyze the advantages and disadvantages of this method, precision and error. Keywords ：NMR; FID; Fluoride Content含氟牙膏在我们护齿健康方面有着不可替代的作用，它是在牙膏中添加少量的含氟物质（CaF 等） ，从而使牙膏中含有少量氟。含氟牙膏是一种局部外用剂，刷牙时有极少氟吸附在牙齿和口腔表面，发挥防龋作用，但是绝不能过量。人体中确实也需要少量的F微量元素，尤其是儿童，但过量摄入氟容
6、易造成中毒。资料表明在氟中毒的3-6岁儿童中，有80%的儿童是因为误吞牙膏导致，这点尤其应该引起重视 1。2009年我国国家标准委员会特别对牙膏的氟含量提出要求，实行牙膏新的国标标准GB8372-2008 规定成人牙膏总氟量在0.05%-0.15%之间；儿童牙膏更严格，氟的含量在0.05%-0.11%之间 2。检测牙膏含氟量的方法很多，传统的主要有硝酸钍滴定法 23和离子色谱法两种。目前国家规定用硝酸钍滴定法来测定水中氟含量，这一方法适用于含氟量大于50mgL -1的废水中氟化物的测定。从原理上来讲，利用氟离子浓度差异导致的颜色差异来测量浓度，其中存在着很大的主观因素。因为溶液颜色的转变以及精
7、确度定位在检测中是至关重要的，但是如何来定性所谓的灰蓝色，这个问题直接影响了整个实验的准确性。基金项目 上海市重点学科（P0502 ） ；上海市高校科技发展基金（07ZZ93）和上 海 市 重 大 科 研 项 目 基 金 （ 06DZ11311） ；上海市科研创新项目（10YZ237）资助项目。作者简介 俞捷（1986） ，男，上海人，汉族，硕士研究生。研究方向：医学影像工程。通讯作者 张学龙（1953） ，男，江苏淮安人，汉族，教授。研究方向：医学物理与工程。作者单位、地址上海市营口路 101 号上海理工大学医疗器械学院， 200093。E-mail：收稿日期2010-04-29 2离子色谱
8、法，在精确度上比硝酸钍滴定法有了很大的提高，与标准色谱的比较只要经过简单的做个比值就能确定F离子的信息，但是整个测量过程中要用到大量的仪器和大量试剂，分析过程也很较长，其中涉及的许多操作过程还需要操作者需要一定的专业知识，这一点和我们所需要的快速便捷的方法相去甚远。本文主要探索应用核磁共振技术，从检测氟的核磁共振信号强度与氟离子含量之间的对比关系完成牙膏氟含量的测量 24 。一 、核磁共振测试技术核磁共振（Nuclear Magnetic Resonance，简称 NMR）是磁矩不 为 零 的原子核在外磁场作用下自旋能级发生塞曼裂分，共振吸收某一定频率的射频辐射的物理过程，也可以理解为于处于静
9、磁场中的原子核在另一交变磁场作用下发生的物理现象。通常人们所说的核磁共振指的是利用核磁共振现象获取分子结构、人体内部结构信息的技术。作为一种对检测物质无任何伤害的安全、快速、准确的检测手段， 核 磁共 振 技 术 已 在 生 物 学 、 医 学 、 化 学 等 领 域 扮 演 着 极 为 重 要 的 角 色 ， 尤 其 是 在 医 学 方 面 ， 核 磁 共 振 更 具有 开 拓 性 。 而 近 年 来 ， 核 磁 共 振 技 术 也 被 逐 渐 更 为 广 泛 的 应 用 于 食 品 农 业 、 地 质 勘 探 、 石 油 化 工 、 医药 材 料 等 领 域 ， 成 为 一 种 越 来 越
10、 重 要 的 无 损 检 测 技 术 。在 静 电 场 中 ， 具 有 磁 性 的 原 子 核 存 在 不 同 的 能 级 。 当 用 一 特 定 频 率 的 电 磁 波 (射 频 场 )照 射 待测 物 时 ， 当 电 磁 波 能 量 等 于 能 极 差 时 ， 原 子 核 吸 收 电 磁 波 发 生 能 及 跃 迁 ， 产 生 共 振 吸 收 信 号 ， 这 即 是核 磁 共 振 的 整 个 过 程 。 通 过 分 析 磁 共 振 基 本 原 理 我 们 可 以 看 出 样 品 的 质 子 数 与 接 收 的 幅 值 信 号 成 正 比关 系 ， 于 是 我 们 得 出 SN式 中 ： N
11、 为 样 品 的 质 子 数 ， S 为 接 收 的 信 号 幅 值 。 本 实 验 便 是 基 于 信 号 幅 值 测 量 的 原 理 来 进 行 试 验 测量 。测量信号的方法有绝对测量与相对测量两种，绝对测量是指对某具体样品直接进行信号获取，通过信号幅值对设备、测量条件及其他因素进行各种修正因子的校正，得到含量结果。对于核磁共振直接测量而言，要考虑到温度变化、谐振点偏移、增益误差等诸多因素，因此复杂且精度低。通常采用相对测量法，即先采用已知含量的标准样品（含量为A sta），在某测量条件下获取核磁共振信号幅值( Ssta)。然后在测量条件完全相同的情况下，对未知样品进行核磁共振扫描，测量
12、信号为(S)，得出未知样品的含量为：A=SAsta/Ssta此式为信号-含量之间的刻度直线 A=kS，为一条过零点的直线，k=A sta/Ssta为含量与信号之间的比例关系，即刻度直线的斜率。实际过程中由于存在着其他因素导致的本底，因此刻度直线不一定为过零的直线，因此需要采用最少2个标准取样的含量与信号之间的比例关系才能确定出刻度曲线：A=kS+b为了减少标准取样较少导致的比例误差，采用多个标准取样，最后通过最小二乘法对多个标准样的含量- 信号关系拟合出最佳的 K值。该方法具有快速、简便、计算方便的特点，一次定标，可快速进行样品检测。相对测量法中，测量条件的完全相同是很关键的，能够创造完全相同
13、的测量条件是精确测量的前提。二、测试设备与参数选择根 据 核 磁 共 振 的 条 件 可 知 ， 质 子 数 与 中 子 数 至 少 有 一 个 奇 数 的 原 子 核 均 可 发 生 核 磁 共 振 现 象 。 氟-18的 质 子 数 为 奇 数 ， 因 此 也 是 共 振 核 。 其 自 旋 量 子 数 为 12， 天 然 丰 度 为 100%， 相 对 灵 敏 度 为0.845， 也 较 高 ， 其 共 振 频 率 约 为 H-1的 94.067%5。 因 此 ， 其 通 量 核 磁 共 振 信 号 与 H-1相 差 不 大 ， 在H核 核 磁 共 振 设 备 上 ， 只 需 调 整 设
14、 备 探 头 的 共 振 频 率 ， 通 过 累 加 扫 描 即 可 得 到 可 观 的 氟 信 号 。 氟虽然是强核（旋磁比高），天然丰度高，但由于牙膏中氟的含量很低，因此采用核磁共振技术检测氟含量3需要在仪器设备和测试技术方面进行相关改进和探索， 1、探头改造本研究采用 NMI20 核磁共振仪作为测量设备。NMI20 核磁共振仪主要用于检测 H 核信号。由于 F和 H 的共振频率不同，为了检测氟信号，需要将设备探头从氢的共振频率调谐到氟的共振频率上。H 核和 F 核的旋磁比比例为 F/H=94.06%。通过软件测得测量用仪器的探头共振频率在恒温 320C 时为21.91MHz，因此相同温度
15、下 F 的共振频率应为 21.54 MHz。接下来通过网络分析仪调节探头的可调谐振电容使探头谐振在 21.54MHz。谐振电容的调节对于探头匹配性能也会造成影响，因此同时还需要调节匹配电容，反复多次调节后使探头能较好的匹配并处于氟信号的谐振点上。同时对于微量氟信号的检测，还需要探头具有较高的品质因素（Q ） ，一般要求 Q30dB。否则对于浓度低的氟含量检测误差较大，甚至难以检测。2、测试参数选择基于幅值的分析测试技术一般采用硬脉冲 FID 序列。它是 90硬脉冲激励后采集 FID 信号的序列。相关主要参数如下 ：P 1(s):射频硬脉冲的宽度；D 3（S ） : 90射频结束到线圈开始采集信
16、号之间的时间间隔； D0（ms）:脉冲序列的重复时间，即 TR；TD:NMR 信号的采样点数（复数点） ；SW(kHZ): 谱宽，采样频率；DFW(kHZ): 数字滤波器的截止频率； SF1(MHZ):射频信号频率 f 的主值；O 1(kHZ):射频信号频率f 的偏移量；NS: 累加次数；DS:过采样倍数；图1、图2中可以观察到对测量信号起到了重要的作用参数：NS、D 3、采样时间T。NS即累加次数，适当的选取能够让设备接受更好的信号，同时减少噪声的影响。在D3的选择中，测量中由于线圈自身外包裹含氟物质加之测量物又为含氟物对于准确的测量造成了不小的问题，所以对D 3的取值变得尤为关键。为了确保
17、测量到的信号为线圈自身氟化物的信号，放入一只载有水的试管于探头中，用同样的方法进行测试。由于氟的共振频率与氢很接近所以最有可能的影响还是来自于氢物质，如果探头会捕捉到氢的信号的话，那之后溶液中大量的水中所含的氢信号将在很大程度上影响测量，实验结果的信号并没有明显变化，可以确认测得的是氟化物的信号，由此可以证明放入仪器中的水并没有影响检测，可以对之后的溶液进行测试。 图 1.硬脉冲 FID 序列的序列形式图2.参数设置 Parameter Value ValueP1(s) 15.0 15.0D3(s) 50 200D0(s) 500 500TD 650 650SW(KHZ) 200.0 200.
18、0DFW(KHZ) 30.0 30.0SF1(MHZ) 21 21O1(KHZ) 567 567RG 4 4NS 256 256DS 10 104当取D 3=50s对线圈进行不载物的空扫时，如下图 3得到曲线即氟的信号曲线，可看到探头的死时间大约在50s左右，由于固体T 2特别短，少于目前NMR探头的死时间，所以可判断图2中的信号为包裹在线圈外的氟化物信号。采样时间t 即(TD/SW)，采集时间对回波信号的采集有影响。采集信号的选取应该根据是否把信号全部采集完作为标准。 图 3.FID 信号因为小的水分子运动很快，对于纯液体，质子难于丢失的能量。处于高能级状态的质子不能把他们的能量迅速传递给周
19、围晶格，只能慢慢地回到原来的低能级水平，这意味着液体有较长的纵向弛豫时间。一般驰豫时间取决于T 1或T 2中小的一个。由于固体T 2特别短少于目前 NMR探头的死时间，所以目前NMR无法检测固体氟化物的信号( 仪器还没来得及检测到核吸收能量的信号，核就已经回到了低能级) ，所以只能测试牙膏中可溶氟的含量。实验运用核磁共振来测量牙膏含氟量，但是不可能通过简单的测量牙膏就获得其含氟量。我们需要有一种标准溶剂即已知其含氟量的物质，通过测量不同数量该物质的信号，最终统计得出含氟量于该物质的信号的曲线，这里称之为含氟量标准曲线。然后再对牙膏进行测量，在获得牙膏氟信号后，与标准曲线上的信号做对比，通过两者
20、的比值就能够清楚的得到牙膏含氟量 6。三、测量过程1.标准样溶液选取与配制在选取标准溶剂时，首先考虑聚四氟乙烯，氟化钙和氟化钠，前两者虽然可以在浓度和安全性方面胜任，但有一个关键问题，两者的溶解性都不佳，在配制溶液时，很可能只能配得少量极稀的溶液，对于信号极低的氟来讲这点是很不利的，所以决定采用溶解性浓度安全性各方面都比较好的 NaF。作为用来定标的曲线，在配制时一定要保证 NaF 完全溶于水以提高精度。为了获得标准曲线，分别取 0.116g，0.232g，0.464g，0.696g，0.927g NaF，用烧杯称量定量水，在容量瓶配制成 100g 溶液，用玻璃棒充分搅拌直至溶液中不再留有悬浮
21、颗粒，最终配得标配的含氟量为 0.05%，0.1% ，0.2%，0.3%，0.4% ，的 NaF 溶液。2.标准样信号测量分别取 D3=50s，D 3=200s 对不同浓度的 NaF 进行测试记录数据，每个样品测量 3 次，取平均值。实验数据如下：表 1 不同浓度 NaF 溶液分别在 D3=50S 和 200S 时的氟信号强度5图 4. D3=50S 时的信号-浓度关系 图 5. D3=200S 时的信号-浓度关系 3.刻度曲线绘制最小二乘法运算最小二乘法代码在 MATLAB 平台下编写。最后拟合出来的刻度直线如图 6 所示。图 6.最小二乘法拟合出的刻度直线如图 6 用最小二乘法拟合得到两曲
22、线的方程如下：D3=50s,Y=18768000X+352000D3=200s, Y=18117000X+57000氟浓度 D3=50S 氟信号强度 D3=200S 氟信号强度0 344500 535500.05% 371500 755000.1% 401000 955000.2% 431000 1325000.3% 454000 1655000.4% 506000 200000D3=50 s时 氟 信 号 强 度344500371500401000431000 45400050600001000002000003000004000005000006000000.0% 0.1% 0.2% 0.
23、3% 0.4% 0.5%氟 浓 度氟信号强度D3=200 s时 氟 信 号 强 度535507550095500 132500165500 2000000500001000001500002000002500000.0% 0.1% 0.2% 0.3% 0.4% 0.5%氟 浓 度氟信号强度6式中 X 为氟的质量，Y 为氟的不同含量所对应的信号值。4.牙膏含氟量测量在测得标准曲线后进而测量牙膏的含氟量，测试方法同测试样本溶液基本相同，需要注意的是，在牙膏放入试管时，不能有牙膏留在试管壁上，这将造成测量数据的不准。D 3=50s 时测试得含氟牙膏信号为360000；D 3=200s 时为 6975
24、0。本次测试样品中含氟化钠约为 0.11%,经计算含氟量约为 0.047%。本样品中使用 98%纯度的 NaF 粉末，根据 Na 和 F 的化学质量通过计算：18/（18+23 ）=44%，可得NaF 的粉末中 F 的浓度为 43.12%。取 2ml 溶液质量为 2g，5 个溶液中含 F 的量分别为：0.116g、0.232g、0.464g、0.696g、0.927g。在称量配制时其中误差采用平均误差的方法进行分担消除。将测得值 Y=69750 代入 Y=18117000X+57000 中可得 D3=200s 时,含氟牙膏的含氟量为 0.045%。用相同方法计算得到当 D3=50s 时,含氟牙
25、膏的含氟量为 0.031%。四、结果分析1.测量的优缺点由于目前 NMR 无法检测固体氟化物的信号，实验可以通过对一种标准溶剂即已知其含氟量的物质，进行测量不同数量该物质的信号，获得含氟量于该物质的信号的曲线，再对牙膏进行测量获得牙膏氟信号，与标准曲线上的信号做对比，通过两者的比值就能够快速的得到牙膏含氟量。较之于传统的方法而言，可以省去实验过程中大量仪器、大量试剂的使用，无需等待太长的分析过程。但由于磁共振成像仪属于高精度实验仪器，它对温度的变化、样品的温度、称量的准确性有着一定要求，存在着一定的误差。2. 测量准确性在 D3 的选取上当我们取 D3=50s 时，由于时间过短少于目前 NMR
26、 探头的死时间导致氟信号强度受到包裹在探头外的含氟物质的影响，这时测得的氟信号强度为探头含氟物质与溶液中的氟共同作用，这样对准确测量出溶液的含氟量带来了不小的困难。此时必须对 D3 进行多次取值找到探头外含氟物质相对影响最小的情况。当 D3=50s，D 3=200s 时得到曲线方程分别为Y=18768000X+352000，Y=18117000 X+57000。从两式可以得到 K1=18768000，K 2= 18117000 , K=(K1+K2)/2=18442500, Relative error=(K1-k2)/K=3.5%,经过多次测量数据的采集和对比通过最终测试出的实验数据认为选取
27、 D3=200s 时对于测量结果较为准确。3误差来源试验中导致误差的因素有很多，如系统的不稳定性导致频率的飘移，使调节中心频率带来一定困难；由于仪器的限制在称量配制溶液时不可避免的会产生误差，在精度极高的测量实验中，这种误差一定极大程度的影响了实验的准确性，此外探头外所包裹的含 F 物质的信号也影响了测量的精确信，它对所测量的含氟物质的影响还不能准确得知；实验过程中仪器随着工作的进行温度会发生变化对测量也会带来误差；同样实验样品在随着室内环境变化，置入测量仪器中温度的变化都会发生变化。五、结论运用核磁共振能够很好的了解核磁共振测牙膏含氟量的原理与方法，运用台式核磁共振仪只需确定一次标准曲线，便
28、可简捷快速地测得牙膏的含氟量，与实际情况较符合。运用台式核磁共振仪对牙膏含氟量的探索同理也可应用于其他样品如农业苗木中、环境监测等方面的氟含量检测研究。参考文献：1 郑宝山. 用-谈中国对数字说话-含氟牙膏的需求,中华口腔医学网 EB/OL.http:/ 张金生,邱竹贤.硝酸钍滴定法分析电解质分子比的研究J.轻金属,1989.12,12(1):28-34.3 王成龙,李刚,杜新丽,王群梅.含氟牙膏中氟离子测定方法的研究J.实用口腔医学杂志,1999,1(15):30.74 容晓文.离子色谱法测定牙膏中的氟C. 第十一届全国离子色谱学术报告会论文集,浙江杭州.2006.11.5 张泽宝主编.医学
29、影像物理学M.人民卫生出版社,2006.6.16 汪红志,张学龙,武杰编著.核磁共振成像技术实验教程M.科学出版社. 2008.1.通讯地址：200093 上海市营口路 101 号医疗器械学院 张学龙 联系电话：021-65502664， 13916239058 传真电话：021-55620106电子邮件：来稿经初审，认为补充以下用“”符号标明的内容后可送专家进一步评审。1 英文标题；2 作者单位英译、邮编；作者汉语拼音；3 中文摘要(200 字左右)；4 英文摘要(200 个实词左右)；5 中文关键词(38 个)；6 英文关键词；7 中图法分类号；8 图题、表题及英译；9 第一作者简介(包括
30、出生年份、性别、民族、学历、职称、职务等)；10 文稿中示意图请采用 Word 绘图工具、AutoCAD 、Excel 或 Origin 绘制；字体：小 6 号宋体；英文和数字：Times NewRoman，图形用 1 磅线条，尺寸线用 0.75 磅线条；示意图请用黑色不要底色，原始示意图复制到 Word 时，请选用“选择性粘贴” ；原始示意图文件同时发至编辑部；11 文 章 中 公 式 请 用 公 式 编 辑 器 编 制 ， 不 要 改 动 大 小 ， 保 持 原 尺 寸 ， 5 号 字 在 公 式 用 编 辑 器 中 的 大 小 为 10.5磅 ；12 在 文 章 段 落 中 的 英 文
31、字 母 和 数 字 不 要 用 公 式 编 辑 器 编 制； 变 量 用 斜 体 表 示 ， 下 标 不 是 变 量 的 用 正 体 表 示 ；13 文稿 5 号字、字体宋体、行距 16 磅、英文和数字 Times New Roman；不用分栏；14 文中单位 t/h 等改成 th-1 的形式；15 补参考文献内容(3 位作者后省略、卷数、期数、页码、出版社、年份等)；16. 内容太少，需增加至 40005000 字左右；17 第 1 作者详细通讯地址、联系手机和电话及 e-mail；18 基金项目编号；19 请 根 据 科 技 期 刊 论 文 的 基 本 格 式 重 新 整 理 ， 提 供
32、Word2000(或 Word2000 以 上 版 本 )书 写 的 电 子 版 文档 。 再次感谢您对本刊的关心和信任。此致敬礼！能源研究与信息编辑部2010 年 11 月 23 日地 址 上海市杨浦区军工路 516 号 邮 编 2000938电 话 021-55270778 021-55272843 传 真 021-55272843Email 附：期刊： 张金铭, 沈 海. 机组负荷分配优化研究J. 汽轮机技术, 2003, 35(2): 75-80.Brauner N. Vapor absorption into falling filmJ. ASME J, 2001, 34(3): 7
33、6-82.书： 施天谟. 计算传热学M. 北京: 科学出版社, 2005.Baehr H D. Heat and mass transferM. Berlin: Springer-Verlag, 2004.论文集： 邱海林, 李 瑛, 汪卫东, 等. 煤燃烧过程中矿物质行为研究A. 李广凡. 上海市工程热物理学会 2000 年年会论文集C. 上海: 上海交通大学出版社, 2002. 59-62.Leiner A. Optical techniques for heat transfer measurementsA. Collins W. Proc IntConf on Energy and E
34、nvironmentC. New York: Academic Press, 2003. 68-71.网络版： 郝 郁. 国际 LNG 贸易合同及定价研究 EB/OL.http:/ 21/content_21/content_214448.htm, 2006-09-22.图 例 ： 字体、线条、坐标表示参考图 例 ： 字体、线条、坐标表示参考工作段计算机电锅炉铂电阻铂电阻仪表测温引风机0.01 0.04 0.07 0.10 0.13 1 3 5 7 9 112005年 监 测 点 12006年 监 测 点 12005年 监 测 点 22006年 监 测 点 20.020.060.100.140.181 3 5 7 9 112005年 监 测 点 12006年 监 测 点 12005年 监 测 点 22006年 监 测 点 2氧份/(mgm-3)PM10/(mgm-3)月 份 月 份
