A new architecture for Spectral Amplitude Coding Optical Code Division Multiple Access (SAC-OCDMA) system based on two Dimensional Multi Diagonal (2D-MD) codes named 2D-MD SAC-OCDMA and utilizing a laser optical source is proposed for Long-Reach Passive Optical Network (LR-PON). In this work, a computer simulator tool is used, for the first time, as a SAC-OCDMA simulation set-up utilizing the unique combination of a coherent laser array and 2D-MD codes. In addition, the system performance is addressed numerically by taking into account Multiple Access Interference (MAI), optical coherent source noise, first, second and third order fiber dispersion, nonlinear effects and photo-detector noise. Simulation results indicate that for a single user (i.e., without considering MAI), the system can operate at a maximum bit rate of 55 Gb/s over 250 km of Single Mode Fiber (SMF), with a Bit Error Rate (BER) below 10−9 (Q-limit = 15.5 dB), when only first order fiber dispersion is considered. However, including the effects of second and third order fiber dispersion as frequency domain parameters, results in a reduction of the maximum bit rate to 40 Gb/s, while maintaining a Q-factor above the Q-limit under the same transmission distance. Furthermore, we demonstrate that the proposed architecture extends the SMF transmission reach up to 600 km and 480 km, when considering linear and nonlinear effects, respectively. Finally, we show that our proposed 2D-MD SAC-OCDMA system outperforms existing solutions presented in the literature for LR-PON configuration, in terms of both aggregate bit rate and transmission reach.