《数字PCR仪中荧光信号的检测及处理研究》-郭成祥硕士论文
武汉工程大学硕士学位论文，机械电子工程专业，2017年9月 导师：陈绪兵 教授、王瑜辉 高工 研究方向：物联网技术及应用
核心内容
本文针对微滴式数字PCR（ddPCR）检测系统中微液滴荧光信号的检测与处理进行了系统研究，分为三大板块：荧光信号数据采集硬件设计、峰值检测算法、荧光补偿算法。基于FPGA+USB2.0搭建了完整的采集系统平台，并对信号处理全流程（数据段截取→基线校准→滤波去噪→峰值提取→荧光补偿）进行了算法对比研究。
第一章 绪论
1.1 研究背景及意义
- PCR技术是基因工程核心环节，广泛应用于医学诊断、食品安全、刑侦检测
- 数字PCR（dPCR）是继传统PCR、qPCR之后的第三代技术，可绝对定量DNA拷贝数
- 全球数字PCR仪市场被4家外企垄断（Bio-Rad、RainDance、Life Technologies、Fluidigm），亟需国产替代
- 本文目标：深入研究数字PCR仪中荧光信号检测与处理的关键技术
1.2 国内外研究现状
PCR技术发展三代：
- 传统PCR — 琼脂糖电泳分析，操作复杂，仅定性
- 实时荧光定量PCR（qPCR） — 相对定量，需标准曲线和参考基因，Ct值易受扩增效率等因素干扰
- 数字PCR（dPCR） — 绝对定量，不依赖Ct值，灵敏度/特异性/精确性更高
数字PCR定量原理：
- 将样品分散到大量微反应单元中，每个单元含0个或多个DNA分子
- PCR扩增后，阳性单元（含模板）发强光，阴性单元发弱光
- 统计阳性数和总单元数，结合泊松分布计算起始浓度
- 公式：C = -ln(1 - p) / v（p为阳性比例，v为单滴体积）
商用数字PCR产品分类：
| 类型 | 代表产品 | 厂商 | 反应单元数 | 荧光通道数 |
|---|---|---|---|---|
| 微反应室式 | Bio-Mark系统 | Fluidigm (2006) | — | — |
| 微反应室式 | QuantStudio 3D | Life Tech (2013) | 20,000/芯片 | 2重 |
| 微液滴式 | QX100/QX200 | Bio-Rad (2011/2013) | 20,000 | 2重 |
| 微液滴式 | Raindrop System | RainDance (2012) | 10,000,000 | 10重 |
1.3 论文章节安排
- 绪论（背景/现状/主要工作）
- 微液滴生成方法与荧光信号检测原理
- 荧光信号采集系统硬件设计与分析
- 荧光信号峰值检测方法研究
- 荧光信号峰值补偿算法研究
- 总结与展望
第二章 微液滴生成方法与荧光信号检测原理
2.1 微液滴生成方法对比
四大类生成技术： 水动力法、气动法、光学控制法、电流激发法 数字PCR系统选用水动力法，又分三种方式：
| 方法 | 原理 | 优点 | 缺点 |
|---|---|---|---|
| T型通道法 | 样品与油在T形交汇处受剪切力生成液滴 | 技术难度低、芯片工艺要求低 | 液滴大小不均一、速度慢 |
| 共流聚焦法 | 油在毛细管中包裹样品形成液滴 | 液滴大小均匀、尺寸更小 | 技术难度高、油品要求高、芯片工艺要求高 |
| 流动聚焦法 | 连续相和分散相在交汇处通过接触力形成液滴 | 生成速度快、大小均匀、技术难度适中 | 油相消耗量大、芯片工艺要求较高 |
本文选用流动聚焦法，综合考虑液滴质量、生成速度、技术难度和成本。 实验装置：PDMS材质微液滴生成芯片 + 电磁阀（控制流向流速） + 下压装置
- 射流孔直径：90~120 μm
- 生成液滴直径：约100 μm
2.2 荧光信号检测系统原理
荧光产生机理： 荧光染料分子吸收激光能量→激发态→跃迁回基态→释放荧光（斯托克斯位移：发射波长 > 激发波长）
激发光源选型： 激光光源（宽光谱光源 vs 激光光源，后者聚焦性好、能量密度高、灵敏度高）
光学检测系统工作原理：
- 激光经滤波/反射后聚焦到检测芯片
- 微液滴流经激光斑点时，PCR扩增后的荧光探针被激发
- 荧光经物镜 + 二向色镜分流至两个PMT探测器（FAM通道和VIC通道）
- PMT将光信号转化为电信号
- 两个光电二极管监控激光稳定性，反馈控制激光电源
- 实验采用荧光染料：FAM 和 VIC
第三章 荧光信号采集系统硬件设计
3.1 系统组成
三大功能模块：
- PMT光电传感器 + A/D转换模块 — 模拟信号调理 + 数字化
- FPGA主控制 + DDR2 SDRAM数据缓存模块 — 采集控制 + 数据缓冲
- USB数据传输模块 — 向上位机传输
总体架构： FPGA为主控制器 → DDR2 SDRAM缓存 → USB 2.0高速传输 → 上位机处理
采集系统特殊要求：
- 两路荧光信号（FAM + VIC）
- 对采集速度和传输速度要求高（微液滴数量庞大）
3.2 主要芯片选型
| 模块 | 芯片型号 | 说明 |
|---|---|---|
| FPGA主控制器 | Altera Cyclone III EP3C40Q | 39,600个逻辑单元，4个PLL，18个嵌入式乘法器，535个I/O |
| USB接口 | Cypress CY7C68013A | USB 2.0高速（480Mbps），增强型8051内核，GPIF可编程接口 |
| PMT探测器 | — | 光电倍增管（型号未明确，FAM/VIC双通道） |
| A/D转换 | — | 通过I2C配置（具体型号未详述） |
| 数据缓存 | DDR2 SDRAM | 高速大容量缓存 |
接口传输速度对比： USB 2.0 (480Mbps) > 并口 > 串口，综合成本和速度选USB 2.0
3.3 核心电路模块设计
- PMT信号调理电路 — 光电倍增管输出信号的放大、滤波等调理
- FPGA控制电路及配置电路 — 系统核心，总体控制 + 数据处理缓存
- USB接口电路 — CY7C68013A与FPGA的连接电路
- 时钟电路 — 为各模块提供时钟源
FPGA主要功能：
- 通过串口调节PMT电源，控制信号采集
- 通过I2C配置A/D芯片，实现模数转换
- 控制荧光信号数据的随机存储（DDR2 SDRAM）
- 通过USB芯片向上位机传输数据
第四章 荧光信号峰值检测方法研究
峰值检测四步法：有效数据段截取 → 基线校准 → 滤波处理 → 有效峰值提取
4.1 有效数据段的截取
- 从连续采集的数据流中定位包含微液滴荧光信号的有效数据段
- 为后续基线校准和峰值提取缩小处理范围
4.2 基线校准
方法一：滑动窗口分段法
- 将信号分段，每段用线性拟合基线
- 简单快速，但精度有限
方法二：基于补偿最小二乘法的基线校准
- 在传统最小二乘基础上引入惩罚项（加权系数）
- 拟合波峰时权重为0，拟合基线时权重为1
- 张志敏等人提出的迭代重加权惩罚最小二乘法进一步改进
- 本文提出了更精确、计算复杂度更低的加权系数计算公式
- 效果：基线拟合比滑动窗口法更准确
4.3 滤波处理
噪声来源：环境光噪声 + PMT元件噪声 + 光学通道折射干扰 + 电路板噪声
方法一：小波包分解法去噪
- 小波包是小波函数的扩展，对高频和低频都进行分解
- 步骤：选择滤波器→小波包变换→选取代价函数找最佳基→处理系数→重构信号
- 分解层数越多，信噪比越大，但误差也越大
方法二：小波阈值法去噪
- 基本思路：大幅值小波系数为有效信号，小幅值为噪声，设阈值分割
- 硬阈值：绝对值小于阈值置零，大于保留（不连续，振荡）
- 软阈值：大于阈值的系数减去阈值（更平滑，效果更好）
- 步骤：分解信号→估计噪声方差σ→计算阈值→软/硬阈值处理→重构
对比结论： 小波软阈值法去噪得到的信号波峰更平滑、峰间拉开效果更好，优于小波包分解法
4.4 有效峰值提取
算法步骤：
- 计算两通道滤波后信号的一阶导数（数字信号用差分近似）
- 求两通道信号平均值
- 结合一阶导数和平均值提取有效峰值
- 满足两个条件：
- 一阶导数过零点（极大值/极小值点）
- 信号值大于通道均值（剔除极小值和异常峰）
- 剔除重叠峰、不规则峰，保留规则平滑的波峰
目的： 为后续荧光补偿计算提供高质量的峰值数据
第五章 荧光信号峰值补偿算法研究
5.1 荧光补偿矩阵
串扰问题： FAM和VIC两种荧光染料的发射光谱有重叠，一个通道的PMT会接收到另一种染料的泄漏荧光，导致峰值数据有偏差。
荧光泄漏模型（2×2矩阵）：
PMT1 = (1-k12)·FL1 + k21·FL2 + a1
PMT2 = k12·FL1 + (1-k22)·FL2 + a2
其中：
- PMT1、PMT2：两个检测器实际测到的信号
- FL1、FL2：FAM和VIC的真实荧光强度（待求）
- k12：FAM泄漏到PMT2的比例；k21：VIC泄漏到PMT1的比例
- a1、a2：自发荧光/背景噪声
核心问题： 已知PMT测量值，求FL真实值 → 需要估计补偿矩阵参数
5.2 递推最小二乘法（RLS）
一般最小二乘法（LS）： 批量处理所有数据，一次性求解最优参数矩阵
递推最小二乘法（RLS）：
- 新估计值 = 旧估计值 + 修正项
- 可边采集边计算，逐次更新参数，直到满足精度
- 适合实时处理场景
- MIMO系统模型下的递推公式组：
- P(k) = P(k-1) - P(k-1)Φ(k-1)T(k-1) （协方差更新）
- K(k) = P(k-1)Φ(k-1) / [1 + Φ^T(k-1)P(k-1)Φ(k-1)] （增益向量）
- θ̂(k) = θ̂(k-1) + P(k)Φ(k-1)[Y(k) - Φ^T(k-1)θ̂(k-1)] （参数更新）
5.3 实验分析
实验设计：
- 背景荧光校准 — 空白微液滴（不含DNA）测背景均值
- 纯染料校准 — 9个梯度浓度FAM/VIC单色染料，各检测10,000个微液滴
- 混合样本检测 — 同时标记两种荧光的微液滴，检测10,000个
- 分别用LS和RLS计算补偿矩阵，对比参数误差
评价指标：
- 参数误差 — 与标准值的偏差
- 变异系数（CV, Coefficient of Variation） — 荧光染料分辨率的标准指标 CV = 标准差 / 均值 × 100%
实验结论：
- 递推最小二乘法（RLS）的CV值更接近标准数据的CV值
- RLS的参数误差更小
- RLS能更实时有效地估算荧光补偿矩阵
第六章 总结与展望
6.1 全文总结
三大成果：
- 硬件设计 — 完成基于FPGA+USB2.0+DDR2的微液滴荧光信号采集系统设计，包括芯片选型和核心电路设计
- 峰值检测四步法 — 有效数据段截取 + 基线校准（补偿最小二乘法优于滑动窗口法） + 滤波去噪（小波软阈值法优于小波包分解法） + 有效峰值提取
- 荧光补偿算法 — 对比一般最小二乘和递推最小二乘，实验证明RLS算法参数误差更小、CV值更优，适合实时荧光补偿矩阵估计
6.2 工作展望
- 硬件方面：仅研究了核心电路，激光控制模块和液滴生成器控制模块待研究
- USB传输：仅完成硬件设计，USB驱动和上位机客户端程序待开发
- 信号处理：仅对比了两类常用方法，可尝试更多信号处理方法
- 系统集成：将各模块整合为完整的数字PCR检测系统
关键技术提炼
| 技术点 | 方法/选型 | 结论 |
|---|---|---|
| 微液滴生成 | 流动聚焦法（水动力） | 速度快、均一性好，综合最优 |
| 激发光源 | 激光光源 | 聚焦性好、能量密度高、灵敏度高 |
| 检测通道 | FAM + VIC 双通道PMT | 双色荧光检测 |
| 主控制器 | FPGA Cyclone III EP3C40Q | 并行处理、可配置、适合高速采集 |
| 数据传输 | USB 2.0 (CY7C68013A) | 480Mbps高速，平衡成本与性能 |
| 基线校准 | 补偿最小二乘法 > 滑动窗口法 | 补偿法基线更准确 |
| 信号去噪 | 小波软阈值法 > 小波包分解法 | 波形更平滑，峰间分离效果好 |
| 荧光补偿 | RLS递推最小二乘法 > LS一般最小二乘法 | 参数误差更小，CV更优，支持实时更新 |
关联笔记
- 数字PCR荧光检测仪数据分析方法的设计与实现-李鑫鑫-硕士论文 — 同方向另一篇硕士论文，侧重软件和k-means聚类
- 荧光溢出补偿-数字PCR多色荧光检测校正技术 — 数字PCR荧光补偿技术文档
- 多重定量PCR系统中多色荧光检测和光谱串扰校正方法-臧留琴 — 磁光开关+四维聚类迭代串扰校正
- YY/T 1173-2010 聚合酶链反应分析仪行业标准 — PCR仪行业标准
- FCS 3.1 流式细胞术数据文件标准 — 流式细胞/数字PCR相关的荧光检测数据标准