随着LHCb的运行,更多的B介子衰变道可以得到更加详细的研究。升级改造后的SuperKEKB也将产生大量的Bs介子,粗略的估计,每年大概能产生5.9×108的Bs事例数。这么多的数据足够我们对Bs介子进行深入的研究和分析。通过对Bs介子衰变过程的分支比,CP对称性破缺,动量分布等的研究,使得标准模型和新物理面临一系列新的更严峻的挑战。本文用pQCD因子化方法研究了W交换的稀有衰变过程Bs→η(J/ψ)D。我们首先简单介绍了B物理研究常用的理论工具和强子矩阵元的计算方法,对pQCD因子化方法进行了详细的介绍。在工作部分,我们用pQCD方法研究了Bs→ηc(J/ψ)D衰变过程,采用了三种Bs介子的波函数,两种ηc(J/ψ)介子的波函数,给出了数值结果,并进行分析和讨论。首先分析说明了这个过程是以硬胶子交换为主的,即pQCD方法对Bs→η(J/ψ)。衰变过程是适用的。我们讨论Bs介子波函数中ΦB+和ΦB-的贡献,ηc(J/ψ)介子中扭度为2和扭度为3的贡献,以及不同拓扑的贡献等。结果表明Bs→ηcD,ηcD,J/ψD,J/ψD的分支比,分别在10,10,10,109的量级。我们也讨论了各种波函数以及硬标度等的不确定性。其中,理论预言的不确定性主要来自于非微扰部分波函数的选取。