论文基于非相对论量子色动力学有效理论（NRQCD）研究了重夸克偶素产生过程γγ?→ηc形状因子的相对论修正。NRQCD是一个非常成功的量子色动力学低能有效场论，它可以用来研究重夸克偶素的产生和湮灭过程。NRQCD因子化定理把一个过程因子化成一系列短程系数和长程矩阵元的乘积。其中短程系数包含了标度大于mQ部分的贡献，而小于此标度的贡献被归入非微扰的长程矩阵元。短程系数可以进行微扰计算，因此可以按耦合常数展开。而各个长程矩阵元是按重夸克偶素系统里面重夸克的特征速度υ展开的。这样NRQCD因子化公式就把一个产生过程的振幅或者截面按耦合常数和特征速度做双重展开。具体计算到哪一阶需要根据精度要求以及展开级数的收敛性来决定。本论文计算了γγ?→ηc形状因子领头阶矩阵元对应的短程系数（F）以及次领头阶矩阵元对应的短程系数(G)。其中该过程的F目前在文献已经计算到了两圈图，并结合实验做了唯象讨论。结果表明仅仅考虑相对论修正的领头阶贡献尚不能跟实验取得一致。为了检验 NRQCD因子化公式的正确性以及预言能力，利用因子化理论领域常用的匹配方法，本论文研究了该过程的次领头阶矩阵元对应的短程系数G，我们把G计算到微扰展开的次领头阶（也就是一圈图阶）。结合F和G并通过适当的选择输入参数做了详尽的唯象分析和讨论。此外论文也研究了γγ?→ηb的形状因子，未来实验可以检验理论计算的正确性，从而检验NRQCD理论。　　论文的计算过程中，利用高能软件包 FeynArts生成费曼图并生产费曼振幅；利用FeynCalc进行Dirac矩阵求积、Lorentz指标缩并以及化简表达式；利用Apart和FIRE软件包对圈图振幅进行张量约化；最后一圈图振幅只包含一些简单的标量圈积分，通过计算这些圈积分就最终得到了过程的振幅表达式。为了保证计算的正确性,我们对一些已有经典过程重新进行了计算，这其中包括：J/ψ→e+e?，ηc→γγ，e+e?→ηc+γ的αs及αsυ2修正，计算结果与已有文献一致。