Ahmad,
To be clear, your goal is not to identify differences across runs at the group level, right? In that case, you can simply average each subject's effect estimate across the available runs, and use 3dttest++, or use 3dMEMA if you average them within 3dDeconvolve/3dREMLfit
Gang