Three-way instrumental data offer the second-order advantage to analysts, a property of great utility in the field of complex sample analysis in the presence of unsuspected components as potential interferents. The available multivariate methodologies for obtaining this advantage are all based on linear models, and hence they are not applicable to spectral information behaving in a non-linear manner with respect to target analyte concentrations. This work describes the combination of a backpropagation artificial neural network model with a technique known as residual bilinearization, applicable to second-order spectral information. The joint model allows one to efficiently extract analyte concentrations from intrinsically non-linear data, even in the presence of unsuspected constituents. Simulations have been performed by mimicking deviations from linearity brought about by: (1) exponential relationship between fluorescence and concentration, (2) kinetic evolution of responsive reaction products and (3) analytes acting as reaction catalysts. In all of these cases, successful prediction of the analyte concentrations was achieved on large test sample sets, which included the presence of overlapping components not included in the training step. The new method not only obtains the second-order advantage, but also correctly retrieves the contribution of the unsuspected components to the total test sample signals. The comparison with a multivariate methodology based on partial least-squares regression with second-order advantage shows that the presently described method displays better predictive ability.