D值法计算